Center for Computational Chemistry at Shanghai

Shanghai, China

Center for Computational Chemistry at Shanghai

Shanghai, China
Time filter
Source Type

Li M.,East China Normal University | Zhang J.Z.H.,East China Normal University | Zhang J.Z.H.,Center for Computational Chemistry at Shanghai | Zhang J.Z.H.,New York University
Physical Chemistry Chemical Physics | Year: 2017

The development of polarizable water models at coarse-grained (CG) levels is of much importance to CG molecular dynamics simulations of large biomolecular systems. In this work, we combined the newly developed two-bead multipole force field (TMFF) for proteins with the two-bead polarizable water models to carry out CG molecular dynamics simulations for benchmark proteins. In our simulations, two different two-bead polarizable water models are employed, the RTPW model representing five water molecules by Riniker et al. and the LTPW model representing four water molecules. The LTPW model is developed in this study based on the Martini three-bead polarizable water model. Our simulation results showed that the combination of TMFF with the LTPW model significantly stabilizes the protein's native structure in CG simulations, while the use of the RTPW model gives better agreement with all-atom simulations in predicting the residue-level fluctuation dynamics. Overall, the TMFF coupled with the two-bead polarizable water models enables one to perform an efficient and reliable CG dynamics study of the structural and functional properties of large biomolecules. © the Owner Societies 2017.

Liu J.,China Pharmaceutical University | He X.,East China Normal University | He X.,Center for Computational Chemistry at Shanghai
Physical Chemistry Chemical Physics | Year: 2017

Accurate prediction of physicochemical properties of ionic liquids (ILs) is of great significance to understand and design novel ILs with unique properties. This study employed the electrostatically embedded generalized molecular fractionation (EE-GMF) method for accurate energy calculation of IL clusters. The accuracy and efficiency of the EE-GMF method are systematically assessed at different ab initio levels (including HF, DFT and MP2) with diverse basis sets. With the fixed charge model for the embedding field, the deviations of the EE-GMF approach from conventional full system calculations are within 2.58 kcal mol-1 for all IL clusters with up to 30 ion pairs (720 atoms), tested in this study. Moreover, this linear-scaling fragment quantum mechanical (QM) method can significantly reduce the total computational cost for post-HF methods. The EE-GMF approach is well-suited for studying the energetic, structural and dynamical properties of ILs using high-level ab initio theories. © 2017 the Owner Societies.

Dupuis R.,DIPC | Benoit M.,CEMES CNRS UPR8011 | Tuckerman M.E.,New York University | Tuckerman M.E.,Courant Institute of Mathematical Sciences | And 2 more authors.
Accounts of Chemical Research | Year: 2017

ConspectusEquilibrium fractionation of stable isotopes is critically important in fields ranging from chemistry, including medicinal chemistry, electrochemistry, geochemistry, and nuclear chemistry, to environmental science. The dearth of reliable estimates of equilibrium fractionation factors, from experiment or from natural observations, has created a need for accurate computational approaches. Because isotope fractionation is a purely quantum mechanical phenomenon, exact calculation of fractionation factors is nontrivial. Consequently, a severe approximation is often made, in which it is assumed that the system can be decomposed into a set of independent harmonic oscillators. Reliance on this often crude approximation is one of the primary reasons that theoretical prediction of isotope fractionation has lagged behind experiment. A class of problems for which one might expect the harmonic approximation to perform most poorly is the isotopic fractionation between solid and solution phases.In order to illustrate the errors associated with the harmonic approximation, we have considered the fractionation of Li isotopes between aqueous solution and phyllosilicate minerals, where we find that the harmonic approximation overestimates isotope fractionation factors by as much as 30% at 25 °C. Lithium is a particularly interesting species to examine, as natural lithium isotope signatures provide information about hydrothermal processes, carbon cycle, and regulation of the Earth's climate by continental alteration. Further, separation of lithium isotopes is of growing interest in the nuclear industry due to a need for pure 6Li and 7Li isotopes. Moving beyond the harmonic approximation entails performing exact quantum calculations, which can be achieved using the Feynman path integral formulation of quantum statistical mechanics. In the path integral approach, a system of quantum particles is represented as a set of classical-like ring-polymer chains, whose interparticle interactions are determined by the rules of quantum mechanics. Because a classical isomorphism exists between the true quantum system and the system of ring-polymers, classical-like methods can be applied. Recent developments of efficient path integral approaches for the exact calculation of isotope fractionation now allow the case of the aforementioned dissolved Li fractionation properties to be studied in detail. Applying this technique, we find that the calculations yield results that are in good agreement with both experimental data and natural observations. Importantly, path integral methods, being fully atomistic, allow us to identify the origins of anharmonic effects and to make reliable predictions at temperatures that are experimentally inaccessible yet are, nevertheless, relevant for natural phenomena. © 2017 American Chemical Society.

