Importance of explicit smeared lone-pairs in anisotropic polarizable molecular mechanics. Torture track angular tests for exchange-repulsion and charge transfer contributions

Development of AMOEBA Polarizable Force Field for Rare-Earth La3+ Interaction with Bioinspired Ligands

Rare-earth metals (REMs) are crucial for many important industries, such as power generation and storage, in addition to cancer treatment and medical imaging. One promising new REM refinement approach involves mimicking the highly selective and efficient binding of REMs observed in relatively recently discovered proteins. However, realizing any such bioinspired approach requires an understanding of the biological recognition mechanisms. Here, we developed a new classical polarizable force field based on the AMOEBA framework for modeling a lanthanum ion (La3+) interacting with water, acetate, and acetamide, which have been found to coordinate the ion in proteins. The parameters were derived by comparing to high-level ab initio quantum mechanical (QM) calculations that include relativistic effects. The AMOEBA model, with advanced atomic multipoles and electronic polarization, is successful in capturing both the QM distance-dependent La3+-ligand interaction energies and experimental hydration free energy. A new scheme for pairwise polarization damping (POLPAIR) was developed to describe the polarization energy in La3+ interactions with both charged and neutral ligands. Simulations of La3+ in water showed water coordination numbers and ion-water distances consistent with previous experimental and theoretical findings. Water residence time analysis revealed both fast and slow kinetics in water exchange around the ion. This new model will allow investigation of fully solvated lanthanum ion-protein systems using GPU-accelerated dynamics simulations to gain insights on binding selectivity, which may be applied to the design of synthetic analogues.

Development of the Quantum-Inspired SIBFA Many-Body Polarizable Force Field: Enabling Condensed-Phase Molecular Dynamics Simulations

We present the extension of the Sum of Interactions Between Fragments Ab initio Computed (SIBFA) many-body polarizable force field to condensed-phase molecular dynamics (MD) simulations. The quantum-inspired SIBFA procedure is grounded on simplified integrals obtained from localized molecular orbital theory and achieves full separability of its intermolecular potential. It embodies long-range multipolar electrostatics (up to quadrupole) coupled to a short-range penetration correction (up to charge-quadrupole), exchange repulsion, many-body polarization, many-body charge transfer/delocalization, exchange dispersion, and dispersion (up to C10). This enables the reproduction of all energy contributions of ab initio symmetry-adapted perturbation theory (SAPT(DFT)) gas-phase reference computations. The SIBFA approach has been integrated within the Tinker-HP massively parallel MD package. To do so, all SIBFA energy gradients have been derived and the approach has been extended to enable periodic boundary conditions simulations using smooth particle mesh Ewald. This novel implementation also notably includes a computationally tractable simplification of the many-body charge transfer/delocalization contribution. As a proof of concept, we perform a first computational experiment defining a water model fitted on a limited set of SAPT(DFT) data. SIBFA is shown to enable a satisfactory reproduction of both gas-phase energetic contributions and condensed-phase properties highlighting the importance of its physically motivated functional form.

Multimolecular complexes of the phosphodiester anion with Zn(II) or Mg(II) and water molecules-Preliminary validations of a polarizable potential by ab initio quantum chemistry

Dimethyl phosphate (DMP^- ) is a model for the phosphodiester backbone of DNA, RNA, and phospholipids. It is central for the binding of divalent cations and water along the backbone of nucleic acids. Significant polarization and charge-transfer contributions and nonadditivity come into play in the multimolecular complexes organized around phosphate. Prior to large-scale molecular dynamics (MD) with advanced polarizable potentials, it is essential to evaluate how well the values and trends of intermolecular interaction energies (ΔE) from ab initio quantum chemistry (QC) and their individual contributions are reproduced in a diversity of such complexes. These differ by the starting binding modes of a divalent cation, Zn(II), namely direct, bi- or mono-dentate to anionic and/or ester oxygens, versus through-water binding. We present first the results from automated refinements of the individual contributions of the SIBFA potential with respect to their QC counterparts using a Zn(II) or a water probe. This is followed by validations on eight relaxed multimolecular complexes of DMP^- with Zn(II) or Mg(II) and seven waters, then on sixteen complexes of DMP^- with Zn(II) and eight waters in arrangements extracted from MD or energy-minimization on a droplet of sixty-four waters. This monitors the compared evolutions of SIBFA and QC ΔE and their individual contributions in the competing arrangements. Some waters, bridging Zn(II) and DMP^- , were found to have exceptionally large dipole moments, of up to 3.8 Debye. The perspectives of extension to a flexible phosphodiester backbone are discussed in the context of the SIBFA potential for DNA and RNA.

Spectrometric and computational studies of the binding of HIV-1 integrase inhibitors to viral DNA extremities

