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
Rowlands Heather Ann
(2009)
Modelling of the sugar processing enzymes chitinase B and Phosphomannomutase/ Phosphoglucomutase
Ross G
(2015)
Water Sites, Networks, And Free Energies with Grand Canonical Monte Carlo
in Journal of the American Chemical Society
Román-Meléndez GD
(2014)
Role of active site residues in promoting cobalt-carbon bond homolysis in adenosylcobalamin-dependent mutases revealed through experiment and computation.
in Biochemistry
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
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
Rhys GG
(2018)
Maintaining and breaking symmetry in homomeric coiled-coil assemblies.
in Nature communications
Ren Q
(2010)
Optimal control design of laser pulses for mode specific vibrational excitation in an enzyme-substrate complex
in Chemical Physics Letters
Rego Campello H
(2018)
Unlocking Nicotinic Selectivity via Direct C?H Functionalization of (-)-Cytisine
in Chem
Reddy T
(2016)
Computational virology: From the inside out
in Biochimica et Biophysica Acta (BBA) - Biomembranes
Raza S
(2019)
Visualizing protein-ligand binding with chemical energy-wise decomposition (CHEWD): application to ligand binding in the kallikrein-8 S1 Site.
in Journal of computer-aided molecular design
Rao S
(2020)
Characterizing Membrane Association and Periplasmic Transfer of Bacterial Lipoproteins through Molecular Dynamics Simulations.
in Structure (London, England : 1993)
Ranaghan KE
(2010)
Computer simulations of quantum tunnelling in enzyme-catalysed hydrogen transfer reactions.
in Interdisciplinary sciences, computational life sciences
Ranaghan KE
(2007)
Analysis of classical and quantum paths for deprotonation of methylamine by methylamine dehydrogenase.
in Chemphyschem : a European journal of chemical physics and physical chemistry
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
Ranaghan KE
(2017)
Ab Initio QM/MM Modeling of the Rate-Limiting Proton Transfer Step in the Deamination of Tryptamine by Aromatic Amine Dehydrogenase.
in The journal of physical chemistry. B
Ranaghan K
(2016)
Simulating Enzyme Reactivity - Computational Methods in Enzyme Catalysis
Ranaghan K
(2014)
A catalytic role for methionine revealed by a combination of computation and experiments on phosphite dehydrogenase
in Chem. Sci.
RANAGHAN K
(2009)
Insights into enzyme catalysis from QM/MM modelling: transition state stabilization in chorismate mutase
in Molecular Physics
Ranaghan K
(2010)
Investigations of enzyme-catalysed reactions with combined quantum mechanics/molecular mechanics (QM/MM) methods
in International Reviews in Physical Chemistry
Punkvang A
(2019)
Simulations of Shikimate Dehydrogenase from Mycobacterium tuberculosis in Complex with 3-Dehydroshikimate and NADPH Suggest Strategies for MtbSDH Inhibition.
in Journal of chemical information and modeling
Pucheta-Martínez E
(2016)
An Allosteric Cross-Talk Between the Activation Loop and the ATP Binding Site Regulates the Activation of Src Kinase.
in Scientific reports
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
Pradon Juliette Jeanne Marie
(2010)
Molecular modelling studies of the relationship between structure and mechanism in (beta)-lactamases and related enzymes
Phintha A
(2021)
Dissecting the low catalytic capability of flavin-dependent halogenases.
in The Journal of biological chemistry
Pentikäinen U
(2008)
Cooperative symmetric to asymmetric conformational transition of the apo-form of scavenger decapping enzyme revealed by simulations.
in Proteins
Pentikäinen U
(2009)
Lennard-Jones Parameters for B3LYP/CHARMM27 QM/MM Modeling of Nucleic Acid Bases.
in Journal of chemical theory and computation
Pennifold RC
(2017)
Correcting density-driven errors in projection-based embedding.
in The Journal of chemical physics
Pasi M
(2014)
µABC: a systematic microsecond molecular dynamics study of tetranucleotide sequence effects in B-DNA.
in Nucleic acids research
Parker JL
(2021)
Cryo-EM structure of PepT2 reveals structural basis for proton-coupled peptide and prodrug transport in mammals.
in Science advances
Paramo T
(2017)
Functional Validation of Heteromeric Kainate Receptor Models.
in Biophysical journal
Palaiokostas M
(2018)
Effects of lipid composition on membrane permeation
in Soft Matter
Pakamwong B
(2022)
Identification of Potent DNA Gyrase Inhibitors Active against Mycobacterium tuberculosis.
in Journal of chemical information and modeling
Oláh J
(2011)
Understanding the determinants of selectivity in drug metabolism through modeling of dextromethorphan oxidation by cytochrome P450.
in Proceedings of the National Academy of Sciences of the United States of America
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)
Oliveira ASF
(2020)
Simulations support the interaction of the SARS-CoV-2 spike protein with nicotinic acetylcholine receptors.
in bioRxiv : the preprint server for biology
Oliveira ASF
(2021)
Dynamical nonequilibrium molecular dynamics reveals the structural basis for allostery and signal propagation in biomolecular systems.
in The European physical journal. B
Oliveira ASF
(2023)
SARS-CoV-2 spike variants differ in their allosteric responses to linoleic acid.
in Journal of molecular cell biology
Oliveira ASF
(2021)
A potential interaction between the SARS-CoV-2 spike protein and nicotinic acetylcholine receptors.
in Biophysical journal
Oliveira A. S. F.
(2019)
Nicotine-induced conformational changes in the a4ß2 nicotinic receptor
in EUROPEAN BIOPHYSICS JOURNAL WITH BIOPHYSICS LETTERS
Oliveira A
(2023)
Fluctuation Relations to Calculate Protein Redox Potentials from Molecular Dynamics Simulations
in Journal of Chemical Theory and Computation
Oliveira A
(2019)
A General Mechanism for Signal Propagation in the Nicotinic Acetylcholine Receptor Family
in Journal of the American Chemical Society
O'Hagan MP
(2020)
Visible-light photoswitching of ligand binding mode suggests G-quadruplex DNA as a target for photopharmacology.
in Chemical communications (Cambridge, England)
O'Hagan MP
(2020)
Enhanced sampling molecular dynamics simulations correctly predict the diverse activities of a series of stiff-stilbene G-quadruplex DNA ligands.
in Chemical science
O'Hagan MP
(2019)
A Photoresponsive Stiff-Stilbene Ligand Fuels the Reversible Unfolding of G-Quadruplex DNA.
in Angewandte Chemie (International ed. in English)
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 |