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.