Combined quantum mechanics/molecular mechanics (QM/MM) Monte Carlo free energy simulations: a feasibility study
Lead Research Organisation:
University of Bristol
Department Name: Chemistry
Abstract
Despite the advances of science, millions of people still die every year from incurable diseases. Unfortunately, the costs of drug development are so high that the focus of medicinal research is into profitable Western diseases. To reduce the costs of developing new medicinal drugs, we would like to be able to use computers to model how a potential drug works within the body, and to use this knowledge to design new and better drugs. Building computational models like this is challenging, requiring a delicate balance between putting enough detail into the model to get realistic behaviour, and making the model as simple as possible so that it doesn't take too long to run the calculations. Until now, the majority of models used have been very simple, modelling the atoms of a drug as balls on springs. By treating the atoms as solid balls, the models neglect the atom's most chemically important part, namely the electrons. This is a severe oversight, as it is the interactions of electrons that determine whether the drug could dissolve in your blood, work its way into your cells, and bind to, and thus neutralize, the proteins of any attacking bacteria or virus. It is possible to model electrons in molecules using quantum mechanics. However, to model the entire protein/drug system using quantum mechanics would be too computationally expensive. We propose to research the use of quantum mechanics to model just the electrons that are part of, and near to, the drug molecule. The rest of the protein can still be treated by simple ball and springs models to make the calculations possible. The new methods we will develop add important extra detail, making them more realistic and better able to model how drugs interact. At the same time, this combined approach should mean that the calculations are practical to do. What makes our planned work different is that it will involve the development of a mixed model specifically tailored for medicinal drug design. Creating a mixed model for this use will require that significant challenges are overcome, and that new ways are developed to handle the interactions between the quantum mechanics part of the model with the ball on springs part.
Organisations
Publications
Van Der Kamp MW
(2008)
Biomolecular simulation and modelling: status, progress and prospects.
in Journal of the Royal Society, Interface
Amaro RE
(2020)
Biomolecular Simulations in the Time of COVID19, and After.
in Computing in science & engineering
Huggins D
(2018)
Biomolecular simulations: From dynamics and mechanisms to computational assays of biological activity
in WIREs Computational Molecular Science
Lester Hedges
(2019)
BioSimSpace: An interoperable Python framework for biomolecular simulation
in Zenodo
Hedges L
(2019)
BioSimSpace: An interoperable Python framework for biomolecular simulation
in Journal of Open Source Software
Lester Hedges
(2019)
BioSimSpace: An interoperable Python framework for biomolecular simulation
in Zenodo
Chudyk Ewa Iwona
(2013)
Calculating free energy profiles for enzyme catalysed reactions
Mulholland A
(1998)
Calculations on the substrates of citrate synthase I. Oxaloacetate
in Journal of Molecular Structure: THEOCHEM
Douglas-Gallardo O
(2022)
Carbon Dioxide Fixation in RuBisCO Is Protonation-State-Dependent and Irreversible
in ACS Catalysis
Mulholland A
(2020)
Catalytic Mechanism of the Colistin Resistance Protein MCR-1
Suardiaz R
(2020)
Catalytic Mechanism of the Colistin Resistance Protein MCR-1
Suardíaz R
(2021)
Catalytic mechanism of the colistin resistance protein MCR-1.
in Organic & biomolecular chemistry
Pucheta-Martinez E
(2016)
Changes in the folding landscape of the WW domain provide a molecular mechanism for an inherited genetic syndrome.
in Scientific reports
Rao S
(2020)
Characterizing Membrane Association and Periplasmic Transfer of Bacterial Lipoproteins through Molecular Dynamics Simulations.
in Structure (London, England : 1993)
Mulholland AJ
(2007)
Chemical accuracy in QM/MM calculations on enzyme-catalysed reactions.
in Chemistry Central journal
Woods C
(2008)
Chemical Modelling - Applications and Theory
Lonsdale R
(2012)
ChemInform Abstract: A Practical Guide to Modelling Enzyme-Catalyzed Reactions
in ChemInform
Lonsdale R
(2010)
ChemInform Abstract: Computational Enzymology
in ChemInform
Van Der Kamp M
(2009)
ChemInform Abstract: Computational Enzymology: Insight into Biological Catalyst from Modelling
in ChemInform
Hedger G
(2019)
Cholesterol Interaction Sites on the Transmembrane Domain of the Hedgehog Signal Transducer and Class F G Protein-Coupled Receptor Smoothened.
in Structure (London, England : 1993)
Elcock A
(2002)
Combined Quantum and Molecular Mechanical Study of DNA Crosslinking by Nitrous Acid
in Journal of the American Chemical Society
Ridder L
(1999)
Combined quantum mechanical and molecular mechanical reaction pathway calculation for aromatic hydroxylation by p-hydroxybenzoate-3-hydroxylase.
in Journal of molecular graphics & modelling
Ainsley J
(2018)
Combined Quantum Mechanics and Molecular Mechanics Studies of Enzymatic Reaction Mechanisms.
in Advances in protein chemistry and structural biology
Van Der Kamp M
(2013)
Combined Quantum Mechanics/Molecular Mechanics (QM/MM) Methods in Computational Enzymology
in Biochemistry
Woods CJ
(2015)
Combined quantum mechanics/molecular mechanics (QM/MM) simulations for protein-ligand complexes: free energies of binding of water molecules in influenza neuraminidase.
in The journal of physical chemistry. B
Lonsdale R
(2011)
Comment on "A stationary-wave model of enzyme catalysis" by Carlo Canepa.
in Journal of computational chemistry
Harvey JN
(2006)
Comment on "Molecular dynamics DFT:B3LYP study of guanosinetriphosphate conversion into guanosinemonophosphate upon Mg2+ chelation of alpha and beta phosphate oxygens of the triphosphate tail" by Alexander A. Tulub, Phys. Chem. Chem. Phys., 2006, 8, 2187.
in Physical chemistry chemical physics : PCCP
Lear A
(2023)
Comment on: "Computer Simulations Reveal an Entirely Entropic Activation Barrier for the Chemical Step in a Designer Enzyme"
in ACS Catalysis
Mlýnský V
(2014)
Comparison of ab Initio, DFT, and Semiempirical QM/MM Approaches for Description of Catalytic Mechanism of Hairpin Ribozyme.
in Journal of chemical theory and computation
Lawan N
(2014)
Comparison of DFT and ab initio QM/MM methods for modelling reaction in chorismate synthase
in Chemical Physics Letters
Rodríguez A
(2007)
Comparison of different quantum mechanical/molecular mechanics boundary treatments in the reaction of the hepatitis C virus NS3 protease with the NS5A/5B substrate.
in The journal of physical chemistry. B
Shaw K
(2009)
Compatibility of Quantum Chemical Methods and Empirical (MM) Water Models in Quantum Mechanics/Molecular Mechanics Liquid Water Simulations
in The Journal of Physical Chemistry Letters
Lonsdale R
(2010)
Compound I reactivity defines alkene oxidation selectivity in cytochrome P450cam.
in The journal of physical chemistry. B
Smith LJ
(2008)
Computational and experimental studies on the catalytic mechanism of biliverdin-IXbeta reductase.
in The Biochemical journal
Woods CJ
(2013)
Computational assay of H7N9 influenza neuraminidase reveals R292K mutation reduces drug binding affinity.
in Scientific reports
Jenkins Jonathan
(2019)
Computational design, construction, and characterisation of artificial peroxidases
Lonsdale R
(2010)
Computational enzymology.
in Chemical communications (Cambridge, England)
Lodola A
(2013)
Computational enzymology.
in Methods in molecular biology (Clifton, N.J.)
Van Der Kamp MW
(2008)
Computational enzymology: insight into biological catalysts from modelling.
in Natural product reports
Mulholland AJ
(2008)
Computational enzymology: modelling the mechanisms of biological catalysts.
in Biochemical Society transactions
Reddy T
(2016)
Computational virology: From the inside out
in Biochimica et Biophysica Acta (BBA) - Biomembranes
Çinaroglu SS
(2023)
Computed Protein-Protein Enthalpy Signatures as a Tool for Identifying Conformation Sampling Problems.
in Journal of chemical information and modeling
Ranaghan KE
(2010)
Computer simulations of quantum tunnelling in enzyme-catalysed hydrogen transfer reactions.
in Interdisciplinary sciences, computational life sciences
Grant Ian M.
(2009)
Conformation and Catalysis in Lysozyme. A computational study
Mcgeagh John David
(2011)
Conformation and cooperativity in homodimeric enzymes investigated by molecular dynamics simulations
Van Der Kamp MW
(2013)
Conformational change and ligand binding in the aristolochene synthase catalytic cycle.
in Biochemistry
Caldwell TA
(2022)
Conformational dynamics of the membrane enzyme LspA upon antibiotic and substrate binding.
in Biophysical journal
Lodola A
(2007)
Conformational effects in enzyme catalysis: reaction via a high energy conformation in fatty acid amide hydrolase.
in Biophysical journal
Christov C
(2013)
Conformational Effects on the pro - S Hydrogen Abstraction Reaction in Cyclooxygenase-1: An Integrated QM/MM and MD Study
in Biophysical Journal
Description | EPSRC |
Amount | £188,950 (GBP) |
Funding ID | E/EP/G007705/1 |
Organisation | Engineering and Physical Sciences Research Council (EPSRC) |
Sector | Public |
Country | United Kingdom |
Start | 10/2013 |
End | 03/2014 |
Title | Sire 2007.1 |
Description | 2007.1 (first official) release of the Sire molecular simulation framework. This included new methods developed to calculate QM/MM free energies. |
Type Of Technology | Software |
Year Produced | 2007 |
Open Source License? | Yes |
Impact | Sire is used in several pharmaceutical companies. This version of the code was used to run the simulations in "An efficient method for the calculation of quantum mechanics/molecular mechanics free energies" Christopher J. Woods, Frederick R. Manby and Adrian J. Mulholland J. Chem. Phys. 128 014109 (2008) doi:10.1063/1.2805379 The combination of quantum mechanics (QM) with molecular mechanics (MM) offers a route to improved accuracy in the study of biological systems, and there is now significant research effort being spent to develop QM/MM methods that can be applied to the calculation of relative free energies. Currently, the computational expense of the QM part of the calculation means that there is no single method that achieves both efficiency and rigor; either the QM/MM free energy method is rigorous and computationally expensive, or the method introduces efficiency-led assumptions that can lead to errors in the result, or a lack of generality of application. In this paper we demonstrate a combined approach to form a single, efficient, and, in principle, exact QM/MM free energy method. We demonstrate the application of this method by using it to explore the difference in hydration of water and methane. We demonstrate that it is possible to calculate highly converged QM/MM relative free energies at the MP2/aug-cc-pVDZ/OPLS level within just two days of computation, using commodity processors, and show how the method allows consistent, high-quality sampling of complex solvent configurational change, both when perturbing hydrophilic water into hydrophobic methane, and also when moving from a MM Hamiltonian to a QM/MM Hamiltonian. The results demonstrate the validity and power of this methodology, and raise important questions regarding the compatibility of MM and QM/MM forcefields, and offer a potential route to improved compatibility. |
URL | http://www.siremol.org/Sire/Home.html |