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
Lawan N
(2022)
QM/MM Molecular Modeling Reveals Mechanism Insights into Flavin Peroxide Formation in Bacterial Luciferase.
in Journal of chemical information and modeling
Lence E
(2018)
QM/MM simulations identify the determinants of catalytic activity differences between type II dehydroquinase enzymes
in Organic & Biomolecular Chemistry
Bowman AL
(2008)
QM/MM simulations predict a covalent intermediate in the hen egg white lysozyme reaction with its natural substrate.
in Chemical communications (Cambridge, England)
Chudyk EI
(2022)
QM/MM Simulations Reveal the Determinants of Carbapenemase Activity in Class A ß-Lactamases.
in ACS infectious diseases
Nutho B
(2019)
QM/MM Study on Cleavage Mechanism Catalyzed by Zika Virus NS2B/NS3 Serine Protease
in Biophysical Journal
Szeto M
(2009)
QM/MM study on the mechanism of peptide hydrolysis by carboxypeptidase A
in Journal of Molecular Structure: THEOCHEM
Szefczyk B
(2007)
Quantum chemical analysis of reaction paths in chorismate mutase: Conformational effects and electrostatic stabilization
in International Journal of Quantum Chemistry
Nutho B
(2019)
Quantum Mechanics/Molecular Mechanics (QM/MM) Calculations Support a Concerted Reaction Mechanism for the Zika Virus NS2B/NS3 Serine Protease with Its Substrate.
in The journal of physical chemistry. B
Jitonnom J
(2017)
Quantum Mechanics/Molecular Mechanics Simulations Identify the Ring-Opening Mechanism of Creatininase.
in Biochemistry
Limb MAL
(2019)
Quantum Mechanics/Molecular Mechanics Simulations Show Saccharide Distortion is Required for Reaction in Hen Egg-White Lysozyme.
in Chemistry (Weinheim an der Bergstrasse, Germany)
Woods CJ
(2014)
Rapid decomposition and visualisation of protein-ligand binding free energies by residue and by water.
in Faraday discussions
Beker W
(2017)
Rapid Estimation of Catalytic Efficiency by Cumulative Atomic Multipole Moments: Application to Ketosteroid Isomerase Mutants.
in Journal of chemical theory and computation
Daniels AD
(2014)
Reaction mechanism of N-acetylneuraminic acid lyase revealed by a combination of crystallography, QM/MM simulation, and mutagenesis.
in ACS chemical biology
Singh J
(2021)
Real-time super-resolution mapping of locally anisotropic grain orientations for ultrasonic non-destructive evaluation of crystalline material
in Neural Computing and Applications
Ansell TB
(2021)
Relative Affinities of Protein-Cholesterol Interactions from Equilibrium Molecular Dynamics Simulations.
in Journal of chemical theory and computation
Lythell E
(2020)
Resistance to the "last resort" antibiotic colistin: a single-zinc mechanism for phosphointermediate formation in MCR enzymes.
in Chemical communications (Cambridge, England)
Hindson S
(2021)
Rigidifying a De Novo Enzyme Increases Activity and Induces a Negative Activation Heat Capacity
in ACS Catalysis
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
O'Connor M
(2018)
Sampling molecular conformations and dynamics in a multiuser virtual reality framework
in Science Advances
Oliveira ASF
(2023)
SARS-CoV-2 spike variants differ in their allosteric responses to linoleic acid.
in Journal of molecular cell biology
Mulholland AJ
(2023)
Science after Brexit: bright spots on the Horizon?
in EMBO reports
Catlow CR
(2020)
Science to enable the circular economy.
in Philosophical transactions. Series A, Mathematical, physical, and engineering sciences
Espejo-Román JM
(2022)
Selective Anticancer Therapy Based on a HA-CD44 Interaction Inhibitor Loaded on Polymeric Nanoparticles.
in Pharmaceutics
Ranaghan K
(2016)
Simulating Enzyme Reactivity - Computational Methods in Enzyme Catalysis
Bowman A
(2007)
Simulating Enzyme-Catalyzed Reactions
in ChemInform
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
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
Schiffrin B
(2016)
Skp is a multivalent chaperone of outer-membrane proteins.
in Nature structural & molecular biology
Hirvonen V
(2020)
Small Changes in Hydration Determine Cephalosporinase Activity of OXA-48 ß-Lactamases
in ACS Catalysis
Larsen AH
(2022)
Specific interactions of peripheral membrane proteins with lipids: what can molecular simulations show us?
in Bioscience reports
Fan Betty
(2017)
Stage IV breast cancer is increased by omitting screening mammography
in ANNALS OF SURGICAL ONCOLOGY
Song W
(2019)
State-dependent Lipid Interactions with the A2a Receptor Revealed by MD Simulations Using In Vivo-Mimetic Membranes.
in Structure (London, England : 1993)
Musgaard M
(2016)
Steered Molecular Dynamics Simulations Predict Conformational Stability of Glutamate Receptors.
in Journal of chemical information and modeling
Van Den Berg B
(2016)
Structural basis for Mep2 ammonium transceptor activation by phosphorylation.
in Nature communications
Karuppiah V
(2017)
Structural Basis of Catalysis in the Bacterial Monoterpene Synthases Linalool Synthase and 1,8-Cineole Synthase.
in ACS catalysis
Lodola A
(2010)
Structural Fluctuations in Enzyme-Catalyzed Reactions: Determinants of Reactivity in Fatty Acid Amide Hydrolase from Multivariate Statistical Analysis of Quantum Mechanics/Molecular Mechanics Paths.
in Journal of chemical theory and computation
Ainsley J
(2018)
Structural Insights from Molecular Dynamics Simulations of Tryptophan 7-Halogenase and Tryptophan 5-Halogenase.
in ACS omega
Gupta K
(2022)
Structural insights in cell-type specific evolution of intra-host diversity by SARS-CoV-2.
in Nature communications
Dawson WM
(2021)
Structural resolution of switchable states of a de novo peptide assembly.
in Nature communications
Wells S
(2015)
Structure and Function in Homodimeric Enzymes: Simulations of Cooperative and Independent Functional Motions
in PLOS ONE
Dunseath O
(2019)
Studies of Black Diamond as an antibacterial surface for Gram Negative bacteria: the interplay between chemical and mechanical bactericidal activity.
in Scientific reports
Hazell G
(2018)
Studies of black silicon and black diamond as materials for antibacterial surfaces.
in Biomaterials science
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 |