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
Fröhlich C
(2024)
Epistasis arises from shifting the rate-limiting step during enzyme evolution of a ß-lactamase
in Nature Catalysis
Çinaroglu S
(2021)
Evaluating the Performance of Water Models with Host-Guest Force Fields in Binding Enthalpy Calculations for Cucurbit[7]uril-Guest Systems
in The Journal of Physical Chemistry B
Bunzel H
(2021)
Evolution of dynamical networks enhances catalysis in a designer enzyme
in Nature Chemistry
Leferink NGH
(2019)
Experiment and Simulation Reveal How Mutations in Functional Plasticity Regions Guide Plant Monoterpene Synthase Product Outcome.
in ACS catalysis
Evans LE
(2019)
Exploitation of Antibiotic Resistance as a Novel Drug Target: Development of a ß-Lactamase-Activated Antibacterial Prodrug.
in Journal of medicinal chemistry
Ashraf S
(2021)
Exploration of the structural requirements of Aurora Kinase B inhibitors by a combined QSAR, modelling and molecular simulation approach.
in Scientific reports
Shannon RJ
(2021)
Exploring human-guided strategies for reaction network exploration: Interactive molecular dynamics in virtual reality as a tool for citizen scientists.
in The Journal of chemical physics
Loeffler HH
(2015)
FESetup: Automating Setup for Alchemical Free Energy Simulations.
in Journal of chemical information and modeling
Oliveira A
(2023)
Fluctuation Relations to Calculate Protein Redox Potentials from Molecular Dynamics Simulations
in Journal of Chemical Theory and Computation
Deeks HM
(2023)
Free energy along drug-protein binding pathways interactively sampled in virtual reality.
in Scientific reports
Toelzer C
(2020)
Free fatty acid binding pocket in the locked structure of SARS-CoV-2 spike protein.
in Science (New York, N.Y.)
Shoemark D
(2021)
Frontispiece: Molecular Simulations suggest Vitamins, Retinoids and Steroids as Ligands of the Free Fatty Acid Pocket of the SARS-CoV-2 Spike Protein
in Angewandte Chemie International Edition
Shoemark DK
(2021)
Frontispiz: Molecular Simulations suggest Vitamins, Retinoids and Steroids as Ligands of the Free Fatty Acid Pocket of the SARS-CoV-2 Spike Protein.
in Angewandte Chemie (Weinheim an der Bergstrasse, Germany)
Trick JL
(2016)
Functional Annotation of Ion Channel Structures by Molecular Simulation.
in Structure (London, England : 1993)
Paramo T
(2017)
Functional Validation of Heteromeric Kainate Receptor Models.
in Biophysical journal
Maingi V
(2015)
Gating-like Motions and Wall Porosity in a DNA Nanopore Scaffold Revealed by Molecular Simulations.
in ACS nano
Lang EJM
(2022)
Generalized Born Implicit Solvent Models Do Not Reproduce Secondary Structures of De Novo Designed Glu/Lys Peptides.
in Journal of chemical theory and computation
Freeman SL
(2023)
Heme binding to the SARS-CoV-2 spike glycoprotein.
in The Journal of biological chemistry
Hermann JC
(2009)
High level QM/MM modeling of the formation of the tetrahedral intermediate in the acylation of wild type and K73A mutant TEM-1 class A beta-lactamase.
in The journal of physical chemistry. A
Claeyssens F
(2006)
High-accuracy computation of reaction barriers in enzymes.
in Angewandte Chemie (International ed. in English)
Kaiyawet N
(2015)
High-level QM/MM calculations support the concerted mechanism for Michael addition and covalent complex formation in thymidylate synthase.
in Journal of chemical theory and computation
Van Der Kamp Marc W.
(2008)
High-level QM/MM modelling predicts an arginine as the acid in the condensation reaction catalysed by citrate synthase
in CHEMICAL COMMUNICATIONS
Van Der Kamp MW
(2008)
High-level QM/MM modelling predicts an arginine as the acid in the condensation reaction catalysed by citrate synthase.
in Chemical communications (Cambridge, England)
Warman H
(2023)
How proton transfer impacts hachimoji DNA.
in RSC advances
Daniels DE
(2023)
Human cellular model systems of ß-thalassemia enable in-depth analysis of disease phenotype.
in Nature communications
Jitonnom Jitrayut
(2011)
Hybrid QM/MM study on the deglycosylation step of chitin hydrolysis catalysed by chitinase B from
Serratia marcescens
in MAEJO INTERNATIONAL JOURNAL OF SCIENCE AND TECHNOLOGY
Sutcliffe MJ
(2006)
Hydrogen tunnelling in enzyme-catalysed H-transfer reactions: flavoprotein and quinoprotein systems.
in Philosophical transactions of the Royal Society of London. Series B, Biological sciences
Corey R
(2021)
Identification and assessment of cardiolipin interactions with E. coli inner membrane proteins
in Science Advances
Pakamwong B
(2022)
Identification of Potent DNA Gyrase Inhibitors Active against Mycobacterium tuberculosis.
in Journal of chemical information and modeling
Lodola A
(2008)
Identification of productive inhibitor binding orientation in fatty acid amide hydrolase (FAAH) by QM/MM mechanistic modelling.
in Chemical communications (Cambridge, England)
Oliveira ASF
(2019)
Identification of the Initial Steps in Signal Transduction in the a4ß2 Nicotinic Receptor: Insights from Equilibrium and Nonequilibrium Simulations.
in Structure (London, England : 1993)
Ge Y
(2017)
Identification of the quinolinedione inhibitor binding site in Cdc25 phosphatase B through docking and molecular dynamics simulations.
in Journal of computer-aided molecular design
Bueren-Calabuig JA
(2016)
Impact of Ser17 Phosphorylation on the Conformational Dynamics of the Oncoprotein MDM2.
in Biochemistry
Martí S
(2022)
Impact of Warhead Modulations on the Covalent Inhibition of SARS-CoV-2 Mpro Explored by QM/MM Simulations.
in ACS catalysis
Gerrard W
(2020)
IMPRESSION - prediction of NMR parameters for 3-dimensional chemical structures using machine learning with near quantum chemical accuracy.
in Chemical science
Gray A
(2015)
In pursuit of an accurate spatial and temporal model of biomolecules at the atomistic level: a perspective on computer simulation.
in Acta crystallographica. Section D, Biological crystallography
Hanwarinroj C
(2022)
In silico design of novel quinazoline-based compounds as potential Mycobacterium tuberculosis PknB inhibitors through 2D and 3D-QSAR, molecular dynamics simulations combined with pharmacokinetic predictions.
in Journal of molecular graphics & modelling
Kamsri P
(2019)
In silico study directed towards identification of the key structural features of GyrB inhibitors targeting MTB DNA gyrase: HQSAR, CoMSIA and molecular dynamics simulations.
in SAR and QSAR in environmental research
Von Kügelgen A
(2020)
In Situ Structure of an Intact Lipopolysaccharide-Bound Bacterial Surface Layer.
in Cell
Buzas D
(2023)
In vitro generated antibodies guide thermostable ADDomer nanoparticle design for nasal vaccination and passive immunization against SARS-CoV-2
in Antibody Therapeutics
Lonsdale R
(2010)
Inclusion of Dispersion Effects Significantly Improves Accuracy of Calculated Reaction Barriers for Cytochrome P450 Catalyzed Reactions
in The Journal of Physical Chemistry Letters
Songsiriritthigul C
(2022)
Inhibition of Mycobacterium tuberculosis InhA by 3-nitropropanoic acid.
in Proteins
Jitonnom J
(2012)
Insights into conformational changes of procarboxypeptidase A and B from simulations: a plausible explanation for different intrinsic activity
in Theoretical Chemistry Accounts
RANAGHAN K
(2009)
Insights into enzyme catalysis from QM/MM modelling: transition state stabilization in chorismate mutase
in Molecular Physics
Sharma V
(2017)
Insights into functions of the H channel of cytochrome c oxidase from atomistic molecular dynamics simulations.
in Proceedings of the National Academy of Sciences of the United States of America
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 |