1
|
Shen C, Jin X, Glover WJ, He X. Accurate Prediction of Absorption Spectral Shifts of Proteorhodopsin Using a Fragment-Based Quantum Mechanical Method. Molecules 2021; 26:4486. [PMID: 34361639 PMCID: PMC8347797 DOI: 10.3390/molecules26154486] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/10/2021] [Revised: 07/19/2021] [Accepted: 07/20/2021] [Indexed: 11/17/2022] Open
Abstract
Many experiments have been carried out to display different colors of Proteorhodopsin (PR) and its mutants, but the mechanism of color tuning of PR was not fully elucidated. In this study, we applied the Electrostatically Embedded Generalized Molecular Fractionation with Conjugate Caps (EE-GMFCC) method to the prediction of excitation energies of PRs. Excitation energies of 10 variants of Blue Proteorhodopsin (BPR-PR105Q) in residue 105GLN were calculated with the EE-GMFCC method at the TD-B3LYP/6-31G* level. The calculated results show good correlation with the experimental values of absorption wavelengths, although the experimental wavelength range among these systems is less than 50 nm. The ensemble-averaged electric fields along the polyene chain of retinal correlated well with EE-GMFCC calculated excitation energies for these 10 PRs, suggesting that electrostatic interactions from nearby residues are responsible for the color tuning. We also utilized the GMFCC method to decompose the excitation energy contribution per residue surrounding the chromophore. Our results show that residues ASP97 and ASP227 have the largest contribution to the absorption spectral shift of PR among the nearby residues of retinal. This work demonstrates that the EE-GMFCC method can be applied to accurately predict the absorption spectral shifts for biomacromolecules.
Collapse
Affiliation(s)
- Chenfei Shen
- Shanghai Engineering Research Center of Molecular Therapeutics and New Drug Development, School of Chemistry and Molecular Engineering, East China Normal University, Shanghai 200062, China; (C.S.); (X.J.)
| | - Xinsheng Jin
- Shanghai Engineering Research Center of Molecular Therapeutics and New Drug Development, School of Chemistry and Molecular Engineering, East China Normal University, Shanghai 200062, China; (C.S.); (X.J.)
| | - William J. Glover
- NYU Shanghai, 1555 Century Avenue, Shanghai 200122, China;
- NYU-ECNU Center for Computational Chemistry at NYU Shanghai, Shanghai 200062, China
- Department of Chemistry, New York University, New York, NY 10003, USA
| | - Xiao He
- Shanghai Engineering Research Center of Molecular Therapeutics and New Drug Development, School of Chemistry and Molecular Engineering, East China Normal University, Shanghai 200062, China; (C.S.); (X.J.)
- NYU-ECNU Center for Computational Chemistry at NYU Shanghai, Shanghai 200062, China
| |
Collapse
|
2
|
Cong Y, Huang K, Li Y, Zhong S, Zhang JZH, Duan L. Entropic effect and residue specific entropic contribution to the cooperativity in streptavidin-biotin binding. NANOSCALE 2020; 12:7134-7145. [PMID: 32191786 DOI: 10.1039/c9nr08380d] [Citation(s) in RCA: 16] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/10/2023]
Abstract
Molecular dynamics (MD) simulations were performed employing the polarized protein-specific charge (PPC) to explore the origin of the cooperativity in streptavidin-biotin systems (wild type, two single mutations and one double-mutation). The results of the experiment found that the existence of cooperativity is mainly the result of the entropic effect. In this study, the entropic contribution to the binding free energy was calculated using the recently developed interaction entropy (IE) method, and computational results are in excellent agreement with the experimental observations and are further verified by the calculation of the thermodynamic integration. Comparison of different force fields in terms of predicted binding strength ordering, cooperativity of energy and the stability of hydrogen bonding suggests that the PPC force field combined IE method is a suitable choice. In addition, the IE method enables us to obtain the residue-specific entropic contributions to the streptavidin-biotin binding affinity and identify ten hot-spot residues providing the dominant contribution to the cooperative binding. Importantly, the overall cooperativity obtained from the ten residues also comes mainly from the entropic effect in our study. The calculation of the potential of mean force shows that the unbinding of streptavidin-biotin is a multi-step process, and each step corresponds to the formation and rupture of the hydrogen bond network. And S45A mutation may increase the rigidity of the linker region, making the flap region relatively difficult to open. The present study provides significant molecular insight into the binding cooperativity of the streptavidin-biotin complex.
Collapse
Affiliation(s)
- Yalong Cong
- School of Physics and Electronics, Shandong Normal University, Jinan 250014, China.
| | | | | | | | | | | |
Collapse
|
3
|
Zhang Z, Sun T, Li L, Zhang D. Effects of PHB and SFC charge models on the side chain-side chain interactions in the simulation of β-hairpins. Chem Phys Lett 2019. [DOI: 10.1016/j.cplett.2019.136801] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
|
4
|
Zhong S, Huang K, Xiao Z, Sheng X, Li Y, Duan L. Binding Mechanism of Thrombin–Ligand Systems Investigated by a Polarized Protein-Specific Charge Force Field and Interaction Entropy Method. J Phys Chem B 2019; 123:8704-8716. [DOI: 10.1021/acs.jpcb.9b08064] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/04/2023]
Affiliation(s)
- Susu Zhong
- School of Physics and Electronics, Shandong Normal University, Jinan 250014, China
| | - Kaifang Huang
- School of Physics and Electronics, Shandong Normal University, Jinan 250014, China
| | - Zhengrong Xiao
- School of Physics and Electronics, Shandong Normal University, Jinan 250014, China
| | - Xiehuang Sheng
- College of Chemistry, Chemical Engineering and Materials Science, Shandong Normal University, Jinan 250014, China
| | - Yuchen Li
- School of Physics and Electronics, Shandong Normal University, Jinan 250014, China
| | - Lili Duan
- School of Physics and Electronics, Shandong Normal University, Jinan 250014, China
| |
Collapse
|
5
|
Li Y, Cong Y, Feng G, Zhong S, Zhang JZH, Sun H, Duan L. The impact of interior dielectric constant and entropic change on HIV-1 complex binding free energy prediction. STRUCTURAL DYNAMICS (MELVILLE, N.Y.) 2018; 5:064101. [PMID: 30868080 PMCID: PMC6404944 DOI: 10.1063/1.5058172] [Citation(s) in RCA: 31] [Impact Index Per Article: 5.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/17/2018] [Accepted: 11/21/2018] [Indexed: 06/01/2023]
Abstract
At present, the calculated binding free energy obtained using the molecular mechanics/Poisson-Boltzmann (Generalized-Born) surface area (MM/PB(GB)SA) method is overestimated due to the lack of knowledge of suitable interior dielectric constants in the simulation on the interaction of Human Immunodeficiency Virus (HIV-1) protease systems with inhibitors. Therefore, the impact of different values of the interior dielectric constant and the entropic contribution when using the MM/PB(GB)SA method to calculate the binding free energy was systemically evaluated. Our results show that the use of higher interior dielectric constants (1.4-2.0) can clearly improve the predictive accuracy of the MM/PBSA and MM/GBSA methods, and computational errors are significantly reduced by including the effects of electronic polarization and using a new highly efficient interaction entropy (IE) method to calculate the entropic contribution. The suitable range for the interior dielectric constant is 1.4-1.6 for the MM/PBSA method; within this range, the correlation coefficient fluctuates around 0.84, and the mean absolute error fluctuates around 2 kcal/mol. Similarly, an interior dielectric constant of 1.8-2.0 produces a correlation coefficient of approximately 0.76 when using the MM/GBSA method. In addition, the entropic contribution of each individual residue was further calculated using the IE method to predict hot-spot residues, and the detailed binding mechanisms underlying the interactions of the HIV-1 protease, its inhibitors, and bridging water molecules were investigated. In this study, the use of a higher interior dielectric constant and the IE method can improve the calculation accuracy of the HIV-1 system.
Collapse
Affiliation(s)
- Yuchen Li
- School of Physics and Electronics, Shandong Normal University, Jinan 250014, China
| | | | - Guoqiang Feng
- School of Physics and Electronics, Shandong Normal University, Jinan 250014, China
| | - Susu Zhong
- School of Physics and Electronics, Shandong Normal University, Jinan 250014, China
| | | | - Huiyong Sun
- Department of Medicinal Chemistry, School of Pharmacy, State Key Laboratory of Natural Medicines, China Pharmaceutical University, Nanjing 210009, China
| | - Lili Duan
- School of Physics and Electronics, Shandong Normal University, Jinan 250014, China
| |
Collapse
|
6
|
Li Y, Wang X, Ren L, Cao X, Ji C, Xia F, Zhang JZH. Electrostatic Polarization Effect on Cooperative Aggregation of Full Length Human Islet Amyloid. J Chem Inf Model 2018; 58:1587-1595. [DOI: 10.1021/acs.jcim.8b00215] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/15/2023]
Affiliation(s)
- Yang Li
- College of Information Science and Engineering, Shandong Agricultural University, Taian 271018, China
- College of Chemistry and Molecular Engineering, East China Normal University, Shanghai 200062, China
| | - Xianwei Wang
- Institute of Laser and Optoelectronic Technology, Zhejiang University of Technology, Hangzhou 310023, China
| | - Longlong Ren
- College of Mechanical and Electronic Engineering, Shandong Agricultural University, Taian 271018, China
| | - Xuecheng Cao
- College of Information Science and Engineering, Shandong Agricultural University, Taian 271018, China
| | - Changge Ji
- 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
| | - Fei Xia
- 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
| | - 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
| |
Collapse
|
7
|
Zhang D, Lazim R, Mun Yip Y. Ab initio folding of mixed-fold FSD-EY protein using formula-based polarizable hydrogen bond (PHB) charge model. Chem Phys Lett 2017. [DOI: 10.1016/j.cplett.2017.06.020] [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]
|
8
|
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
|
9
|
Collins MA. Can Systematic Molecular Fragmentation Be Applied to Direct Ab Initio Molecular Dynamics? J Phys Chem A 2016; 120:9281-9291. [DOI: 10.1021/acs.jpca.6b08739] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Michael A. Collins
- Research School of Chemistry, Australian National University, Canberra, ACT 0200, Australia
| |
Collapse
|
10
|
Large-scale molecular dynamics simulation: Effect of polarization on thrombin-ligand binding energy. Sci Rep 2016; 6:31488. [PMID: 27507430 PMCID: PMC4979035 DOI: 10.1038/srep31488] [Citation(s) in RCA: 46] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/02/2016] [Accepted: 07/19/2016] [Indexed: 01/17/2023] Open
Abstract
Molecular dynamics (MD) simulations lasting 500 ns were performed in explicit water to investigate the effect of polarization on the binding of ligands to human α-thrombin based on the standard nonpolarizable AMBER force field and the quantum-derived polarized protein-specific charge (PPC). The PPC includes the electronic polarization effect of the thrombin-ligand complex, which is absent in the standard force field. A detailed analysis and comparison of the results of the MD simulation with experimental data provided strong evidence that intra-protein, protein-ligand hydrogen bonds and the root-mean-square deviation of backbone atoms were significantly stabilized through electronic polarization. Specifically, two critical hydrogen bonds between thrombin and the ligand were broken at approximately 190 ns when AMBER force field was used and the number of intra-protein backbone hydrogen bonds was higher under PPC than under AMBER. The thrombin-ligand binding energy was computed using the molecular mechanics Poisson-Boltzmann surface area (MM/PBSA) method, and the results were consistent with the experimental value obtained using PPC. Because hydrogen bonds were unstable, it was failed to predict the binding affinity under the AMBER force field. Furthermore, the results of the present study revealed that differences in the binding free energy between AMBER and PPC almost comes from the electrostatic interaction. Thus, this study provides evidence that protein polarization is critical to accurately describe protein-ligand binding.
Collapse
|
11
|
Abstract
The specific electric field generated by a protease at its active site is considered as an important source of the catalytic power. Accurate calculation of electric field at the active site of an enzyme has both fundamental and practical importance. Measuring site-specific changes of electric field at internal sites of proteins due to, eg, mutation, has been realized by using molecular probes with CO or CN groups in the context of vibrational Stark effect. However, theoretical prediction of change in electric field inside a protein based on a conventional force field, such as AMBER or OPLS, is often inadequate. For such calculation, quantum chemical approach or quantum-based polarizable or polarized force field is highly preferable. Compared with the result from conventional force field, significant improvement is found in predicting experimentally measured mutation-induced electric field change using quantum-based methods, indicating that quantum effect such as polarization plays an important role in accurate description of electric field inside proteins. In comparison, the best theoretical prediction comes from fully quantum mechanical calculation in which both polarization and inter-residue charge transfer effects are included for accurate prediction of electrostatics in proteins.
Collapse
Affiliation(s)
- X Wang
- Center for Optics & Optoelectronics Research, College of Science, Zhejiang University of Technology, Hangzhou, Zhejiang, China; School of Chemistry and Molecular Engineering, East China Normal University, Shanghai, China
| | - X He
- School of Chemistry and Molecular Engineering, East China Normal University, Shanghai, China; NYU-ECNU Center for Computational Chemistry at NYU Shanghai, Shanghai, China.
| | - J Z H Zhang
- School of Chemistry and Molecular Engineering, East China Normal University, Shanghai, China; NYU-ECNU Center for Computational Chemistry at NYU Shanghai, Shanghai, China; New York University, New York, NY, United States.
| |
Collapse
|
12
|
Chen J, Wang J, Zhang Q, Chen K, Zhu W. Probing Origin of Binding Difference of inhibitors to MDM2 and MDMX by Polarizable Molecular Dynamics Simulation and QM/MM-GBSA Calculation. Sci Rep 2015; 5:17421. [PMID: 26616018 PMCID: PMC4663504 DOI: 10.1038/srep17421] [Citation(s) in RCA: 36] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/05/2015] [Accepted: 10/29/2015] [Indexed: 12/13/2022] Open
Abstract
Binding abilities of current inhibitors to MDMX are weaker than to MDM2. Polarizable molecular dynamics simulations (MD) followed by Quantum mechanics/molecular mechanics generalized Born surface area (QM//MM-GBSA) calculations were performed to investigate the binding difference of inhibitors to MDM2 and MDMX. The predicted binding free energies not only agree well with the experimental results, but also show that the decrease in van der Walls interactions of inhibitors with MDMX relative to MDM2 is a main factor of weaker bindings of inhibitors to MDMX. The analyses of dihedral angles based on MD trajectories suggest that the closed conformation formed by the residues M53 and Y99 in MDMX leads to a potential steric clash with inhibitors and prevents inhibitors from arriving in the deep of MDMX binding cleft, which reduces the van der Waals contacts of inhibitors with M53, V92, P95 and L98. The calculated results using the residue-based free energy decomposition method further prove that the interaction strength of inhibitors with M53, V92, P95 and L98 from MDMX are obviously reduced compared to MDM2. We expect that this study can provide significant theoretical guidance for designs of potent dual inhibitors to block the p53-MDM2/MDMX interactions.
Collapse
Affiliation(s)
- Jianzhong Chen
- School of Science, Shandong Jiaotong University, Jinan, 250014, China
| | - Jinan Wang
- Discovery and Design Center, CAS Key Laboratory of Receptor Research, Shanghai Institute of Materia Medica, Chinese Academy of Sciences, 555 Zuchongzhi Road, Shanghai, 201203, China
| | - Qinggang Zhang
- College of Physics and Electronics, Shandong Normal University, Jinan, 250014, China
| | - Kaixian Chen
- Discovery and Design Center, CAS Key Laboratory of Receptor Research, Shanghai Institute of Materia Medica, Chinese Academy of Sciences, 555 Zuchongzhi Road, Shanghai, 201203, China
| | - Weiliang Zhu
- Discovery and Design Center, CAS Key Laboratory of Receptor Research, Shanghai Institute of Materia Medica, Chinese Academy of Sciences, 555 Zuchongzhi Road, Shanghai, 201203, China
| |
Collapse
|
13
|
Cheng W, Liang Z, Wang W, Yi C, Li H, Zhang S, Zhang Q. Insight into binding modes of p53 and inhibitors to MDM2 based on molecular dynamic simulations and principal component analysis. Mol Phys 2015. [DOI: 10.1080/00268976.2015.1087598] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
|
14
|
Chen J, Wang X, Zhu T, Zhang Q, Zhang JZH. A Comparative Insight into Amprenavir Resistance of Mutations V32I, G48V, I50V, I54V, and I84V in HIV-1 Protease Based on Thermodynamic Integration and MM-PBSA Methods. J Chem Inf Model 2015; 55:1903-13. [DOI: 10.1021/acs.jcim.5b00173] [Citation(s) in RCA: 98] [Impact Index Per Article: 10.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022]
Affiliation(s)
- Jianzhong Chen
- School
of Science, Shandong Jiaotong University, Jinan 250357 China
| | - Xingyu Wang
- NYU−ECNU
Center for Computational Chemistry at NYU Shanghai, Shanghai 200062, China
| | - Tong Zhu
- State
Key Laboratory of Precision Spectroscopy, Institute of Theoretical
and Computational Science, East China Normal University, Shanghai 200062, China
| | - Qinggang Zhang
- Collage
of Physics and Electronic Science, Shandong Normal University, Jinan 250014, China
| | - John Z. H. Zhang
- NYU−ECNU
Center for Computational Chemistry at NYU Shanghai, Shanghai 200062, China
- State
Key Laboratory of Precision Spectroscopy, Institute of Theoretical
and Computational Science, East China Normal University, Shanghai 200062, China
| |
Collapse
|
15
|
Duan R, Lazim R, Zhang D. Understanding the basis of I50V-induced affinity decrease in HIV-1 protease via molecular dynamics simulations using polarized force field. J Comput Chem 2015. [PMID: 26198456 DOI: 10.1002/jcc.24020] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Abstract
Human immunodeficiency virus (HIV)-1 protease is one of the most promising drug target commonly utilized to combat Acquired Immune Deficiency Syndrome (AIDS). However, with the emergence of drug resistance arising from mutations, the efficiency of protease inhibitors (PIs) as a viable treatment for AIDS has been greatly reduced. I50V mutation as one of the most significant mutations occurring in HIV-1 protease will be investigated in this study. Molecular dynamics (MD) simulation was utilized to examine the effect of I50V mutation on the binding of two PIs namely indinavir and amprenavir to HIV-1 protease. Prior to the simulations conducted, the electron density distributions of the PI and each residue in HIV-1 protease are derived by combining quantum fragmentation approach molecular fractionation with conjugate caps and Poisson-Boltzmann solvation model based on polarized protein-specific charge scheme. The atomic charges of the binding complex are subsequently fitted using delta restrained electrostatic potential (delta-RESP) method to overcome the poor charge determination of buried atom. This way, both intraprotease polarization and the polarization between protease and the PI are incorporated into partial atomic charges. Through this study, the mutation-induced affinity variations were calculated and significant agreement between experiments and MD simulations conducted was observed for both HIV-1 protease-drug complexes. In addition, the mechanism governing the decrease in the binding affinity of PI in the presence of I50V mutation was also explored to provide insights pertaining to the design of the next generation of anti-HIV drugs.
Collapse
Affiliation(s)
- Rui Duan
- Division of Chemistry and Biological Chemistry, School of Physical and Mathematical Sciences, Nanyang Technological University, Singapore, 637371, Singapore
| | - Raudah Lazim
- Division of Chemistry and Biological Chemistry, School of Physical and Mathematical Sciences, Nanyang Technological University, Singapore, 637371, Singapore
| | - Dawei Zhang
- Division of Chemistry and Biological Chemistry, School of Physical and Mathematical Sciences, Nanyang Technological University, Singapore, 637371, Singapore
| |
Collapse
|
16
|
Collins MA. Molecular forces, geometries, and frequencies by systematic molecular fragmentation including embedded charges. J Chem Phys 2015; 141:094108. [PMID: 25194365 DOI: 10.1063/1.4894185] [Citation(s) in RCA: 37] [Impact Index Per Article: 4.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The accuracy of energies, energy gradients, and hessians evaluated by systematic molecular fragmentation is examined for a wide range of neutral molecules, zwitterions, and ions. A protocol is established that may employ embedded charges in conjunction with fragmentation to provide accurate evaluation of minimum energy geometries and vibrational frequencies in an automated procedure.
Collapse
Affiliation(s)
- Michael A Collins
- Research School of Chemistry, Australian National University, Canberra, ACT, Australia
| |
Collapse
|
17
|
Zhu T, Zhang JZH, He X. Correction of erroneously packed protein's side chains in the NMR structure based on ab initio chemical shift calculations. Phys Chem Chem Phys 2015; 16:18163-9. [PMID: 25052367 DOI: 10.1039/c4cp02553a] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
In this work, protein side chain (1)H chemical shifts are used as probes to detect and correct side-chain packing errors in protein's NMR structures through structural refinement. By applying the automated fragmentation quantum mechanics/molecular mechanics (AF-QM/MM) method for ab initio calculation of chemical shifts, incorrect side chain packing was detected in the NMR structures of the Pin1 WW domain. The NMR structure is then refined by using molecular dynamics simulation and the polarized protein-specific charge (PPC) model. The computationally refined structure of the Pin1 WW domain is in excellent agreement with the corresponding X-ray structure. In particular, the use of the PPC model yields a more accurate structure than that using the standard (nonpolarizable) force field. For comparison, some of the widely used empirical models for chemical shift calculations are unable to correctly describe the relationship between the particular proton chemical shift and protein structures. The AF-QM/MM method can be used as a powerful tool for protein NMR structure validation and structural flaw detection.
Collapse
Affiliation(s)
- Tong Zhu
- State Key Laboratory of Precision Spectroscopy, Institute of Theoretical and Computational Science, East China Normal University, Shanghai, 200062, China.
| | | | | |
Collapse
|
18
|
Chen J, Wang J, Zhang Q, Chen K, Zhu W. A comparative study of trypsin specificity based on QM/MM molecular dynamics simulation and QM/MM GBSA calculation. J Biomol Struct Dyn 2015; 33:2606-18. [PMID: 25562613 DOI: 10.1080/07391102.2014.1003146] [Citation(s) in RCA: 23] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/30/2022]
Abstract
Hydrogen bonding and polar interactions play a key role in identification of protein-inhibitor binding specificity. Quantum mechanics/molecular mechanics molecular dynamics (QM/MM MD) simulations combined with DFT and semi-empirical Hamiltonian (AM1d, RM1, PM3, and PM6) methods were performed to study the hydrogen bonding and polar interactions of two inhibitors BEN and BEN1 with trypsin. The results show that the accuracy of treating the hydrogen bonding and polar interactions using QM/MM MD simulation of PM6 can reach the one obtained by the DFT QM/MM MD simulation. Quantum mechanics/molecular mechanics generalized Born surface area (QM/MM-GBSA) method was applied to calculate binding affinities of inhibitors to trypsin and the results suggest that the accuracy of binding affinity prediction can be significantly affected by the accurate treatment of the hydrogen bonding and polar interactions. In addition, the calculated results also reveal the binding specificity of trypsin: (1) the amidinium groups of two inhibitors generate favorable salt bridge interaction with Asp189 and form hydrogen bonding interactions with Ser190 and Gly214, (2) the phenyl of inhibitors can produce favorable van der Waals interactions with the residues His58, Cys191, Gln192, Trp211, Gly212, and Cys215. This systematic and comparative study can provide guidance for the choice of QM/MM MD methods and the designs of new potent inhibitors targeting trypsin.
Collapse
Affiliation(s)
- Jianzhong Chen
- a School of Science , Shandong Jiaotong University , Jinan , 250014 , China
| | - Jinan Wang
- b Discovery and Design Center , CAS Key Laboratory of Receptor Research, Shanghai Institute of Materia Medica, Chinese Academy of Sciences , 555 Zuchongzhi Road, Shanghai , 201203 , China
| | - Qinggang Zhang
- c College of Physics and Electronics , Shandong Normal University , Jinan , 250014 , China
| | - Kaixian Chen
- b Discovery and Design Center , CAS Key Laboratory of Receptor Research, Shanghai Institute of Materia Medica, Chinese Academy of Sciences , 555 Zuchongzhi Road, Shanghai , 201203 , China
| | - Weiliang Zhu
- b Discovery and Design Center , CAS Key Laboratory of Receptor Research, Shanghai Institute of Materia Medica, Chinese Academy of Sciences , 555 Zuchongzhi Road, Shanghai , 201203 , China
| |
Collapse
|
19
|
Duan LL, Gao Y, Ji CG, Mei Y, Zhang QG, Tang B, Zhang JZH. Energetics of protein backbone hydrogen bonds and their local electrostatic environment. Sci China Chem 2014. [DOI: 10.1007/s11426-014-5246-0] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/24/2022]
|
20
|
Chen J, Liang Z, Wang W, Yi C, Zhang S, Zhang Q. Revealing origin of decrease in potency of darunavir and amprenavir against HIV-2 relative to HIV-1 protease by molecular dynamics simulations. Sci Rep 2014; 4:6872. [PMID: 25362963 PMCID: PMC4217091 DOI: 10.1038/srep06872] [Citation(s) in RCA: 63] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/28/2014] [Accepted: 10/13/2014] [Indexed: 12/30/2022] Open
Abstract
Clinical inhibitors Darunavir (DRV) and Amprenavir (APV) are less effective on HIV-2 protease (PR2) than on HIV-1 protease (PR1). To identify molecular basis associated with the lower inhibition, molecular dynamics (MD) simulations and molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) calculations were performed to investigate the effectiveness of the PR1 inhibitors DRV and APV against PR1/PR2. The rank of predicted binding free energies agrees with the experimental determined one. Moreover, our results show that two inhibitors bind less strongly to PR2 than to PR1, again in agreement with the experimental findings. The decrease in binding free energies for PR2 relative to PR1 is found to arise from the reduction of the van der Waals interactions induced by the structural adjustment of the triple mutant V32I, I47V and V82I. This result is further supported by the difference between the van der Waals interactions of inhibitors with each residue in PR2 and in PR1. The results from the principle component analysis suggest that inhibitor binding tends to make the flaps of PR2 close and the one of PR1 open. We expect that this study can theoretically provide significant guidance and dynamics information for the design of potent dual inhibitors targeting PR1/PR2.
Collapse
Affiliation(s)
- Jianzhong Chen
- School of Science, Shandong Jiaotong University, Jinan 250357, China
| | - Zhiqiang Liang
- School of Science, Shandong Jiaotong University, Jinan 250357, China
| | - Wei Wang
- School of Science, Shandong Jiaotong University, Jinan 250357, China
| | - Changhong Yi
- School of Science, Shandong Jiaotong University, Jinan 250357, China
| | - Shaolong Zhang
- College of Physics and Electronics, Shandong Normal University, Jinan, 250014, China
| | - Qinggang Zhang
- College of Physics and Electronics, Shandong Normal University, Jinan, 250014, China
| |
Collapse
|
21
|
Cui Q, Elstner M. Density functional tight binding: values of semi-empirical methods in an ab initio era. Phys Chem Chem Phys 2014; 16:14368-77. [PMID: 24850383 PMCID: PMC4836871 DOI: 10.1039/c4cp00908h] [Citation(s) in RCA: 102] [Impact Index Per Article: 10.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
Abstract
Semi-empirical (SE) methods are derived from Hartree-Fock (HF) or Density Functional Theory (DFT) by neglect and approximation of electronic integrals. Thereby, parameters are introduced which have to be determined from reference calculations and/or by fitting to available experimental data. This leads to computational methods that are about 2-3 orders of magnitude faster than the standard HF/DFT methods using medium sized basis sets while being about 3 orders of magnitude slower than empirical force field methods (Molecular Mechanics: MM). Therefore, SE methods are most appropriate for a specific range of applications. These include the study of systems that contain a large number of atoms and therefore being too large for ab initio or DFT methods and also problems where dynamic or entropic effects are particularly important. In the latter case, the errors made by considering a very limited number of molecular structures or neglecting entropic contributions can be much larger than the accuracy lost due to the use of SE methods. Another area where SE methods are attractive concerns the analysis of systems for which reliable MM models are not readily available. Therefore, even in an era when rapid progress is being made in ab initio methods, there is considerable interest in further developing SE methods. We illustrate this point by focusing on the discussion of recent development and application of the Density Functional Tight Binding method.
Collapse
|
22
|
Duan LL, Mei Y, Zhang QG, Tang B, Zhang JZH. Protein's native structure is dynamically stabilized by electronic polarization. JOURNAL OF THEORETICAL & COMPUTATIONAL CHEMISTRY 2014. [DOI: 10.1142/s0219633614400057] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/18/2022]
Abstract
In this paper, molecular dynamics (MD) simulations were performed for a number of benchmark proteins using both the standard assisted model building with energy refinement (AMBER) charge and the dynamically adjusted polarized protein-specific charge (DPPC) from quantum fragment calculations to provide accurate electrostatic interactions. Our result shows that proteins' dynamic structures drifted away from the native structures in simulations under standard (nonpolarizable) AMBER force field. For comparison, proteins' native structures were dynamically stable after a long time simulation under DPPC. The free energy landscape reveals that the native structure is the lowest energy conformation under DPPC, while it is not under standard AMBER charge. To further investigate the polarization effect on the stability of native structures of proteins, we restarted from some decoy structures generated from simulations using standard AMBER charges and then carried out further MD simulation using DPPC to refine those structures. Our study shows that the native structures from these decoy structures can be mostly recovered using DPPC and that the dynamic structures with the highest population in cluster analysis are in close agreement with the corresponding native structures. The current study demonstrates the importance of electronic polarization of protein in stabilizing the native structure.
Collapse
Affiliation(s)
- Li L. Duan
- College of Physics and Electronics, Shandong Normal University, Jinan 250014, P. R. China
- College of Chemistry, Chemical Engineering and Materials Science, Synergetic Innovation Center of Chemical, Imaging Functionalized Probes, Key Laboratory of Molecular and Nano Probes, Ministry of Education, Shandong Normal University, Jinan 250014, P. R. China
| | - Ye Mei
- 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
| | - Qing G. Zhang
- College of Physics and Electronics, Shandong Normal University, Jinan 250014, P. R. China
| | - Bo Tang
- College of Chemistry, Chemical Engineering and Materials Science, Synergetic Innovation Center of Chemical, Imaging Functionalized Probes, Key Laboratory of Molecular and Nano Probes, Ministry of Education, Shandong Normal University, Jinan 250014, P. R. China
| | - John Z. H. Zhang
- 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
|
23
|
Wei C, Lazim R, Zhang D. Importance of polarization effect in the study of metalloproteins: application of polarized protein specific charge scheme in predicting the reduction potential of azurin. Proteins 2014; 82:2209-19. [PMID: 24753270 DOI: 10.1002/prot.24584] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/04/2013] [Revised: 03/07/2014] [Accepted: 04/12/2014] [Indexed: 11/08/2022]
Abstract
Molecular dynamics (MD) simulation is commonly used in the study of protein dynamics, and in recent years, the extension of MD simulation to the study of metalloproteins is gaining much interest. Choice of force field is crucial in MD studies, and the inclusion of metal centers complicates the process of accurately describing the electrostatic environment that surrounds the redox centre. Herein, we would like to explore the importance of including electrostatic contribution from both protein and solvent in the study of metalloproteins. MD simulations with the implementation of thermodynamic integration will be conducted to model the reduction process of azurin from Pseudomonas aeruginosa. Three charge schemes will be used to derive the partial charges of azurin. These charge schemes differ in terms of the amount of immediate environment, respective to copper, considered during charge fitting, which ranges from the inclusion of copper and residues in the first coordination sphere during density functional theory charge fitting to the comprehensive inclusion of protein and solvent effect surrounding the metal centre using polarized protein-specific charge scheme. From the simulations conducted, the relative reduction potential of the mutated azurins respective to that of wild-type azurin (ΔEcal) were calculated and compared with experimental values. The ΔEcal approached experimental value with increasing consideration of environmental effect hence substantiating the importance of polarization effect in the study of metalloproteins. This study also attests the practicality of polarized protein-specific charge as a computational tool capable of incorporating both protein environment and solvent effect into MD simulations.
Collapse
Affiliation(s)
- Caiyi Wei
- Division of Chemistry and Biological Chemistry, School of Physical and Mathematical Sciences, Nanyang Technological University, Singapore
| | | | | |
Collapse
|
24
|
Duan LL, Zhu T, Zhang QG, Tang B, Zhang JZH. Electronic polarization stabilizes tertiary structure prediction of HP-36. J Mol Model 2014; 20:2195. [PMID: 24715046 PMCID: PMC3996369 DOI: 10.1007/s00894-014-2195-7] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/04/2013] [Accepted: 03/02/2014] [Indexed: 01/10/2023]
Abstract
Molecular dynamic (MD) simulations with both implicit and explicit solvent models have been carried out to study the folding dynamics of HP-36 protein. Starting from the extended conformation, the secondary structure of all three helices in HP-36 was formed in about 50 ns and remained stable in the remaining simulation. However, the formation of the tertiary structure was difficult. Although some intermediates were close to the native structure, the overall conformation was not stable. Further analysis revealed that the large structure fluctuation of loop and hydrophobic core regions was devoted mostly to the instability of the structure during MD simulation. The backbone root-mean-square deviation (RMSD) of the loop and hydrophobic core regions showed strong correlation with the backbone RMSD of the whole protein. The free energy landscape indicated that the distribution of main chain torsions in loop and turn regions was far away from the native state. Starting from an intermediate structure extracted from the initial AMBER simulation, HP-36 was found to generally fold to the native state under the dynamically adjusted polarized protein-specific charge (DPPC) simulation, while the peptide did not fold into the native structure when AMBER force filed was used. The two best folded structures were extracted and taken into further simulations in water employing AMBER03 charge and DPPC for 25 ns. Result showed that introducing polarization effect into interacting potential could stabilize the near-native protein structure.
Collapse
Affiliation(s)
- Li L Duan
- College of Physics and Electronics, Shandong Normal University, Jinan, 250014, China
| | | | | | | | | |
Collapse
|
25
|
Coupling between inter-helical hydrogen bonding and water dynamics in a proton transporter. J Struct Biol 2014; 186:95-111. [DOI: 10.1016/j.jsb.2014.02.010] [Citation(s) in RCA: 30] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/08/2013] [Revised: 02/13/2014] [Accepted: 02/15/2014] [Indexed: 12/20/2022]
|
26
|
Chapman DE, Steck JK, Nerenberg PS. Optimizing Protein–Protein van der Waals Interactions for the AMBER ff9x/ff12 Force Field. J Chem Theory Comput 2013; 10:273-81. [DOI: 10.1021/ct400610x] [Citation(s) in RCA: 23] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/14/2023]
Affiliation(s)
- Dail E. Chapman
- W. M. Keck Science Department, Claremont McKenna, Pitzer, and Scripps Colleges, Claremont, California 91711-5916, United States
| | - Jonathan K. Steck
- W. M. Keck Science Department, Claremont McKenna, Pitzer, and Scripps Colleges, Claremont, California 91711-5916, United States
| | - Paul S. Nerenberg
- W. M. Keck Science Department, Claremont McKenna, Pitzer, and Scripps Colleges, Claremont, California 91711-5916, United States
| |
Collapse
|
27
|
Reid DM, Collins MA. Molecular electrostatic potentials by systematic molecular fragmentation. J Chem Phys 2013; 139:184117. [DOI: 10.1063/1.4827020] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
28
|
Zeng J, Jia X, Zhang JZH, Mei Y. The F130L mutation in streptavidin reduces its binding affinity to biotin through electronic polarization effect. J Comput Chem 2013; 34:2677-86. [PMID: 24000160 DOI: 10.1002/jcc.23421] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/20/2013] [Revised: 07/30/2013] [Accepted: 08/04/2013] [Indexed: 11/06/2022]
Abstract
Recently, Baugh et al. discovered that a distal point mutation (F130L) in streptavidin causes no distinct variation to the structure of the binding pocket but a 1000-fold reduction in biotin binding affinity. In this work, we carry out molecular dynamics simulations and apply an end-state free energy method to calculate the binding free energies of biotin to wild type streptavidin and its F130L mutant. The absolute binding affinities based on AMBER charge are repulsive, and the mutation induced binding loss is underestimated. When using the polarized protein-specific charge, the absolute binding affinities are significantly enhanced. In particular, both the absolute and relative binding affinities are in line with the experimental measurements. Further investigation indicates that polarization effect is indispensable in both the generation of structural ensembles and the calculation of interaction energies. This work verifies Baugh's conjecture that electrostatic polarization effect plays an essential role in modulating the binding affinity of biotin to the streptavidin through F130L mutation.
Collapse
Affiliation(s)
- Juan Zeng
- State Key Laboratory of Precision Spectroscopy, Department of Physics and Institute of Theoretical and Computational Science, East China Normal University, Shanghai, 200062, China
| | | | | | | |
Collapse
|
29
|
Lazim R, Wei C, Sun T, Zhang D. Ab initio folding of extended α-helix: a theoretical study about the role of electrostatic polarization in the folding of helical structures. Proteins 2013; 81:1610-20. [PMID: 23670702 DOI: 10.1002/prot.24319] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/27/2012] [Revised: 03/22/2013] [Accepted: 04/19/2013] [Indexed: 11/06/2022]
Abstract
In this work, we report the ab initio folding of three different extended helical peptides namely 2khk, N36, and C34 through conventional molecular dynamics simulation at room temperature using implicit solvation model. Employing adaptive hydrogen bond specific charge (AHBC) scheme to account for the polarization effect of hydrogen bonds established during the simulation, the effective folding of the three extended helices were observed with best backbone RMSDs in comparison to the experimental structures over the helical region determined to be 1.30 Å for 2khk, 0.73 Å for N36 and 0.72 Å for C34. In this study, 2khk will be used as a benchmark case serving as a means to compare the ability of polarized (AHBC) and nonpolarized force field in the folding of an extended helix. Analyses conducted revealed the ability of the AHBC scheme in effectively folding the extended helix by promoting helix growth through the stabilization of backbone hydrogen bonds upon formation during the folding process. Similar observations were also noted when AHBC scheme was employed during the folding of C34 and N36. However, under Amber03 force field, helical structures formed during the folding of 2khk was not accompanied by stabilization thus highlighting the importance of electrostatic polarization in the folding of helical structures.
Collapse
Affiliation(s)
- Raudah Lazim
- Division of Chemistry and Biological Chemistry, School of Physical and Mathematical Sciences, Nanyang Technological University, Singapore, 637371, Singapore
| | | | | | | |
Collapse
|
30
|
Sun T, Wei C, Neo NWC, Zhang D. Misfolding of a polyalanine variant due to lack of electrostatic polarization effects. Theor Chem Acc 2013. [DOI: 10.1007/s00214-013-1354-8] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
|
31
|
Zhu T, Zhang JZH, He X. Automated Fragmentation QM/MM Calculation of Amide Proton Chemical Shifts in Proteins with Explicit Solvent Model. J Chem Theory Comput 2013; 9:2104-14. [PMID: 26583557 DOI: 10.1021/ct300999w] [Citation(s) in RCA: 64] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Abstract
We have performed a density functional theory (DFT) calculation of the amide proton NMR chemical shift in proteins using a recently developed automated fragmentation quantum mechanics/molecular mechanics (AF-QM/MM) approach. Systematic investigation was carried out to examine the influence of explicit solvent molecules, cooperative hydrogen bonding effects, density functionals, size of the basis sets, and the local geometry of proteins on calculated chemical shifts. Our result demonstrates that the predicted amide proton ((1)HN) NMR chemical shift in explicit solvent shows remarkable improvement over that calculated with the implicit solvation model. The cooperative hydrogen bonding effect is also shown to improve the accuracy of (1)HN chemical shifts. Furthermore, we found that the OPBE exchange-correlation functional is the best density functional for the prediction of protein (1)HN chemical shifts among a selective set of DFT methods (namely, B3LYP, B3PW91, M062X, M06L, mPW1PW91, OB98, OPBE), and the locally dense basis set of 6-311++G**/4-31G* is shown to be sufficient for (1)HN chemical shift calculation. By taking ensemble averaging into account, (1)HN chemical shifts calculated by the AF-QM/MM approach can be used to validate the performance of various force fields. Our study underscores that the electronic polarization of protein is of critical importance to stabilizing hydrogen bonding, and the AF-QM/MM method is able to describe the local chemical environment in proteins more accurately than most widely used empirical models.
Collapse
Affiliation(s)
- Tong Zhu
- State Key Laboratory of Precision Spectroscopy and Department of Physics, Institute of Theoretical and Computational Science, East China Normal University, Shanghai, China 200062
| | - John Z H Zhang
- State Key Laboratory of Precision Spectroscopy and Department of Physics, Institute of Theoretical and Computational Science, East China Normal University, Shanghai, China 200062.,Department of Chemistry, New York University, New York, New York 10003, United States
| | - Xiao He
- State Key Laboratory of Precision Spectroscopy and Department of Physics, Institute of Theoretical and Computational Science, East China Normal University, Shanghai, China 200062
| |
Collapse
|
32
|
Hartono YD, Mun YY, Zhang D. Adsorption and folding dynamics of MPER of HIV-1 gp41 in the presence of dpc micelle. Proteins 2013; 81:933-44. [DOI: 10.1002/prot.24256] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/18/2012] [Revised: 12/07/2012] [Accepted: 01/04/2013] [Indexed: 01/20/2023]
|
33
|
Abstract
Molecular dynamics simulations of biomolecules have matured into powerful tools of structural biology. In addition to the commonly used empirical force field potentials, quantum mechanical descriptions are gaining popularity for structure optimization and dynamic simulations of peptides and proteins. In this chapter, we introduce methodological developments such as the QM/MM framework and linear-scaling QM that make efficient calculations on large biomolecules possible. We identify the most common scenarios in which quantum descriptions of peptides and proteins are employed, such as structural refinement, force field development, treatment of unusual residues, and predicting spectroscopic and exited state properties. The benefits and shortcomings of QM potentials, in comparison to classical force fields, are discussed, with special emphasis on the sampling problems of protein conformational space. Finally, recent examples of QM/MM calculations in light-sensitive membrane proteins illustrate typical applications of the reviewed methods.
Collapse
Affiliation(s)
- Thomas Steinbrecher
- Institute of Physical Chemistry, Karlsruhe Institute of Technology, Karlsruhe, Germany
| | | |
Collapse
|
34
|
Yao XX, Ji CG, Xie DQ, Zhang JZ. Molecular dynamics study of DNA binding by INT-DBD under a polarized force field. J Comput Chem 2013; 34:1136-42. [DOI: 10.1002/jcc.23244] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/30/2012] [Revised: 11/08/2012] [Accepted: 12/30/2012] [Indexed: 11/06/2022]
|
35
|
Li Y, Ji C, Xu W, Zhang JZH. Dynamical stability and assembly cooperativity of β-sheet amyloid oligomers--effect of polarization. J Phys Chem B 2012; 116:13368-73. [PMID: 23101885 DOI: 10.1021/jp3086599] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/26/2023]
Abstract
The soluble intermediate oligomers of amyloidogenic proteins are suspected to be more cytotoxic than the mature fibrils in neurodegenerative disorders. Here, the dynamic stability and assembly cooperativity of a model oligomer of human islet amyloid polypeptide (hIAPP) segments were explored by means of all-atom molecular dynamics (MD) simulations under different force fields including AMBER99SB, OPLS, and polarized protein-specific charge (PPC) model. Simulation results show that the dynamic stability of β-sheet oligomers is seriously impacted by electrostatic polarization. Without inclusion of polarization (simulation under standard AMBER and OPLS force field), the β-sheet oligomers are dynamically unstable during MD simulation. For comparison, simulation results under PPC give significantly more stable dynamical structures of the oligomers. Furthermore, calculation of electrostatic interaction energy between the neighboring β strands with an approximate polarizable method produces energetic evidence for cooperative assembly of β-strand oligomers. This result supports a picture of downhill-like cooperative assembly of β strands during fibrillation process. The present study demonstrates the critical role of polarization in dynamic stability and assembly cooperativity of β-sheet-rich amyloid oligomers.
Collapse
Affiliation(s)
- Yang Li
- State Key Laboratory of Precision Spectroscopy and Department of Physics, Institute of Theoretical and Computational Science, East China Normal University, Shanghai, China
| | | | | | | |
Collapse
|
36
|
Gao Y, Guo M, Mei Y, Zhang JZ. Protein–water hydrogen bonds are stabilized by electrostatic polarization. Mol Phys 2012. [DOI: 10.1080/00268976.2012.668965] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/28/2022]
|
37
|
Ji CG, Zhang JZH. Effect of interprotein polarization on protein-protein binding energy. J Comput Chem 2012; 33:1416-20. [PMID: 22495971 DOI: 10.1002/jcc.22969] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/17/2012] [Accepted: 03/02/2012] [Indexed: 12/26/2022]
Abstract
Molecular dynamics simulation in explicit water for the binding of the benchmark barnase-barstar complex was carried out to investigate the effect polarization of interprotein hydrogen bonds on its binding free energy. Our study is based on the AMBER force field but with polarized atomic charges derived from fragment quantum mechanical calculation for the protein complex. The quantum-derived atomic charges include the effect of polarization of interprotein hydrogen bonds, which was absent in the standard force fields that were used in previous theoretical calculations of barnase-barstar binding energy. This study shows that this polarization effect impacts both the static (electronic) and dynamic interprotein electrostatic interactions and significantly lowers the free energy of the barnase-barstar complex.
Collapse
Affiliation(s)
- Chang G Ji
- State Key Laboratory of Precision Spectroscopy, Department of Physics, Institute of Theoretical and Computational Science, East China Normal University, Shanghai 200062, China
| | | |
Collapse
|
38
|
Mei Y, Li YL, Zeng J, Zhang JZH. Electrostatic polarization is critical for the strong binding in streptavidin-biotin system. J Comput Chem 2012; 33:1374-82. [PMID: 22467070 DOI: 10.1002/jcc.22970] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/28/2011] [Revised: 01/29/2012] [Accepted: 03/02/2012] [Indexed: 11/10/2022]
Abstract
The origin of strong affinity of biotin and its analogs binding to (strept)avidin is still the subject of an ongoing controversy. In this work, thermodynamic integration is carried out to study of the difference of binding free energies between biotin and iminobiotin to streptavidin. Three atomic charge schemes are implemented and compared. One is the traditional AMBER charge, and the other two, termed the polarized protein-specific charge, are based on a linear scaling quantum mechanical method and a continuous solvation model and have polarization effect partially or fully included. The result indicates that when nonpolarized AMBER force field is applied, the result is much underestimated. When electronic polarization is gradually included, the difference of binding affinity increases along with it. Using the linear-response approximation to eliminate the error in self-charging process, the corrected binding affinity agrees well with the experimental observation. This study is direct evidence indicating that polarization effect is critical for the strong binding in streptavidin-biotin system.
Collapse
Affiliation(s)
- Ye Mei
- State Key Laboratory of Precision Spectroscopy, Department of Physics, Institute of Theoretical and Computational Science, East China Normal University, Shanghai 200062, China.
| | | | | | | |
Collapse
|
39
|
Frähmcke JS, Wanko M, Elstner M. Building a model of the blue cone pigment based on the wild type rhodopsin structure with QM/MM methods. J Phys Chem B 2012; 116:3313-21. [PMID: 22332756 DOI: 10.1021/jp2086472] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Understanding the mechanism of color tuning of the retinal chromophore by its host protein became one of the key issues in the research on rhodopsins. While early mutation studies addressed its genetic origin, recent studies advanced to investigate its structural origin, based on X-ray crystallographic structures. For the human cone pigments, no crystal structures have been produced, and homology models were employed to elucidate the origin of its blue-shifted absorption. In this theoretical study, we take a different route to establish a structural model for human blue. Starting from the well-resolved structure of bovine rhodopsin, we derive multiple mutant models by stepwise mutation and equilibration using molecular dynamics simulations in a hybrid quantum mechanics/molecular mechanics framework. Our 30fold mutant reproduces the experimental UV-vis absorption shift of 0.45 eV and provides new insights about both structural and genetic factors that affect the excitation energy. Electrostatic effects of individual amino acids and collaborative structural effects are analyzed using semiempirical (OM2/MRCI) and ab initio (SORCI) multireference approaches.
Collapse
Affiliation(s)
- Jan S Frähmcke
- Institute for Physical and Theoretical Chemistry, TU Braunschweig, Hans-Sommer-Strasse 10, D-38106 Braunschweig, Germany
| | | | | |
Collapse
|
40
|
Duan LL, Gao Y, Mei Y, Zhang QG, Tang B, Zhang JZH. Folding of a helix is critically stabilized by polarization of backbone hydrogen bonds: study in explicit water. J Phys Chem B 2012; 116:3430-5. [PMID: 22369598 DOI: 10.1021/jp212516g] [Citation(s) in RCA: 25] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/05/2023]
Abstract
Multiple single-trajectory molecular dynamics (MD) simulation at room temperature (300 K) in explicit water was carried out to study the folding dynamics of an α-helix (PDB 2I9M ) using a polarized charge scheme that includes electronic polarization of backbone hydrogen bonds. Starting from an extended conformation, the 17-residue peptide was successfully folded into the native structure (α-helix) between 80 and 130 ns with a root-mean-square deviation of ~1.0 Å. Analysis of the time-dependent trajectories revealed that helix formation of the peptide started at the terminals and progressed toward the center of the peptide. For comparison, MD trajectories generated under various versions of standard AMBER force fields failed to show any significant or stable helix formation in our simulation. Our result shows clear evidence that the electronic polarization of backbone hydrogen bonds energetically stabilizes the helix formation and is critical to the stable folding of the short helix structure.
Collapse
Affiliation(s)
- Li L Duan
- State Key Laboratory of Precision Spectroscopy, Department of Physics, and Institute of Theoretical and Computational Science, East China Normal University, Shanghai 200062, China
| | | | | | | | | | | |
Collapse
|
41
|
Collins MA. Systematic fragmentation of large molecules by annihilation. Phys Chem Chem Phys 2012; 14:7744-51. [DOI: 10.1039/c2cp23832b] [Citation(s) in RCA: 72] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
42
|
Gao Y, Lu X, Duan LL, Zhang JZH, Mei Y. Polarization of intraprotein hydrogen bond is critical to thermal stability of short helix. J Phys Chem B 2011; 116:549-54. [PMID: 22126129 DOI: 10.1021/jp208953x] [Citation(s) in RCA: 36] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/06/2023]
Abstract
Simulation result for protein folding/unfolding is highly dependent on the accuracy of the force field employed. Even for the simplest structure of protein such as a short helix, simulations using the existing force fields often fail to produce the correct structural/thermodynamic properties of the protein. Recent research indicated that lack of polarization is at least partially responsible for the failure to successfully fold a short helix. In this work, we develop a simple formula-based atomic charge polarization model for intraprotein (backbone) hydrogen bonding based on the existing AMBER force field to study the thermal stability of a short helix (2I9M) by replica exchange molecular dynamics simulation. By comparison of the simulation results with those obtained by employing the standard AMBER03 force field, the formula-based atomic charge polarization model gave the helix melting curve in close agreement with the NMR experiment. However, in simulations using the standard AMBER force field, the helix was thermally unstable at the temperature of the NMR experiment, with a melting temperature almost below the freezing point. The difference in observed thermal stability from these two simulations is the effect of backbone intraprotein polarization, which was included in the formula-based atomic charge polarization model. The polarization of backbone hydrogen bonding thus plays a critical role in the thermal stability of helix or more general protein structures.
Collapse
Affiliation(s)
- Ya Gao
- State Key Laboratory of Precision Spectroscopy and Department of Physics, Institute of Theoretical and Computational Science, East China Normal University, Shanghai, China
| | | | | | | | | |
Collapse
|
43
|
Folding of EK peptide and its dependence on salt concentration and pH: A computational study. Sci China Chem 2011. [DOI: 10.1007/s11426-011-4399-3] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/15/2022]
|
44
|
Gordon MS, Fedorov DG, Pruitt SR, Slipchenko LV. Fragmentation Methods: A Route to Accurate Calculations on Large Systems. Chem Rev 2011; 112:632-72. [DOI: 10.1021/cr200093j] [Citation(s) in RCA: 836] [Impact Index Per Article: 64.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Affiliation(s)
- Mark S. Gordon
- Department of Chemistry and Ames Laboratory, Iowa State University, Ames Iowa 50011, United States
| | - Dmitri G. Fedorov
- Nanosystem Research Institute, National Institute of Advanced Industrial Science and Technology (AIST), 1-1-1 Umezono, Tsukuba, Ibaraki 305-8568, Japan
| | - Spencer R. Pruitt
- Department of Chemistry and Ames Laboratory, Iowa State University, Ames Iowa 50011, United States
| | - Lyudmila V. Slipchenko
- Department of Chemistry, Purdue University, West Lafayette, Indiana 47907, United States
| |
Collapse
|
45
|
Li YL, Mei Y, Zhang DW, Xie DQ, Zhang JZH. Structure and dynamics of a dizinc metalloprotein: effect of charge transfer and polarization. J Phys Chem B 2011; 115:10154-62. [PMID: 21766867 DOI: 10.1021/jp203505v] [Citation(s) in RCA: 22] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Structures and dynamics of a recently designed dizinc metalloprotein (DFsc) (J. Mol. Biol. 2003, 334, 1101) are studied by molecular dynamics simulation using a dynamically adapted polarized force field derived from fragment quantum calculation for protein in solvent. To properly describe the effect of charge transfer and polarization in the present approach, quantum chemistry calculation of the zinc-binding group is periodically performed (on-the-fly) to update the atomic charges of the zinc-binding group during the MD simulation. Comparison of the present result with those obtained from simulations under standard AMBER force field reveals that charge transfer and polarization are critical to maintaining the correct asymmetric metal coordination in the DFsc. Detailed analysis of the result also shows that dynamic fluctuation of the zinc-binding group facilitates solvent interaction with the zinc ions. In particular, the dynamic fluctuation of the zinc-zinc distance is shown to be an important feature of the catalytic function of the di-ion zinc-binding group. Our study demonstrates that the dynamically adapted polarization approach is computationally practical and can be used to study other metalloprotein systems.
Collapse
Affiliation(s)
- Yong L Li
- Key Laboratory of Mesoscopic Chemistry of MOE, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210093, China
| | | | | | | | | |
Collapse
|
46
|
Wei C, Tung D, Yip YM, Mei Y, Zhang D. Communication: The electrostatic polarization is essential to differentiate the helical propensity in polyalanine mutants. J Chem Phys 2011; 134:171101. [DOI: 10.1063/1.3581888] [Citation(s) in RCA: 21] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
47
|
|
48
|
Duan LL, Mei Y, Zhang D, Zhang QG, Zhang JZH. Folding of a helix at room temperature is critically aided by electrostatic polarization of intraprotein hydrogen bonds. J Am Chem Soc 2010; 132:11159-64. [PMID: 20698682 DOI: 10.1021/ja102735g] [Citation(s) in RCA: 81] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
We report direct folding of a 17-residue helix protein (pdb:2I9M) by standard molecular dynamics simulation (single trajectory) at room temperature with implicit solvent. Starting from a fully extended structure, 2I9M successfully folds into the native conformation within 16 ns using adaptive hydrogen bond-specific charges to take into account the electrostatic polarization effect. Cluster analysis shows that conformations in the native state cluster have the highest population (78.4%) among all sampled conformations. Folding snapshots and the secondary structure analysis demonstrate that the folding of 2I9M begins at terminals and progresses toward the center. A plot of the free energy landscape indicates that there is no significant free energy barrier during folding, which explains the observed fast folding speed. For comparison, exactly the same molecular dynamics simulation but carried out under existing AMBER charges failed to fold 2I9M into native-like structures. The current study demonstrates that electrostatic polarization of intraprotein hydrogen bonding, which stabilizes the helix, is critical to the successful folding of 2I9m.
Collapse
Affiliation(s)
- Li L Duan
- College of Physics and Electronics, Shandong Normal University, Jinan 250014, China
| | | | | | | | | |
Collapse
|
49
|
Lu Y, Mei Y, Zhang JZH, Zhang D. Communications: Electron polarization critically stabilizes the Mg2+ complex in the catalytic core domain of HIV-1 integrase. J Chem Phys 2010; 132:131101. [PMID: 20387913 DOI: 10.1063/1.3360769] [Citation(s) in RCA: 25] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
In this paper, we present a detailed dynamics study of the catalytic core domain (CCD) of HIV-1 integrase using both polarized and nonpolarized force fields. The numerical results reveal the critical role of protein polarization in stabilizing Mg(2+) coordination complex in CCD. Specifically, when nonpolarized force field is used, a remarkable drift of the Mg(2+) complex away from its equilibrium position is observed, which causes the binding site blocked by the Mg(2+) complex. In contrast, when polarized force field is employed in MD simulation, HIV-1 integrase CCD structure is stabilized and both the position of the Mg(2+) complex and the binding site are well preserved. The detailed analysis shows the transition of alpha-helix to 3(10)-helix adjacent to the catalytic loop (residues 139-147), which correlates with the dislocation of the Mg(2+) complex. The current study demonstrates the importance of electronic polarization of protein in stabilizing the metal complex in the catalytic core domain of HIV-1 integrase.
Collapse
Affiliation(s)
- Yunpeng Lu
- Division of Chemistry and Biological Chemistry, School of Physical and Mathematical Sciences, Nanyang Technological University, Singapore 637371, Singapore
| | | | | | | |
Collapse
|
50
|
|