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
Kiani YS
(2019)
Molecular Dynamics Simulation Framework to Probe the Binding Hypothesis of CYP3A4 Inhibitors.
in International journal of molecular sciences
Tooke CL
(2019)
Molecular Basis of Class A ß-Lactamase Inhibition by Relebactam.
in Antimicrobial agents and chemotherapy
Bunzel HA
(2019)
Emergence of a Negative Activation Heat Capacity during Evolution of a Designed Enzyme.
in Journal of the American Chemical Society
Lester Hedges
(2019)
BioSimSpace: An interoperable Python framework for biomolecular simulation
in Zenodo
Evans LE
(2019)
Exploitation of Antibiotic Resistance as a Novel Drug Target: Development of a ß-Lactamase-Activated Antibacterial Prodrug.
in Journal of medicinal chemistry
Tooke Catherine L.
(2019)
Interactions of ß-lactamases with antibiotics and inhibitors
in JOURNAL OF BIOMOLECULAR STRUCTURE & DYNAMICS
Hirvonen VHA
(2019)
An Efficient Computational Assay for ß-Lactam Antibiotic Breakdown by Class A ß-Lactamases.
in Journal of chemical information and modeling
Jenkins Jonathan
(2019)
Computational design, construction, and characterisation of artificial peroxidases
Hedges L
(2019)
BioSimSpace: An interoperable Python framework for biomolecular simulation
in Journal of Open Source Software
Meletiou A
(2019)
Tios: The Internet of Simulations. Turning Molecular Dynamics into a Data Streaming Web Application.
in Journal of chemical information and modeling
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)
Sartori GR
(2019)
Ligand-induced conformational selection predicts the selectivity of cysteine protease inhibitors.
in PloS one
Lester Hedges
(2019)
BioSimSpace: An interoperable Python framework for biomolecular simulation
in Zenodo
Oliveira A. S. F.
(2019)
Nicotine-induced conformational changes in the a4ß2 nicotinic receptor
in EUROPEAN BIOPHYSICS JOURNAL WITH BIOPHYSICS LETTERS
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)
Newport TD
(2019)
The MemProtMD database: a resource for membrane-embedded protein structures and their lipid interactions.
in Nucleic acids research
Nutho B
(2019)
The reaction mechanism of Zika virus NS2B/NS3 serine protease inhibition by dipeptidyl aldehyde: a QM/MM study.
in Physical chemistry chemical physics : PCCP
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
O'Hagan M
(2019)
A Photoresponsive Stiff-Stilbene Ligand Fuels the Reversible Unfolding of G-Quadruplex DNA
in Angewandte Chemie
Leferink NGH
(2019)
Experiment and Simulation Reveal How Mutations in Functional Plasticity Regions Guide Plant Monoterpene Synthase Product Outcome.
in ACS catalysis
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
Voice A
(2019)
Limitations of Ligand-Only Approaches for Predicting the Reactivity of Covalent Inhibitors.
in Journal of chemical information and modeling
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
Ali HS
(2019)
Entropy of Simulated Liquids Using Multiscale Cell Correlation.
in Entropy (Basel, Switzerland)
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)
Bennie S
(2019)
Teaching Enzyme Catalysis Using Interactive Molecular Dynamics in Virtual Reality
in Journal of Chemical Education
Oliveira A
(2019)
A General Mechanism for Signal Propagation in the Nicotinic Acetylcholine Receptor Family
in Journal of the American Chemical Society
Shahane G
(2019)
Interaction of Antimicrobial Lipopeptides with Bacterial Lipid Bilayers
in The Journal of Membrane Biology
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
O'Hagan MP
(2019)
A Photoresponsive Stiff-Stilbene Ligand Fuels the Reversible Unfolding of G-Quadruplex DNA.
in Angewandte Chemie (International ed. in English)
De Simone A
(2019)
A computationally designed binding mode flip leads to a novel class of potent tri-vector cyclophilin inhibitors.
in Chemical science
Wang L
(2020)
Mixing and matching genes of marine and terrestrial origin in the biosynthesis of the mupirocin antibiotics.
in Chemical science
Rao S
(2020)
Characterizing Membrane Association and Periplasmic Transfer of Bacterial Lipoproteins through Molecular Dynamics Simulations.
in Structure (London, England : 1993)
Arcus VL
(2020)
Enzyme evolution and the temperature dependence of enzyme catalysis.
in Current opinion in structural biology
Mulholland A
(2020)
Catalytic Mechanism of the Colistin Resistance Protein MCR-1
Tooke Catherine L.
(2020)
Substrate and inhibitor interactions of class A ß-lactamases
Kwon H
(2020)
Visualizing the protons in a metalloenzyme electron proton transfer pathway.
in Proceedings of the National Academy of Sciences of the United States of America
Yamamoto E
(2020)
Multiple lipid binding sites determine the affinity of PH domains for phosphoinositide-containing membranes.
in Science advances
Douglas-Gallardo OA
(2020)
Electronic structure benchmark calculations of CO2 fixing elementary chemical steps in RuBisCO using the projector-based embedding approach.
in Journal of computational chemistry
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)
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 |