Abstract:Conspectus Electrostatic interaction plays a significant role in determining many properties of biomolecules, which exist and function in aqueous solution, a highly polar environment. For example, proteins are composed of amino acids with charged, polar, and nonpolar side chains and their specific electrostatic properties are fundamental to the structure and function of proteins. An important issue that arises in computational study of biomolecular interaction and dynamics based on classical force field is lack of polarization. Polarization is a phenomenon in which the charge distribution of an isolated molecule will be distorted when interacting with another molecule or presented in an external electric field. The distortion of charge distribution is intended to lower the overall energy of the molecular system, which is counter balanced by the increased internal energy of individual molecules due to the distorted charge distributions. The amount of the charge redistribution, which characterizes the polarizability of a molecule, is determined by the level of the charge distortion. Polarization is inherently quantum mechanical, and therefore classical force fields with fixed atomic charges are incapable of capturing this important effect. As a result, simulation studies based on popular force fields, AMBER, CHARMM, etc., lack the polarization effect, which is a widely known deficiency in most computational studies of biomolecules today. Many efforts have been devoted to remedy this deficiency, such as adding additional movable charge on the atom, allowing atomic charges to fluctuate, or including induced multipoles. Although various successes have been achieved and progress at various levels has been reported over the past decades, the issue of lacking polarization in force field based simulations is far from over. For example, some of these methods do not always give converged results, and other methods require huge computational cost. This Account reviews recent work on developing polarized and polarizable force fields based on fragment quantum mechanical calculations for proteins. The methods described here are based on quantum mechanical calculations of proteins in solution, but with a different level of rigor and different computational efficiency for the molecular dynamics applications. In the general approach, a fragment quantum mechanical calculation for protein with implicit solvation is carried out to derive a polarized protein-specific charge (PPC) for any given protein structure. The PPC correctly reflects the polarization state of the protein in a given conformation, and it can also be dynamically changed as the protein changes conformation in dynamics simulations. Another approach that is computationally more efficient is the effective polarizable bond method in which only polar bonds or groups can be polarized and their polarizabilities are predetermined from quantum mechanical calculations of these groups in external electric fields. Both methods can be employed for applications in various situations by taking advantage of their unique features.

Development of a Polarizable Force Field Using Multiple Fluctuating Charges per Atom

Parametrizing A Polarizable Force Field from Ab Initio Data. I. the Fluctuating Point Charge Model

Development of an Effective Polarizable Bond Method for Biomolecular Simulation.

A Polarizable Atomic Multipole-Based Force Field for Molecular Dynamics Simulations of Anionic Lipids

Polarizable Force Fields Based On Physical Models And Quantum Chemical Calculations

A fast, convenient, polarizable electrostatic model for molecular dynamics

Energetic Molecule Specific Polarizable Force Field

Development of an anisotropic polarizable model for the all-atom AMOEBA force field

Construction of atomistic force field model for perfluorosulfonic acid polymer

Fast Evaluation Of Polarizable Forces

Polarizable Multipole-Based Force Field for Dimethyl and Trimethyl Phosphate

From bonds to interactions: comprehensive molecular characterization via polarizable bond-dipole approach

Higher Accuracy Achieved in the Simulations of Protein Structure Refinement, Protein Folding, and Intrinsically Disordered Proteins Using Polarizable Force Fields

Constant-pH Simulations with the Polarizable Atomic Multipole AMOEBA Force Field

Development of an effective polarizable bond method for biomolecular simulation

Drude Polarizable Force Field for Molecular Dynamics Simulations of Saturated and Unsaturated Zwitterionic Lipids.

Balanced polarizable Drude force field parameters for molecular anions: phosphates, sulfates, sulfamates, and oxides

Some practical approaches to treating electrostatic polarization of proteins.

Coarse-Grained Force Field for Polyethylene Oxide and Polyethylene Glycol Aqueous Solutions Based on a Polarizable Water Model.

Validation of Polarizable Force Field Parameters for Nucleic Acids by Inter-Molecular Interactions

Development of Accurate Coarse-Grained Force Fields for Weakly Polar Groups by an Indirect Parameterization Strategy.