Zhou J.,Sun Yat Sen University | Li M.,Sun Yat Sen University | Chen N.,Sun Yat Sen University | Wang S.,New York University | And 4 more authors.
ACS Chemical Biology | Year: 2015

Development of isoform-selective histone deacetylase (HDAC) inhibitors is of great biological and medical interest. Among 11 zinc-dependent HDAC isoforms, it is particularly challenging to achieve isoform inhibition selectivity between HDAC1 and HDAC2 due to their very high structural similarities. In this work, by developing and applying a novel de novo reaction-mechanism-based inhibitor design strategy to exploit the reactivity difference, we have discovered the first HDAC2-selective inhibitor, β-hydroxymethyl chalcone. Our bioassay experiments show that this new compound has a unique time-dependent selective inhibition on HDAC2, leading to about 20-fold isoform-selectivity against HDAC1. Furthermore, our ab initio QM/MM molecular dynamics simulations, a state-of-the-art approach to study reactions in biological systems, have elucidated how the β-hydroxymethyl chalcone can achieve the distinct time-dependent inhibition toward HDAC2. © 2014 American Chemical Society.

Margul D.T.,New York University | Tuckerman M.E.,New York University | Tuckerman M.E.,Courant Institute of Mathematical Sciences | Tuckerman M.E.,Center for Computational Chemistry at Shanghai
Journal of Chemical Theory and Computation | Year: 2016

Molecular dynamics remains one of the most widely used computational tools in the theoretical molecular sciences to sample an equilibrium ensemble distribution and/or to study the dynamical properties of a system. The efficiency of a molecular dynamics calculation is limited by the size of the time step that can be employed, which is dictated by the highest frequencies in the system. However, many properties of interest are connected to low-frequency, long time-scale phenomena, requiring many small time steps to capture. This ubiquitous problem can be ameliorated by employing multiple time-step algorithms, which assign different time steps to forces acting on different time scales. In such a scheme, fast forces are evaluated more frequently than slow forces, and as the former are often computationally much cheaper to evaluate, the savings can be significant. Standard multiple time-step approaches are limited, however, by resonance phenomena, wherein motion on the fastest time scales limits the step sizes that can be chosen for the slower time scales. In atomistic models of biomolecular systems, for example, the largest time step is typically limited to around 5 fs. Previously, we introduced an isokinetic extended phase-space algorithm (Minary et al. Phys. Rev. Lett. 2004, 93, 150201) and its stochastic analog (Leimkuhler et al. Mol. Phys. 2013, 111, 3579) that eliminate resonance phenomena through a set of kinetic energy constraints. In simulations of a fixed-charge flexible model of liquid water, for example, the time step that could be assigned to the slow forces approached 100 fs. In this paper, we develop a stochastic isokinetic algorithm for multiple time-step molecular dynamics calculations using a polarizable model based on fluctuating dipoles. The scheme developed here employs two sets of induced dipole moments, specifically, those associated with short-range interactions and those associated with a full set of interactions. The scheme is demonstrated on the polarizable AMOEBA water model. As was seen with fixed-charge models, we are able to obtain large time steps exceeding 100 fs, allowing calculations to be performed 10 to 20 times faster than standard thermostated molecular dynamics. © 2016 American Chemical Society.

Xiao X.,East China Normal University | Zhu T.,East China Normal University | Ji C.G.,East China Normal University | Ji C.G.,Center for Computational Chemistry at Shanghai | And 2 more authors.
Journal of Physical Chemistry B | Year: 2013

An effective polarizable bond (EPB) model has been developed for computer simulation of proteins. In this partial polarizable approach, all polar groups of amino acids are treated as polarizable, and the relevant polarizable parameters were determined by fitting to quantum calculated electrostatic properties of these polar groups. Extensive numerical tests on a diverse set of proteins (including 1IEP, 1MWE, 1NLJ, 4COX, 1PGB, 1K4C, 1MHN, 1UBQ, 1IGD) showed that this EPB model is robust in MD simulation and can correctly describe the structure and dynamics of proteins (both soluble and membrane proteins). Comparison of the computed hydrogen bond properties and dynamics of proteins with experimental data and with results obtained from the nonpolarizable force field clearly demonstrated that EPB can produce results in much better agreement with experiment. The averaged deviation of the simulated backbone N-H order parameter of the B3 immunoglobobulin-binding domain of streptococcal protein G from experimental observation is 0.0811 and 0.0332 for Amber99SB and EPB, respectively. This new model inherited the effective character of the classic force field and the fluctuating feature of previous polarizable models. Different from other polarizable models, the polarization cost energy is implicitly included in the present method. As a result, the present method avoids the problem of over polarization and is numerically stable and efficient for dynamics simulation. Finally, compared to the traditional fixed AMBER charge model, the present method only adds about 5% additional computational time and is therefore highly efficient for practical applications. © 2013 American Chemical Society.

