Realistic modeling of biomolecular systems requires an accurate treatment of electrostatics, including electronic polarization. Due to recent advances in physical models, simulation algorithms, and computing hardware, biomolecular simulations with advanced force fields at biologically relevant timescales are becoming increasingly promising. These advancements have not only led to new biophysical insights but also afforded opportunities to advance our understanding of fundamental intermolecular forces. This article describes the recent advances and applications, as well as future directions, of polarizable force fields in biomolecular simulations.

Tinker-HP is massively parallel software dedicated to polarizable molecular dynamics.

The AMOEBA polarizable atomic multipole force field for nucleic acids is presented. Valence and electrostatic parameters were determined from high-level quantum mechanical data, including structures, conformational energy, and electrostatic potentials, of nucleotide model compounds. Previously derived parameters for the phosphate group and nucleobases were incorporated. A total of over 35 μs of condensed-phase molecular dynamics simulations of DNA and RNA molecules in aqueous solution and crystal lattice were performed to validate and refine the force field. The solution and/or crystal structures of DNA B-form duplexes, RNA duplexes, and hairpins were captured with an average root-mean-squared deviation from NMR structures below or around 2.0 Å. Structural details, such as base pairing and stacking, sugar puckering, backbone and χ-torsion angles, groove geometries, and crystal packing interfaces, agreed well with NMR and/or X-ray. The interconversion between A- and B-form DNAs was observed in ethanol-water mixtures at 328 K. Crystal lattices of B- and Z-form DNA and A-form RNA were examined with simulations. For the RNA tetraloop, single strand tetramers, and HIV TAR with 29 residues, the simulated conformational states, J-coupling, nuclear Overhauser effect, and residual dipolar coupling data were compared with NMR results. Starting from a totally unstacked/unfolding state, the rCAAU tetranucleotide was folded into A-form-like structures during ∼1 μs molecular dynamics simulations.

One of the challenges in material science has been to prepare macro-or mesoporous zeolite. Although examples of their synthesis exist, there is a need for a facile yet versatile approach to such hierarchical structures. Here we report a concept for designing a single quaternary ammonium head amphiphilic template with strong ordered self-assembling ability through p-p stacking in hydrophobic side, which stabilizes the mesostructure to form singlecrystalline mesostructured zeolite nanosheets. The concept is demonstrated for the formation of a new type of MFI (zeolite framework code by International Zeolite Association) nanosheets joined with a 90°rotational boundary, which results in a mesoporous zeolite with highly specific surface area even after calcination. Low binding energies for this selfassembling system are supported by a theoretical analysis. A geometrical matching between the arrangement of aromatic groups and the zeolitic framework is speculated for the formation of single-crystalline MFI nanosheets.

Calcium ion is a versatile messenger in many cell-signaling processes. To achieve their functions, calcium-binding proteins selectively bind Ca against a background of competing ions such as Mg The high specificity of calcium-binding proteins has been intriguing since Mg has a higher charge density than Ca and is expected to bind more tightly to the carboxylate groups in calcium-binding pockets. Here, we showed that the specificity for Ca is dictated by the many-body polarization effect, which is an energetic cost arising from the dense packing of multiple residues around the metal ion. Since polarization has stronger distance dependence compared with permanent electrostatics, the cost associated with the smaller Mg is much higher than that with Ca and outweighs the electrostatic attraction favorable for Mg With the AMOEBA (atomic multipole optimized energetics for biomolecular simulation) polarizable force field, our simulations captured the relative binding free energy between Ca and Mg for proteins with various types of binding pockets and explained the nonmonotonic size dependence of the binding free energy in EF-hand proteins. Without electronic polarization, the smaller ions are always favored over larger ions and the relative binding free energy is roughly proportional to the net charge of the pocket. The many-body effect depends on both the number and the arrangement of charged residues. Fine-tuning of the ion selectivity could be achieved by combining the many-body effect and geometric constraint.

The extensibility of force field is a key to solve the missing parameter problem commonly found in force field applications. The extensibility of conventional force fields is traditionally managed in the parameterization procedure, which becomes impractical as the coverage of the force field increases above a threshold. A hierarchical atom-type definition (HAD) scheme is proposed to make extensible atom type definitions, which ensures that the force field developed based on the definitions are extensible. To demonstrate how HAD works and to prepare a foundation for future developments, two general force fields based on AMBER and DFF functional forms are parameterized for common organic molecules. The force field parameters are derived from the same set of quantum mechanical data and experimental liquid data using an automated parameterization tool, and validated by calculating molecular and liquid properties. The hydration free energies are calculated successfully by introducing a polarization scaling factor to the dispersion term between the solvent and solute molecules. © 2015 Wiley Periodicals, Inc.

The interactions between metal ions and proteins are ubiquitous in biology. The selective binding of metal ions has a variety of regulatory functions. Therefore, there is a need to understand the mechanism of protein-ion binding. The interactions involving metal ions are complicated in nature, where short-range charge-penetration, charge transfer, polarization, and many-body effects all contribute significantly, and a quantitative description of all these interactions is lacking. In addition, it is unclear how well current polarizable force fields can capture these energy terms and whether these polarization models are good enough to describe the many-body effects. In this work, two energy decomposition methods, absolutely localized molecular orbitals and symmetry-adapted perturbation theory, were utilized to study the interactions between Mg/Ca and model compounds for amino acids. Comparison of individual interaction components revealed that while there are significant charge-penetration and charge-transfer effects in Ca complexes, these effects can be captured by the van der Waals (vdW) term in the AMOEBA force field. The electrostatic interaction in Mg complexes is well described by AMOEBA since the charge penetration is small, but the distance-dependent polarization energy is problematic. Many-body effects were shown to be important for protein-ion binding. In the absence of many-body effects, highly charged binding pockets will be over-stabilized, and the pockets will always favor Mg and thus lose selectivity. Therefore, many-body effects must be incorporated in the force field in order to predict the structure and energetics of metalloproteins. Also, the many-body effects of charge transfer in Ca complexes were found to be non-negligible. The absorption of charge-transfer energy into the additive vdW term was a main source of error for the AMOEBA many-body interaction energies.

Molecular recognition is of paramount interest in many applications. Here we investigate a series of host-guest systems previously used in the SAMPL4 blind challenge by using molecular simualtions and the AMOEBA polarizable force field. The free energy results computed by Bennett’s acceptance ratio (BAR) method using the AMOEBA polarizable force field ranked favorably among the entries submitted to the SAMPL4 host-guest competition [Muddana et al., J. Comput.-Aided Mol. Des., 2014, 28, 305-317.]. In this work we conduct an in-depth analysis of the AMOEBA force-field host-guest binding thermodynamics by using both BAR and the orthogonal space random walk (OSRW) methods. The binding entropy-enthalpy contributions are analyzed for each host-guest system. For systems of inordinate binding entropy-enthalpy values, we further examine the hydrogen bonding patterns and configurational entropy contribution. The binding mechanism of this series of host-guest systems varies from ligand to ligand, driven by enthalpy and/or entropy changes. Convergence of BAR and OSRW binding free energy methods are discussed. Ultimately, this work illustrates the value of molecular modelling and advanced force fields for the exploration and interpretation of binding thermodynamics.