1
|
Pederson JP, McDaniel JG. PyDFT-QMMM: A modular, extensible software framework for DFT-based QM/MM molecular dynamics. J Chem Phys 2024; 161:034103. [PMID: 39007371 DOI: 10.1063/5.0219851] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/20/2024] [Accepted: 06/24/2024] [Indexed: 07/16/2024] Open
Abstract
PyDFT-QMMM is a Python-based package for performing hybrid quantum mechanics/molecular mechanics (QM/MM) simulations at the density functional level of theory. The program is designed to treat short-range and long-range interactions through user-specified combinations of electrostatic and mechanical embedding procedures within periodic simulation domains, providing necessary interfaces to external quantum chemistry and molecular dynamics software. To enable direct embedding of long-range electrostatics in periodic systems, we have derived and implemented force terms for our previously described QM/MM/PME approach [Pederson and McDaniel, J. Chem. Phys. 156, 174105 (2022)]. Communication with external software packages Psi4 and OpenMM is facilitated through Python application programming interfaces (APIs). The core library contains basic utilities for running QM/MM molecular dynamics simulations, and plug-in entry-points are provided for users to implement custom energy/force calculation and integration routines, within an extensible architecture. The user interacts with PyDFT-QMMM primarily through its Python API, allowing for complex workflow development with Python scripting, for example, interfacing with PLUMED for free energy simulations. We provide benchmarks of forces and energy conservation for the QM/MM/PME and alternative QM/MM electrostatic embedding approaches. We further demonstrate a simple example use case for water solute in a water solvent system, for which radial distribution functions are computed from 100 ps QM/MM simulations; in this example, we highlight how the solvation structure is sensitive to different basis-set choices due to under- or over-polarization of the QM water molecule's electron density.
Collapse
Affiliation(s)
- John P Pederson
- School of Chemistry and Biochemistry, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, USA
| | - Jesse G McDaniel
- School of Chemistry and Biochemistry, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, USA
| |
Collapse
|
2
|
Cerutti DS, Wiewiora R, Boothroyd S, Sherman W. STORMM: Structure and topology replica molecular mechanics for chemical simulations. J Chem Phys 2024; 161:032501. [PMID: 39007368 DOI: 10.1063/5.0211032] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/28/2024] [Accepted: 06/26/2024] [Indexed: 07/16/2024] Open
Abstract
The Structure and TOpology Replica Molecular Mechanics (STORMM) code is a next-generation molecular simulation engine and associated libraries optimized for performance on fast, vectorized central processor units and graphics processing units (GPUs) with independent memory and tens of thousands of threads. STORMM is built to run thousands of independent molecular mechanical calculations on a single GPU with novel implementations that tune numerical precision, mathematical operations, and scarce on-chip memory resources to optimize throughput. The libraries are built around accessible classes with detailed documentation, supporting fine-grained parallelism and algorithm development as well as copying or swapping groups of systems on and off of the GPU. A primary intention of the STORMM libraries is to provide developers of atomic simulation methods with access to a high-performance molecular mechanics engine with extensive facilities to prototype and develop bespoke tools aimed toward drug discovery applications. In its present state, STORMM delivers molecular dynamics simulations of small molecules and small proteins in implicit solvent with tens to hundreds of times the throughput of conventional codes. The engineering paradigm transforms two of the most memory bandwidth-intensive aspects of condensed-phase dynamics, particle-mesh mapping, and valence interactions, into compute-bound problems for several times the scalability of existing programs. Numerical methods for compressing and streamlining the information present in stored coordinates and lookup tables are also presented, delivering improved accuracy over methods implemented in other molecular dynamics engines. The open-source code is released under the MIT license.
Collapse
Affiliation(s)
| | | | | | - Woody Sherman
- Psivant Therapeutics, Boston, Massachusetts 02210, USA
| |
Collapse
|
3
|
La Force H, Freindorf M, Kraka E. Ligand Characterization and DNA Intercalation of Ru(II) Polypyridyl Complexes: A Local Vibrational Mode Study. J Phys Chem A 2024. [PMID: 38990174 DOI: 10.1021/acs.jpca.4c02954] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/12/2024]
Abstract
We investigated in this work ruthenium-ligand bonding across the RuN framework in 12 Ru(II) polypyridyl complexes in the gas phase and solution for both singlet and triplet states, in addition to their affinity for DNA binding through π-π stacking interactions with DNA nucleobases. As a tool to assess the intrinsic strength of the ruthenium-ligand bonds, we determined local vibrational force constants via our local vibrational mode analysis software. We introduced a novel local force constant that directly accounts for the intrinsic strength of the π-π stacking interaction between DNA and the intercalated Ru(II) complex. According to our findings, [Ru(phen)2(dppz)]2+ and [Ru(phen)2(11-CN-dppz)]2+ provide an intriguing trade-off between photoinduced complex excitation and the strength of the subsequent π-π stacking interaction with DNA. [Ru(phen)2(dppz)]2+ displays a small singlet-triplet splitting and a strong π-π stacking interaction in its singlet state, suggesting a favorable photoexcitation but potentially weaker interaction with DNA in the excited state. Conversely, [Ru(phen)2(11-CN-dppz)]2+ exhibits a larger singlet-triplet splitting and a stronger π-π stacking interaction with DNA in its triplet state, indicating a less favorable photoinduced transition but a stronger interaction with DNA postexcitation. We hope our study will inspire future experimental and computational work aimed at the design of novel Ru-polypyridyl drug candidates and that our new quantitative measure of π-π stacking interactions in DNA will find a general application in the field.
Collapse
Affiliation(s)
- Hunter La Force
- Computational and Theoretical Chemistry Group (CATCO), Department of Chemistry, Southern Methodist University, 3215 Daniel Avenue, Dallas, Texas 75275-0314, United States
| | - Marek Freindorf
- Computational and Theoretical Chemistry Group (CATCO), Department of Chemistry, Southern Methodist University, 3215 Daniel Avenue, Dallas, Texas 75275-0314, United States
| | - Elfi Kraka
- Computational and Theoretical Chemistry Group (CATCO), Department of Chemistry, Southern Methodist University, 3215 Daniel Avenue, Dallas, Texas 75275-0314, United States
| |
Collapse
|
4
|
Huang H, Peng J, Zhang Y, Gu FL, Lan Z, Xu C. The development of the QM/MM interface and its application for the on-the-fly QM/MM nonadiabatic dynamics in JADE package: Theory, implementation, and applications. J Chem Phys 2024; 160:234101. [PMID: 38884395 DOI: 10.1063/5.0215036] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/22/2024] [Accepted: 05/15/2024] [Indexed: 06/18/2024] Open
Abstract
Understanding the nonadiabatic dynamics of complex systems is a challenging task in computational photochemistry. Herein, we present an efficient and user-friendly quantum mechanics/molecular mechanics (QM/MM) interface to run on-the-fly nonadiabatic dynamics. Currently, this interface consists of an independent set of codes designed for general-purpose use. Herein, we demonstrate the ability and feasibility of the QM/MM interface by integrating it with our long-term developed JADE package. Tailored to handle nonadiabatic processes in various complex systems, especially condensed phases and protein environments, we delve into the theories, implementations, and applications of on-the-fly QM/MM nonadiabatic dynamics. The QM/MM approach is established within the framework of the additive QM/MM scheme, employing electrostatic embedding, link-atom inclusion, and charge-redistribution schemes to treat the QM/MM boundary. Trajectory surface-hopping dynamics are facilitated using the fewest switches algorithm, encompassing classical and quantum treatments for nuclear and electronic motions, respectively. Finally, we report simulations of nonadiabatic dynamics for two typical systems: azomethane in water and the retinal chromophore PSB3 in a protein environment. Our results not only illustrate the power of the QM/MM program but also reveal the important roles of environmental factors in nonadiabatic processes.
Collapse
Affiliation(s)
- Haiyi Huang
- MOE Key Laboratory of Environmental Theoretical Chemistry and Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety, SCNU Environmental Research Institute, School of Environment, South China Normal University, Guangzhou 510006, China
- Center for Advanced Materials Research, Beijing Normal University, Zhuhai 519087, China
- MOE Key Laboratory of Theoretical and Computational Photochemistry, College of Chemistry, Beijing Normal University, Beijing 100875, China
| | - Jiawei Peng
- MOE Key Laboratory of Environmental Theoretical Chemistry and Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety, SCNU Environmental Research Institute, School of Environment, South China Normal University, Guangzhou 510006, China
| | - Yulin Zhang
- MOE Key Laboratory of Environmental Theoretical Chemistry and Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety, SCNU Environmental Research Institute, School of Environment, South China Normal University, Guangzhou 510006, China
| | - Feng Long Gu
- MOE Key Laboratory of Environmental Theoretical Chemistry and Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety, SCNU Environmental Research Institute, School of Environment, South China Normal University, Guangzhou 510006, China
| | - Zhenggang Lan
- MOE Key Laboratory of Environmental Theoretical Chemistry and Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety, SCNU Environmental Research Institute, School of Environment, South China Normal University, Guangzhou 510006, China
| | - Chao Xu
- MOE Key Laboratory of Environmental Theoretical Chemistry and Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety, SCNU Environmental Research Institute, School of Environment, South China Normal University, Guangzhou 510006, China
| |
Collapse
|
5
|
Khatun S, Bhagat RP, Amin SA, Jha T, Gayen S. Density functional theory (DFT) studies in HDAC-based chemotherapeutics: Current findings, case studies and future perspectives. Comput Biol Med 2024; 175:108468. [PMID: 38657469 DOI: 10.1016/j.compbiomed.2024.108468] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/02/2023] [Revised: 04/08/2024] [Accepted: 04/09/2024] [Indexed: 04/26/2024]
Abstract
Density Functional Theory (DFT) is a quantum chemical computational method used to predict and analyze the electronic properties of atoms, molecules, and solids based on the density of electrons rather than wavefunctions. It provides insights into the structure, bonding, and behavior of different molecules, including those involved in the development of chemotherapeutic agents, such as histone deacetylase inhibitors (HDACis). HDACs are a wide group of metalloenzymes that facilitate the removal of acetyl groups from acetyl-lysine residues situated in the N-terminal tail of histones. Abnormal HDAC recruitment has been linked to several human diseases, especially cancer. Therefore, it has been recognized as a prospective target for accelerating the development of anticancer therapies. Researchers have studied HDACs and its inhibitors extensively using a combination of experimental methods and diverse in-silico approaches such as machine learning and quantitative structure-activity relationship (QSAR) methods, molecular docking, molecular dynamics, pharmacophore mapping, and more. In this context, DFT studies can make significant contribution by shedding light on the molecular properties, interactions, reaction pathways, transition states, reactivity and mechanisms involved in the development of HDACis. This review attempted to elucidate the scope in which DFT methodologies may be used to enhance our comprehension of the molecular aspects of HDAC inhibitors, aiding in the rational design and optimization of these compounds for therapeutic applications in cancer and other ailments. The insights gained can guide experimental efforts toward developing more potent and selective HDAC inhibitors.
Collapse
Affiliation(s)
- Samima Khatun
- Laboratory of Drug Design and Discovery, Department of Pharmaceutical Technology, Jadavpur University, Kolkata, 700032, India
| | - Rinki Prasad Bhagat
- Laboratory of Drug Design and Discovery, Department of Pharmaceutical Technology, Jadavpur University, Kolkata, 700032, India
| | - Sk Abdul Amin
- Department of Pharmaceutical Technology, JIS University, 81, Nilgunj Road, Agarpara, Kolkata, West Bengal, India
| | - Tarun Jha
- Natural Science Laboratory, Division of Medicinal and Pharmaceutical Chemistry, Department of Pharmaceutical Technology, Jadavpur University, Kolkata 700032, India
| | - Shovanlal Gayen
- Laboratory of Drug Design and Discovery, Department of Pharmaceutical Technology, Jadavpur University, Kolkata, 700032, India.
| |
Collapse
|
6
|
Lu J, Lv X, Yu W, Zhang J, Lu J, Liu Y, Li J, Du G, Chen J, Liu L. Reshaping Phosphatase Substrate Preference for Controlled Biosynthesis Using a "Design-Build-Test-Learn" Framework. ADVANCED SCIENCE (WEINHEIM, BADEN-WURTTEMBERG, GERMANY) 2024; 11:e2309852. [PMID: 38504470 PMCID: PMC11165480 DOI: 10.1002/advs.202309852] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/15/2023] [Revised: 02/29/2024] [Indexed: 03/21/2024]
Abstract
Biosynthesis is the application of enzymes in microbial cell factories and has emerged as a promising alternative to chemical synthesis. However, natural enzymes with limited catalytic performance often need to be engineered to meet specific needs through a time-consuming trial-and-error process. This study presents a quantum mechanics (QM)-incorporated design-build-test-learn (DBTL) framework to rationally design phosphatase BT4131, an enzyme with an ambiguous substrate spectrum involved in N-acetylglucosamine (GlcNAc) biosynthesis. First, mutant M1 (L129Q) is designed using force field-based methods, resulting in a 1.4-fold increase in substrate preference (kcat/Km) toward GlcNAc-6-phosphate (GlcNAc6P). QM calculations indicate that the shift in substrate preference is caused by a 13.59 kcal mol-1 reduction in activation energy. Furthermore, an iterative computer-aided design is conducted to stabilize the transition state. As a result, mutant M4 (I49Q/L129Q/G172L) with a 9.5-fold increase in kcat-GlcNAc6P/Km-GlcNAc6P and a 59% decrease in kcat-Glc6P/Km-Glc6P is highly desirable compared to the wild type in the GlcNAc-producing chassis. The GlcNAc titer increases to 217.3 g L-1 with a yield of 0.597 g (g glucose)-1 in a 50-L bioreactor, representing the highest reported level. Collectively, this DBTL framework provides an easy yet fascinating approach to the rational design of enzymes for industrially viable biocatalysts.
Collapse
Affiliation(s)
- Jiangong Lu
- Key Laboratory of Carbohydrate Chemistry and BiotechnologyMinistry of EducationJiangnan UniversityWuxi214122China
- Science Center for Future FoodsJiangnan UniversityWuxi214122China
| | - Xueqin Lv
- Key Laboratory of Carbohydrate Chemistry and BiotechnologyMinistry of EducationJiangnan UniversityWuxi214122China
- Science Center for Future FoodsJiangnan UniversityWuxi214122China
| | - Wenwen Yu
- Key Laboratory of Carbohydrate Chemistry and BiotechnologyMinistry of EducationJiangnan UniversityWuxi214122China
- Science Center for Future FoodsJiangnan UniversityWuxi214122China
| | - Jianing Zhang
- Key Laboratory of Carbohydrate Chemistry and BiotechnologyMinistry of EducationJiangnan UniversityWuxi214122China
- Science Center for Future FoodsJiangnan UniversityWuxi214122China
| | - Jianxing Lu
- Shandong Runde Biotechnology co., LTDTaian271200China
| | - Yanfeng Liu
- Key Laboratory of Carbohydrate Chemistry and BiotechnologyMinistry of EducationJiangnan UniversityWuxi214122China
- Science Center for Future FoodsJiangnan UniversityWuxi214122China
| | - Jianghua Li
- Science Center for Future FoodsJiangnan UniversityWuxi214122China
| | - Guocheng Du
- Key Laboratory of Carbohydrate Chemistry and BiotechnologyMinistry of EducationJiangnan UniversityWuxi214122China
- Science Center for Future FoodsJiangnan UniversityWuxi214122China
| | - Jian Chen
- Science Center for Future FoodsJiangnan UniversityWuxi214122China
| | - Long Liu
- Key Laboratory of Carbohydrate Chemistry and BiotechnologyMinistry of EducationJiangnan UniversityWuxi214122China
- Science Center for Future FoodsJiangnan UniversityWuxi214122China
| |
Collapse
|
7
|
Mirgaux M, Leherte L, Wouters J. Human indoleamine-2,3-dioxygenase 2 cofactor lability and low substrate affinity explained by homology modeling, molecular dynamics and molecular docking. J Biomol Struct Dyn 2024; 42:4475-4488. [PMID: 37301605 DOI: 10.1080/07391102.2023.2220830] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/07/2023] [Accepted: 05/28/2023] [Indexed: 06/12/2023]
Abstract
The human indoleamine-2,3-dioxygenase 2 (hIDO2) protein is growing of interest as it is increasingly implicated in multiple diseases (cancer, autoimmune diseases, COVID-19). However, it is only poorly reported in the literature. Its mode of action remains unknown because it does not seem to catalyze the reaction for which it is attributed: the degradation of the L-Tryptophan into N-formyl-kynurenine. This contrasts with its paralog, the human indoleamine-2,3-dioxygenase 1 (hIDO1), which has been extensively studied in the literature and for which several inhibitors are already in clinical trials. Yet, the recent failure of one of the most advanced hIDO1 inhibitors, the Epacadostat, could be caused by a still unknown interaction between hIDO1 and hIDO2. In order to better understand the mechanism of hIDO2, and in the absence of experimental structural data, a computational study mixing homology modeling, Molecular Dynamics, and molecular docking was conducted. The present article highlights an exacerbated lability of the cofactor as well as an inadequate positioning of the substrate in the active site of hIDO2, which might bring part of an answer to its lack of activity.Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Manon Mirgaux
- Department of Chemistry, Laboratoire de Chimie Biologique Structurale (CBS), Namur Institute of Structured Matter (NISM), Namur Research Institute for Life Science (NARILIS), University of Namur (UNamur), Namur, Belgium
| | - Laurence Leherte
- Department of Chemistry, Laboratoire de Chimie Biologique Structurale (CBS), Namur Institute of Structured Matter (NISM), Namur Research Institute for Life Science (NARILIS), University of Namur (UNamur), Namur, Belgium
| | - Johan Wouters
- Department of Chemistry, Laboratoire de Chimie Biologique Structurale (CBS), Namur Institute of Structured Matter (NISM), Namur Research Institute for Life Science (NARILIS), University of Namur (UNamur), Namur, Belgium
| |
Collapse
|
8
|
Carrasco-Busturia D, Ippoliti E, Meloni S, Rothlisberger U, Olsen JMH. Multiscale biomolecular simulations in the exascale era. Curr Opin Struct Biol 2024; 86:102821. [PMID: 38688076 DOI: 10.1016/j.sbi.2024.102821] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/21/2024] [Revised: 04/03/2024] [Accepted: 04/04/2024] [Indexed: 05/02/2024]
Abstract
The complexity of biological systems and processes, spanning molecular to macroscopic scales, necessitates the use of multiscale simulations to get a comprehensive understanding. Quantum mechanics/molecular mechanics (QM/MM) molecular dynamics (MD) simulations are crucial for capturing processes beyond the reach of classical MD simulations. The advent of exascale computing offers unprecedented opportunities for scientific exploration, not least within life sciences, where simulations are essential to unravel intricate molecular mechanisms underlying biological processes. However, leveraging the immense computational power of exascale computing requires innovative algorithms and software designs. In this context, we discuss the current status and future prospects of multiscale biomolecular simulations on exascale supercomputers with a focus on QM/MM MD. We highlight our own efforts in developing a versatile and high-performance multiscale simulation framework with the aim of efficient utilization of state-of-the-art supercomputers. We showcase its application in uncovering complex biological mechanisms and its potential for leveraging exascale computing.
Collapse
Affiliation(s)
- David Carrasco-Busturia
- DTU Chemistry, Technical University of Denmark (DTU), Kongens Lyngby, DK-2800, Denmark. https://twitter.com/@DavidCdeB
| | - Emiliano Ippoliti
- Computational Biomedicine, Institute of Advanced Simulations IAS-5/Institute for Neuroscience and Medicine INM-9, Forschungszentrum Jülich GmbH, Jülich, DE-52428, Germany
| | - Simone Meloni
- Dipartimento di Scienze Chimiche, Farmaceutiche ed Agrarie (DOCPAS), Università degli Studi di Ferrara (Unife), Ferrara, I-44121, Italy. https://twitter.com/@smeloni99
| | - Ursula Rothlisberger
- Laboratory of Computational Chemistry and Biochemistry, Institute of Chemical Sciences and Engineering, École Polytechnique Fédérale de Lausanne (EPFL), Lausanne, CH-1015, Switzerland. https://twitter.com/@lcbc_epfl
| | | |
Collapse
|
9
|
Joy C, Mandal B, Bostan D, Dubernet ML, Babikov D. Mixed quantum/classical theory (MQCT) approach to the dynamics of molecule-molecule collisions in complex systems. Faraday Discuss 2024. [PMID: 38770664 DOI: 10.1039/d3fd00166k] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/22/2024]
Abstract
We developed a general theoretical approach and a user-ready computer code that permit study of the dynamics of collisional energy transfer and ro-vibrational energy exchange in complex molecule-molecule collisions. The method is a mixture of classical and quantum mechanics. The internal ro-vibrational motion of collision partners is treated quantum mechanically using a time-dependent Schrödinger equation that captures many quantum phenomena including state quantization and zero-point energy, propensity and selection rules for state-to-state transitions, quantum symmetry and interference phenomena. A significant numerical speed up is obtained by describing the translational motion of collision partners classically, using the Ehrenfest mean-field trajectory approach. Within this framework a family of approximate methods for collision dynamics is developed. Several benchmark studies for diatomic and triatomic molecules, such as H2O and ND3 collided with He, H2 and D2, show that the results of MQCT are in good agreement with full-quantum calculations in a broad range of energies, especially at high collision energies where they become nearly identical to the full quantum results. Numerical efficiency of the method and massive parallelism of the MQCT code permit us to embrace some of the most complicated collisional systems ever studied, such as C6H6 + He, CH3COOH + He and H2O + H2O. Application of MQCT to the collisions of chiral molecules such as CH3CHCH2O + He, and to molecule-surface collisions is also possible and will be pursued in the future.
Collapse
Affiliation(s)
- Carolin Joy
- Chemistry Department, Wehr Chemistry Building, Marquette University, Milwaukee, Wisconsin 53201-1881, USA.
| | - Bikramaditya Mandal
- Chemistry Department, Wehr Chemistry Building, Marquette University, Milwaukee, Wisconsin 53201-1881, USA.
| | - Dulat Bostan
- Chemistry Department, Wehr Chemistry Building, Marquette University, Milwaukee, Wisconsin 53201-1881, USA.
| | - Marie-Lise Dubernet
- Observatoire de Paris, PSL University, Sorbonne Universite, CNRS, SYRTE, Paris, France
| | - Dmitri Babikov
- Chemistry Department, Wehr Chemistry Building, Marquette University, Milwaukee, Wisconsin 53201-1881, USA.
| |
Collapse
|
10
|
Nochebuena J, Liu S, Cisneros GA. Relative cooperativity in neutral and charged molecular clusters using QM/MM calculations. J Chem Phys 2024; 160:134301. [PMID: 38557841 DOI: 10.1063/5.0203020] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/08/2024] [Accepted: 03/13/2024] [Indexed: 04/04/2024] Open
Abstract
QM/MM methods have been used to study electronic structure properties and chemical reactivity in complex molecular systems where direct electronic structure calculations are not feasible. In our previous work, we showed that non-polarizable force fields, by design, describe intermolecular interactions through pairwise interactions, overlooking many-body interactions involving three or more particles. In contrast, polarizable force fields account partially for many-body effects through polarization, but still handle van der Waals and permanent electrostatic interactions pairwise. We showed that despite those limitations, polarizable and non-polarizable force fields can reproduce relative cooperativity achieved using density functional theory due to error compensation mechanisms. In this contribution, we assess the performance of QM/MM methods in reproducing these phenomena. Our study highlights the significance of the QM region size and force field choice in QM/MM calculations, emphasizing the importance of parameter validation to obtain accurate interaction energy predictions.
Collapse
Affiliation(s)
- Jorge Nochebuena
- Department of Physics, University of Texas at Dallas, Richardson, Texas 75080, USA
| | - Shubin Liu
- Research Computing Center, University of North Carolina, Chapel Hill, North Carolina 27599, USA
- Department of Chemistry, University of North Carolina, Chapel Hill, North Carolina 27599, USA
| | - G Andrés Cisneros
- Department of Physics, University of Texas at Dallas, Richardson, Texas 75080, USA
- Department of Chemistry and Biochemistry, University of Texas at Dallas, Richardson, Texas 75080, USA
| |
Collapse
|
11
|
Csóka J, Hégely B, Nagy PR, Kállay M. Development of analytic gradients for the Huzinaga quantum embedding method and its applications to large-scale hybrid and double hybrid DFT forces. J Chem Phys 2024; 160:124113. [PMID: 38530010 DOI: 10.1063/5.0194463] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/27/2023] [Accepted: 03/06/2024] [Indexed: 03/27/2024] Open
Abstract
The theory of analytic gradients is presented for the projector-based density functional theory (DFT) embedding approach utilizing the Huzinaga-equation. The advantages of the Huzinaga-equation-based formulation are demonstrated. In particular, it is shown that the projector employed does not appear in the Lagrangian, and the potential risk of numerical problems is avoided at the evaluation of the gradients. The efficient implementation of the analytic gradient theory is presented for approaches where hybrid DFT, second-order Møller-Plesset perturbation theory, or double hybrid DFT are embedded in lower-level DFT environments. To demonstrate the applicability of the method and to gain insight into its accuracy, it is applied to equilibrium geometry optimizations, transition state searches, and potential energy surface scans. Our results show that bond lengths and angles converge rapidly with the size of the embedded system. While providing structural parameters close to high-level quality for the embedded atoms, the embedding approach has the potential to relax the coordinates of the environment as well. Our demonstrations on a 171-atom zeolite and a 570-atom protein system show that the Huzinaga-equation-based embedding can accelerate (double) hybrid gradient computations by an order of magnitude with sufficient active regions and enables affordable force evaluations or geometry optimizations for molecules of hundreds of atoms.
Collapse
Affiliation(s)
- József Csóka
- Department of Physical Chemistry and Materials Science, Faculty of Chemical Technology and Biotechnology, Budapest University of Technology and Economics, Műegyetem rkp. 3., H-1111 Budapest, Hungary
- HUN-REN-BME Quantum Chemistry Research Group, Műegyetem rkp. 3., H-1111 Budapest, Hungary
- MTA-BME Lendület Quantum Chemistry Research Group, Műegyetem rkp. 3., H-1111 Budapest, Hungary
| | - Bence Hégely
- Department of Physical Chemistry and Materials Science, Faculty of Chemical Technology and Biotechnology, Budapest University of Technology and Economics, Műegyetem rkp. 3., H-1111 Budapest, Hungary
- HUN-REN-BME Quantum Chemistry Research Group, Műegyetem rkp. 3., H-1111 Budapest, Hungary
- MTA-BME Lendület Quantum Chemistry Research Group, Műegyetem rkp. 3., H-1111 Budapest, Hungary
| | - Péter R Nagy
- Department of Physical Chemistry and Materials Science, Faculty of Chemical Technology and Biotechnology, Budapest University of Technology and Economics, Műegyetem rkp. 3., H-1111 Budapest, Hungary
- HUN-REN-BME Quantum Chemistry Research Group, Műegyetem rkp. 3., H-1111 Budapest, Hungary
- MTA-BME Lendület Quantum Chemistry Research Group, Műegyetem rkp. 3., H-1111 Budapest, Hungary
| | - Mihály Kállay
- Department of Physical Chemistry and Materials Science, Faculty of Chemical Technology and Biotechnology, Budapest University of Technology and Economics, Műegyetem rkp. 3., H-1111 Budapest, Hungary
- HUN-REN-BME Quantum Chemistry Research Group, Műegyetem rkp. 3., H-1111 Budapest, Hungary
- MTA-BME Lendület Quantum Chemistry Research Group, Műegyetem rkp. 3., H-1111 Budapest, Hungary
| |
Collapse
|
12
|
Song P, Wang S, Han R, Wang H, Hu B, Luan J, Zhang H, Wang Z, Ma C, Wang J. Insights into the selective mechanism of PDE2/9a inhibitors from silico aspects. J Biomol Struct Dyn 2024:1-18. [PMID: 38525932 DOI: 10.1080/07391102.2024.2331098] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/17/2024] [Accepted: 03/10/2024] [Indexed: 03/26/2024]
Abstract
The selective design of competitive enzyme inhibitors is an extremely difficult task but necessary work for certain types of systems, such as the phosphodiesterase (PDE) system addressed in this article. In the PDE family, PDE2A and PDE9 respectively target the central nervous system and heart failure, and share many conserved amino acids at their binding sites. Therefore, gaining a deep understanding of the selective mechanisms of PDE2A/9A is crucial for designing highly selective drugs. In this study, various computer-aided drug design (CADD) methods, including molecular docking, molecular dynamics simulations (MD), and binding free energy calculations, are employed to explore the selective mechanisms of PDE2A/9A. Overall, our research results indicate a selective design strategy for PDE2A, which involves incorporating hydrophobic or aromatic moieties into the molecular structure to better accommodate the hydrophobic pocket of PDE2A. Additionally, it is recommended to introduce functional groups capable of forming connections with selective residues, such as Phe830 and Gln812 for PDE2A, or Ala452 and Tyr424 for PDE9A, to enhance the selectivity of inhibitors targeting PDE2A/9A. This achievement is anticipated to pave the way for the development of innovative and selective small molecules targeting PDE2A/9A.Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Pengfei Song
- Key Laboratory of Structure-Based Drug Design &Discovery of Ministry of Education, Shenyang Pharmaceutical University, Shenyang, People's Republic of China
- Key Laboratory of Intelligent Drug Design and New Drug Discovery of Liaoning Province, Shenyang Pharmaceutical University, Shenyang, China
- School of Pharmaceutical Engineering, Shenyang Pharmaceutical University, Shenyang, People's Republic of China
| | - Shizhun Wang
- Key Laboratory of Structure-Based Drug Design &Discovery of Ministry of Education, Shenyang Pharmaceutical University, Shenyang, People's Republic of China
- Key Laboratory of Intelligent Drug Design and New Drug Discovery of Liaoning Province, Shenyang Pharmaceutical University, Shenyang, China
- School of Pharmaceutical Engineering, Shenyang Pharmaceutical University, Shenyang, People's Republic of China
| | - Ruiheng Han
- Pharmaceutical department, Avanc Pharmaceutical Co., Ltd., China
| | - Hanxun Wang
- Key Laboratory of Structure-Based Drug Design &Discovery of Ministry of Education, Shenyang Pharmaceutical University, Shenyang, People's Republic of China
- Key Laboratory of Intelligent Drug Design and New Drug Discovery of Liaoning Province, Shenyang Pharmaceutical University, Shenyang, China
- School of Pharmaceutical Engineering, Shenyang Pharmaceutical University, Shenyang, People's Republic of China
| | - Baichun Hu
- Key Laboratory of Structure-Based Drug Design &Discovery of Ministry of Education, Shenyang Pharmaceutical University, Shenyang, People's Republic of China
- Key Laboratory of Intelligent Drug Design and New Drug Discovery of Liaoning Province, Shenyang Pharmaceutical University, Shenyang, China
- School of Pharmaceutical Engineering, Shenyang Pharmaceutical University, Shenyang, People's Republic of China
| | - Jiasi Luan
- Key Laboratory of Structure-Based Drug Design &Discovery of Ministry of Education, Shenyang Pharmaceutical University, Shenyang, People's Republic of China
- Key Laboratory of Intelligent Drug Design and New Drug Discovery of Liaoning Province, Shenyang Pharmaceutical University, Shenyang, China
- School of Medical Devices, Shenyang Pharmaceutical University, Shenyang, People's Republic of China
| | - Haoyu Zhang
- Key Laboratory of Structure-Based Drug Design &Discovery of Ministry of Education, Shenyang Pharmaceutical University, Shenyang, People's Republic of China
- Key Laboratory of Intelligent Drug Design and New Drug Discovery of Liaoning Province, Shenyang Pharmaceutical University, Shenyang, China
- School of Pharmaceutical Engineering, Shenyang Pharmaceutical University, Shenyang, People's Republic of China
| | - Zhijian Wang
- Key Laboratory of Structure-Based Drug Design &Discovery of Ministry of Education, Shenyang Pharmaceutical University, Shenyang, People's Republic of China
- Key Laboratory of Intelligent Drug Design and New Drug Discovery of Liaoning Province, Shenyang Pharmaceutical University, Shenyang, China
- School of Pharmaceutical Engineering, Shenyang Pharmaceutical University, Shenyang, People's Republic of China
| | - Chao Ma
- Key Laboratory of Structure-Based Drug Design &Discovery of Ministry of Education, Shenyang Pharmaceutical University, Shenyang, People's Republic of China
- Key Laboratory of Intelligent Drug Design and New Drug Discovery of Liaoning Province, Shenyang Pharmaceutical University, Shenyang, China
- School of Pharmaceutical Engineering, Shenyang Pharmaceutical University, Shenyang, People's Republic of China
| | - Jian Wang
- Key Laboratory of Structure-Based Drug Design &Discovery of Ministry of Education, Shenyang Pharmaceutical University, Shenyang, People's Republic of China
- Key Laboratory of Intelligent Drug Design and New Drug Discovery of Liaoning Province, Shenyang Pharmaceutical University, Shenyang, China
- School of Pharmaceutical Engineering, Shenyang Pharmaceutical University, Shenyang, People's Republic of China
| |
Collapse
|
13
|
Quintano M, Moura RT, Kraka E. Exploring Jahn-Teller distortions: a local vibrational mode perspective. J Mol Model 2024; 30:102. [PMID: 38478107 DOI: 10.1007/s00894-024-05882-8] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/10/2023] [Accepted: 02/20/2024] [Indexed: 04/06/2024]
Abstract
The characterization of normal mode (CNM) procedure coupled with an adiabatic connection scheme (ACS) between local and normal vibrational modes, both being a part of the Local Vibrational Mode theory developed in our group, can identify spectral changes as structural fingerprints that monitor symmetry alterations, such as those caused by Jahn-Teller (JT) distortions. Employing the PBE0/Def2-TZVP level of theory, we investigated in this proof-of-concept study the hexaaquachromium cation case,[ Cr ( OH 2 ) 6 ] 3 + /[ Cr ( OH 2 ) 6 ] 2 + , as a commonly known example for a JT distortion, followed by the more difficult ferrous and ferric hexacyanide anion case,[ Fe ( CN ) 6 ] 4 - /[ Fe ( CN ) 6 ] 3 - . We found that in both cases CNM of the characteristic normal vibrational modes reflects delocalization consistent with high symmetry and ACS confirms symmetry breaking, as evidenced by the separation of axial and equatorial group frequencies. As underlined by the Cremer-Kraka criterion for covalent bonding, from[ Cr ( OH 2 ) 6 ] 3 + to[ Cr ( OH 2 ) 6 ] 2 + there is an increase in axial covalency whereas the equatorial bonds shift toward electrostatic character. From[ Fe ( CN ) 6 ] 4 - to[ Fe ( CN ) 6 ] 3 - we observed an increase in covalency without altering the bond nature. Distinct π back-donation disparity could be confirmed by comparison with the isolated CN- system. In summary, our study positions the CNM/ACS protocol as a robust tool for investigating less-explored JT distortions, paving the way for future applications.
Collapse
Affiliation(s)
- Mateus Quintano
- Computational and Theoretical Chemistry Group (CATCO), Department of Chemistry, Southern Methodist University, 3215 Daniel Ave, Dallas, TX, 75275-0314, USA
| | - Renaldo T Moura
- Computational and Theoretical Chemistry Group (CATCO), Department of Chemistry, Southern Methodist University, 3215 Daniel Ave, Dallas, TX, 75275-0314, USA
- Department of Chemistry and Physics, Center of Agrarian Sciences, Federal University of Paraiba, Areia, PB, 58397-000, Brazil
| | - Elfi Kraka
- Computational and Theoretical Chemistry Group (CATCO), Department of Chemistry, Southern Methodist University, 3215 Daniel Ave, Dallas, TX, 75275-0314, USA.
| |
Collapse
|
14
|
Jafari S, Ryde U, Irani M. QM/MM study of the catalytic reaction of aphid myrosinase. Int J Biol Macromol 2024; 262:130089. [PMID: 38360236 DOI: 10.1016/j.ijbiomac.2024.130089] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/29/2023] [Revised: 01/07/2024] [Accepted: 02/08/2024] [Indexed: 02/17/2024]
Abstract
Brevicoryne brassicae, an aphid species, exclusively consumes plants from the Brassicaceae family and employs a sophisticated defense mechanism involving a myrosinase enzyme that breaks down glucosinolates obtained from its host plants. In this work, we employed combined quantum mechanical and molecular mechanical (QM/MM) calculations and molecular dynamics (MD) simulations to study the catalytic reaction of aphid myrosinase. A proper QM region to study the myrosinase reaction should contain the whole substrate, models of Gln-19, His-122, Asp-124, Asn-166, Glu-167, Lys-173, Tyr-180, Val-228, Tyr-309, Tyr-346, Ile-347, Glu-374, Glu-423, Trp-424, and a water molecule. The calculations show that Asp-124 and Glu-423 must be charged, His-122 must be protonated on NE2, and Glu-167 must be protonated on OE2. Our model reproduces the anomeric retaining characteristic of myrosinase and indicates that the deglycosylation reaction is the rate-determining step of the reaction. Based on the calculations, we propose a reaction mechanism for aphid myrosinase-mediated hydrolysis of glucosinolates with an overall barrier of 15.2 kcal/mol. According to the results, removing a proton from Arg-312 or altering it to valine or methionine increases glycosylation barriers but decreases the deglycosylation barrier.
Collapse
Affiliation(s)
- Sonia Jafari
- Department of Chemistry, University of Kurdistan, P.O. Box 66175-416, Sanandaj, Iran
| | - Ulf Ryde
- Department of Theoretical Chemistry, Lund University, P.O. Box 124, SE-221 00 Lund, Sweden
| | - Mehdi Irani
- Department of Chemistry, University of Kurdistan, P.O. Box 66175-416, Sanandaj, Iran.
| |
Collapse
|
15
|
McFarlane NR, Harvey JN. Exploration of biochemical reactivity with a QM/MM growing string method. Phys Chem Chem Phys 2024; 26:5999-6007. [PMID: 38293892 DOI: 10.1039/d3cp05772k] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/01/2024]
Abstract
In this work, we have implemented the single-ended growing string method using a hybrid internal/Cartesian coordinate scheme within our in-house QM/MM package, QoMMMa, representing the first implementation of the growing string method in the QM/MM framework. The goal of the implementation was to facilitate generation of QM/MM reaction pathways with minimal user input, and also to improve the quality of the pathways generated as compared to the widely used adiabatic mapping approach. We have validated the algorithm against a reaction which has been studied extensively in previous computational investigations - the Claisen rearrangement catalysed by chorismate mutase. The nature of the transition state and the height of the barrier was predicted well using our algorithm, where more than 88% of the pathways generated were deemed to be of production quality. Directly compared to using adiabatic mapping, we found that while our QM/MM single-ended growing string method is slightly less efficient, it readily produces reaction pathways with fewer discontinuites and thus minimises the need for involved remapping of unsatisfactory energy profiles.
Collapse
Affiliation(s)
- Neil R McFarlane
- Department of Chemistry, KU Leuven, B-3001 Leuven, Celestijnenlaan 200f, 2404, Belgium.
| | - Jeremy N Harvey
- Department of Chemistry, KU Leuven, B-3001 Leuven, Celestijnenlaan 200f, 2404, Belgium.
| |
Collapse
|
16
|
Brooks CL, MacKerell AD, Post CB, Nilsson L. Biomolecular dynamics in the 21st century. Biochim Biophys Acta Gen Subj 2024; 1868:130534. [PMID: 38065235 PMCID: PMC10842176 DOI: 10.1016/j.bbagen.2023.130534] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/26/2023] [Revised: 11/28/2023] [Accepted: 11/29/2023] [Indexed: 01/03/2024]
Abstract
The relevance of motions in biological macromolecules has been clear since the early structural analyses of proteins by X-ray crystallography. Computer simulations have been applied to provide a deeper understanding of the dynamics of biological macromolecules since 1976, and are now a standard tool in many labs working on the structure and function of biomolecules. In this mini-review we highlight some areas of current interest and active development for simulations, in particular all-atom molecular dynamics simulations.
Collapse
Affiliation(s)
- Charles L Brooks
- University of Michigan, Department of Chemistry, Ann Arbor, MI 48109, USA.
| | | | - Carol B Post
- Purdue University, Department of Medicinal Chemistry and Molecular Pharmacology, West Lafayette, IN 47907-2091, USA.
| | - Lennart Nilsson
- Karolinska Institutet, Department of Biosciences and Nutrition, SE-14183 Huddinge, Sweden.
| |
Collapse
|
17
|
Ding Y, Huang J. DP/MM: A Hybrid Model for Zinc-Protein Interactions in Molecular Dynamics. J Phys Chem Lett 2024; 15:616-627. [PMID: 38198685 DOI: 10.1021/acs.jpclett.3c03158] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/12/2024]
Abstract
Zinc-containing proteins are vital for many biological processes, yet accurately modeling them using classical force fields is hindered by complicated polarization and charge transfer effects. This study introduces DP/MM, a hybrid force field scheme that utilizes a deep potential model to correct the atomic forces of zinc ions and their coordinated atoms, elevating them from MM to QM levels of accuracy. Trained on the difference between MM and QM atomic forces across diverse zinc coordination groups, the DP/MM model faithfully reproduces structural characteristics of zinc coordination during simulations, such as the tetrahedral coordination of Cys4 and Cys3His1 groups. Furthermore, DP/MM allows water exchange in the zinc coordination environment. With its unique blend of accuracy, efficiency, flexibility, and transferability, DP/MM serves as a valuable tool for studying structures and dynamics of zinc-containing proteins and also represents a pioneering approach in the evolving landscape of machine learning potentials for molecular modeling.
Collapse
Affiliation(s)
- Ye Ding
- College of Life Sciences, Zhejiang University, Hangzhou, Zhejiang 310027, China
- School of Life Sciences, Westlake University, Hangzhou, Zhejiang 310024, China
- Westlake AI Therapeutics Lab, Westlake Laboratory of Life Sciences and Biomedicine, Hangzhou, Zhejiang 310024, China
| | - Jing Huang
- School of Life Sciences, Westlake University, Hangzhou, Zhejiang 310024, China
- Westlake AI Therapeutics Lab, Westlake Laboratory of Life Sciences and Biomedicine, Hangzhou, Zhejiang 310024, China
| |
Collapse
|
18
|
Almeida NMS, Bali SK, James D, Wang C, Wilson AK. Binding of Per- and Polyfluoroalkyl Substances (PFAS) to the PPARγ/RXRα-DNA Complex. J Chem Inf Model 2023; 63:7423-7443. [PMID: 37990410 DOI: 10.1021/acs.jcim.3c01384] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/23/2023]
Abstract
Nuclear receptors are the fundamental building blocks of gene expression regulation and the focus of many drug targets. While binding to DNA, nuclear receptors act as transcription factors, governing a multitude of functions in the human body. Peroxisome proliferator-activator receptor γ (PPARγ) and the retinoid X receptor α (RXRα) form heterodimers with unique properties and have a primordial role in insulin sensitization. This PPARγ/RXRα heterodimer has been shown to be impacted by per- and polyfluoroalkyl substances (PFAS) and linked to a variety of significant health conditions in humans. Herein, a selection of the most common PFAS (legacy and emerging) was studied utilizing molecular dynamics simulations for PPARγ/RXRα. The local and global structural effects of PFAS binding on the known ligand binding pockets of PPARγ and RXRα as well as the DNA binding domain (DBD) of RXRα were inspected. The binding free energies were predicted computationally and were compared between the different binding pockets. In addition, two electronic structure approaches were utilized to model the interaction of PFAS within the DNA binding domain, density functional theory (DFT) and domain-based pair natural orbital coupled cluster with perturbative triples (DLPNO-CCSD(T)) approaches, with implicit solvation. Residue decomposition and hydrogen-bonding analysis were also performed, detailing the role of prominent residues in molecular recognition. The role of l-carnitine is explored as a potential in vivo remediation strategy for PFAS interaction with the PPARγ/RXRα heterodimer. In this work, it was found that PFAS can bind and act as agonists for all of the investigated pockets. For the first time in the literature, PFAS are postulated to bind to the DNA binding domain in a nonspecific manner. In addition, for the PPARγ ligand binding domain, l-carnitine shows promise in replacing smaller PFAS from the pocket.
Collapse
Affiliation(s)
- Nuno M S Almeida
- Department of Chemistry, Michigan State University, East Lansing, Michigan 48864, United States
| | - Semiha Kevser Bali
- Department of Chemistry, Michigan State University, East Lansing, Michigan 48864, United States
| | - Deepak James
- Department of Chemistry, Michigan State University, East Lansing, Michigan 48864, United States
| | - Cong Wang
- Department of Chemistry, Michigan State University, East Lansing, Michigan 48864, United States
| | - Angela K Wilson
- Department of Chemistry, Michigan State University, East Lansing, Michigan 48864, United States
| |
Collapse
|
19
|
Mysovsky AS, Bogdanov AI. Seamless Multilayer─A Novel Total Energy Partition Scheme for Embedded and Hybrid Calculations. J Chem Theory Comput 2023. [PMID: 37973151 DOI: 10.1021/acs.jctc.3c00666] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/19/2023]
Abstract
In this paper, we provide general formulation of a multilayer approach, covering both additive and subtractive quantum mechanics/molecular mechanics (QM/MM) as special cases. After that, we suggest a novel definition of QM/MM total energy based on the consideration of a system divided into three layers. In a simplified form, it is E = E Q M ( 1 + 2 ) - E Q M ( 2 ) + E M M ( 2 + 3 ) , where layers 1, 2, and 3 represent inner QM, outer QM, and classical MM regions, respectively. The novel formulation is also not limited by only QM/MM combination of methods─in fact, any computational methods can be combined in a hybrid calculation. In this paper, we call the new approach seamless multilayer. Test calculations performed for silica and boric oxide show that the new approach requires no QM/MM interface parameterization as well as no or very simple correction terms for boundary atoms. This can greatly facilitate QM/MM studies of covalent inorganic solids. However, test calculations of α-Al2O3 show that for ionic compounds, the new method requires some additional development.
Collapse
Affiliation(s)
- Andrey S Mysovsky
- A.P. Vinogradov Institute of Geochemistry SB RAS, 1a Favorsky Street, 664033 Irkutsk, Russia
- Institute of Quantum Physics, Irkutsk National Research Technical University, 83 Lermontov Street, 664074 Irkutsk, Russia
| | - Alexander I Bogdanov
- A.P. Vinogradov Institute of Geochemistry SB RAS, 1a Favorsky Street, 664033 Irkutsk, Russia
| |
Collapse
|
20
|
Freindorf M, Antonio J, Kraka E. Hydrogen Sulfide Ligation in Hemoglobin I of Lucina pectinata─A QM/MM and Local Mode Study. J Phys Chem A 2023; 127:8316-8329. [PMID: 37774120 DOI: 10.1021/acs.jpca.3c04399] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/01/2023]
Abstract
In this study, we investigated the interaction between the H2S ligand and the heme pocket of hemoglobin I (HbI) of Lucina pectinata for the wild-type protein; three known mutations where distal glutamine is replaced by hydrophobic valine (Gln64Val) and hydrophilic histidine in both protonation forms (Gln64Hisϵ and Gln64Hisδ); five known mutations of the so-called phenyl cage, replacing the hydrophobic phenylalanines Phe29 and Phe43 with tyrosine (Tyr), valine (Val), or leucine (Leu); and two additional mutations, Phe68Tyr and Phe68Val, in order to complement previous studies with new insights about the binding mechanism at the molecular level. A particular focus was on the intrinsic strengths of the chemical bonds involved, utilizing local vibrational force constants based on combined quantum mechanical-molecular mechanical calculations. Wild-type protein and mutations clustered into two distinct groups: Group 1 protein systems with a proton acceptor in the distal protein pocket, close to one of the H2S bonds, and Group 2 protein systems without a hydrogen acceptor close by in the active site of the protein. According to our results, the interactions between H2S and HbI of Lucina pectinata involve two important elements, namely, binding of H2S to Fe of the heme group, followed by the proton transfer from the HS bond to the distal residue. The distal residue is additionally stabilized by a second proton transfer from the distal residue to COO- of the propionate group in heme. We could identify the FeS bond as a key player and discovered that the strength of this bond depends on two mutual factors, namely, the strength of the HS bond involved in the proton transfer and the electrostatic field of the protein pocket qualifying the FeS bond as a sensitive probe for monitoring changes in H2S ligation upon protein mutations. We hope our study will inspire and guide future experimental studies, targeting new promising mutations such as Phe68Tyr, Phe68Val, or Phe43Tyr/Phe68Val.
Collapse
Affiliation(s)
- Marek Freindorf
- Computational and Theoretical Chemistry Group (CATCO), Department of Chemistry, Southern Methodist University, 3215 Daniel Avenue, Dallas, Texas 75275-0314, United States
| | - Juliana Antonio
- Computational and Theoretical Chemistry Group (CATCO), Department of Chemistry, Southern Methodist University, 3215 Daniel Avenue, Dallas, Texas 75275-0314, United States
| | - Elfi Kraka
- Computational and Theoretical Chemistry Group (CATCO), Department of Chemistry, Southern Methodist University, 3215 Daniel Avenue, Dallas, Texas 75275-0314, United States
| |
Collapse
|
21
|
Zhang F, Zeng T, Wu R. QM/MM Modeling Aided Enzyme Engineering in Natural Products Biosynthesis. J Chem Inf Model 2023; 63:5018-5034. [PMID: 37556841 DOI: 10.1021/acs.jcim.3c00779] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/11/2023]
Abstract
Natural products and their derivatives are widely used across various industries, particularly pharmaceuticals. Modern engineered biosynthesis provides an alternative way of producing and meeting the growing need for diverse natural products. Natural enzymes, on the other hand, often exhibit unsatisfactory catalytic characteristics and necessitate further enzyme engineering modifications. QM/MM, as a powerful and extensively used computational tool in the field of enzyme catalysis, has been increasingly applied in rational enzyme engineering over the past decade. In this review, we summarize recent advances in QM/MM computational investigation on enzyme catalysis and enzyme engineering for natural product biosynthesis. The challenges and perspectives for future QM/MM applications aided enzyme engineering in natural product biosynthesis will also be discussed.
Collapse
Affiliation(s)
- Fan Zhang
- Guangdong Provincial Key Laboratory of Translational Cancer Research of Chinese Medicines, Joint International Research Laboratory of Translational Cancer Research of Chinese Medicines, International Institute for Translational Chinese Medicine, School of Pharmaceutical Sciences, Guangzhou University of Chinese Medicine, Guangzhou 510006, P. R. China
- School of Pharmaceutical Sciences, Guangdong Provincial Key Laboratory of New Drug Design and Evaluation, Sun Yat-sen University, Guangzhou 510006, P. R. China
| | - Tao Zeng
- School of Pharmaceutical Sciences, Guangdong Provincial Key Laboratory of New Drug Design and Evaluation, Sun Yat-sen University, Guangzhou 510006, P. R. China
| | - Ruibo Wu
- School of Pharmaceutical Sciences, Guangdong Provincial Key Laboratory of New Drug Design and Evaluation, Sun Yat-sen University, Guangzhou 510006, P. R. China
| |
Collapse
|
22
|
Wang T, He X, Li M, Shao B, Liu TY. AIMD-Chig: Exploring the conformational space of a 166-atom protein Chignolin with ab initio molecular dynamics. Sci Data 2023; 10:549. [PMID: 37607915 PMCID: PMC10444755 DOI: 10.1038/s41597-023-02465-9] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/13/2023] [Accepted: 08/11/2023] [Indexed: 08/24/2023] Open
Abstract
Molecular dynamics (MD) simulations have revolutionized the modeling of biomolecular conformations and provided unprecedented insight into molecular interactions. Due to the prohibitive computational overheads of ab initio simulation for large biomolecules, dynamic modeling for proteins is generally constrained on force field with molecular mechanics, which suffers from low accuracy as well as ignores the electronic effects. Here, we report AIMD-Chig, an MD dataset including 2 million conformations of 166-atom protein Chignolin sampled at the density functional theory (DFT) level with 7,763,146 CPU hours. 10,000 conformations were initialized covering the whole conformational space of Chignolin, including folded, unfolded, and metastable states. Ab initio simulations were driven by M06-2X/6-31 G* with a Berendsen thermostat at 340 K. We reported coordinates, energies, and forces for each conformation. AIMD-Chig brings the DFT level conformational space exploration from small organic molecules to real-world proteins. It can serve as the benchmark for developing machine learning potentials for proteins and facilitate the exploration of protein dynamics with ab initio accuracy.
Collapse
Affiliation(s)
- Tong Wang
- Microsoft Research AI4Science, Beijing, China.
| | - Xinheng He
- Microsoft Research AI4Science, Beijing, China
- Work done during an internship at Microsoft Research AI4Science, Beijing, China
- State Key Laboratory of Drug Research and CAS Key Laboratory of Receptor Research and, Shanghai Institute of Materia Medica, Chinese Academy of Sciences, Shanghai, China
- University of Chinese Academy of Sciences, Beijing, China
| | - Mingyu Li
- Microsoft Research AI4Science, Beijing, China
- Work done during an internship at Microsoft Research AI4Science, Beijing, China
- Department of Pathophysiology, Key Laboratory of Cell Differentiation and Apoptosis of Chinese Ministry of Education, Shanghai Jiao Tong University, School of Medicine, Shanghai, China
| | - Bin Shao
- Microsoft Research AI4Science, Beijing, China.
| | - Tie-Yan Liu
- Microsoft Research AI4Science, Beijing, China
| |
Collapse
|
23
|
Plett C, Katbashev A, Ehlert S, Grimme S, Bursch M. ONIOM meets xtb: efficient, accurate, and robust multi-layer simulations across the periodic table. Phys Chem Chem Phys 2023. [PMID: 37378957 DOI: 10.1039/d3cp02178e] [Citation(s) in RCA: 5] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/29/2023]
Abstract
The computational treatment of large molecular structures is of increasing interest in fields of modern chemistry. Accordingly, efficient quantum chemical approaches are needed to perform sophisticated investigations on such systems. This engaged the development of the well-established "Our own N-layered integrated molecular orbital and molecular mechanics" (ONIOM) multi-layer scheme [L. W. Chung et al., Chem. Rev., 2015, 115, 5678-5796]. In this work, we present the specific implementation of the ONIOM scheme into the xtb semi-empirical extended tight-binding program package and its application to challenging transition-metal complexes. The efficient and broadly applicable GFNn-xTB and -FF methods are applied in the ONIOM framework to elucidate reaction energies, geometry optimizations, and explicit solvation effects for metal-organic systems with up to several hundreds of atoms. It is shown that an ONIOM-based combination of density functional theory, semi-empirical, and force-field methods can be used to drastically reduce the computational costs and thus enable the investigation of huge systems at almost no significant loss in accuracy.
Collapse
Affiliation(s)
- Christoph Plett
- Mulliken Center for Theoretical Chemistry, Universität Bonn, Beringstr. 4, 53115 Bonn, Germany.
| | - Abylay Katbashev
- Mulliken Center for Theoretical Chemistry, Universität Bonn, Beringstr. 4, 53115 Bonn, Germany.
| | - Sebastian Ehlert
- Microsoft Research AI4Science, Evert van de Beekstraat 254, 1118 CZ Schiphol, The Netherlands
| | - Stefan Grimme
- Mulliken Center for Theoretical Chemistry, Universität Bonn, Beringstr. 4, 53115 Bonn, Germany.
| | - Markus Bursch
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, 45470 Mülheim an der Ruhr, Germany.
| |
Collapse
|
24
|
Bassani D, Moro S. Past, Present, and Future Perspectives on Computer-Aided Drug Design Methodologies. Molecules 2023; 28:molecules28093906. [PMID: 37175316 PMCID: PMC10180087 DOI: 10.3390/molecules28093906] [Citation(s) in RCA: 19] [Impact Index Per Article: 19.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/22/2023] [Revised: 04/28/2023] [Accepted: 05/02/2023] [Indexed: 05/15/2023] Open
Abstract
The application of computational approaches in drug discovery has been consolidated in the last decades. These families of techniques are usually grouped under the common name of "computer-aided drug design" (CADD), and they now constitute one of the pillars in the pharmaceutical discovery pipelines in many academic and industrial environments. Their implementation has been demonstrated to tremendously improve the speed of the early discovery steps, allowing for the proficient and rational choice of proper compounds for a desired therapeutic need among the extreme vastness of the drug-like chemical space. Moreover, the application of CADD approaches allows the rationalization of biochemical and interactive processes of pharmaceutical interest at the molecular level. Because of this, computational tools are now extensively used also in the field of rational 3D design and optimization of chemical entities starting from the structural information of the targets, which can be experimentally resolved or can also be obtained with other computer-based techniques. In this work, we revised the state-of-the-art computer-aided drug design methods, focusing on their application in different scenarios of pharmaceutical and biological interest, not only highlighting their great potential and their benefits, but also discussing their actual limitations and eventual weaknesses. This work can be considered a brief overview of computational methods for drug discovery.
Collapse
Affiliation(s)
- Davide Bassani
- Pharmaceutical Research & Early Development, Roche Innovation Center Basel, F. Hoffmann-La Roche Ltd., 4070 Basel, Switzerland
- Molecular Modeling Section (MMS), Department of Pharmaceutical and Pharmacological Sciences, University of Padova, Via Marzolo 5, 35131 Padova, Italy
| | - Stefano Moro
- Molecular Modeling Section (MMS), Department of Pharmaceutical and Pharmacological Sciences, University of Padova, Via Marzolo 5, 35131 Padova, Italy
| |
Collapse
|
25
|
Clemente CM, Capece L, Martí MA. Best Practices on QM/MM Simulations of Biological Systems. J Chem Inf Model 2023; 63:2609-2627. [PMID: 37100031 DOI: 10.1021/acs.jcim.2c01522] [Citation(s) in RCA: 10] [Impact Index Per Article: 10.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/28/2023]
Abstract
During the second half of the 20th century, following structural biology hallmark works on DNA and proteins, biochemists shifted their questions from "what does this molecule look like?" to "how does this process work?". Prompted by the theoretical and practical developments in computational chemistry, this led to the emergence of biomolecular simulations and, along with the 2013 Nobel Prize in Chemistry, to the development of hybrid QM/MM methods. QM/MM methods are necessary whenever the problem we want to address involves chemical reactivity and/or a change in the system's electronic structure, with archetypal examples being the studies of an enzyme's reaction mechanism and a metalloprotein's active site. In the last decades QM/MM methods have seen an increasing adoption driven by their incorporation in widely used biomolecular simulation software. However, properly setting up a QM/MM simulation is not an easy task, and several issues need to be properly addressed to obtain meaningful results. In the present work, we describe both the theoretical concepts and practical issues that need to be considered when performing QM/MM simulations. We start with a brief historical perspective on the development of these methods and describe when and why QM/MM methods are mandatory. Then we show how to properly select and analyze the performance of the QM level of theory, the QM system size, and the position and type of the boundaries. We show the relevance of performing prior QM model system (or QM cluster) calculations in a vacuum and how to use the corresponding results to adequately calibrate those derived from QM/MM. We also discuss how to prepare the starting structure and how to select an adequate simulation strategy, including those based on geometry optimizations as well as free energy methods. In particular, we focus on the determination of free energy profiles using multiple steered molecular dynamics (MSMD) combined with Jarzynski's equation. Finally, we describe the results for two illustrative and complementary examples: the reaction performed by chorismate mutase and the study of ligand binding to hemoglobins. Overall, we provide many practical recommendations (or shortcuts) together with important conceptualizations that we hope will encourage more and more researchers to incorporate QM/MM studies into their research projects.
Collapse
Affiliation(s)
- Camila M Clemente
- Departamento de Química Biológica, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires (FCEyN-UBA) e Instituto de Química Biológica de la Facultad de Ciencias Exactas y Naturales (IQUIBICEN) CONICET, Pabellón 2 de Ciudad Universitaria, Ciudad de Buenos Aires C1428EHA, Argentina
| | - Luciana Capece
- Departamento de Química Inorgánica Analítica y Química Física, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires (FCEyN-UBA) e Instituto de Química de los Materiales, Ambiente y Energía (INQUIMAE) CONICET, Pabellòn 2 de Ciudad Universitaria, Ciudad de Buenos Aires C1428EHA, Argentina
| | - Marcelo A Martí
- Departamento de Química Biológica, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires (FCEyN-UBA) e Instituto de Química Biológica de la Facultad de Ciencias Exactas y Naturales (IQUIBICEN) CONICET, Pabellón 2 de Ciudad Universitaria, Ciudad de Buenos Aires C1428EHA, Argentina
| |
Collapse
|
26
|
Bramley GA, Beynon OT, Stishenko PV, Logsdail AJ. The application of QM/MM simulations in heterogeneous catalysis. Phys Chem Chem Phys 2023; 25:6562-6585. [PMID: 36810655 DOI: 10.1039/d2cp04537k] [Citation(s) in RCA: 7] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/10/2023]
Abstract
The QM/MM simulation method is provenly efficient for the simulation of biological systems, where an interplay of extensive environment and delicate local interactions drives a process of interest through a funnel on a complex energy landscape. Recent advances in quantum chemistry and force-field methods present opportunities for the adoption of QM/MM to simulate heterogeneous catalytic processes, and their related systems, where similar intricacies exist on the energy landscape. Herein, the fundamental theoretical considerations for performing QM/MM simulations, and the practical considerations for setting up QM/MM simulations of catalytic systems, are introduced; then, areas of heterogeneous catalysis are explored where QM/MM methods have been most fruitfully applied. The discussion includes simulations performed for adsorption processes in solvent at metallic interfaces, reaction mechanisms within zeolitic systems, nanoparticles, and defect chemistry within ionic solids. We conclude with a perspective on the current state of the field and areas where future opportunities for development and application exist.
Collapse
Affiliation(s)
- Gabriel Adrian Bramley
- Cardiff Catalysis Institute, School of Chemistry, Cardiff University, Park Place, CF10 3AT, UK.
| | - Owain Tomos Beynon
- Cardiff Catalysis Institute, School of Chemistry, Cardiff University, Park Place, CF10 3AT, UK.
| | | | - Andrew James Logsdail
- Cardiff Catalysis Institute, School of Chemistry, Cardiff University, Park Place, CF10 3AT, UK.
| |
Collapse
|
27
|
Manathunga M, Aktulga HM, Götz AW, Merz KM. Quantum Mechanics/Molecular Mechanics Simulations on NVIDIA and AMD Graphics Processing Units. J Chem Inf Model 2023; 63:711-717. [PMID: 36720086 DOI: 10.1021/acs.jcim.2c01505] [Citation(s) in RCA: 14] [Impact Index Per Article: 14.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/02/2023]
Abstract
We have ported and optimized the graphics processing unit (GPU)-accelerated QUICK and AMBER-based ab initio quantum mechanics/molecular mechanics (QM/MM) implementation on AMD GPUs. This encompasses the entire Fock matrix build and force calculation in QUICK including one-electron integrals, two-electron repulsion integrals, exchange-correlation quadrature, and linear algebra operations. General performance improvements to the QUICK GPU code are also presented. Benchmarks carried out on NVIDIA V100 and AMD MI100 cards display similar performance on both hardware for standalone HF/DFT calculations with QUICK and QM/MM molecular dynamics simulations with QUICK/AMBER. Furthermore, with respect to the QUICK/AMBER release version 21, significant speedups are observed for QM/MM molecular dynamics simulations. This significantly increases the range of scientific problems that can be addressed with open-source QM/MM software on state-of-the-art computer hardware.
Collapse
Affiliation(s)
- Madushanka Manathunga
- Department of Chemistry and Department of Biochemistry and Molecular Biology, Michigan State University, East Lansing, Michigan48824-1322, United States
| | - Hasan Metin Aktulga
- Department of Computer Science and Engineering, Michigan State University, East Lansing, Michigan48824-1322, United States
| | - Andreas W Götz
- San Diego Supercomputer Center, University of California San Diego, La Jolla, California92093-0505, United States
| | - Kenneth M Merz
- Department of Chemistry and Department of Biochemistry and Molecular Biology, Michigan State University, East Lansing, Michigan48824-1322, United States
| |
Collapse
|
28
|
Quantum Mechanical Study of Oxygen Ligands Protonation for the Stable States of the Laccase Active Site. Int J Mol Sci 2023; 24:ijms24032990. [PMID: 36769314 PMCID: PMC9917769 DOI: 10.3390/ijms24032990] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/14/2023] [Revised: 01/29/2023] [Accepted: 01/31/2023] [Indexed: 02/05/2023] Open
Abstract
Laccases are enzymes catalyzing the oxidation of a wide range of organic and inorganic substrates accompanied by molecular oxygen reduction to water. Recently, oxygen reduction by laccases has been studied by single-crystal serial X-ray crystallography with increasing absorption doses at subatomic resolution. There were two determined structures corresponding to the reduced and oxidized stable states of the laccase active site. However, the protonation of the oxygen ligands involved cannot be determined even at subatomic resolution. In the present work, the protonation of oxygen ligands in the active site of laccase for the two stable states determined in the X-ray study was explored using quantum mechanical and continuum-electrostatics calculations. This is important for understanding the reaction of the oxygen reduction mechanism in laccases. The high precision of X-ray data at subatomic resolutions allowed us to optimize the quantum mechanical calculations.
Collapse
|
29
|
Liebschner D, Moriarty NW, Poon BK, Adams PD. In situ ligand restraints from quantum-mechanical methods. Acta Crystallogr D Struct Biol 2023; 79:100-110. [PMID: 36762856 PMCID: PMC9912925 DOI: 10.1107/s2059798323000025] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/05/2022] [Accepted: 01/02/2023] [Indexed: 01/21/2023] Open
Abstract
In macromolecular crystallographic structure refinement, ligands present challenges for the generation of geometric restraints due to their large chemical variability, their possible novel nature and their specific interaction with the binding pocket of the protein. Quantum-mechanical approaches are useful for providing accurate ligand geometries, but can be plagued by the number of minima in flexible molecules. In an effort to avoid these issues, the Quantum Mechanical Restraints (QMR) procedure optimizes the ligand geometry in situ, thus accounting for the influence of the macromolecule on the local energy minima of the ligand. The optimized ligand geometry is used to generate target values for geometric restraints during the crystallographic refinement. As demonstrated using a sample of >2330 ligand instances in >1700 protein-ligand models, QMR restraints generally result in lower deviations from the target stereochemistry compared with conventionally generated restraints. In particular, the QMR approach provides accurate torsion restraints for ligands and other entities.
Collapse
Affiliation(s)
- Dorothee Liebschner
- Molecular Biosciences and Integrated Bioimaging Division, Lawrence Berkeley National Laboratory, Berkeley, CA 94720, USA
| | - Nigel W. Moriarty
- Molecular Biosciences and Integrated Bioimaging Division, Lawrence Berkeley National Laboratory, Berkeley, CA 94720, USA
| | - Billy K. Poon
- Molecular Biosciences and Integrated Bioimaging Division, Lawrence Berkeley National Laboratory, Berkeley, CA 94720, USA
| | - Paul D. Adams
- Molecular Biosciences and Integrated Bioimaging Division, Lawrence Berkeley National Laboratory, Berkeley, CA 94720, USA
- Department of Bioengineering, University of California, Berkeley, CA 94720, USA
| |
Collapse
|
30
|
Wang K. GPDOCK: highly accurate docking strategy for metalloproteins based on geometric probability. Brief Bioinform 2023; 24:6987821. [PMID: 36642411 DOI: 10.1093/bib/bbac620] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/20/2022] [Revised: 12/09/2022] [Accepted: 12/17/2022] [Indexed: 01/17/2023] Open
Abstract
Accurately predicting the interaction modes for metalloproteins remains extremely challenging in structure-based drug design and mechanism analysis of enzymatic catalysis due to the complexity of metal coordination in metalloproteins. Here, we report a docking method for metalloproteins based on geometric probability (GPDOCK) with unprecedented accuracy. The docking tests of 10 common metal ions with 9360 metalloprotein-ligand complexes demonstrate that GPDOCK has an accuracy of 94.3% in predicting binding pose. What is more, it can accurately realize the docking of metalloproteins with ligand when one or two water molecules are engaged in the metal ion coordination. Since GPDOCK only depends on the three-dimensional structure of metalloprotein and ligand, structure-based machine learning model is employed for the scoring of binding poses, which significantly improves computational efficiency. The proposed docking strategy can be an effective and efficient tool for drug design and further study of binding mechanism of metalloproteins. The manual of GPDOCK and the code for the logistical regression model used to re-rank the docking results are available at https://github.com/wangkai-zhku/GPDOCK.git.
Collapse
Affiliation(s)
- Kai Wang
- School of Agriculture and Biology, Zhongkai University of Agriculture and Engineering, Guangzhou 510225, P. R. China.,Abinitio Technology Company, Ltd, Guangzhou 510640, P. R. China
| |
Collapse
|