Hujun Shen

Fast and accurate computation schemes for evaluating vibrational entropy of proteins

Standard normal mode analysis (NMA) method is able to calculate vibrational entropy of proteins, but it is computationally intensive, especially for large proteins. To evaluate vibrational entropy efficiently and accurately, we, here, propose computation schemes based on coarse‐grained NMA methods. This can be achieved by rescaling coarse‐grained results with a specific factor that is derived on the basis of the linear correlation of protein vibrational entropy between standard NMA and coarse‐grained NMA. Our coarse‐grained NMA computation schemes can repeat correctly and efficiently the results of standard NMA for large proteins.

Anisotropic Coarse-Grained Model for Proteins Based On Gay−Berne and Electric Multipole Potentials

Gay–Berne anisotropic potential has been widely used to evaluate the nonbonded interactions between coarse-grained particles being described as elliptical rigid bodies. In this paper, we are presenting a coarse-grained model for twenty kinds of amino acids and proteins, based on the anisotropic Gay–Berne and point electric multipole (EMP) potentials. We demonstrate that the anisotropic coarse-grained model, namely GBEMP model, is able to reproduce many key features observed from experimental protein structures (Dunbrack Library), as well as from atomistic force field simulations (using AMOEBA, AMBER, and CHARMM force fields), while saving the computational cost by a factor of about 10–200 depending on specific cases and atomistic models. More importantly, unlike other coarse-grained approaches, our framework is based on the fundamental intermolecular forces with explicit treatment of electrostatic and repulsion-dispersion forces. As a result, the coarse-grained protein model presented an accurate description of nonbonded interactions (particularly electrostatic component) between hetero/homodimers (such as peptide–peptide, peptide–water). In addition, the encouraging performance of the model was reflected by the excellent correlation between GBEMP and AMOEBA models in the calculations of the dipole moment of peptides. In brief, the GBEMP model given here is general and transferable, suitable for simulating complex biomolecular systems.

Exploring the proton conductance and drug resistance of BM2 channel through molecular dynamics simulations and free energy calculations at different pH conditions.

BM2 channel plays an important role in the replication of influenza virus B. However, few studies attempt to investigate the mechanism of the proton conductance in BM2 channel, as well as the drug resistance of the BM2 channel. The first experimental structure of the BM2 protein channel has recently been solved, enabling us to theoretically study BM2 systems with different protonation states of histidine. By performing molecular dynamics simulations on the BM2 systems with different protonation states of four His19 residues, we provided our understanding of the structure, dynamics, and drug resistance of the BM2 channel. In general, the results of our study and other investigations both have demonstrated that whether the BM2 channel adopts an open or a closed form depends on the protonation state of His19. Meanwhile, we discovered that a drug (amantadine) was unable to enter into the center of the BM2 channel even at a low pH condition probably due to the number of hydrophilic residues of the BM2 channel. Finally, potentials of mean force (PMF) calculations were performed for the drug binding BM2 channel, energetically explaining why the BM2 channel exhibited drug resistance to two inhibitors of the AM2 channel, amantadine and rimantadine.

Coarse-Grained Modeling of Nucleic Acids Using Anisotropic Gay-Berne and Electric Multipole Potentials.

In this work, we attempt to apply a coarse-grained (CG) model, which is based on anisotropic Gay-Berne and electric multipole (EMP) potentials, to the modeling of nucleic acids. First, a comparison has been made between the CG and atomistic models (AMBER point-charge model) in the modeling of DNA and RNA hairpin structures. The CG results have demonstrated a good quality in maintaining the nucleic acid hairpin structures, in reproducing the dynamics of backbone atoms of nucleic acids, and in describing the hydrogen-bonding interactions between nucleic acid base pairs. Second, the CG and atomistic AMBER models yield comparable results in modeling double-stranded DNA and RNA molecules. It is encouraging that our CG model is capable of reproducing many elastic features of nucleic acid base pairs in terms of the distributions of the interbase pair step parameters (such as shift, slide, tilt, and twist) and the intrabase pair parameters (such as buckle, propeller, shear, and stretch). Finally, The GBEMP model has shown a promising ability to predict the melting temperatures of DNA duplexes with different lengths.

An anisotropic coarse‐grained model based on Gay–Berne and electric multipole potentials and its application to simulate a DMPC bilayer in an implicit solvent model

In this work, we aim at optimizing the performance of the anisotropic GBEMP model, which adopts a framework by combining a Gay–Berne (GB) anisotropic potential with an electric multipole (EMP) potential, in simulating a DMPC lipid bilayer in an implicit solvent model. First, the Gay–Berne parameters were initially obtained by fitting to atomistic profiles of van der Waals interactions between homodimers of molecular fragments while EMP parameters was directly derived from the expansion of point multipoles at predefined EMP sites. Second, the GB and EMP parameters for DMPC molecule were carefully optimized to be comparable to AMBER atomistic model in the calculations of the dipole moments of DMPC monomers adopting different conformations as well as the nonbonded interactions between two DMPC molecules adopting different conformations and separated at various distances. Finally, the GB parameters for DMPC were slightly adjusted in simulating a 72 DMPC bilayer system so that our GBEMP model would be able to reproduce a few important structural properties, namely, thickness ( DHH ), area per lipid ( AL ) and volume per lipid ( VL ). Meanwhile, the atomistic and experimental results for electron density profiles and order parameters were reproduced reasonably well by the GBEMP model, demonstrating the promising feature of GBEMP model in modeling lipid systems. Finally, we have shown that current GBEMP model is more efficient by a factor of about 25 than AMBER atomistic point charge model.