Three integrase strand transfer inhibitors are in intensive clinical use, raltegravir (RAL), elvitegravir (EVG) and dolutegravir (DTG). The onset of integrase resistance mutations limits their therapeutic efficiency. As put forth earlier, the drug affinity for the intasome could be improved by targeting preferentially the retroviral nucleobases, which are little, if at all, mutation-prone. We report experimental results of anisotropy fluorescence titrations of viral DNA by these three drugs. These show the DTG > EVG > RAL ranking of their inhibitory activities of the intasome to correspond to that of their free energies of binding, ∆Gs, to retroviral DNA, and that such a ranking is only governed by the binding enthalpies, ∆H, the entropy undergoing marginal variations. We sought whether this ranking might be reproduced through quantum chemistry (QC) Density Functional Theory calculations of intermolecular interaction energies between simplified models consisting of sole halobenzene ring and the highly conserved retroviral nucleobases G4 and C16. These calculations showed that binding of EVG has a small preference over DTG, while RAL ranked third. This indicates that additional interactions of the diketoacid parts of the drugs with DNA could be necessary to further enable preferential binding of DTG. The corresponding ∆Etotvalues computed with a polarizable molecular mechanics/dynamics procedure, Sum of Interactions Between Fragments Ab initio computed (SIBFA), showed good correlations with this ∆E(QC) ranking. These validations are an important step toward the use of polarizable molecular dynamics simulations on DTG or EVG derivatives in their complexes with the complete intasome, an application now motivated and enabled by the advent of currently developed and improved massively parallel software.

Describing Molecular Polarizability by a Bond Capacity Model

We propose a bond capacity model for describing molecular polarization in force field energy functions at the charge-only level. Atomic charges are calculated by allowing charge to flow between atom pairs according to a bond capacity and a difference in electrostatic potential. The bond capacity is closely related to the bond order and decays to zero as the bond distance is increased. The electrostatic potential is composed of an intrinsic potential, identified as the electronegativity, and a screened Coulomb potential from all other charges. The bond capacity model leads to integer fragment charges upon bond dissociation and displays linear scaling of the polarizability with system size. Bond capacity parameters can be derived from reference molecular polarizabilities, while electronegativity parameters can be derived from reference atomic charges or a reference molecular electrostatic potential. Out-of-plane polarization for planar systems is modeled by off-nuclei charge sites. The model is shown to be able to reproduce anisotropic reference molecular polarizabilities with an accuracy of ∼10% using a limited set of bond capacity parameters and can describe both inter- and intramolecular polarization.

A physically grounded damped dispersion model with particle mesh Ewald summation

Accurate modeling of dispersion is critical to the goal of predictive biomolecular simulations. To achieve this accuracy, a model must be able to correctly capture both the short-range and asymptotic behavior of dispersion interactions. We present here a damped dispersion model based on the overlap of charge densities that correctly captures both regimes. The overlap damped dispersion model represents a classical physical interpretation of dispersion: the interaction between the instantaneous induced dipoles of two distinct charge distributions. This model is shown to be an excellent fit with symmetry adapted perturbation theory dispersion energy calculations, yielding an RMS error on the S101x7 database of 0.5 kcal/mol. Moreover, the damping function used in this model is wholly derived and parameterized from the electrostatic dipole-dipole interaction, making it not only physically grounded but transferable as well.

On the Computational Characterization of Charge-Transfer Effects in Noncovalently Bound Molecular Complexes

Charge-transfer (CT) is an important binding force in the formation of intermolecular complexes, and there have been a variety of theoretical models proposed to quantify this effect. These approaches, which typically rely on a definition of a "CT-free" state based on a partition of the system, sometimes yield significantly different results for a given intermolecular complex. Two widely used definitions of the "CT-free" state, the absolutely localized molecular orbitals (ALMO) method (where only on-fragment orbital mixings are permitted) and the constrained density functional theory (CDFT) approach (where fragment electron populations are fixed), are carefully examined in this work. Natural bond orbital (NBO) and the regularized symmetry-adapted perturbation theory (SAPT) are also briefly considered. Results for the ALMO and CDFT definitions of CT are compared on a broad range of model systems, including hydrogen-bonding systems, borane complexes, metal-carbonyl complexes, and complexes formed by water and metal cations. For most of these systems, CDFT yields a much smaller equilibrium CT energy compared to that given by the ALMO-based definition. This is mainly because the CDFT population constraint does not fully inhibit CT, which means that the CDFT "CT-free" state is in fact CT-contaminated. Examples of this contamination include (i) matching forward and backward donation (e.g., formic acid dimer) and (ii) unidirectional CT without changing fragment populations. The magnitude of the latter effect is quantified in systems such as the water dimer by employing a 3-space density constraint in addition to the orbital constraint. Furthermore, by means of the adiabatic EDA, it is shown that several observable effects of CT, such as the "pyramidalization" of the planar BH₃ molecule upon the complexation with Lewis bases, already appear on the "CT-free" CDFT surface. These results reveal the essential distinctions between the ALMO and CDFT definitions of CT and suggest that the former is more consistent with accepted understanding of the role of CT in intermolecular binding.