Quantum Mechanics/Fluctuating Charge Protocol to Compute Solvatochromic Shifts.
Ontology highlight
ABSTRACT:
SUBMITTER: Ambrosetti M
PROVIDER: S-EPMC8582258 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:The coordination of the magnesium ion in proteins by triphosphates plays an important role in catalytic hydrolysis of GTP or ATP, either in signal transduction or energy conversion. For example, in Ras the magnesium ion contributes to the catalysis of GTP hydrolysis. The cleavage of GTP to GDP and P(i) in Ras switches off cellular signaling. We analyzed GTP hydrolysis in water, Ras, and Ras·Ras-GTPase-activating protein using quantum mechanics/molecular mechanics simulations. By comparison of the theoretical IR-difference spectra for magnesium ion coordinated triphosphate to experimental ones, the simulations are validated. We elucidated thereby how the magnesium ion contributes to catalysis. It provides a temporary storage for the electrons taken from the triphosphate and it returns them after bond cleavage and P(i) release back to the diphosphate. Furthermore, the Ras·Mg(2+) complex forces the triphosphate into a stretched conformation in which the ?- and ?-phosphates are coordinated in a bidentate manner. In this conformation, the triphosphate elongates the bond, which has to be cleaved during hydrolysis. Furthermore, the ?-phosphate adopts a more planar structure, driving the conformation of the molecule closer to the hydrolysis transition state. GTPase-activating protein enhances these changes in GTP conformation and charge distribution via the intruding arginine finger.
| S-EPMC3400779 | biostudies-literature
Project description:We present the ProCS method for the rapid and accurate prediction of protein backbone amide proton chemical shifts--sensitive probes of the geometry of key hydrogen bonds that determine protein structure. ProCS is parameterized against quantum mechanical (QM) calculations and reproduces high level QM results obtained for a small protein with an RMSD of 0.25 ppm (r = 0.94). ProCS is interfaced with the PHAISTOS protein simulation program and is used to infer statistical protein ensembles that reflect experimentally measured amide proton chemical shift values. Such chemical shift-based structural refinements, starting from high-resolution X-ray structures of Protein G, ubiquitin, and SMN Tudor Domain, result in average chemical shifts, hydrogen bond geometries, and trans-hydrogen bond ((h3)J(NC')) spin-spin coupling constants that are in excellent agreement with experiment. We show that the structural sensitivity of the QM-based amide proton chemical shift predictions is needed to obtain this agreement. The ProCS method thus offers a powerful new tool for refining the structures of hydrogen bonding networks to high accuracy with many potential applications such as protein flexibility in ligand binding.
| S-EPMC3877219 | biostudies-literature
Project description:Selection of residues and other molecular fragments for inclusion in the quantum mechanics (QM) region for QM/molecular mechanics (MM) simulations is an important step for these calculations. Here, we present an approach that combines protein sequence/structure evolution and electron localization function (ELF) analyses. The combination of these two analyses allows the determination of whether a residue needs to be included in the QM subsystem or can be represented by the MM environment. We have applied this approach on two systems previously investigated by QM/MM simulations, 4-oxalocrotonate tautomerase (4OT) and ten-eleven translocation-2 (TET2), that provide examples where fragments may or may not need to be included in the QM subsystem. Subsequently, we present the use of this approach to determine the appropriate QM subsystem to calculate the minimum energy path (MEP) for the reaction catalyzed by human DNA polymerase λ (Polλ) with a third cation in the active site. Our results suggest that the combination of protein evolutionary and ELF analyses provides insights into residue/molecular fragment selection for QM/MM simulations.
| S-EPMC8477969 | biostudies-literature
Project description:Magic angle spinning NMR spectroscopy is uniquely suited to probe the structure and dynamics of insoluble proteins and protein assemblies at atomic resolution, with NMR chemical shifts containing rich information about biomolecular structure. Access to this information, however, is problematic, since accurate quantum mechanical calculation of chemical shifts in proteins remains challenging, particularly for 15NH. Here we report on isotropic chemical shift predictions for the carbohydrate recognition domain of microcrystalline galectin-3, obtained from using hybrid quantum mechanics/molecular mechanics (QM/MM) calculations, implemented using an automated fragmentation approach, and using very high resolution (0.86 Å lactose-bound and 1.25 Å apo form) X-ray crystal structures. The resolution of the X-ray crystal structure used as an input into the AF-NMR program did not affect the accuracy of the chemical shift calculations to any significant extent. Excellent agreement between experimental and computed shifts is obtained for 13C?, while larger scatter is observed for 15NH chemical shifts, which are influenced to a greater extent by electrostatic interactions, hydrogen bonding, and solvation.
| S-EPMC5892201 | biostudies-literature
Project description:The solvent Stark effect on the spectral shifts of anthracene is studied with temperature-dependent solvatochromic measurements. The Stark contribution ?vStark to the absorption shift ?vp in polar solvents is measured to be ?vStark =(53±35)?cm-1 , in reasonable agreement with dielectric continuum theory estimate of 28?cm-1 , whereas the major shift ?vp ?300?cm-1 presumably originates from the solute quadrupole. We pay attention to the accurate correction of ?vp for the nonpolar contribution that is crucial when the shifts are modest in magnitude.
| S-EPMC7187296 | biostudies-literature
Project description:In quantum mechanics, the theory of quantum transitions is grounded on the convergence of a series of time-dependent perturbation theory. In nuclear and atomic physics, this series converges because the dynamics of quantum transitions (quantum jumps) are absent by definition. In molecular and chemical physics, on the contrary, the dynamics of "quantum" transitions, being determined by the joint motion of a light electron (or electrons) and very heavy nuclei, are present by definition, and the series of time-dependent perturbation theory becomes singular. An exception is the dynamic problem for stationary states in the Born-Oppenheimer adiabatic approximation, when the electronic subsystem turns out to be "off" from the general dynamic process and therefore is not dynamically full-fledged: it only forms an electric potential in which the nuclei oscillate. Removing the aforementioned singularity can be accomplished in two ways. The first method was consisted of introducing an additional postulate in the form of the Franck-Condon principle into molecular quantum mechanics, in which the adiabatic approximation is used. The second method was proposed by the author and consisted of damping the singular dynamics of the joint motion of an electron and nuclei in the intermediate (transient) state of molecular "quantum" transitions by introducing chaos. This chaos arises only during molecular quantum transitions and is called dozy chaos. Formally, the damping is carried out by replacing an infinitely small imaginary addition in the spectral representation of the complete Green's function of the system with its finite quantity. The damping chaos (dozy chaos) leads to the continuity of the energy spectrum in the molecular transient state, which is a sign of classical mechanics. Meanwhile, the initial and final states of the molecule obey quantum mechanics in the adiabatic approximation. Molecular quantum mechanics, which takes into account the chaotic dynamics of the transient state of molecular "quantum" transitions, can be called quantum-classical (dozy-chaos) mechanics. The efficacy of the damping for the aforementioned singularity was previously shown by dozy-chaos mechanics of elementary electron transfers in condensed matter, which is the simplest case of dozy-chaos mechanics, and its applications to a whole number of problems, especially to the optical spectra in polymethine dyes and their aggregates. This paper provides a regular exposition of this dozy-chaos (quantum-classical) mechanics of the elementary electron transfers. The main results of its applications presented in the introduction are also described.
| S-EPMC6909101 | biostudies-literature
Project description:An expression of unexpected simplicity is derived for the shift in optical transition energies of solute molecules in nonpolar solvents. The expression reveals a new spectroscopic rule that says: The higher the excited state of the solute, the larger the solvatochromic red shift. A puzzle formulated >50 years ago by Bayliss is solved. Bayliss, based on arguments from classical physics, assumed that the shift scales with the oscillator strength of the solute transition, but noted strong quantitative deviations from this rule in experiments. As the present expression shows, the shift does not depend on the oscillator strength of the transition, but reflects the change in dispersive solute-solvent interactions between the ground and excited states of the solute, that are determined by the anisotropy of intramolecular electron correlation. The theory is applied to explain the solvatochromic shifts of the two lowest electronic excitations of bacteriochlorophyll a and bacteriopheophytin a.
| S-EPMC2533174 | biostudies-other
Project description:We present a toolkit that allows for the preparation of QM/MM input files from a conformational ensemble of molecular geometries. The package is currently compatible with trajectory and topology files in Amber, CHARMM, GROMACS and NAMD formats, and has the possibility to generate QM/MM input files for Gaussian (09 and 16), Orca (≥4.0), NWChem and (Open)Molcas. The toolkit can be used in command line, so that no programming experience is required, although it presents some features that can also be employed as a python application programming interface. We apply the toolkit in four situations in which different electronic-structure properties of organic molecules in the presence of a solvent or a complex biological environment are computed: the reduction potential of the nucleobases in acetonitrile, an energy decomposition analysis of tyrosine interacting with water, the absorption spectrum of an azobenzene derivative integrated into a voltage-gated ion channel, and the absorption and emission spectra of the luciferine/luciferase complex. These examples show that the toolkit can be employed in a manifold of situations for both the electronic ground state and electronically excited states. It also allows for the automatic correction of the active space in the case of CASSCF calculations on an ensemble of geometries, as it is shown for the azobenzene derivative photoswitch case.
| S-EPMC10107847 | biostudies-literature
Project description:Non-relativistic quantum mechanics is reformulated here based on the idea that relational properties among quantum systems, instead of the independent properties of a quantum system, are the most fundamental elements to construct quantum mechanics. This idea, combining with the emphasis that measurement of a quantum system is a bidirectional interaction process, leads to a new framework to calculate the probability of an outcome when measuring a quantum system. In this framework, the most basic variable is the relational probability amplitude. Probability is calculated as summation of weights from the alternative measurement configurations. The properties of quantum systems, such as superposition and entanglement, are manifested through the rules of counting the alternatives. Wave function and reduced density matrix are derived from the relational probability amplitude matrix. They are found to be secondary mathematical tools that equivalently describe a quantum system without explicitly calling out the reference system. Schrödinger Equation is obtained when there is no entanglement in the relational probability amplitude matrix. Feynman Path Integral is used to calculate the relational probability amplitude, and is further generalized to formulate the reduced density matrix. In essence, quantum mechanics is reformulated as a theory that describes physical systems in terms of relational properties.
| S-EPMC6127378 | biostudies-literature
Project description:We raise fundamental questions about the very meaning of conservation laws in quantum mechanics, and we argue that the standard way of defining conservation laws, while perfectly valid as far as it goes, misses essential features of nature and has to be revisited and extended.
| S-EPMC7817196 | biostudies-literature