Binding Modes and Interaction Mechanism Between Different Base Pairs and Methylene Blue Trihydrate: A Quantum Mechanics Study

Different quantum mechanic methods have been evaluated for the calculation of binding modes and interactions between intercalators with different DNA base pairs by comparing with the results of MP2, which is very expensive, indicating that WB97XD method under 6-311+G* basis set is able to efficiently reproduce MP2 results. We discovered that the methylene blue trihydrate intercalated into the DNA base pairs, and DNA intercalation increased the distance between DNA base pairs, depending on the types of DNA bases. According to the binding energy results, it was found that the intercalation of methylene blue trihydrate into AA-TT base pair was more favorable in the orientation of nitrogen than other directions and intercalation, and the electric charge was transferred from methylene blue trihydrate to the AA-TT base pair. The analysis of change in the charge density shows that changes often take place in the heavy atom in the middle of the system which the charge density changes most remarkable.

Mechanistic insight into the functional transition of the enzyme guanylate kinase induced by a single mutation

Dramatic functional changes of enzyme usually require scores of alterations in amino acid sequence. However, in the case of guanylate kinase (GK), the functional novelty is induced by a single (S→P) mutation, leading to the functional transition of the enzyme from a phosphoryl transfer kinase into a phosphorprotein interaction domain. Here, by using molecular dynamic (MD) and metadynamics simulations, we provide a comprehensive description of the conformational transitions of the enzyme after mutating serine to proline. Our results suggest that the serine plays a crucial role in maintaining the closed conformation of wild-type GK and the GMP recognition. On the contrary, the S→P mutant exhibits a stable open conformation and loses the ability of ligand binding, which explains its functional transition from the GK enzyme to the GK domain. Furthermore, the free energy profiles (FEPs) obtained by metadymanics clearly demonstrate that the open-closed conformational transition in WT GK is positive correlated with the process of GMP binding, indicating the GMP-induced closing motion of GK enzyme, which is not observed in the mutant. In addition, the FEPs show that the S→P mutation can also leads to the mis-recognition of GMP, explaining the vanishing of catalytic activity of the mutant.

Theoretical Elucidation of the Origin for Assembly of the DAP12 Dimer with Only One NKG2C in the Lipid Membrane

In this work, we have investigated in details the origin of the assembly of the DAP12 dimer with only one NKG2C in the activating immunoreceptor complex from thew two aspects of electronic properties and dynamic structures by performing density functional theory (DFT) calculations and molecular dynamics (MD) simulations. In the DFT calculations, we studied the aggregation ability of the NKG2C(TM) with the DAP12(TM) dimer and the DAP12(TM)-DAP12(TM)-NKG2C(TM) complex by analyzing the electrostatic potentials and frontier molecular orbitals (FMOs), and in the MD simulations we mainly investigated the dynamic structures of the DAP12(TM)-DAP12(TM)-NKG2C(TM) complex and its mutants, as well as the tetramer complex consisting of two DAP12(TM) and two NKG2C(TM) helixes without any restriction. Through the studies of the electrostatic potential, the FMOs, and the dynamic structures, we have provided reasonable explanations to some extent for the experimental observation that only one NKG2C can associate with the DAP12 homodimer. The present theoretical results are expected to give valuable information for further studying the assembly between receptors and signaling subunits.

Theoretical Studies on the Folding Mechanisms for Different DNA G-quadruplexes

The G-quadruplex DNA formed by the stack of guanines in human telomere sequence is a promising anticancer target. In this study we used the energy landscape theory to elucidate the folding mechanisms for the thrombin aptamer, Form 1 and Form 3 G-quadruplexes. The three G-quadruplexes were simulated with all-atom Gō-model. Results show that, the three G-quadruplexes fold through a two-state mechanism. In the initial stage of the folding process, the compact structures are formed. The G-quadruplexes need to form G-triplex structures on the basis of the compact structures before folding to the native states. The folding free energy barrier of Form 3 G-quadruplex is higher than thrombin aptamer and Form 1, which shows that the structure of Form 3 G-quadruplex has more stability than the other two G-quadruplexes.

Folding mechanisms of Trefoil Knot proteins studied by molecular dynamics simulations and Go-model.

Most proteins need to avoid the complex topologies when folding into the native structures, but some proteins with nontrivial topologies have been found in nature. Here we used protein unfolding simulations under high temperature and all-atom Gō-model to investigate the folding mechanisms for two trefoil knot proteins. Results show that, the contacts in β-sheet are important to the formation of knot protein, and if these contacts disappeared, the knot protein would be easy to untie. In the Gō-model simulations, the folding processes of the two knot proteins are similar. The compact structures of the two knot proteins with the native contacts in β-sheet are formed in transition state, and the intermediate state has loose C-terminal. This model also reveals the detailed folding mechanisms for the two proteins.