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
Morando MA
(2016)
Conformational Selection and Induced Fit Mechanisms in the Binding of an Anticancer Drug to the c-Src Kinase.
in Scientific reports
Scott AJ
(2021)
Constructing ion channels from water-soluble a-helical barrels.
in Nature chemistry
Watkins DW
(2017)
Construction and in vivo assembly of a catalytically proficient and hyperthermostable de novo enzyme.
in Nature communications
Pentikäinen U
(2008)
Cooperative symmetric to asymmetric conformational transition of the apo-form of scavenger decapping enzyme revealed by simulations.
in Proteins
Pennifold RC
(2017)
Correcting density-driven errors in projection-based embedding.
in The Journal of chemical physics
Mulholland AJ
(2020)
COVID19 - Computational Chemists Meet the Moment.
in Journal of chemical information and modeling
Parker JL
(2021)
Cryo-EM structure of PepT2 reveals structural basis for proton-coupled peptide and prodrug transport in mammals.
in Science advances
Twidale RM
(2021)
Crystallography and QM/MM Simulations Identify Preferential Binding of Hydrolyzed Carbapenem and Penem Antibiotics to the L1 Metallo-ß-Lactamase in the Imine Form.
in Journal of chemical information and modeling
Tooke CL
(2020)
Cyclic boronates as versatile scaffolds for KPC-2 ß-lactamase inhibition.
in RSC medicinal chemistry
Lonsdale Richard
(2009)
Cytochrome P450 reactivity and specificity from QM/MM modelling
Thomas F
(2018)
De Novo-Designed a-Helical Barrels as Receptors for Small Molecules.
in ACS synthetic biology
Castelli M
(2023)
Decrypting Allostery in Membrane-Bound K-Ras4B Using Complementary In Silico Approaches Based on Unbiased Molecular Dynamics Simulations
in Journal of the American Chemical Society
Bunzel HA
(2021)
Designing better enzymes: Insights from directed evolution.
in Current opinion in structural biology
Lonsdale R
(2012)
Determinants of reactivity and selectivity in soluble epoxide hydrolase from quantum mechanics/molecular mechanics modeling.
in Biochemistry
Ingram M
(2020)
Determination of Bubble Size Distribution Using Ultrasound Array Imaging.
in IEEE transactions on ultrasonics, ferroelectrics, and frequency control
May PW
(2016)
Diamond-coated 'black silicon' as a promising material for high-surface-area electrochemical electrodes and antibacterial surfaces.
in Journal of materials chemistry. B
Miles BT
(2017)
Direct Evidence of Lack of Colocalisation of Fluorescently Labelled Gold Labels Used in Correlative Light Electron Microscopy.
in Scientific reports
Steventon Jack
(2020)
Directed evolution of artificial oxidoreductases
Kamsri P
(2020)
Discovery of New and Potent InhA Inhibitors as Antituberculosis Agents: Structure-Based Virtual Screening Validated by Biological Assays and X-ray Crystallography.
in Journal of chemical information and modeling
Hanwarinroj C
(2022)
Discovery of novel and potent InhA inhibitors by an in silico screening and pharmacokinetic prediction
in Future Medicinal Chemistry
Chan HTH
(2021)
Discovery of SARS-CoV-2 Mpro peptide inhibitors from modelling substrate and ligand binding.
in Chemical science
Mulholland A
(2016)
Dispelling the effects of a sorceress in enzyme catalysis
in Proceedings of the National Academy of Sciences
Phintha A
(2021)
Dissecting the low catalytic capability of flavin-dependent halogenases.
in The Journal of biological chemistry
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
Van Der Kamp M
(2017)
Dynamical origins of heat capacity changes in enzyme catalysed reactions
Van Der Kamp MW
(2018)
Dynamical origins of heat capacity changes in enzyme-catalysed reactions.
in Nature communications
Lonsdale R
(2012)
Effects of Dispersion in Density Functional Based Quantum Mechanical/Molecular Mechanical Calculations on Cytochrome P450 Catalyzed Reactions.
in Journal of chemical theory and computation
Ding W
(2017)
Effects of High Pressure on Phospholipid Bilayers.
in The journal of physical chemistry. B
Ding W
(2015)
Effects of Lipid Composition on Bilayer Membranes Quantified by All-Atom Molecular Dynamics.
in The journal of physical chemistry. B
Palaiokostas M
(2018)
Effects of lipid composition on membrane permeation
in Soft Matter
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
Calabrò G
(2016)
Elucidation of Nonadditive Effects in Protein-Ligand Binding Energies: Thrombin as a Case Study.
in The journal of physical chemistry. B
Bunzel HA
(2019)
Emergence of a Negative Activation Heat Capacity during Evolution of a Designed Enzyme.
in Journal of the American Chemical Society
Mujika J
(2012)
Encyclopedia of Inorganic and Bioinorganic Chemistry
Mujika J
(2005)
Encyclopedia of Inorganic Chemistry
Mendoza-Martinez C
(2022)
Energetics of a protein disorder-order transition in small molecule recognition
in Chemical Science
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
Hanpaibool C
(2023)
Enhancement by pyrazolones of colistin efficacy against mcr-1-expressing E. coli: an in silico and in vitro investigation.
in Journal of computer-aided molecular design
Zinovjev K
(2020)
Enlighten2: molecular dynamics simulations of protein-ligand systems made accessible.
in Bioinformatics (Oxford, England)
Ali HS
(2019)
Entropy of Simulated Liquids Using Multiscale Cell Correlation.
in Entropy (Basel, Switzerland)
Dooher T
(2022)
Environmental stress cracking of polymers: Case studies from industry (ABS and LDPE)
in Engineering Failure Analysis
Mulholland A
(2012)
Enzyme dynamics and catalysis in the mechanism of DNA polymerase
in Theoretical Chemistry Accounts
Arcus VL
(2020)
Enzyme evolution and the temperature dependence of enzyme catalysis.
in Current opinion in structural biology
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 |