1
|
Ibrahim MAA. Performance Assessment of Semiempirical Molecular Orbital Methods in Describing Halogen Bonding: Quantum Mechanical and Quantum Mechanical/Molecular Mechanical-Molecular Dynamics Study. J Chem Inf Model 2011; 51:2549-59. [DOI: 10.1021/ci2002582] [Citation(s) in RCA: 32] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Mahmoud A. A. Ibrahim
- School of Chemistry, University of Manchester, Oxford Road, Manchester M139PL, United Kingdom
- Chemistry Department, Faculty of Science, Minia University, Minia 61519, Egypt
| |
Collapse
|
2
|
Cummins PL, Gready JE. The Influence of Starting Coordinates in Free Energy Simulations of Ligand Binding to Dihydrofolate Reductase. MOLECULAR SIMULATION 2006. [DOI: 10.1080/08927029508024052] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/23/2022]
|
3
|
Monajjemi M, Heshmat M, Haeri HH, Kaveh F. Theoretical study of vitamin properties from combined QM-MM methods: Comparison of chemical shifts and energy. RUSSIAN JOURNAL OF PHYSICAL CHEMISTRY A 2006. [DOI: 10.1134/s0036024406070119] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/07/2022]
|
4
|
Sagarik K, Chaiyapongs S. Structures and stability of salt-bridge in aqueous solution. Biophys Chem 2005; 117:119-40. [PMID: 15935545 DOI: 10.1016/j.bpc.2005.04.012] [Citation(s) in RCA: 14] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/16/2005] [Revised: 04/22/2005] [Accepted: 04/21/2005] [Indexed: 11/29/2022]
Abstract
Structures and stability of salt-bridges in aqueous solutions were investigated using a complex formed from the guanidinium (Gdm+) and formate (FmO-) ions as a model system. The Test-particle model (T-model) potentials to describe the interactions in the Gdm+-H2O, FmO(-)-H2O and Gdm+-FmO- complexes were constructed, tested and applied in molecular dynamics (MD) simulations of the aqueous solutions at 298 K. The three-dimensional structures and energetic of the hydrogen bond (H-bond) networks of water in the first hydration shells of the Gdm+ and FmO- ions, as well as the Gdm+-FmO- complex, were visualized and analyzed using various probability distribution (PD) maps. The structures of the average potential energy landscapes at the H-bond networks were employed to characterize the stability and dynamic behavior of water molecules in the first hydration shells of the solutes. It was observed that water molecules in the first hydration shell of the close-contact Gdm+-FmO- complex form associated H-bond networks, which introduce a net stabilization effect to the ion-pair, whereas those in the interstitial H-bond network destabilize and break the solvent-separated Gdm+-FmO- complex. The present results showed that, in order to provide complete insights into the structures and stability of ion-pairs in aqueous solutions, explicit water molecules have to be included in the model calculations.
Collapse
Affiliation(s)
- Kritsana Sagarik
- School of Chemistry, Institute of Science, Suranaree University of Technology, Nakhon Ratchasima 30000, Thailand.
| | | |
Collapse
|
5
|
Gustiananda M, Liggins JR, Cummins PL, Gready JE. Conformation of prion protein repeat peptides probed by FRET measurements and molecular dynamics simulations. Biophys J 2004; 86:2467-83. [PMID: 15041684 PMCID: PMC1304095 DOI: 10.1016/s0006-3495(04)74303-9] [Citation(s) in RCA: 54] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022] Open
Abstract
We report the combined use of steady-state fluorescence resonance energy transfer (FRET) experiments and molecular dynamics (MD) simulations to investigate conformational distributions of the prion protein (PrP) repeat system. FRET was used for the first time to probe the distance, as a function of temperature and pH, between a donor Trp residue and an acceptor dansyl group attached to the N-terminus in seven model peptides containing one to three repeats of the second decarepeat of PrP from marsupial possum (PHPGGSNWGQ)nG, and one and two human PrP consensus octarepeats (PHGGGWGQ)nG. In multirepeat peptides, single-Trp mutants were made by replacing other Trp(s) with Phe. As previous work has shown PrP repeats do not adopt a single preferred stable conformation, the FRET values are averages reflecting heterogeneity in the donor-acceptor distances. The T-dependence of the conformational distributions, and derived average dansyl-Trp distances, were obtained directly from MD simulation of the marsupial dansyl-PHPGGSNWGQG peptide. The results show excellent agreement between the FRET and MD T-dependent distances, and demonstrate the remarkable sensitivity and reproducibility of the FRET method in this first-time use for a set of disordered peptides. Based on the results, we propose a model involving cation-pi or pi-pi His-Trp interactions to explain the T- (5-85 degrees C) and pH- (6.0, 7.2) dependencies on distance, with HW i, i + 4 or WH i, i + 4 separations in sequence being more stable than HW i, i + 6 or WH i, i + 6 separations. The model has peptides adopting loosely folded conformations, with dansyl-Trp distances very much less than estimates for fully extended conformations, for example, approximately 16 vs. 33, approximately 21 vs. 69, and approximately 22 vs. 106 A for 1-3 decarepeats, and approximately 14 vs. 25 and approximately 19 vs. 54 A for 1-2 octarepeats, respectively. The study demonstrates the usefulness of combining FRET with MD, a combination reported only once previously. Initial "mapping" of the conformational distribution of flexible peptides by simulation can assist in designing and interpreting experiments using steady-state intensity methods, and indicating how time-resolved or anisotropy methods might be used.
Collapse
Affiliation(s)
- Marsia Gustiananda
- Computational Proteomics Group, John Curtin School of Medical Research, Australian National University, Canberra ACT 2601, Australia
| | | | | | | |
Collapse
|
6
|
Sherwood P, de Vries AH, Guest MF, Schreckenbach G, Catlow CA, French SA, Sokol AA, Bromley ST, Thiel W, Turner AJ, Billeter S, Terstegen F, Thiel S, Kendrick J, Rogers SC, Casci J, Watson M, King F, Karlsen E, Sjøvoll M, Fahmi A, Schäfer A, Lennartz C. QUASI: A general purpose implementation of the QM/MM approach and its application to problems in catalysis. ACTA ACUST UNITED AC 2003. [DOI: 10.1016/s0166-1280(03)00285-9] [Citation(s) in RCA: 642] [Impact Index Per Article: 30.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/27/2022]
|
7
|
Titmuss SJ, Cummins PL, Rendell AP, Bliznyuk AA, Gready JE. Comparison of linear-scaling semiempirical methods and combined quantum mechanical/molecular mechanical methods for enzymic reactions. II. An energy decomposition analysis. J Comput Chem 2002; 23:1314-22. [PMID: 12214314 DOI: 10.1002/jcc.10122] [Citation(s) in RCA: 25] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
Abstract
QM/MM methods have been developed as a computationally feasible solution to QM simulation of chemical processes, such as enzyme-catalyzed reactions, within a more approximate MM representation of the condensed-phase environment. However, there has been no independent method for checking the quality of this representation, especially for highly nonisotropic protein environments such as those surrounding enzyme active sites. Hence, the validity of QM/MM methods is largely untested. Here we use the possibility of performing all-QM calculations at the semiempirical PM3 level with a linear-scaling method (MOZYME) to assess the performance of a QM/MM method (PM3/AMBER94 force field). Using two model pathways for the hydride-ion transfer reaction of the enzyme dihydrofolate reductase studied previously (Titmuss et al., Chem Phys Lett 2000, 320, 169-176), we have analyzed the reaction energy contributions (QM, QM/MM, and MM) from the QM/MM results and compared them with analogous-region components calculated via an energy partitioning scheme implemented into MOZYME. This analysis further divided the MOZYME components into Coulomb, resonance and exchange energy terms. For the model in which the MM coordinates are kept fixed during the reaction, we find that the MOZYME and QM/MM total energy profiles agree very well, but that there are significant differences in the energy components. Most significantly there is a large change (approximately 16 kcal/mol) in the MOZYME MM component due to polarization of the MM region surrounding the active site, and which arises mostly from MM atoms close to (<10 A) the active-site QM region, which is not modelled explicitly by our QM/MM method. However, for the model where the MM coordinates are allowed to vary during the reaction, we find large differences in the MOZYME and QM/MM total energy profiles, with a discrepancy of 52 kcal/mol between the relative reaction (product-reactant) energies. This is largely due to a difference in the MM energies of 58 kcal/mol, of which we can attribute approximately 40 kcal/mol to geometry effects in the MM region and the remainder, as before, to MM region polarization. Contrary to the fixed-geometry model, there is no correlation of the MM energy changes with distance from the QM region, nor are they contributed by only a few residues. Overall, the results suggest that merely extending the size of the QM region in the QM/MM calculation is not a universal solution to the MOZYME- and QM/MM-method differences. They also suggest that attaching physical significance to MOZYME Coulomb, resonance and exchange components is problematic. Although we conclude that it would be possible to reparameterize the QM/MM force field to reproduce MOZYME energies, a better way to account for both the effects of the protein environment and known deficiencies in semiempirical methods would be to parameterize the force field based on data from DFT or ab initio QM linear-scaling calculations. Such a force field could be used efficiently in MD simulations to calculate free energies.
Collapse
Affiliation(s)
- Stephen J Titmuss
- Computational Proteomics and Therapy Design Group, Division of Biochemistry and Molecular Biology, John Curtin School of Medical Research, Australian National University, P.O. Box 334, Canberra, ACT 2601, Australia
| | | | | | | | | |
Collapse
|
8
|
Cummins PL, Titmuss SJ, Jayatilaka D, Bliznyuk AA, Rendell AP, Gready JE. Comparison of semiempirical and ab initio QM decomposition analyses for the interaction energy between molecules. Chem Phys Lett 2002. [DOI: 10.1016/s0009-2614(01)01417-8] [Citation(s) in RCA: 14] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/27/2022]
|
9
|
Orozco M, Luque FJ. Theoretical Methods for the Description of the Solvent Effect in Biomolecular Systems. Chem Rev 2000; 100:4187-4226. [PMID: 11749344 DOI: 10.1021/cr990052a] [Citation(s) in RCA: 454] [Impact Index Per Article: 18.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/13/2023]
Affiliation(s)
- Modesto Orozco
- Departament de Bioquímica i Biologia Molecular, Facultat de Química, Universitat de Barcelona, Martí i Franqués 1, E-08028 Barcelona, Spain, and Departament de Fisicoquímica, Facultat de Farmàcia, Universitat de Barcelona, Avgda. Diagonal s/n, E-08028 Barcelona, Spain
| | | |
Collapse
|
10
|
Luque FJ, Reuter N, Cartier A, Ruiz-López MF. Calibration of the Quantum/Classical Hamiltonian in Semiempirical QM/MM AM1 and PM3 Methods. J Phys Chem A 2000. [DOI: 10.1021/jp001974g] [Citation(s) in RCA: 48] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- F. J. Luque
- Departament de Fisicoquímica, Facultat de Farmàcia, Universitat de Barcelona, Av. Diagonal s/n, 08028 Barcelona, Spain
| | - N. Reuter
- Unité de Recherche CNRS-UHP 7565, Laboratoire de Chimie Théorique, Université Henri Poincaré-Nancy I, BP 239, 54506 Vandoeuvre-lès-Nancy, France
| | - A. Cartier
- Unité de Recherche CNRS-UHP 7565, Laboratoire de Chimie Théorique, Université Henri Poincaré-Nancy I, BP 239, 54506 Vandoeuvre-lès-Nancy, France
| | - M. F. Ruiz-López
- Unité de Recherche CNRS-UHP 7565, Laboratoire de Chimie Théorique, Université Henri Poincaré-Nancy I, BP 239, 54506 Vandoeuvre-lès-Nancy, France
| |
Collapse
|
11
|
Svanberg M, Pettersson JBC, Bolton K. Coupled QM/MM Molecular Dynamics Simulations of HCl Interacting with Ice Surfaces and Water Clusters − Evidence of Rapid Ionization. J Phys Chem A 2000. [DOI: 10.1021/jp0012698] [Citation(s) in RCA: 96] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Marcus Svanberg
- Department of Chemistry, Physical Chemistry, Göteborg University, SE-412 96 Göteborg, Sweden, School of Engineering, University of Borås, SE-501 90 Borås, Sweden, and School of Environmental Sciences, Göteborg University, SE-412 96 Göteborg, Sweden
| | - Jan B. C. Pettersson
- Department of Chemistry, Physical Chemistry, Göteborg University, SE-412 96 Göteborg, Sweden, School of Engineering, University of Borås, SE-501 90 Borås, Sweden, and School of Environmental Sciences, Göteborg University, SE-412 96 Göteborg, Sweden
| | - Kim Bolton
- Department of Chemistry, Physical Chemistry, Göteborg University, SE-412 96 Göteborg, Sweden, School of Engineering, University of Borås, SE-501 90 Borås, Sweden, and School of Environmental Sciences, Göteborg University, SE-412 96 Göteborg, Sweden
| |
Collapse
|
12
|
Titmuss SJ, Cummins PL, Bliznyuk AA, Rendell AP, Gready JE. Comparison of linear-scaling semiempirical methods and combined quantum mechanical/molecular mechanical methods applied to enzyme reactions. Chem Phys Lett 2000. [DOI: 10.1016/s0009-2614(00)00215-3] [Citation(s) in RCA: 30] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/27/2022]
|
13
|
Cummins PL, Gready JE. Coupled semiempirical quantum mechanics and molecular mechanics (QM/MM) calculations on the aqueous solvation free energies of ionized molecules. J Comput Chem 1999. [DOI: 10.1002/(sici)1096-987x(19990730)20:10<1028::aid-jcc5>3.0.co;2-7] [Citation(s) in RCA: 23] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022]
|
14
|
Cummins PL, Gready JE. Coupled semiempirical molecular orbital and molecular mechanics model (QM/MM) for organic molecules in aqueous solution. J Comput Chem 1997. [DOI: 10.1002/(sici)1096-987x(199709)18:12<1496::aid-jcc7>3.0.co;2-e] [Citation(s) in RCA: 62] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/17/2023]
|
15
|
Ranganathan S, Gready JE. Hybrid Quantum and Molecular Mechanical (QM/MM) Studies on the Pyruvate to l-Lactate Interconversion in l-Lactate Dehydrogenase. J Phys Chem B 1997. [DOI: 10.1021/jp971051u] [Citation(s) in RCA: 54] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/27/2023]
Affiliation(s)
- Shoba Ranganathan
- Division of Biochemistry and Molecular Biology, John Curtin School of Medical Research, Australian National University, PO Box 334, Canberra ACT 2601, Australia
| | - Jill E. Gready
- Division of Biochemistry and Molecular Biology, John Curtin School of Medical Research, Australian National University, PO Box 334, Canberra ACT 2601, Australia
| |
Collapse
|
16
|
|
17
|
|
18
|
|