1
|
Oliveira MP, Gonçalves YMH, Ol Gheta SK, Rieder SR, Horta BAC, Hünenberger PH. Comparison of the United- and All-Atom Representations of (Halo)alkanes Based on Two Condensed-Phase Force Fields Optimized against the Same Experimental Data Set. J Chem Theory Comput 2022; 18:6757-6778. [PMID: 36190354 DOI: 10.1021/acs.jctc.2c00524] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The level of accuracy that can be achieved by a force field is influenced by choices made in the interaction-function representation and in the relevant simulation parameters. These choices, referred to here as functional-form variants (FFVs), include for example the model resolution, the charge-derivation procedure, the van der Waals combination rules, the cutoff distance, and the treatment of the long-range interactions. Ideally, assessing the effect of a given FFV on the intrinsic accuracy of the force-field representation requires that only the specific FFV is changed and that this change is performed at an optimal level of parametrization, a requirement that may prove extremely challenging to achieve in practice. Here, we present a first attempt at such a comparison for one specific FFV, namely the choice of a united-atom (UA) versus an all-atom (AA) resolution in a force field for saturated acyclic (halo)alkanes. Two force-field versions (UA vs AA) are optimized in an automated way using the CombiFF approach against 961 experimental values for the pure-liquid densities ρliq and vaporization enthalpies ΔHvap of 591 compounds. For the AA force field, the torsional and third-neighbor Lennard-Jones parameters are also refined based on quantum-mechanical rotational-energy profiles. The comparison between the UA and AA resolutions is also extended to properties that have not been included as parameterization targets, namely the surface-tension coefficient γ, the isothermal compressibility κT, the isobaric thermal-expansion coefficient αP, the isobaric heat capacity cP, the static relative dielectric permittivity ϵ, the self-diffusion coefficient D, the shear viscosity η, the hydration free energy ΔGwat, and the free energy of solvation ΔGche in cyclohexane. For the target properties ρliq and ΔHvap, the UA and AA resolutions reach very similar levels of accuracy after optimization. For the nine other properties, the AA representation leads to more accurate results in terms of η; comparably accurate results in terms of γ, κT, αP, ϵ, D, and ΔGche; and less accurate results in terms of cP and ΔGwat. This work also represents a first step toward the calibration of a GROMOS-compatible force field at the AA resolution.
Collapse
Affiliation(s)
- Marina P Oliveira
- Laboratorium für Physikalische Chemie, ETH Zürich, ETH-Hönggerberg, HCI, CH-8093 Zürich, Switzerland
| | - Yan M H Gonçalves
- Laboratorium für Physikalische Chemie, ETH Zürich, ETH-Hönggerberg, HCI, CH-8093 Zürich, Switzerland
| | - S Kashef Ol Gheta
- Laboratorium für Physikalische Chemie, ETH Zürich, ETH-Hönggerberg, HCI, CH-8093 Zürich, Switzerland
| | - Salomé R Rieder
- Laboratorium für Physikalische Chemie, ETH Zürich, ETH-Hönggerberg, HCI, CH-8093 Zürich, Switzerland
| | - Bruno A C Horta
- Laboratorium für Physikalische Chemie, ETH Zürich, ETH-Hönggerberg, HCI, CH-8093 Zürich, Switzerland
| | - Philippe H Hünenberger
- Laboratorium für Physikalische Chemie, ETH Zürich, ETH-Hönggerberg, HCI, CH-8093 Zürich, Switzerland
| |
Collapse
|
2
|
Kang W, Jiang F, Wu YD. Universal Implementation of a Residue-Specific Force Field Based on CMAP Potentials and Free Energy Decomposition. J Chem Theory Comput 2018; 14:4474-4486. [PMID: 29906395 DOI: 10.1021/acs.jctc.8b00285] [Citation(s) in RCA: 22] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The coupling between neighboring backbone ϕ and ψ dihedral angles (torsions) has been well appreciated in protein force field development, as in correction map (CMAP) potentials. However, although preferences of backbone torsions are significantly affected by side-chain conformation, there has been no easy way to optimize this coupling. Herein, we prove that the three-dimensional (3D) free energy hypersurface of joint (ϕ, ψ, χ1) torsions can be decomposed into three separated 2D surfaces. Thus, each of the 2D torsional surfaces can be efficiently and automatically optimized using a CMAP potential. This strategy is then used to reparameterize an AMBER force field such that the resulting χ1-dependent backbone conformational preference can agree excellently with the reference protein coil library statistics. In various validation simulations (including the folding of seven peptides/proteins, backbone dynamics of three folded proteins, and two intrinsically disordered peptides), the new RSFF2C (residue-specific force field with CMAP potentials) force field gives similar or better performance compared with RSFF2. This strategy can be used to implement our RSFF force fields into a variety of molecular dynamics packages easily.
Collapse
Affiliation(s)
- Wei Kang
- Laboratory of Computational Chemistry and Drug Design, State Key Laboratory of Chemical Oncogenomics , Peking University Shenzhen Graduate School , Shenzhen 518055 , China.,College of Chemistry and Molecular Engineering , Peking University , Beijing 100871 , China
| | - Fan Jiang
- Laboratory of Computational Chemistry and Drug Design, State Key Laboratory of Chemical Oncogenomics , Peking University Shenzhen Graduate School , Shenzhen 518055 , China
| | - Yun-Dong Wu
- Laboratory of Computational Chemistry and Drug Design, State Key Laboratory of Chemical Oncogenomics , Peking University Shenzhen Graduate School , Shenzhen 518055 , China.,College of Chemistry and Molecular Engineering , Peking University , Beijing 100871 , China
| |
Collapse
|
3
|
Duan LL, Zhu T, Li YC, Zhang QG, Zhang JZH. Effect of polarization on HIV-1protease and fluoro-substituted inhibitors binding energies by large scale molecular dynamics simulations. Sci Rep 2017; 7:42223. [PMID: 28155907 PMCID: PMC5290483 DOI: 10.1038/srep42223] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/18/2016] [Accepted: 01/06/2017] [Indexed: 01/01/2023] Open
Abstract
Molecular dynamics simulations in explicit water are carried out to study the binding of six inhibitors to HIV-1 protease (PR) for up to 700 ns using the standard AMBER force field and polarized protein-specific charge (PPC). PPC is derived from quantum mechanical calculation for protein in solution and therefore it includes electronic polarization effect. Our results show that in all six systems, the bridging water W301 drifts away from the binding pocket in AMBER simulation. However, it is very stable in all six complexes systems using PPC. Especially, intra-protease, protease-inhibitor hydrogen bonds are dynamic stabilized in MD simulation. The computed binding free energies of six complexes have a significantly linear correlation with those experiment values and the correlation coefficient is found to be 0.91 in PPC simulation. However, the result from AMBER simulation shows a weaker correlation with the correlation coefficient of −0.51 due to the lack of polarization effect. Detailed binding interactions of W301, inhibitors with PR are further analyzed and discussed. The present study provides important information to quantitative understanding the interaction mechanism of PR-inhibitor and PR-W301 and these data also emphasizes the importance of both the electronic polarization and the bridging water molecule in predicting precisely binding affinities.
Collapse
Affiliation(s)
- Li L Duan
- School of Physics and Electronics, Shandong Normal University, Jinan 250014, China
| | - T Zhu
- Department of Chemistry, East China Normal University, Shanghai 200062, China.,NYU-ECNU Center for Computational Chemistry at NYU Shanghai, Shanghai 200062, China
| | - Yu C Li
- School of Physics and Electronics, Shandong Normal University, Jinan 250014, China
| | - Qing G Zhang
- School of Physics and Electronics, Shandong Normal University, Jinan 250014, China
| | - John Z H Zhang
- Department of Chemistry, East China Normal University, Shanghai 200062, China.,NYU-ECNU Center for Computational Chemistry at NYU Shanghai, Shanghai 200062, China
| |
Collapse
|
4
|
Gao Y, Zhang C, Zhang JZH, Mei Y. Evaluation of the Coupled Two-Dimensional Main Chain Torsional Potential in Modeling Intrinsically Disordered Proteins. J Chem Inf Model 2017; 57:267-274. [PMID: 28095698 DOI: 10.1021/acs.jcim.6b00589] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023]
Abstract
Intrinsically disordered proteins (IDPs) carry out crucial biological functions in essential biological processes of life. Because of the highly dynamic and conformationally heterogeneous nature of the disordered states of IDPs, molecular dynamics simulations are becoming an indispensable tool for the investigation of the conformational ensembles and dynamic properties of IDPs. Nevertheless, there is still no consensus on the most reliable force field in molecular dynamics simulations for IDPs hitherto. In this work, the recently proposed AMBER99SB2D force field is evaluated in modeling some disordered polypeptides and proteins by checking its ability to reproduce experimental NMR data. The results highlight that when the ildn side-chain corrections are included, AMBER99SB2D-ildn exhibits reliable results that agree with experiments compared with its predecessors, the AMBER14SB, AMBER99SB, AMBER99SB-ildn, and AMBER99SB2D force fields, and that decreasing the overall magnitude of protein-protein interactions in favor of protein-water interactions is a key ingredient behind the improvement.
Collapse
Affiliation(s)
- Ya Gao
- College of Fundamental Studies, Shanghai University of Engineering Science , Shanghai 201620, China
| | - Chaomin Zhang
- College of Fundamental Studies, Shanghai University of Engineering Science , Shanghai 201620, China
| | - John Z H Zhang
- College of Chemistry and Molecular Engineering, East China Normal University , Shanghai 200062, China.,NYU-ECNU Center for Computational Chemistry at NYU Shanghai , Shanghai 200062, China.,Collaborative Innovation Center of Extreme Optics, Shanxi University , Taiyuan, Shanxi 030006, China
| | - Ye Mei
- NYU-ECNU Center for Computational Chemistry at NYU Shanghai , Shanghai 200062, China.,Collaborative Innovation Center of Extreme Optics, Shanxi University , Taiyuan, Shanxi 030006, China.,State Key Laboratory of Precision Spectroscopy, School of Physics and Materials Science, East China Normal University , Shanghai 200062, China
| |
Collapse
|
5
|
Duan L, Feng G, Wang X, Wang L, Zhang Q. Effect of electrostatic polarization and bridging water on CDK2–ligand binding affinities calculated using a highly efficient interaction entropy method. Phys Chem Chem Phys 2017; 19:10140-10152. [DOI: 10.1039/c7cp00841d] [Citation(s) in RCA: 46] [Impact Index Per Article: 6.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
A new highly efficient interaction entropy (IE) method combined with the polarized protein-specific charge (PPC) force field is employed to investigate the interaction mechanism of CDK2–ligand binding and the effect of the bridging water.
Collapse
Affiliation(s)
- Lili Duan
- School of Physics and Electronics
- Shandong Normal University
- Jinan 250014
- China
| | - Guoqiang Feng
- School of Physics and Electronics
- Shandong Normal University
- Jinan 250014
- China
| | - Xianwei Wang
- Center for Optics & Optoelectronics Research
- College of Science
- Zhejiang University of Technology
- Hangzhou 310023
- China
| | - Lizhi Wang
- School of Physics
- Ludong University
- Yantai 264025
- China
| | - Qinggang Zhang
- School of Physics and Electronics
- Shandong Normal University
- Jinan 250014
- China
| |
Collapse
|
6
|
Liu J, Zhu T, Wang X, He X, Zhang JZH. Quantum Fragment Based ab Initio Molecular Dynamics for Proteins. J Chem Theory Comput 2015; 11:5897-905. [PMID: 26642993 DOI: 10.1021/acs.jctc.5b00558] [Citation(s) in RCA: 51] [Impact Index Per Article: 5.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Developing ab initio molecular dynamics (AIMD) methods for practical application in protein dynamics is of significant interest. Due to the large size of biomolecules, applying standard quantum chemical methods to compute energies for dynamic simulation is computationally prohibitive. In this work, a fragment based ab initio molecular dynamics approach is presented for practical application in protein dynamics study. In this approach, the energy and forces of the protein are calculated by a recently developed electrostatically embedded generalized molecular fractionation with conjugate caps (EE-GMFCC) method. For simulation in explicit solvent, mechanical embedding is introduced to treat protein interaction with explicit water molecules. This AIMD approach has been applied to MD simulations of a small benchmark protein Trpcage (with 20 residues and 304 atoms) in both the gas phase and in solution. Comparison to the simulation result using the AMBER force field shows that the AIMD gives a more stable protein structure in the simulation, indicating that quantum chemical energy is more reliable. Importantly, the present fragment-based AIMD simulation captures quantum effects including electrostatic polarization and charge transfer that are missing in standard classical MD simulations. The current approach is linear-scaling, trivially parallel, and applicable to performing the AIMD simulation of proteins with a large size.
Collapse
Affiliation(s)
- Jinfeng Liu
- State Key Laboratory of Precision Spectroscopy, Institute of Theoretical and Computational Science, East China Normal University , Shanghai 200062, China
| | - Tong Zhu
- State Key Laboratory of Precision Spectroscopy, Institute of Theoretical and Computational Science, East China Normal University , Shanghai 200062, China.,NYU-ECNU Center for Computational Chemistry at NYU Shanghai, Shanghai 200062, China
| | - Xianwei Wang
- Center for Optics & Optoelectronics Research, College of Science, Zhejiang University of Technology , Hangzhou, Zhejiang 310023, China
| | - Xiao He
- State Key Laboratory of Precision Spectroscopy, Institute of Theoretical and Computational Science, East China Normal University , Shanghai 200062, China.,NYU-ECNU Center for Computational Chemistry at NYU Shanghai, Shanghai 200062, China
| | - John Z H Zhang
- State Key Laboratory of Precision Spectroscopy, Institute of Theoretical and Computational Science, East China Normal University , Shanghai 200062, China.,NYU-ECNU Center for Computational Chemistry at NYU Shanghai, Shanghai 200062, China.,Department of Chemistry, New York University , New York, New York 10003, United States
| |
Collapse
|
7
|
Correct folding of an α-helix and a β-hairpin using a polarized 2D torsional potential. Sci Rep 2015; 5:10359. [PMID: 26039188 PMCID: PMC5380191 DOI: 10.1038/srep10359] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/29/2014] [Accepted: 04/09/2015] [Indexed: 01/12/2023] Open
Abstract
A new modification to the AMBER force field that incorporates the coupled two-dimensional main chain torsion energy has been evaluated for the balanced representation of secondary structures. In this modified AMBER force field (AMBER03(2D)), the main chain torsion energy is represented by 2-dimensional Fourier expansions with parameters fitted to the potential energy surface generated by high-level quantum mechanical calculations of small peptides in solution. Molecular dynamics simulations are performed to study the folding of two model peptides adopting either α-helix or β-hairpin structures. Both peptides are successfully folded into their native structures using an AMBER03(2D) force field with the implementation of a polarization scheme (AMBER03(2D)p). For comparison, simulations using a standard AMBER03 force field with and without polarization, as well as AMBER03(2D) without polarization, fail to fold both peptides successfully. The correction to secondary structure propensity in the AMBER03 force field and the polarization effect are critical to folding Trpzip2; without these factors, a helical structure is obtained. This study strongly suggests that this new force field is capable of providing a more balanced preference for helical and extended conformations. The electrostatic polarization effect is shown to be indispensable to the growth of secondary structures.
Collapse
|
8
|
Gao Y, Li Y, Mou L, Hu W, Zheng J, Zhang JZH, Mei Y. Coupled Two-Dimensional Main-Chain Torsional Potential for Protein Dynamics II: Performance and Validation. J Phys Chem B 2015; 119:4188-93. [DOI: 10.1021/jp510215c] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/29/2023]
Affiliation(s)
- Ya Gao
- College
of Fundamental Studies, Shanghai University of Engineering Science, Shanghai 201620, China
| | - Yongxiu Li
- College
of Fundamental Studies, Shanghai University of Engineering Science, Shanghai 201620, China
- Key
Laboratory of Catalysis and Materials Science of the State Ethnic
Affairs Commission and Ministry of Education, Hubei Province, South-Central University for Nationalities, Wuhan 430074, China
| | - Lirong Mou
- Institutes
for Advanced Interdisciplinary Research, East China Normal University, Shanghai 200062, China
| | - Wenxin Hu
- Computing Center, School of Information Science & Technology, East China Normal University, Shanghai 200062, China
| | - Jun Zheng
- Computing Center, School of Information Science & Technology, East China Normal University, Shanghai 200062, China
| | - John Z. H. Zhang
- College
of Fundamental Studies, Shanghai University of Engineering Science, Shanghai 201620, China
- NYU-ECNU Center
for Computational Chemistry at NYU Shanghai, Shanghai 200062, China
| | - Ye Mei
- Center
for Laser and Computational Biophysics, State Key Laboratory of Precision
Spectroscopy, Department of Physics and Institute of Theoretical and
Computational Science, East China Normal University, Shanghai 200062, China
- NYU-ECNU Center
for Computational Chemistry at NYU Shanghai, Shanghai 200062, China
| |
Collapse
|
9
|
Mou L, Jia X, Gao Y, Li Y, Zhang JZH, Mei Y. Folding simulation of Trp-cage utilizing a new AMBER compatible force field with coupled main chain torsions. JOURNAL OF THEORETICAL & COMPUTATIONAL CHEMISTRY 2014. [DOI: 10.1142/s0219633614500266] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/18/2022]
Abstract
A newly developed AMBER compatible force field with coupled backbone torsion potential terms (AMBER032D) is utilized in a folding simulation of a mini-protein Trp-cage. Through replica exchange and direct molecular dynamics (MD) simulations, a multi-step folding mechanism with a synergetic folding of the hydrophobic core (HPC) and the α-helix in the final stage is suggested. The native structure has the lowest free energy and the melting temperature predicted from the specific heat capacity Cvis only 12 K higher than the experimental measurement. This study, together with our previous study, shows that AMBER032Dis an accurate force field that can be used for protein folding simulations.
Collapse
Affiliation(s)
- Lirong Mou
- Institute for Advanced Interdisciplinary Research, East China Normal University, Shanghai 200062, P. R. China
| | - Xiangyu Jia
- Center for Laser and Computational Biophysics, State Key Laboratory of Precision Spectroscopy and Department of Physics, Institute of Theoretical and Computational Science, East China Normal University, Shanghai 200062, P. R. China
| | - Ya Gao
- Center for Laser and Computational Biophysics, State Key Laboratory of Precision Spectroscopy and Department of Physics, Institute of Theoretical and Computational Science, East China Normal University, Shanghai 200062, P. R. China
| | - Yongxiu Li
- Center for Laser and Computational Biophysics, State Key Laboratory of Precision Spectroscopy and Department of Physics, Institute of Theoretical and Computational Science, East China Normal University, Shanghai 200062, P. R. China
| | - John Z. H. Zhang
- Center for Laser and Computational Biophysics, State Key Laboratory of Precision Spectroscopy and Department of Physics, Institute of Theoretical and Computational Science, East China Normal University, Shanghai 200062, P. R. China
- NYU-ECNU Center for Computational Chemistry at NYU Shanghai, Shanghai 200062, P. R. China
| | - Ye Mei
- Center for Laser and Computational Biophysics, State Key Laboratory of Precision Spectroscopy and Department of Physics, Institute of Theoretical and Computational Science, East China Normal University, Shanghai 200062, P. R. China
- NYU-ECNU Center for Computational Chemistry at NYU Shanghai, Shanghai 200062, P. R. China
| |
Collapse
|