The Born–Oppenheimer approximation underlies much of chemical simulation and provides the framework defining the potential energy surfaces that are used for much of our pictorial understanding of ...chemical phenomena. However, this approximation breaks down when the dynamics of molecules in excited electronic states are considered. Describing dynamics when the Born–Oppenheimer approximation breaks down requires a quantum mechanical description of the nuclei. Chemical reaction dynamics on excited electronic states is critical for many applications in renewable energy, chemical synthesis, and bioimaging. Furthermore, it is necessary in order to connect with many ultrafast pump–probe spectroscopic experiments. In this review, we provide an overview of methods that can describe nonadiabatic dynamics, with emphasis on those that are able to simultaneously address the quantum mechanics of both electrons and nuclei. Such ab initio quantum molecular dynamics methods solve the electronic Schrödinger equation alongside the nuclear dynamics and thereby avoid the need for precalculation of potential energy surfaces and nonadiabatic coupling matrix elements. Two main families of methods are commonly employed to simulate nonadiabatic dynamics in molecules: full quantum dynamics, such as the multiconfigurational time-dependent Hartree method, and classical trajectory-based approaches, such as trajectory surface hopping. In this review, we describe a third class of methods that is intermediate between the two: Gaussian basis set expansions built around trajectories.
The increasing availability of high-quality experimental data and first-principles calculations creates opportunities for developing more accurate empirical force fields for simulation of proteins. ...We developed the AMBER-FB15 protein force field by building a high-quality quantum chemical data set consisting of comprehensive potential energy scans and employing the ForceBalance software package for parameter optimization. The optimized potential surface allows for more significant thermodynamic fluctuations away from local minima. In validation studies where simulation results are compared to experimental measurements, AMBER-FB15 in combination with the updated TIP3P-FB water model predicts equilibrium properties with equivalent accuracy, and temperature dependent properties with significantly improved accuracy, in comparison with published models. We also discuss the effect of changing the protein force field and water model on the simulation results.
We develop an extension of the variational quantum eigensolver (VQE) algorithm-multistate contracted VQE (MC-VQE)-that allows for the efficient computation of the transition energies between the ...ground state and several low-lying excited states of a molecule, as well as the oscillator strengths associated with these transitions. We numerically simulate MC-VQE by computing the absorption spectrum of an ab initio exciton model of an 18-chromophore light-harvesting complex from purple photosynthetic bacteria.
Few would dispute that theoretical chemistry tools can now provide keen insights into chemical phenomena. Yet the holy grail of efficient and reliable prediction of complex reactivity has remained ...elusive. Fortunately, recent advances in electronic structure theory based on the concepts of both element- and rank-sparsity, coupled with the emergence of new highly parallel computer architectures, have led to a significant increase in the time and length scales which can be simulated using first principles molecular dynamics. This opens the possibility of new discovery-based approaches to chemical reactivity, such as the recently proposed ab initio nanoreactor. We argue that due to these and other recent advances, the holy grail of computational discovery for complex chemical reactivity is rapidly coming within our reach.
Hybrid quantum mechanical–molecular mechanical (QM/MM) simulations are widely used in studies of enzymatic catalysis. Until recently, it has been cost prohibitive to determine the asymptotic limit of ...key energetic and structural properties with respect to increasingly large QM regions. Leveraging recent advances in electronic structure efficiency and accuracy, we investigate catalytic properties in catechol O-methyltransferase, a prototypical methyltransferase critical to human health. Using QM regions ranging in size from reactants-only (64 atoms) to nearly one-third of the entire protein (940 atoms), we show that properties such as the activation energy approach within chemical accuracy of the large-QM asymptotic limits rather slowly, requiring approximately 500–600 atoms if the QM residues are chosen simply by distance from the substrate. This slow approach to asymptotic limit is due to charge transfer from protein residues to the reacting substrates. Our large QM/MM calculations enable identification of charge separation for fragments in the transition state as a key component of enzymatic methyl transfer rate enhancement. We introduce charge shift analysis that reveals the minimum number of protein residues (approximately 11–16 residues or 200–300 atoms for COMT) needed for quantitative agreement with large-QM simulations. The identified residues are not those that would be typically selected using criteria such as chemical intuition or proximity. These results provide a recipe for a more careful determination of QM region sizes in future QM/MM studies of enzymes.
Repetitive DNA sequences are abundant in a broad range of species, from bacteria to mammals, and they cover nearly half of the human genome. Repeats have always presented technical challenges for ...sequence alignment and assembly programs. Next-generation sequencing projects, with their short read lengths and high data volumes, have made these challenges more difficult. From a computational perspective, repeats create ambiguities in alignment and assembly, which, in turn, can produce biases and errors when interpreting results. Simply ignoring repeats is not an option, as this creates problems of its own and may mean that important biological phenomena are missed. We discuss the computational problems surrounding repeats and describe strategies used by current bioinformatics systems to solve them.
Celotno besedilo
Dostopno za:
DOBA, IJS, IZUM, KILJ, NUK, PILJ, PNG, SAZU, UILJ, UKNU, UL, UM, UPUK
The spinal dorsal horn receives input from primary afferent axons, which terminate in a modality-specific fashion in different laminae. The incoming somatosensory information is processed through ...complex synaptic circuits involving excitatory and inhibitory interneurons, before being transmitted to the brain via projection neurons for conscious perception. The dorsal horn is important, firstly because changes in this region contribute to chronic pain states, and secondly because it contains potential targets for the development of new treatments for pain. However, at present, we have only a limited understanding of the neuronal circuitry within this region, and this is largely because of the difficulty in defining functional populations among the excitatory and inhibitory interneurons. The recent discovery of specific neurochemically defined interneuron populations, together with the development of molecular genetic techniques for altering neuronal function in vivo, are resulting in a dramatic improvement in our understanding of somatosensory processing at the spinal level.