Mechanical/molecular Mechanical (molecular + mechanical)

Distribution by Scientific Domains

Kinds of Mechanical/molecular Mechanical

  • quantum molecular mechanical

  • Terms modified by Mechanical/molecular Mechanical

  • molecular mechanical methods

  • Selected Abstracts


    QM/MM calculation of solvent effects on absorption spectra of guanine

    JOURNAL OF COMPUTATIONAL CHEMISTRY, Issue 1 2010
    Maja Parac
    Abstract Electronic spectra of guanine in the gas phase and in water were studied by quantum mechanical/molecular mechanical (QM/MM) methods. Geometries for the excited-state calculations were extracted from ground-state molecular dynamics (MD) simulations using the self-consistent-charge density functional tight binding (SCC-DFTB) method for the QM region and the TIP3P force field for the water environment. Theoretical absorption spectra were generated from excitation energies and oscillator strengths calculated for 50 to 500 MD snapshots of guanine in the gas phase (QM) and in solution (QM/MM). The excited-state calculations used time-dependent density functional theory (TDDFT) and the DFT-based multireference configuration interaction (DFT/MRCI) method of Grimme and Waletzke, in combination with two basis sets. Our investigation covered keto-N7H and keto-N9H guanine, with particular focus on solvent effects in the low-energy spectrum of the keto-N9H tautomer. When compared with the vertical excitation energies of gas-phase guanine at the optimized DFT (B3LYP/TZVP) geometry, the maxima in the computed solution spectra are shifted by several tenths of an eV. Three effects contribute: the use of SCC-DFTB-based rather than B3LYP-based geometries in the MD snapshots (red shift of ca. 0.1 eV), explicit inclusion of nuclear motion through the MD snapshots (red shift of ca. 0.1 eV), and intrinsic solvent effects (differences in the absorption maxima in the computed gas-phase and solution spectra, typically ca. 0.1,0.3 eV). A detailed analysis of the results indicates that the intrinsic solvent effects arise both from solvent-induced structural changes and from electrostatic solute,solvent interactions, the latter being dominant. © 2009 Wiley Periodicals, Inc. J Comput Chem 2010 [source]


    Dynamic structures of phosphodiesterase-5 active site by combined molecular dynamics simulations and hybrid quantum mechanical/molecular mechanical calculations

    JOURNAL OF COMPUTATIONAL CHEMISTRY, Issue 8 2008
    Ying Xiong
    Abstract Various quantum mechanical/molecular mechanical (QM/MM) geometry optimizations starting from an x-ray crystal structure and from the snapshot structures of constrained molecular dynamics (MD) simulations have been performed to characterize two dynamically stable active site structures of phosphodiesterase-5 (PDE5) in solution. The only difference between the two PDE5 structures exists in the catalytic, second bridging ligand (BL2) which is HO, or H2O. It has been shown that, whereas BL2 (i.e. HO,) in the PDE5(BL2 = HO,) structure can really bridge the two positively charged metal ions (Zn2+ and Mg2+), BL2 (i.e. H2O) in the PDE5(BL2 = H2O) structure can only coordinate Mg2+. It has been demonstrated that the results of the QM/MM geometry optimizations are remarkably affected by the solvent water molecules, the dynamics of the protein environment, and the electronic embedding charges of the MM region in the QM part of the QMM/MM calculation. The PDE5(BL2 = H2O) geometries optimized by using the QM/MM method in different ways show strong couplings between these important factors. It is interesting to note that the PDE5(BL2 = HO - ) and PDE5(BL2 = H2O) geometries determined by the QM/MM calculations neglecting these three factors are all consistent with the corresponding geometries determined by the QM/MM calculations that account for all of these three factors. These results suggest the overall effects of these three important factors on the optimized geometries can roughly cancel out. However, the QM/MM calculations that only account for some of these factors could lead to considerably different geometries. These results might be useful also in guiding future QM/MM geometry optimizations on other enzymes. © 2007 Wiley Periodicals, Inc. J Comput Chem, 2008 [source]


    Quantum chemical modeling of the reduction of cis -diammineplatinum(IV) tetrachloride [Pt(NH3)2Cl4] by methyl thiolate anion

    JOURNAL OF COMPUTATIONAL CHEMISTRY, Issue 9 2005
    Ia. I. Dobrogorskaia, Méreau
    Abstract We present here both an ab initio and quantum mechanical/molecular mechanical (QM/MM) study of the cis -[Pt(NH3)2Cl4] complex reduction by methyl thiolate anion, SCH, which is used as a model of glutathione. Geometry and electronic structure of cis -[Pt(NH3)2Cl4] are determined without and in aqueous medium. The mechanism of the reaction of reduction is characterized. The calculated activation energy of the reaction compares remarkably well with the experimental value. © 2005 Wiley Periodicals, Inc. J Comput Chem 26: 865,870, 2005 [source]


    Quantum Mechanical/Molecular Mechanical Studies on Spectral Tuning Mechanisms of Visual Pigments and Other Photoactive Proteins,

    PHOTOCHEMISTRY & PHOTOBIOLOGY, Issue 4 2008
    Ahmet Altun
    The protein environments surrounding the retinal tune electronic absorption maximum from 350 to 630 nm. Hybrid quantum mechanical/molecular mechanical (QM/MM) methods can be used in calculating excitation energies of retinal in its native protein environments and in studying the molecular basis of spectral tuning. We hereby review recent QM/MM results on the phototransduction of bovine rhodopsin, bacteriorhodopsin, sensory rhodopsin II, nonretinal photoactive yellow protein and their mutants. [source]


    The Jahn,Teller Effect of the TiIII Ion in Aqueous Solution: Extended Ab Initio QM/MM Molecular Dynamics Simulations,

    CHEMPHYSCHEM, Issue 10 2004
    Chinapong Kritayakornupong Dr.
    Abstract Combined ab initio quantum mechanical/molecular mechanical (QM/MM) molecular dynamics simulations, including only the first and the first and second hydration shells in the QM region, were performed for TiIIIin aqueous solution. The hydration structure of TiIIIis discussed in terms of radial distribution functions, coordination-number distributions and several angle distributions. Dynamical properties, such as librational and vibrational motions and TiIIIO vibrations, were evaluated. A fast dynamical Jahn,Teller effect of TiIII(aq) was observed in the QM/MM simulations, in particular when the second hydration shell was included into the QM region. The results justify the computational effort required for the inclusion of the second hydration shell into the QM region and show the importance of this effort for obtaining accurate hydration-shell geometries, dynamical properties, and details of the Jahn,Teller effect. [source]


    Modeling dioxygen binding to the non-heme iron-containing enzymes

    INTERNATIONAL JOURNAL OF QUANTUM CHEMISTRY, Issue 10 2006
    A. V. Nemukhin
    Abstract The structures and properties of the complexes formed upon binding the oxygen molecule to the iron sites in non-heme 2-oxoglutarate-dependent enzymes are characterized by QM(CASSCF)/MM and density functional theory (DFT) calculations. Molecular models for the calculations are constructed following the crystal structure of hypoxia-inducible factor asparaginyl hydroxylase (FIH-1). DFT calculations for the 37-atomic cluster have been carried out at the B3LYP(LANL2DZdp) level. The flexible effective fragment potential method is used as a combined quantum mechanical,molecular mechanical (QM/MM) technique to characterize the fragment of the enzymatic system, including 1,758 atoms in the MM part and 27 atoms in the QM part. In these calculations, the CASSCF(LANL2DZdp) approach is applied in the QM subsystem, and AMBER force field parameters are used in the MM subsystem. With both approaches, equilibrium geometry configurations have been located for different spin states of the system. In DFT calculations, the order of the states is as follows: septet, triplet (+7.7 kcal/mol), quintet (+10.7 kcal/mol). Geometry configurations correspond to the end-on structures with no evidences of electron transfer from Fe(II) to molecular oxygen. In contrast, QM(CASSCF)/MM calculations predict the quintet state as the lowest one, while the septet structure has slightly (<2 kcal/mol) higher energy, and the triplet state is considerably more energetic. In QM/MM calculations, in both quintet and septet states, the electronic configurations show considerable electron charge transfer from iron to oxygen, and the oxidation state of iron in the metal binding site can be characterized as Fe(III). © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2006 [source]


    Prediction of pKa shifts in proteins using a combination of molecular mechanical and continuum solvent calculations

    JOURNAL OF COMPUTATIONAL CHEMISTRY, Issue 15 2004
    Bernd Kuhn
    Abstract The prediction of pKa shifts of ionizable groups in proteins is of great relevance for a number of important biological phenomena. We present an implementation of the MM-GBSA approach, which combines molecular mechanical (MM) and generalized Born (GB) continuum solvent energy terms, to the calculation of pKa values of a panel of nine proteins, including 69 individual comparisons with experiment. While applied so far mainly to the calculation of biomolecular binding free energies, we show that this method can also be used for the estimation of protein pKa shifts, with an accuracy around 1 pKa unit, even for strongly shifted residues. Our analysis reveals that the nonelectrostatic terms that are part of the MM-GBSA free energy expression are important contributors to improved prediction accuracy. This suggests that most of the previous approaches that focus only on electrostatic interactions could be improved by adding other nonpolar energy terms to their free energy expression. Interestingly, our method yields best accuracy at protein dielectric constants of ,int = 2,4, which is in contrast to previous approaches that peak at higher ,int , 8. An important component of our procedure is an intermediate minimization step of each protonation state involving different rotamers and tautomers as a way to explicitly model protein relaxation upon (de)protonation. © 2004 Wiley Periodicals, Inc. J Comput Chem 25: 1865,1872, 2004 [source]