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
Hutchings G
(2023)
Modern Developments in Catalysis - Volume 2
Song W
(2021)
Modulation of adenosine A2a receptor oligomerization by receptor activation and PIP2 interactions.
in Structure (London, England : 1993)
Tooke CL
(2019)
Molecular Basis of Class A ß-Lactamase Inhibition by Relebactam.
in Antimicrobial agents and chemotherapy
Kiani YS
(2019)
Molecular Dynamics Simulation Framework to Probe the Binding Hypothesis of CYP3A4 Inhibitors.
in International journal of molecular sciences
Kattnig D
(2018)
Molecular dynamics simulations disclose early stages of the photo-activation of cryptochrome 4
in New Journal of Physics
Naafs B
(2021)
Molecular dynamics simulations support the hypothesis that the brGDGT paleothermometer is based on homeoviscous adaptation
in Geochimica et Cosmochimica Acta
Mulholland A
(2007)
Molecular Materials with Specific Interactions - Modeling and Design
Hermann JC
(2006)
Molecular mechanisms of antibiotic resistance: QM/MM modelling of deacylation in a class A beta-lactamase.
in Organic & biomolecular chemistry
Pradon Juliette Jeanne Marie
(2010)
Molecular modelling studies of the relationship between structure and mechanism in (beta)-lactamases and related enzymes
Shoemark D
(2021)
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
Shoemark DK
(2021)
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 ed. in English)
Woods CJ
(2009)
Multicore Parallelization of Kohn-Sham Theory.
in Journal of chemical theory and computation
Yamamoto E
(2020)
Multiple lipid binding sites determine the affinity of PH domains for phosphoinositide-containing membranes.
in Science advances
Amaro RE
(2018)
Multiscale Methods in Drug Design Bridge Chemical and Biological Complexity in the Search for Cures.
in Nature reviews. Chemistry
Jagger BR
(2020)
Multiscale simulation approaches to modeling drug-protein binding.
in Current opinion in structural biology
Hirvonen VHA
(2022)
Multiscale Simulations Identify Origins of Differential Carbapenem Hydrolysis by the OXA-48 ß-Lactamase.
in ACS catalysis
Yang Z
(2021)
Multiscale Workflow for Modeling Ligand Complexes of Zinc Metalloproteins.
in Journal of chemical information and modeling
Tooke C
(2021)
Natural variants modify Klebsiella pneumoniae carbapenemase (KPC) acyl-enzyme conformational dynamics to extend antibiotic resistance
in Journal of Biological Chemistry
Angulo G
(2016)
New methods: general discussion.
in Faraday discussions
Oliveira A. S. F.
(2019)
Nicotine-induced conformational changes in the a4ß2 nicotinic receptor
in EUROPEAN BIOPHYSICS JOURNAL WITH BIOPHYSICS LETTERS
Althorpe SC
(2016)
Non-adiabatic reactions: general discussion.
in Faraday discussions
Chudyk E
(2013)
Nonempirical Energetic Analysis of Reactivity and Covalent Inhibition of Fatty Acid Amide Hydrolase
in The Journal of Physical Chemistry B
Minguez Teresa
(2020)
Novel determinants of agonist selectivity in nicotinic ACh receptors
in BRITISH JOURNAL OF PHARMACOLOGY
Jefferys E
(2014)
NRas slows the rate at which a model lipid bilayer phase separates.
in Faraday discussions
Samsudin F
(2016)
OmpA: A Flexible Clamp for Bacterial Cell Wall Attachment.
in Structure (London, England : 1993)
Sampson C
(2019)
On the magnetosensitivity of lipid peroxidation: two- versus three-radical dynamics
in Physical Chemistry Chemical Physics
Arcus VL
(2016)
On the Temperature Dependence of Enzyme-Catalyzed Rates.
in Biochemistry
Ren Q
(2010)
Optimal control design of laser pulses for mode specific vibrational excitation in an enzyme-substrate complex
in Chemical Physics Letters
Heiblum Robles A
(2018)
Phase transitions in stigmergic territorial systems
in Physical Review E
Jambrina P
(2015)
Phosphorylation of RAF Kinase Dimers Drives Conformational Changes that Facilitate Transactivation
in Angewandte Chemie
Bunzel H
(2022)
Photovoltaic enzymes by design and evolution
Mulholland Adrian J.
(2009)
PHYS 47-Biomolecular simulations of enzymatic reactions
in ABSTRACTS OF PAPERS OF THE AMERICAN CHEMICAL SOCIETY
Shahane G
(2019)
Physical properties of model biological lipid bilayers: insights from all-atom molecular dynamics simulations.
in Journal of molecular modeling
Matheson A
(2018)
Phytosterol-based edible oleogels: A novel way of replacing saturated fat in food.
in Nutrition bulletin
Ranaghan KE
(2019)
Projector-Based Embedding Eliminates Density Functional Dependence for QM/MM Calculations of Reactions in Enzymes and Solution.
in Journal of chemical information and modeling
McGeagh JD
(2011)
Protein dynamics and enzyme catalysis: insights from simulations.
in Biochimica et biophysica acta
Glowacki DR
(2012)
Protein dynamics and enzyme catalysis: the ghost in the machine?
in Biochemical Society transactions
Wu Z
(2019)
Proton Control of Transitions in an Amino Acid Transporter.
in Biophysical journal
Slocombe L
(2022)
Proton transfer during DNA strand separation as a source of mutagenic guanine-cytosine tautomers
in Communications Chemistry
Song W
(2022)
PyLipID: A Python Package for Analysis of Protein-Lipid Interactions from Molecular Dynamics Simulations.
in Journal of chemical theory and computation
Hanpaibool C
(2023)
Pyrazolones Potentiate Colistin Activity against MCR-1-Producing Resistant Bacteria: Computational and Microbiological Study
in ACS Omega
Harvey J
(2006)
QM and QM/MM studies of selectivity in organic and bioorganic chemistry
in Journal of Physical Organic Chemistry
Jitonnom J
(2014)
QM/MM free-energy simulations of reaction in Serratia marcescens Chitinase B reveal the protonation state of Asp142 and the critical role of Tyr214.
in The journal of physical chemistry. B
Bathelt CM
(2008)
QM/MM modeling of benzene hydroxylation in human cytochrome P450 2C9.
in The journal of physical chemistry. A
Harvey JN
(2006)
QM/MM modeling of compound I active species in cytochrome P450, cytochrome C peroxidase, and ascorbate peroxidase.
in Journal of computational chemistry
Karabencheva-Christova Tatyana Georgieva
(2009)
QM/MM modelling of mechanisms of flavin-containing enzymes
Lawan Narin
(2010)
QM/MM modelling of the reaction mechanism of chorismate synthase and mutase
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 | 09/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 |