Chen J.,Shandong Jiaotong University | Wang X.,Center for Computational Chemistry at Shanghai | Zhu T.,East China Normal University | Zhang Q.,Shandong Normal University | And 2 more authors.
Journal of Chemical Information and Modeling | Year: 2015

Drug resistance of mutations V32I, G48V, I50V, I54V, and I84V in HIV-1 protease (PR) was found in clinical treatment of HIV patients with the drug amprenavir (APV). In order to elucidate the molecular mechanism of drug resistance associated with these mutations, the thermodynamic integration (TI) and molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) methods were applied to calculate binding free energies of APV to wild-type PR and these mutated PRs. The relative binding free energy differences from the TI calculations reveal that the decrease in van der Waals interactions of APV with mutated PRs relative to the wild-type PR mainly drives the drug resistance. This result is in good agreement with the previous experimental results and is also consistent with the results from MM-PBSA calculations. Analyses based on molecular dynamics trajectories show that these mutations can adjust the shape and conformation of the binding pocket, which provides main contributions to the decrease in the van der Waals interactions of APV with mutated PRs. The present study could provide important guidance for the design of new potent inhibitors that could alleviate drug resistance of PR due to mutations. © 2015 American Chemical Society.

Schneider E.,New York University | Vogt L.,New York University | Tuckerman M.E.,New York University | Tuckerman M.E.,Courant Institute of Mathematical Sciences | Tuckerman M.E.,Center for Computational Chemistry at Shanghai
Acta Crystallographica Section B: Structural Science, Crystal Engineering and Materials | Year: 2016

Prediction and exploration of possible polymorphism in organic crystal compounds are of great importance for industries ranging from organic electronics to pharmaceuticals to high-energy materials. Here we apply our crystal structure prediction procedure and the enhanced molecular dynamics based sampling approach called the Crystal-Adiabatic Free Energy Dynamics (Crystal-AFED) method to benzene and naphthalene. Crystal-AFED allows the free energy landscape of structures to be explored efficiently at any desired temperature and pressure. For each system, we successfully predict the most stable crystal structures at atmospheric pressure and explore the relative Gibbs free energies of predicted polymorphs at high pressures. Using Crystal-AFED sampling, we find that mixed structures, which typically cannot be discovered by standard crystal structure prediction methods, are prevalent in the solid forms of these compounds at high pressure.The crystal structures of benzene and naphthalene are successfully predicted for atmospheric and high-pressure conditions. Using enhanced molecular dynamics based sampling, we find that mixed structures, which typically cannot be discovered by standard crystal structure prediction methods, are prevalent in the solid forms of these compounds at high pressure. © International Union of Crystallography, 2016.

Xu M.,New York University | Jimenez-Ruiz M.,Laue Langevin Institute | Johnson M.R.,Laue Langevin Institute | Rols S.,Laue Langevin Institute | And 7 more authors.
Physical Review Letters | Year: 2014

We report an inelastic neutron scattering (INS) study of a H2 molecule encapsulated inside the fullerene C60 which confirms the recently predicted selection rule, the first to be established for the INS spectroscopy of aperiodic, discrete molecular compounds. Several transitions from the ground state of para-H2 to certain excited translation-rotation states, forbidden according to the selection rule, are systematically absent from the INS spectra, thus validating the selection rule with a high degree of confidence. Its confirmation sets a precedent, as it runs counter to the widely held view that the INS spectroscopy of molecular compounds is not subject to any selection rules. © 2014 American Physical Society.

Lei J.,Nanjing University | Zhou Y.,Nanjing University | Xie D.,Nanjing University | Xie D.,Anhui University of Science and Technology | And 2 more authors.
Journal of the American Chemical Society | Year: 2015

Aspirin, one of the oldest and most common anti-inflammatory agents, has recently been shown to reduce cancer risks. The principal pharmacological effects of aspirin are known to arise from its covalent modification of cyclooxygenase-2 (COX-2) through acetylation of Ser530, but the detailed mechanism of its biochemical action and specificity remains to be elucidated. In this work, we have filled this gap by employing a state-of-the-art computational approach, Born-Oppenheimer molecular dynamics simulations with ab initio quantum mechanical/molecular mechanical potential and umbrella sampling. Our studies have characterized a substrate-assisted inhibition mechanism for aspirin acetylating COX: it proceeds in two successive stages with a metastable tetrahedral intermediate, in which the carboxyl group of aspirin serves as the general base. The computational results confirmed that aspirin would be 10-100 times more potent against COX-1 than against COX-2, and revealed that this inhibition specificity between the two COX isoforms can be attributed mainly to the difference in kinetics rate of the covalent inhibition reaction, not the aspirin-binding step. The structural origin of this differential inhibition of the COX enzymes by aspirin has also been elucidated. © 2014 American Chemical Society.

Loading Center for Computational Chemistry at Shanghai collaborators
Loading Center for Computational Chemistry at Shanghai collaborators