1
|
Abstract
Glycoscience assembles all the scientific disciplines involved in studying various molecules and macromolecules containing carbohydrates and complex glycans. Such an ensemble involves one of the most extensive sets of molecules in quantity and occurrence since they occur in all microorganisms and higher organisms. Once the compositions and sequences of these molecules are established, the determination of their three-dimensional structural and dynamical features is a step toward understanding the molecular basis underlying their properties and functions. The range of the relevant computational methods capable of addressing such issues is anchored by the specificity of stereoelectronic effects from quantum chemistry to mesoscale modeling throughout molecular dynamics and mechanics and coarse-grained and docking calculations. The Review leads the reader through the detailed presentations of the applications of computational modeling. The illustrations cover carbohydrate-carbohydrate interactions, glycolipids, and N- and O-linked glycans, emphasizing their role in SARS-CoV-2. The presentation continues with the structure of polysaccharides in solution and solid-state and lipopolysaccharides in membranes. The full range of protein-carbohydrate interactions is presented, as exemplified by carbohydrate-active enzymes, transporters, lectins, antibodies, and glycosaminoglycan binding proteins. A final section features a list of 150 tools and databases to help address the many issues of structural glycobioinformatics.
Collapse
Affiliation(s)
- Serge Perez
- Centre de Recherche sur les Macromolecules Vegetales, University of Grenoble-Alpes, Centre National de la Recherche Scientifique, Grenoble F-38041, France
| | - Olga Makshakova
- FRC Kazan Scientific Center of Russian Academy of Sciences, Kazan Institute of Biochemistry and Biophysics, Kazan 420111, Russia
| |
Collapse
|
2
|
Taguchi M, Oyama R, Kaneso M, Hayashi S. Hybrid QM/MM Free-Energy Evaluation of Drug-Resistant Mutational Effect on the Binding of an Inhibitor Indinavir to HIV-1 Protease. J Chem Inf Model 2022; 62:1328-1344. [PMID: 35212226 DOI: 10.1021/acs.jcim.1c01193] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022]
Abstract
A human immunodeficiency virus-1 (HIV-1) protease is a homodimeric aspartic protease essential for the replication of HIV. The HIV-1 protease is a target protein in drug discovery for antiretroviral therapy, and various inhibitor molecules of transition state analogues have been developed. However, serious drug-resistant mutants have emerged. For understanding the molecular mechanism of the drug resistance, an accurate examination of the impacts of the mutations on ligand binding and enzymatic activity is necessary. Here, we present a molecular simulation study on the ligand binding of indinavir, a potent transition state analogue inhibitor, to the wild-type protein and a V82T/I84V drug-resistant mutant of the HIV-1 protease. We employed a hybrid ab initio quantum mechanical/molecular mechanical (QM/MM) free-energy optimization technique which combines a highly accurate QM description of the ligand molecule and its interaction with statistically ample conformational sampling of the MM protein environment by long-time molecular dynamics simulations. Through the free-energy calculations of protonation states of catalytic groups at the binding pocket and of the ligand-binding affinity changes upon the mutations, we successfully reproduced the experimentally observed significant reduction of the binding affinity upon the drug-resistant mutations and elucidated the underlying molecular mechanism. The present study opens the way for understanding the molecular mechanism of drug resistance through the direct quantitative comparison of ligand binding and enzymatic reaction with the same accuracy.
Collapse
Affiliation(s)
- Masahiko Taguchi
- Department of Chemistry, Graduate School of Science, Kyoto University, Kyoto 606-8502, Japan.,Institute for Quantum Life Science, National Institutes for Quantum Science and Technology, Kizugawa, Kyoto 619-0215, Japan
| | - Ryo Oyama
- Department of Chemistry, Graduate School of Science, Kyoto University, Kyoto 606-8502, Japan
| | - Masahiro Kaneso
- Department of Chemistry, Graduate School of Science, Kyoto University, Kyoto 606-8502, Japan
| | - Shigehiko Hayashi
- Department of Chemistry, Graduate School of Science, Kyoto University, Kyoto 606-8502, Japan
| |
Collapse
|
3
|
Jia X, Li P. Solvation Free Energy Calculation Using a Fixed-Charge Model: Implicit and Explicit Treatments of the Polarization Effect. J Phys Chem B 2019; 123:1139-1148. [PMID: 30628452 DOI: 10.1021/acs.jpcb.8b10479] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
In this work, IPolQ-Mod charges and the reference potential scheme are used to calculate the solvation free energies of a set of organic molecules. Both methods could capture the phase transfer of a solute with accompanying polarization cost utilizing a fixed-charge model. The IPolQ-Mod charges, which are the average of two charge sets fitted in a vacuum state and a condensed phase, take account of the polarization effect implicitly. For the reference potential method, the quantum mechanics polarization corrections are calculated explicitly by thermodynamic perturbation. The polarization effect captured by the IPolQ-Mod charges is an approximation to that of the reference potential method theoretically. In the present study, the reference potential method shows a slight improvement over the classical restrained electrostatic potential (RESP) charges, which perform pretty well in predicting the solvation free energy. However, IPolQ-Mod(MP2) shows a poor agreement with the experimental data. Compared with IPolQ-Mod(MP2), IPolQ-Mod(M06-2X) or IPolQ-Mod(ωB97X) is found to give more appropriate prediction of the molecule's dipole and the solvation free energies calculated by IPolQ-Mod(M06-2X) or IPolQ-Mod(ωB97X) are more compatible with those of the RESP charges. If the other force field parameters remain unchanged, M06-2X or ωB97X is recommended to derive the IPolQ-Mod charges.
Collapse
Affiliation(s)
- Xiangyu Jia
- NYU Shanghai , 1555 Century Avenue , Shanghai 200122 , China.,NYU-ECNU Center for Computational Chemistry at NYU Shanghai , 3663 Zhongshan Road North , Shanghai 200127 , China
| | - Pengfei Li
- State Key Laboratory of Precision Spectroscopy and Department of Physics and Institute of Theoretical and Computational Science , East China Normal University , Shanghai 200062 , China
| |
Collapse
|
4
|
Wang M, Mei Y, Ryde U. Predicting Relative Binding Affinity Using Nonequilibrium QM/MM Simulations. J Chem Theory Comput 2018; 14:6613-6622. [DOI: 10.1021/acs.jctc.8b00685] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Affiliation(s)
- Meiting Wang
- State Key Laboratory of Precision Spectroscopy, School of Physics and Materials Science, East China Normal University, Shanghai 200062, China
- Department of Theoretical Chemistry, Lund University, Chemical Centre, P.O. Box 124, SE-221 00 Lund, Sweden
| | - Ye Mei
- State Key Laboratory of Precision Spectroscopy, School of Physics and Materials Science, East China Normal University, Shanghai 200062, China
- NYU−ECNU Center for Computational Chemistry at NYU Shanghai, Shanghai 200062, China
| | - Ulf Ryde
- Department of Theoretical Chemistry, Lund University, Chemical Centre, P.O. Box 124, SE-221 00 Lund, Sweden
| |
Collapse
|
5
|
Efficient Computation of Free Energy Surfaces of Diels⁻Alder Reactions in Explicit Solvent at Ab Initio QM/MM Level. Molecules 2018; 23:molecules23102487. [PMID: 30274188 PMCID: PMC6222833 DOI: 10.3390/molecules23102487] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/05/2018] [Revised: 09/25/2018] [Accepted: 09/26/2018] [Indexed: 11/17/2022] Open
Abstract
For Diels–Alder (DA) reactions in solution, an accurate and converged free energy (FE) surface at ab initio (ai) quantum mechanical/molecular mechanical (QM/MM) level is imperative for the understanding of reaction mechanism. However, this computation is still far too expensive. In a previous work, we proposed a new method termed MBAR+wTP, with which the computation of the ai FE profile can be accelerated by several orders of magnitude via a three-step procedure: (I) an umbrella sampling (US) using a semi-empirical (SE) QM/MM Hamiltonian is performed; (II) the FE profile is generated using the Multistate Bennett Acceptance Ratio (MBAR) analysis; and (III) a weighted Thermodynamic Perturbation (wTP) from the SE Hamiltonian to the ai Hamiltonian is performed to obtain the ai QM/MM FE profile using weight factors from the MBAR analysis. In this work, this method is extended to the calculations of two-dimensional FE surfaces of two Diels–Alder reactions of cyclopentadiene with either acrylonitrile or 1-4-naphthoquinone at ai QM/MM level. The accurate activation free energies at the ai QM/MM level, which are much closer to the experimental measurements than those calculated by other methods, indicate that this MBAR+wTP method can be applied in the studies of complex reactions in condensed phase with much-enhanced efficiency.
Collapse
|
6
|
Mishra SK, Koča J. Assessing the Performance of MM/PBSA, MM/GBSA, and QM-MM/GBSA Approaches on Protein/Carbohydrate Complexes: Effect of Implicit Solvent Models, QM Methods, and Entropic Contributions. J Phys Chem B 2018; 122:8113-8121. [PMID: 30084252 DOI: 10.1021/acs.jpcb.8b03655] [Citation(s) in RCA: 47] [Impact Index Per Article: 7.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/22/2023]
Abstract
Rapid and accurate binding affinity prediction of protein-carbohydrate complexes is a major challenge in glycomimetics design. Among the existing computational techniques, end-point methods have received considerable interest because of their low computational cost. However, significant obstacles remain when such methods are applied to protein-glycan complexes. This article reports the performance of end-point free-energy calculation methods: molecular mechanics/Poisson-Boltzmann surface area (MM/PBSA), MM/generalized Born surface area (MM/GBSA), and quantum mechanics-MM/GBSA (QM-MM/GBSA) on monosaccharides bound to RSL lectin from Ralstonia solanacearum. A careful investigation of the molecular dynamics simulation length, van der Waals radii sets, GB models, QM Hamiltonians, and entropic compensation has been made, and the results are compared with the experimental binding free energies from isothermal titration calorimetry/surface plasmon resonance measurements. The binding free energies using implicit solvent methods are found to be sensitive to the simulation length, radii set, GB model, and QM Hamiltonian. A simulation length of 10 ns using the radii set mbondi provides the best agreement with the experimental values ( r2 = 0.96) by MM/PBSA. The GBHCT model is in accord with the experimental values in MM/GBSA ( r2 = 0.91) or in combination with parameterized model number 6 (PM6) ( r2 = 0.98) in QM-MM/GBSA. Out of 12 QM Hamiltonians tested, PM6, density functional theory-based tight binding (DFTB), and their variants proved to be more efficient than other semiempirical methods. These methods perform equally well in predicting both absolute and relative binding free energies.
Collapse
Affiliation(s)
- Sushil K Mishra
- Central European Institute of Technology (CEITEC), and National Centre for Biomolecular Research, Faculty of Science , Masaryk University , Kamenice-5, 625 00 Brno , Czech Republic
| | - Jaroslav Koča
- Central European Institute of Technology (CEITEC), and National Centre for Biomolecular Research, Faculty of Science , Masaryk University , Kamenice-5, 625 00 Brno , Czech Republic
| |
Collapse
|
7
|
Li PF, Jia XY, Wang MT, Mei Y. Comparison of Accuracy and Convergence Rate between Equilibrium and Nonequilibrium Alchemical Transformations for Calculation of Relative Binding Free Energy. CHINESE J CHEM PHYS 2017. [DOI: 10.1063/1674-0068/30/cjcp1711204] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/18/2023]
Affiliation(s)
- Peng-fei Li
- State Key Laboratory of Precision Spectroscopy, School of Physics and Materials Science, East China Normal University, Shanghai 200062, China
| | - Xiang-yu Jia
- State Key Laboratory of Precision Spectroscopy, School of Physics and Materials Science, East China Normal University, Shanghai 200062, China
- NYU-ECNU Center for Computational Chemistry at NYU Shanghai, Shanghai 200062, China
| | - Mei-ting Wang
- State Key Laboratory of Precision Spectroscopy, School of Physics and Materials Science, East China Normal University, Shanghai 200062, China
| | - Ye Mei
- State Key Laboratory of Precision Spectroscopy, School of Physics and Materials Science, East China Normal University, Shanghai 200062, China
- NYU-ECNU Center for Computational Chemistry at NYU Shanghai, Shanghai 200062, China
- Department of Chemistry and Biochemistry, University of Oklahoma, Norman OK 73019, USA
| |
Collapse
|