1
|
Csizi K, Reiher M. Universal
QM
/
MM
approaches for general nanoscale applications. WIRES COMPUTATIONAL MOLECULAR SCIENCE 2023. [DOI: 10.1002/wcms.1656] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 02/04/2023]
Affiliation(s)
| | - Markus Reiher
- Laboratorium für Physikalische Chemie ETH Zürich Zürich Switzerland
| |
Collapse
|
2
|
Design, synthesis and evaluation of RGD peptidomimetic – Gold nanostar conjugates as M21 cell adhesion inhibitors. Bioorg Chem 2022; 126:105873. [DOI: 10.1016/j.bioorg.2022.105873] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/04/2021] [Revised: 04/12/2022] [Accepted: 05/10/2022] [Indexed: 11/20/2022]
|
3
|
Baldini L, Lenci E, Bianchini F, Trabocchi A. Identification of a Common Pharmacophore for Binding to MMP2 and RGD Integrin: Towards a Multitarget Approach to Inhibit Cancer Angiogenesis and Metastasis. Molecules 2022; 27:molecules27041249. [PMID: 35209039 PMCID: PMC8879803 DOI: 10.3390/molecules27041249] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/04/2022] [Revised: 02/08/2022] [Accepted: 02/10/2022] [Indexed: 12/10/2022] Open
Abstract
During tumor angiogenesis different growth factors, cytokines and other molecules interact closely with each other to facilitate tumor cell invasion and metastatic diffusion. The most intensively studied as molecular targets in anti-angiogenic therapies are vascular endothelial growth factor (VEGF) and related receptors, integrin receptors and matrix metalloproteinases (MMPs). Considering the poor efficacy of cancer angiogenesis monotherapies, we reasoned combining the inhibition of αvβ3 and MMP2 as a multitarget approach to deliver a synergistic blockade of tumor cell migration, invasion and metastasis. Accordingly, we identified a common pharmacophore in the binding cavity of MMP2 and αvβ3, demonstrating such approach with the design, synthesis and bioassays of tyrosine-derived peptidomimetics carrying the necessary functional groups to bind to key pharmacophoric elements of MMP2 and αvβ3 RGD integrin.
Collapse
Affiliation(s)
- Lorenzo Baldini
- Department of Chemistry “Ugo Schiff”, University of Florence, 50019 Sesto Fiorentino, Italy; (L.B.); (E.L.)
| | - Elena Lenci
- Department of Chemistry “Ugo Schiff”, University of Florence, 50019 Sesto Fiorentino, Italy; (L.B.); (E.L.)
| | - Francesca Bianchini
- Department of Biomedical, Experimental and Clinical Sciences “Mario Serio”, University of Florence, 50134 Florence, Italy;
| | - Andrea Trabocchi
- Department of Chemistry “Ugo Schiff”, University of Florence, 50019 Sesto Fiorentino, Italy; (L.B.); (E.L.)
- Correspondence:
| |
Collapse
|
4
|
Naidoo KJ, Bruce-Chwatt T, Senapathi T, Hillebrand M. Multidimensional Free Energy and Accelerated Quantum Library Methods Provide a Gateway to Glycoenzyme Conformational, Electronic, and Reaction Mechanisms. Acc Chem Res 2021; 54:4120-4130. [PMID: 34726899 DOI: 10.1021/acs.accounts.1c00477] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/17/2023]
Abstract
Enzyme reactions are complex to simulate accurately, and none more so than glycoenzymes (glycosyltransferase and glycosidases). A rigorous sampling of the protein frame and the conformationally plural carbohydrate substrate coupled with an unbiased treatment of the electron dynamics is needed to discover the true reaction landscapes. Here, we demonstrate the effectiveness of two computational methods ported in libraries that we have developed. The first is a flat histogram free energy method called FEARCF capable of multidimensional sampling and rapidly converging to a complete coverage of phase space. The second, the Quantum Supercharger Library (QSL), is a method that accelerates the computation of the ab initio electronic wave function as well as the integral derivatives on graphical processing units (GPUs). These QSL accelerated computations form the core components needed for direct quantum dynamics and QM/MM dynamics when coupled with legacy codes such as GAMESS and NWCHEM, making state of the art hyper-parallel electronic computations in chemistry and chemical biology possible. The combination of QSL (acceleration of ab initio QM computation) and FEARCF (multidimensional hyper-parallel reaction dynamics) makes the simulation of ab initio QM/MM reaction dynamics of enzyme catalysis feasible. Enzymes that process carbohydrates pose an added challenge as their pyranose ring substrates span multidimensional conformational space whose sampling is an intimate function of the catalytic mechanism. Here, we use the pairing of FEARCF and QSL to simulate the catalytic effect of the glycoenzyme β-N-acetylglucosamine transferase (OGT). The reaction mechanism is discovered from a variable three bond reaction surface using SCCDFTB. The role of the OGT in distorting the pyranose ring of β-N-acetylglucosamine (GlcNAc) away from the equilibrium 4C1 chair conformation toward the E3 envelope needed for the transition state is discovered from its pucker free energy hypersurfaces (or free energy volume, FEV). A complete GlcNAc ring pucker HF 6-31g FEV is constructed from ab initio QM dynamics in vacuum and ab initio QM/MM dynamics in the OGT catalytic domain. The OGT is shown to clearly lower the pathway toward the transition state E3 ring conformer as well as stabilize it by 1.63 kcal/mol. Illustrated here is the use of QSL accelerated ab initio QM/MM dynamics that thoroughly explores carbohydrate catalyzed reactions through a FEARCF multidimensional sampling of the interdependence between reaction and conformational space. This demonstrates how experimentally inaccessible molecular and electronic mechanisms that underpin enzyme catalysis can be discovered by directly modeling the dynamics of these complex reactions.
Collapse
Affiliation(s)
- Kevin J Naidoo
- Scientific Computing Research Unit and Department of Chemistry, University of Cape Town, Rondebosch 7701, South Africa
- Institute of Infectious Disease and Molecular Medicine, Faculty of Health Science, University of Cape Town, Rondebosch 7701, South Africa
| | - Tomás Bruce-Chwatt
- Scientific Computing Research Unit and Department of Chemistry, University of Cape Town, Rondebosch 7701, South Africa
| | - Tharindu Senapathi
- Scientific Computing Research Unit and Department of Chemistry, University of Cape Town, Rondebosch 7701, South Africa
| | - Malcolm Hillebrand
- Scientific Computing Research Unit and Department of Chemistry, University of Cape Town, Rondebosch 7701, South Africa
- Nonlinear Dynamics and Chaos Group, Department of Mathematics, University of Cape Town, Rondebosch 7701, South Africa
| |
Collapse
|
5
|
Modeling Catalysis in Allosteric Enzymes: Capturing Conformational Consequences. Top Catal 2021; 65:165-186. [DOI: 10.1007/s11244-021-01521-1] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/19/2022]
|
6
|
Nochebuena J, Naseem-Khan S, Cisneros GA. Development and application of quantum mechanics/molecular mechanics methods with advanced polarizable potentials. WILEY INTERDISCIPLINARY REVIEWS. COMPUTATIONAL MOLECULAR SCIENCE 2021; 11:e1515. [PMID: 34367343 PMCID: PMC8341087 DOI: 10.1002/wcms.1515] [Citation(s) in RCA: 12] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/28/2020] [Accepted: 12/19/2020] [Indexed: 01/02/2023]
Abstract
Quantum mechanics/molecular mechanics (QM/MM) simulations are a popular approach to study various features of large systems. A common application of QM/MM calculations is in the investigation of reaction mechanisms in condensed-phase and biological systems. The combination of QM and MM methods to represent a system gives rise to several challenges that need to be addressed. The increase in computational speed has allowed the expanded use of more complicated and accurate methods for both QM and MM simulations. Here, we review some approaches that address several common challenges encountered in QM/MM simulations with advanced polarizable potentials, from methods to account for boundary across covalent bonds and long-range effects, to polarization and advanced embedding potentials.
Collapse
Affiliation(s)
- Jorge Nochebuena
- Department of Chemistry, University of North Texas, Denton, Texas, USA
| | - Sehr Naseem-Khan
- Department of Chemistry, University of North Texas, Denton, Texas, USA
| | - G Andrés Cisneros
- Department of Chemistry, University of North Texas, Denton, Texas, USA
| |
Collapse
|
7
|
Scherbinina SI, Toukach PV. Three-Dimensional Structures of Carbohydrates and Where to Find Them. Int J Mol Sci 2020; 21:E7702. [PMID: 33081008 PMCID: PMC7593929 DOI: 10.3390/ijms21207702] [Citation(s) in RCA: 17] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/26/2020] [Revised: 10/15/2020] [Accepted: 10/16/2020] [Indexed: 02/06/2023] Open
Abstract
Analysis and systematization of accumulated data on carbohydrate structural diversity is a subject of great interest for structural glycobiology. Despite being a challenging task, development of computational methods for efficient treatment and management of spatial (3D) structural features of carbohydrates breaks new ground in modern glycoscience. This review is dedicated to approaches of chemo- and glyco-informatics towards 3D structural data generation, deposition and processing in regard to carbohydrates and their derivatives. Databases, molecular modeling and experimental data validation services, and structure visualization facilities developed for last five years are reviewed.
Collapse
Affiliation(s)
- Sofya I. Scherbinina
- N.D. Zelinsky Institute of Organic Chemistry, Russian Academy of Science, Leninsky prospect 47, 119991 Moscow, Russia
- Higher Chemical College, D. Mendeleev University of Chemical Technology of Russia, Miusskaya Square 9, 125047 Moscow, Russia
| | - Philip V. Toukach
- N.D. Zelinsky Institute of Organic Chemistry, Russian Academy of Science, Leninsky prospect 47, 119991 Moscow, Russia
| |
Collapse
|
8
|
Khan T, Lawrence AJ, Azad I, Raza S, Joshi S, Khan AR. Computational Drug Designing and Prediction Of Important Parameters Using in silico Methods- A Review. Curr Comput Aided Drug Des 2020; 15:384-397. [PMID: 30914032 DOI: 10.2174/1573399815666190326120006] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/21/2018] [Revised: 01/28/2019] [Accepted: 02/25/2019] [Indexed: 01/02/2023]
Abstract
BACKGROUND Computational or in silico studies are undertaken to assess the drug like properties of lead compounds. These studies help in fast prediction of relevant properties. OBJECTIVE Through this review, an effort is made to encapsulate some of the important parameters which should be met by a compound for it to be considered as a potential drug candidate along with an overview of automated softwares which can be used for making various predictions. METHODS Drug uptake, its absorption, evacuation and associated hazardous effects are important factors for consideration in drug designing and should be known in early stages of drug development. Several important physicochemical properties like molecular weight, polar surface area (PSA), molecular flexibility etc. have to be taken into consideration in drug designing. Toxicological assessment is another important aspect of drug discovery which predicts the safety and adverse effects of a drug. RESULTS Additionally, bioactivity scores of probable drug leads against various human receptors can also be predicted to evaluate the probability of them to act as a potential drug candidate. The in vivo biological targets of a molecule can also be efficiently predicted by molecular docking studies. CONCLUSION Some important software like iGEMDOCK, AutoDock, OSIRIS property explorer, Molinspiration, MetaPrint2D, admetSAR and their working methodology and principle of working have been summarized in this review.
Collapse
Affiliation(s)
- Tahmeena Khan
- Department of Chemistry, Integral University, Lucknow, 226026, U.P., India
| | - Alfred J Lawrence
- Department of Chemistry, Isabella Thoburn College, Lucknow, 226007, U.P., India
| | - Iqbal Azad
- Department of Chemistry, Integral University, Lucknow, 226026, U.P., India
| | - Saman Raza
- Department of Chemistry, Isabella Thoburn College, Lucknow, 226007, U.P., India
| | - Seema Joshi
- Department of Chemistry, Isabella Thoburn College, Lucknow, 226007, U.P., India
| | - Abdul Rahman Khan
- Department of Chemistry, Integral University, Lucknow, 226026, U.P., India
| |
Collapse
|
9
|
Innocenti R, Lenci E, Menchi G, Trabocchi A. Combination of multicomponent KA 2 and Pauson-Khand reactions: short synthesis of spirocyclic pyrrolocyclopentenones. Beilstein J Org Chem 2020; 16:200-211. [PMID: 32117477 PMCID: PMC7034245 DOI: 10.3762/bjoc.16.23] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/18/2019] [Accepted: 01/23/2020] [Indexed: 12/23/2022] Open
Abstract
The Cu-catalyzed multicomponent ketone–amine–alkyne (KA2) reaction was combined with a Pauson–Khand cycloaddition to give access of unprecedented constrained spirocyclic pyrrolocyclopentenone derivatives following a DOS couple-pair approach. The polyfunctional molecular scaffolds were tested on the cyclopentenone reactivity to further expand the skeletal diversity, demonstrating the utility of this combined approach in generating novel spiro compounds as starting material for the generation of chemical libraries. The chemoinformatics characterization of the newly-synthesized molecules gave evidence about structural and physicochemical properties with respect to a set of blockbuster drugs, and showed that such scaffolds are drug-like but more spherical and three-dimensional in character than the drugs.
Collapse
Affiliation(s)
- Riccardo Innocenti
- Department of Chemistry "Ugo Schiff", University of Florence, Via della Lastruccia 13, 50019 Sesto Fiorentino, Florence, Italy
| | - Elena Lenci
- Department of Chemistry "Ugo Schiff", University of Florence, Via della Lastruccia 13, 50019 Sesto Fiorentino, Florence, Italy
| | - Gloria Menchi
- Department of Chemistry "Ugo Schiff", University of Florence, Via della Lastruccia 13, 50019 Sesto Fiorentino, Florence, Italy
| | - Andrea Trabocchi
- Department of Chemistry "Ugo Schiff", University of Florence, Via della Lastruccia 13, 50019 Sesto Fiorentino, Florence, Italy.,Interdepartmental Center for Preclinical Development of Molecular Imaging (CISPIM), University of Florence, Viale Morgagni 85, 50134 Florence, Italy
| |
Collapse
|
10
|
Yosberto CM. Theoretical study on binding interactions of laccase-enzyme from Ganoderma weberianum with multiples ligand substrates with environmental impact. ACTA ACUST UNITED AC 2019. [DOI: 10.29328/journal.apb.1001007] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
|
11
|
Hahn DF, Milić JV, Hünenberger PH. Vase
‐
Kite
Equilibrium of Resorcin[4]arene Cavitands Investigated Using Molecular Dynamics Simulations with Ball‐and‐Stick Local Elevation Umbrella Sampling. Helv Chim Acta 2019. [DOI: 10.1002/hlca.201900060] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/06/2022]
Affiliation(s)
- David F. Hahn
- Laboratory of Physical Chemistry, Department of Chemistry and Applied BiosciencesETH Zürich Vladimir-Prelog-Weg 2 CH-8093 Zürich Switzerland
| | - Jovana V. Milić
- Laboratory of Photonics and InterfacesÉcole Polytechnique Fédérale de Lausanne, EPFL SB ISIC LPI, Station 6 CH-1015 Lausanne Switzerland
| | - Philippe H. Hünenberger
- Laboratory of Physical Chemistry, Department of Chemistry and Applied BiosciencesETH Zürich Vladimir-Prelog-Weg 2 CH-8093 Zürich Switzerland
| |
Collapse
|
12
|
Marion A, Gokcan H, Monard G. Semi-Empirical Born-Oppenheimer Molecular Dynamics (SEBOMD) within the Amber Biomolecular Package. J Chem Inf Model 2019; 59:206-214. [PMID: 30433776 DOI: 10.1021/acs.jcim.8b00605] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/05/2023]
Abstract
Semi-empirical quantum methods from the neglect of differential diatomic overlap (NDDO) family such as MNDO, AM1, or PM3 are fast albeit approximate quantum methods. By combining them with linear scaling methods like the divide & conquer (D&C) method, it is possible to quickly evaluate the energy of systems containing hundreds to thousands of atoms. We here present our implementation in the Amber biomolecular package of a SEBOMD module that provides a way to run semi-empirical Born-Oppenheimer molecular dynamics. At each step of a SEBOMD, a fully converged self-consistent field (SCF) calculation is performed to obtain the semiempirical quantum potential energy of a molecular system encaged or not in periodic boundary conditions. We describe the implementation and the features of our SEBOMD implementation. We show the requirements to conserve the total energy in NVE simulations, and how to accelerate SCF convergence through density matrix extrapolation. Specific ways of handling periodic boundary conditions using mechanical embedding or electrostatic embedding through a tailored quantum Ewald summation is developed. The parallel performance of SEBOMD simulations using the D&C scheme are presented for liquid water systems of various sizes, and a comparison between the traditional full diagonalization scheme and the D&C approach for the reproduction of the structure of liquid water illustrates the potentiality of SEBOMD to simulate molecular systems containing several hundreds of atoms for hundreds of picoseconds with a quantum mechanical potential in a reasonable amount of CPU time.
Collapse
Affiliation(s)
- Antoine Marion
- Université de Lorraine, CNRS, LPCT , F-54000 Nancy , France.,Department of Chemistry , Middle East Technical University , 06800 , Ankara , Turkey
| | - Hatice Gokcan
- Université de Lorraine, CNRS, LPCT , F-54000 Nancy , France.,Department of Chemistry , University of North Texas , Denton , Texas 76201 , United States
| | - Gerald Monard
- Université de Lorraine, CNRS, LPCT , F-54000 Nancy , France
| |
Collapse
|
13
|
Plazinski W, Plazinska A. Molecular dynamics simulations of hexopyranose ring distortion in different force fields. PURE APPL CHEM 2017. [DOI: 10.1515/pac-2016-0922] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022]
Abstract
Abstract:The four classical, biomolecular force fields designed to study hexopyranose-based carbohydrates (GROMOS 56a6CARBO/56a6CARBO_R, GROMOS 53a6GLYC, CHARMM and GLYCAM06) have been tested in the context of ring-inversion properties. These properties were evaluated for both unfunctionalized monomers of all hexopyranoses of the d series and for residues in a chain composed of uniform units connected by α(1→4) and β(1→4) glycosidic linkages. The results indicate that the tested force fields differ in their predictions of the ring-inversion properties of both monomers and residues in a chain. The comparison with the available experimental data and with the semi-empirical Angyal scheme reveals that, at the level of monomers, GROMOS 56a6CARBO, GROMOS 53a6GLYC and CHARMM correctly reproduce the ring-inversion free energies. However, due to the lack of analogous reference data we cannot state which force field is more or less accurate in the context of ring distortion of residues in a chain. Therefore, the use of ab initio potentials is recommended in the prospective, quantitative studies on the related subject.
Collapse
Affiliation(s)
- Wojciech Plazinski
- Jerzy Haber Institute of Catalysis and Surface Chemistry, Polish Academy of Sciences, ul. Niezapominajek 8, 30-239 Cracow, Poland, Tel.: +48815375685, Fax: +48815375685, e-mail:
| | - Anita Plazinska
- Department of Biopharmacy, Faculty of Pharmacy, Medical University of Lublin, W. Chodzki Str., 4a, 20-093 Lublin, Poland
| |
Collapse
|
14
|
Rogers IL, Naidoo KJ. Producing DFT/MM enzyme reaction trajectories from SCC-DFTB/MM driving forces to probe the underlying electronics of a glycosyltransferase reaction. J Comput Chem 2017; 38:1789-1798. [PMID: 28488320 DOI: 10.1002/jcc.24820] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/05/2016] [Revised: 03/14/2017] [Accepted: 04/03/2017] [Indexed: 02/05/2023]
Abstract
The SCC-DFTB/MIO/CHARMM free energy surface for a glycosyltransferase, TcTS, is benchmarked against a DFT/MM reaction trajectory using the same CHARMM MM force field ported to the NWChem package. The popular B3LYP functional, against which the MIO parameter set was parameterized is used to optimize TS structures and run DFT reaction dynamics. A novel approach was used to generate reaction forces from a SCC-DFTB/MIO/CHARMM reaction surface to drive B3LYP/6-31G/MM and B3LYP/6-31G(d)/MM reaction trajectories. Although TS structures compare favorably, differences stemming primarily from a minimal basis set approximation prevented a successful 6-31G(d) FEARCF reaction dynamics trajectory. None the less, the dynamic evolution of the B3LYP/6-31G/MM-computed electron density provided an opportunity to perform NBO analysis along the reaction trajectory. Here, we illustrate that a successful ab initio reaction trajectory is computationally accessible when the underlying potential energy function of the semi-empirical method used to produce driving forces is sufficiently close to the ab initio potential. © 2017 Wiley Periodicals, Inc.
Collapse
Affiliation(s)
- Ian L Rogers
- Scientific Computing Research Unit and Department of Chemistry, University of Cape Town, Rondebosch, 7701, South Africa
| | - Kevin J Naidoo
- Scientific Computing Research Unit and Department of Chemistry, University of Cape Town, Rondebosch, 7701, South Africa
| |
Collapse
|
15
|
Ojeda-May P, Nam K. Acceleration of Semiempirical QM/MM Methods through Message Passage Interface (MPI), Hybrid MPI/Open Multiprocessing, and Self-Consistent Field Accelerator Implementations. J Chem Theory Comput 2017. [PMID: 28628742 DOI: 10.1021/acs.jctc.7b00322] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
The strategy and implementation of scalable and efficient semiempirical (SE) QM/MM methods in CHARMM are described. The serial version of the code was first profiled to identify routines that required parallelization. Afterward, the code was parallelized and accelerated with three approaches. The first approach was the parallelization of the entire QM/MM routines, including the Fock matrix diagonalization routines, using the CHARMM message passage interface (MPI) machinery. In the second approach, two different self-consistent field (SCF) energy convergence accelerators were implemented using density and Fock matrices as targets for their extrapolations in the SCF procedure. In the third approach, the entire QM/MM and MM energy routines were accelerated by implementing the hybrid MPI/open multiprocessing (OpenMP) model in which both the task- and loop-level parallelization strategies were adopted to balance loads between different OpenMP threads. The present implementation was tested on two solvated enzyme systems (including <100 QM atoms) and an SN2 symmetric reaction in water. The MPI version exceeded existing SE QM methods in CHARMM, which include the SCC-DFTB and SQUANTUM methods, by at least 4-fold. The use of SCF convergence accelerators further accelerated the code by ∼12-35% depending on the size of the QM region and the number of CPU cores used. Although the MPI version displayed good scalability, the performance was diminished for large numbers of MPI processes due to the overhead associated with MPI communications between nodes. This issue was partially overcome by the hybrid MPI/OpenMP approach which displayed a better scalability for a larger number of CPU cores (up to 64 CPUs in the tested systems).
Collapse
Affiliation(s)
| | - Kwangho Nam
- Department of Chemistry and Biochemistry, University of Texas at Arlington , Arlington, Texas 76019-0065, United States
| |
Collapse
|
16
|
Montgomery AP, Xiao K, Wang X, Skropeta D, Yu H. Computational Glycobiology: Mechanistic Studies of Carbohydrate-Active Enzymes and Implication for Inhibitor Design. STRUCTURAL AND MECHANISTIC ENZYMOLOGY 2017; 109:25-76. [DOI: 10.1016/bs.apcsb.2017.04.003] [Citation(s) in RCA: 23] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
|
17
|
Gaweda K, Plazinski W. Pyranose ring conformations in mono- and oligosaccharides: a combined MD and DFT approach. Phys Chem Chem Phys 2017; 19:20760-20772. [DOI: 10.1039/c7cp02920a] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
Abstract
A two-step computational protocol is proposed to efficiently study the conformational properties of hexopyranoses with a special emphasis on their ring-inversion-properties. By applying it, the errors resulting from overestimating the contribution of the hydrogen bond-rich, low-energy structures that are not abundant in aqueous solutions are avoided.
Collapse
Affiliation(s)
- Karolina Gaweda
- Jerzy Haber Institute of Catalysis and Surface Chemistry
- Polish Academy of Sciences
- 30-239 Cracow
- Poland
| | - Wojciech Plazinski
- Jerzy Haber Institute of Catalysis and Surface Chemistry
- Polish Academy of Sciences
- 30-239 Cracow
- Poland
| |
Collapse
|
18
|
Janoš P, Trnka T, Kozmon S, Tvaroška I, Koča J. Different QM/MM Approaches To Elucidate Enzymatic Reactions: Case Study on ppGalNAcT2. J Chem Theory Comput 2016; 12:6062-6076. [PMID: 27787999 DOI: 10.1021/acs.jctc.6b00531] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Hybrid QM/MM computational studies can provide invaluable insight into the mechanisms of enzymatic reactions that can be exploited for rational drug design. Various approaches are available for such studies. However, their strengths and weaknesses may not be immediately apparent. Using the retaining glycosyltransferase ppGalNAcT2 as a case study, we compare different methodologies used to obtain reaction paths and transition state information. Ab Initio MD using CPMD coupled with the String Method is used to derive the minimum free energy reaction path. The geometrical features of the free energy path, especially around the transition state, agree with the minimum potential energy path obtained by the much less computationally expensive Nudged Elastic Band method. The barrier energy, however, differs by 8 kcal/mol. The free energy surface generated by metadynamics provides a rough overview of the reaction and can confirm the physical relevance of optimized paths or provide an initial guess for path optimization methods. Calculations of enzymatic reactions are usually performed at best at the DFT level of theory. A comparison of widely used functionals with high-level DLPNO-CCSD(T)/CBS data on the potential energy profile serves as a validation of the usability of DFT for this type of enzymatic reaction. The M06-2X meta-hybrid functional in particular matches the DLPNO-CCSD(T)/CBS reference extremely well with errors within 1 kcal/mol. However, even pure-GGA functional OPBE provides sufficient accuracy for this type of study.
Collapse
Affiliation(s)
| | | | - Stanislav Kozmon
- Institute of Chemistry, Slovak Academy of Sciences , 845 38 Bratislava, Slovak Republic
| | - Igor Tvaroška
- Institute of Chemistry, Slovak Academy of Sciences , 845 38 Bratislava, Slovak Republic
| | | |
Collapse
|
19
|
Pshetitsky Y, Eitan R, Verner G, Kohen A, Major DT. Improved Sugar Puckering Profiles for Nicotinamide Ribonucleoside for Hybrid QM/MM Simulations. J Chem Theory Comput 2016; 12:5179-5189. [PMID: 27490188 DOI: 10.1021/acs.jctc.6b00401] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
Abstract
The coenzyme nicotinamide adenine dinucleotide (NAD+) and its reduced form (NADH) play ubiquitous roles as oxidizing and reducing agents in nature. The binding, and possibly the chemical redox step, of NAD+/NADH may be influenced by the cofactor conformational distribution and, in particular, by the ribose puckering of its nicotinamide-ribonucleoside (NR) moiety. In many hybrid quantum mechanics-molecular mechanics (QM/MM) studies of NAD+/NADH dependent enzymes, the QM region is treated by semiempirical (SE) methods. Recent work suggests that SE methods do not adequately describe the ring puckering in sugar molecules. In the present work we adopt an efficient and practical strategy to correct for this deficiency for NAD+/NADH. We have implemented a cost-effective correction to a SE Hamiltonian by adding a correction potential, which is defined as the difference between an accurate benchmark density functional theory (DFT) potential energy surface (PES) and the SE PES. In practice, this is implemented via a B-spline interpolation scheme for the grid-based potential energy difference surface. We find that the puckering population distributions obtained from free energy QM(SE)/MM simulations are in good agreement with DFT and in fair accord with experimental results. The corrected PES should facilitate a more accurate description of the ribose puckering in the NAD+/NADH cofactor in simulations of biological systems.
Collapse
Affiliation(s)
- Yaron Pshetitsky
- Department of Chemistry and the Lise Meitner-Minerva Center of Computational Quantum Chemistry, Bar-Ilan University , Ramat-Gan 52900, Israel
| | - Reuven Eitan
- Department of Chemistry and the Lise Meitner-Minerva Center of Computational Quantum Chemistry, Bar-Ilan University , Ramat-Gan 52900, Israel
| | - Gilit Verner
- Department of Chemistry and the Lise Meitner-Minerva Center of Computational Quantum Chemistry, Bar-Ilan University , Ramat-Gan 52900, Israel
| | - Amnon Kohen
- Department of Chemistry, University of Iowa , Iowa City, Iowa 52242, United States
| | - Dan Thomas Major
- Department of Chemistry and the Lise Meitner-Minerva Center of Computational Quantum Chemistry, Bar-Ilan University , Ramat-Gan 52900, Israel
| |
Collapse
|
20
|
Christensen A, Kubař T, Cui Q, Elstner M. Semiempirical Quantum Mechanical Methods for Noncovalent Interactions for Chemical and Biochemical Applications. Chem Rev 2016; 116:5301-37. [PMID: 27074247 PMCID: PMC4867870 DOI: 10.1021/acs.chemrev.5b00584] [Citation(s) in RCA: 246] [Impact Index Per Article: 30.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/03/2015] [Indexed: 12/28/2022]
Abstract
Semiempirical (SE) methods can be derived from either Hartree-Fock or density functional theory by applying systematic approximations, leading to efficient computational schemes that are several orders of magnitude faster than ab initio calculations. Such numerical efficiency, in combination with modern computational facilities and linear scaling algorithms, allows application of SE methods to very large molecular systems with extensive conformational sampling. To reliably model the structure, dynamics, and reactivity of biological and other soft matter systems, however, good accuracy for the description of noncovalent interactions is required. In this review, we analyze popular SE approaches in terms of their ability to model noncovalent interactions, especially in the context of describing biomolecules, water solution, and organic materials. We discuss the most significant errors and proposed correction schemes, and we review their performance using standard test sets of molecular systems for quantum chemical methods and several recent applications. The general goal is to highlight both the value and limitations of SE methods and stimulate further developments that allow them to effectively complement ab initio methods in the analysis of complex molecular systems.
Collapse
Affiliation(s)
- Anders
S. Christensen
- Department
of Chemistry and Theoretical Chemistry Institute, University of Wisconsin—Madison, 1101 University Avenue, Madison, Wisconsin 53706, United States
| | - Tomáš Kubař
- Institute of Physical
Chemistry & Center for Functional Nanostructures and Institute of Physical
Chemistry, Karlsruhe Institute of Technology, Kaiserstrasse 12, 76131 Karlsruhe, Germany
| | - Qiang Cui
- Department
of Chemistry and Theoretical Chemistry Institute, University of Wisconsin—Madison, 1101 University Avenue, Madison, Wisconsin 53706, United States
| | - Marcus Elstner
- Institute of Physical
Chemistry & Center for Functional Nanostructures and Institute of Physical
Chemistry, Karlsruhe Institute of Technology, Kaiserstrasse 12, 76131 Karlsruhe, Germany
| |
Collapse
|
21
|
Dixit M, Das S, Mhashal A, Eitan R, Major D. Practical Aspects of Multiscale Classical and Quantum Simulations of Enzyme Reactions. Methods Enzymol 2016; 577:251-86. [DOI: 10.1016/bs.mie.2016.05.046] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/11/2023]
|
22
|
Plazinski W, Lonardi A, Hünenberger PH. Revision of the GROMOS 56A6(CARBO) force field: Improving the description of ring-conformational equilibria in hexopyranose-based carbohydrates chains. J Comput Chem 2015; 37:354-65. [PMID: 26525424 DOI: 10.1002/jcc.24229] [Citation(s) in RCA: 71] [Impact Index Per Article: 7.9] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/17/2015] [Revised: 10/03/2015] [Accepted: 10/06/2015] [Indexed: 01/19/2023]
Abstract
This article describes a revised version 56A6(CARBO_R) of the GROMOS 56A6(CARBO) force field for hexopyranose-based carbohydrates. The simulated properties of unfunctionalized hexopyranoses are unaltered with respect to 56A6CARBO . In the context of both O1 -alkylated hexopyranoses and oligosaccharides, the revision stabilizes the regular (4) C1 chair for α-anomers, with the opposite effect for β-anomers. As a result, spurious ring inversions observed in α(1→4)-linked chains when using the original 56A6(CARBO) force field are alleviated. The (4) C1 chair is now the most stable conformation for all d-hexopyranose residues, irrespective of the linkage type and anomery, and of the position of the residue along the chain. The methylation of a d-hexopyranose leads to a systematic shift in the ring-inversion free energy ((4) C1 to (1) C4 ) by 7-8 kJ mol(-1), positive for the α-anomers and negative for the β-anomers, which is qualitatively compatible with the expected enhancement of the anomeric effect upon methylation at O1. The ring-inversion free energies for residues within chains are typically smaller in magnitude compared to those of the monomers, and correlate rather poorly with the latter. This suggests that the crowding of ring substituents upon chain formation alters the ring flexibility in a nonsystematic fashion. In general, the description of carbohydrate chains afforded by 56A6(CARBO_R) suggests a significant extent of ring flexibility, i.e., small but often non-negligible equilibrium populations of inverted chairs, and challenges the "textbook" picture of conformationally locked carbohydrate rings.
Collapse
Affiliation(s)
- Wojciech Plazinski
- Laboratory of Physical Chemistry, ETH Hönggerberg, HCI, Zürich, CH-8093, Switzerland.,J. Haber Institute of Catalysis and Surface Chemistry, Polish Academy of Sciences, Ul. Niezapominajek 8, Cracow, 30-239, Poland
| | - Alice Lonardi
- Laboratory of Physical Chemistry, ETH Hönggerberg, HCI, Zürich, CH-8093, Switzerland
| | | |
Collapse
|
23
|
Fernandes KD, Renison CA, Naidoo KJ. Quantum supercharger library: Hyper-parallelism of the Hartree-Fock method. J Comput Chem 2015; 36:1399-409. [DOI: 10.1002/jcc.23936] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/11/2015] [Revised: 04/02/2015] [Accepted: 04/16/2015] [Indexed: 01/05/2023]
Affiliation(s)
- Kyle D. Fernandes
- Scientific Computing Research Unit and Department of Chemistry; University of Cape Town; Rondebosch 7701 South Africa
| | - C. Alicia Renison
- Scientific Computing Research Unit and Department of Chemistry; University of Cape Town; Rondebosch 7701 South Africa
| | - Kevin J. Naidoo
- Scientific Computing Research Unit and Department of Chemistry; University of Cape Town; Rondebosch 7701 South Africa
| |
Collapse
|
24
|
Renison CA, Fernandes KD, Naidoo KJ. Quantum supercharger library: Hyper-parallel integral derivatives algorithms forab initioQM/MM dynamics. J Comput Chem 2015; 36:1410-9. [DOI: 10.1002/jcc.23938] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/19/2015] [Revised: 04/15/2015] [Accepted: 04/16/2015] [Indexed: 01/23/2023]
Affiliation(s)
- C. Alicia Renison
- Scientific Computing Research Unit and Department of Chemistry; University of Cape Town; Rondebosch, Cape Town 7701 South Africa
| | - Kyle D. Fernandes
- Scientific Computing Research Unit and Department of Chemistry; University of Cape Town; Rondebosch, Cape Town 7701 South Africa
| | - Kevin J. Naidoo
- Scientific Computing Research Unit and Department of Chemistry; University of Cape Town; Rondebosch, Cape Town 7701 South Africa
| |
Collapse
|
25
|
Yilmazer ND, Korth M. Enhanced semiempirical QM methods for biomolecular interactions. Comput Struct Biotechnol J 2015; 13:169-75. [PMID: 25848495 PMCID: PMC4372622 DOI: 10.1016/j.csbj.2015.02.004] [Citation(s) in RCA: 55] [Impact Index Per Article: 6.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/16/2015] [Revised: 02/17/2015] [Accepted: 02/19/2015] [Indexed: 12/21/2022] Open
Abstract
Recent successes and failures of the application of 'enhanced' semiempirical QM (SQM) methods are reviewed in the light of the benefits and backdraws of adding dispersion (D) and hydrogen-bond (H) correction terms. We find that the accuracy of SQM-DH methods for non-covalent interactions is very often reported to be comparable to dispersion-corrected density functional theory (DFT-D), while computation times are about three orders of magnitude lower. SQM-DH methods thus open up a possibility to simulate realistically large model systems for problems both in life and materials science with comparably high accuracy.
Collapse
Affiliation(s)
| | - Martin Korth
- Institute of Theoretical Chemistry, Ulm University, D-89069 Ulm, Germany
| |
Collapse
|
26
|
Govender KK, Naidoo KJ. Evaluating AM1/d-CB1 for Chemical Glycobiology QM/MM Simulations. J Chem Theory Comput 2014; 10:4708-17. [DOI: 10.1021/ct500373p] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/20/2023]
Affiliation(s)
- Krishna K. Govender
- Scientific Computing
Research Unit and Department
of Chemistry, University of Cape Town, Rondebosch 7701, South Africa
| | - Kevin J. Naidoo
- Scientific Computing
Research Unit and Department
of Chemistry, University of Cape Town, Rondebosch 7701, South Africa
| |
Collapse
|