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.

An effective all-atom potential for proteins

A Discontinuous Potential Model for Protein-Protein Interactions.

Iterative derivation of effective potentials to sample the conformational space of proteins at atomistic scale

Monte Carlo Sampling of Protein Folding by Combining an All-Atom Physics-Based Model with a Native State Bias

Development of an Effective Polarizable Bond Method for Biomolecular Simulation.

A New-Protein Folding Recognition Potential Function

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

A Peptide Potential Based on a Bond Dipole Representation of Electrostatics

A Force Balanced Fragmentation Method for Ab Initio Molecular Dynamic Simulation of Protein

Significantly Improved Protein Folding Thermodynamics Using a Dispersion-Corrected Water Model and a New Residue-Specific Force Field

Dissipative Particle Dynamics Simulation of Protein Folding in Explicit and Implicit Solvents: Coarse-Grained Model for Atomic Resolution

A quantum mechanical computational method for modeling electrostatic and solvation effects of protein

Quantum Fragment Based ab Initio Molecular Dynamics for Proteins.

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

Monte Carlo Protein Folding: Simulations of Met-Enkephalin with Solvent-Accessible Area Parameterizations

Can a continuum solvent model reproduce the free energy landscape of a β-hairpin folding in water?

A Test of Implicit Solvent Models on the Folding Simulation of the GB1 Peptide.

Learning effective amino acid interactions through iterative stochastic techniques

Transferable coarse-grained potential for $\textit{de novo}$ protein folding and design

Further optimization of a hybrid united-atom and coarse-grained force field for folding simulations: Improved backbone hydration and interactions between charged side chains.

Pair potentials for protein folding: Choice of reference states and sensitivity of predicted native states to variations in the interaction schemes