Abstract: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.

The Orientational Order of the Α-Helical Backbone of Poly(γ-Benzyl L-Glutamate) – a Molecular Dynamics Simulation

The Orientational Orders of Poly(β-Phenethyl L -Aspartate) in Two Opposite Α-Helical Form: a Molecular Dynamic Simulation

A Molecular Dynamics Model Of The Amorphous Regions Of Polyethylene Oxide

A Study of Molecular Orientation of Polymethylene Chains

Electronic Polarization Stabilizes Tertiary Structure Prediction of HP-36.

Development of an Effective Polarizable Bond Method for Biomolecular Simulation.

A Test of AMBER Force Fields in Predicting the Secondary Structure of Α-Helical and Β-Hairpin Peptides

Direct Folding Simulation of Helical Proteins Using an Effective Polarizable Bond Force Field

Exploration of Multi-State Conformational Dynamics and Underlying Global Functional Landscape of Maltose Binding Protein

Intrinsic Conformational Dynamics of Glycine and Alanine in Polarizable Molecular Dynamics Force Fields: Comparison to Spectroscopic Data

Understanding the role of intra- and intermolecular interactions in the formation of single- and double-helical structures of aromatic oligoamides: a computational study.

Side chain and backbone ordering in a polypeptide

Conformation Of Azobenzene-Modified Poly(Alpha-L-Glutamate) (Azoplga) In Thin Films: Solid State Nmr Studies

Construction of an Intermediate-Resolution Lattice Model and Re-Examination of the Helix-Coil Transition: a Dynamic Monte Carlo Simulation.

Ab initio molecular dynamics with intramolecular noncovalent interactions for unsolvated polypeptides

All-atom Replica Exchange Molecular Simulation of Protein BBL

Multiscale Simulation of Polyglutamine and the Effect of Neighboring Amino Acids on Oligomerization

Exploring the Helix-Coil Transition via All-Atom Equilibrium Ensemble Simulations

Stability and folding dynamics of polyglutamic acid

Pseudopotentials for coarse-grained cross-link-assisted modeling of protein structures

Some Insight into Stability of Amorphous Poly(Ethylene Glycol) Dimethyl Ether Polymers Based on Molecular Dynamics Simulations