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.

Molecular basis for the initiation of DNA primer synthesis.

During the initiation of DNA replication, oligonucleotide primers are synthesized de novo by primases and are subsequently extended by replicative polymerases to complete genome duplication. The primase-polymerase (Prim-Pol) superfamily is a diverse grouping of primases, which includes replicative primases and CRISPR-associated primase-polymerases (CAPPs) involved in adaptive immunity1-3. Although much is known about the activities of these enzymes, the precise mechanism used by primases to initiate primer synthesis has not been elucidated. Here we identify the molecular bases for the initiation of primer synthesis by CAPP and show that this mechanism is also conserved in replicative primases. The crystal structure of a primer initiation complex reveals how the incoming nucleotides are positioned within the active site, adjacent to metal cofactors and paired to the templating single-stranded DNA strand, before synthesis of the first phosphodiester bond. Furthermore, the structure of a Prim-Pol complex with double-stranded DNA shows how the enzyme subsequently extends primers in a processive polymerase mode. The structural and mechanistic studies presented here establish how Prim-Pol proteins instigate primer synthesis, revealing the requisite molecular determinants for primer synthesis within the catalytic domain. This work also establishes that the catalytic domain of Prim-Pol enzymes, including replicative primases, is sufficient to catalyse primer formation.

Impact of Remdesivir Incorporation along the Primer Strand on SARS-CoV-2 RNA-Dependent RNA Polymerase.

Remdesivir was the first antiviral drug that received emergency use authorization from the United States Food and Drug Administration and is now formally approved to treat COVID-19. Remdesivir is a nucleotide analogue that targets the RNA-dependent RNA polymerase (RdRp) of coronaviruses, including SARS-CoV-2. The solution of multiple RdRp structures has been one of the main axes of research in the race against the SARS-CoV-2 virus. Several hypotheses of the mechanism of inhibition of RdRp by remdesivir have been proposed, although open questions remain. This work uses molecular dynamics simulations to explore the impact of remdesivir and two analogues as incoming nucleotides and of up to four incorporations of remdesivir along the primer strand on RdRp. The simulation results suggest that the overall structure and the dynamical behavior of RdRp are destabilized by remdesivir and the two analogues in the incoming position. The incorporation of remdesivir along the primer strand impacts specific non-bonded interactions between the nascent RNA and the polymerase subunit, as well as the overall dynamical networks on RdRp. The strongest impact on the structure and dynamics are observed after three incorporations, when remdesivir is located at position -A3, in agreement with previously reported experimental and computational results. Our results provide atomic-level details of the role played by remdesivir on the disruption of RNA synthesis by RdRp and the main drivers of these disruptions.

Improvement of the Gaussian Electrostatic Model by separate fitting of Coulomb and exchange-repulsion densities and implementation of a new dispersion term.

The description of each separable contribution of the intermolecular interaction is a useful approach to develop polarizable force fields (polFFs). The Gaussian Electrostatic Model (GEM) is based on this approach, coupled with the use of density fitting techniques. In this work, we present the implementation and testing of two improvements of GEM: the Coulomb and exchange-repulsion energies are now computed with separate frozen molecular densities and a new dispersion formulation inspired by the Sum of Interactions Between Fragments Ab initio Computed polFF, which has been implemented to describe the dispersion and charge-transfer interactions. Thanks to the combination of GEM characteristics and these new features, we demonstrate a better agreement of the computed structural and condensed properties for water with experimental results, as well as binding energies in the gas phase with the ab initio reference compared with the previous GEM* potential. This work provides further improvements to GEM and the items that remain to be improved and the importance of the accurate reproduction for each separate contribution.

Development and application of quantum mechanics/molecular mechanics methods with advanced polarizable potentials.

Quantum mechanics/molecular mechanics (QM/MM) simulations are a popular approach to study various features of large systems. A common application of QM/MM calculations is in the investigation of reaction mechanisms in condensed-phase and biological systems. The combination of QM and MM methods to represent a system gives rise to several challenges that need to be addressed. The increase in computational speed has allowed the expanded use of more complicated and accurate methods for both QM and MM simulations. Here, we review some approaches that address several common challenges encountered in QM/MM simulations with advanced polarizable potentials, from methods to account for boundary across covalent bonds and long-range effects, to polarization and advanced embedding potentials.

Assessment of SAPT and Supermolecular EDA Approaches for the Development of Separable and Polarizable Force Fields.

Which is the best reference quantum chemical approach to decipher the energy components of the total interaction energy: Symmetry-Adapted Perturbation Theory (SAPT) or Supermolecular Energy Decomposition Analysis (EDA) methods? With the rise of physically motivated polarizable force fields (polFF) grounded on these procedures, the need to answer such a question becomes critical. We report a systematic and detailed assessment of three variants of SAPT (namely SAPT2, SAPT2+3, and SAPT(DFT)) and three supermolecular EDA approaches (ALMO, CSOV, and RVS). A set of challenging, strongly bound water complexes, (H2O)2, Zn2+. . .H2O, and F-. . .H2O, is used as "stress tests" for these electronic structure methods. We have developed a procedure to separate the induction energy into the polarization and charge-delocalization using an infinite-order strategy based on SAPT(DFT). This paper aims to provide not only an overview of the capabilities and limitations but also similarities of SAPT and supermolecular EDA approaches for polFF developments. Our results show that SAPT(DFT)/noS2 and ωB97X-D∥ALMO are the most accurate and reliable techniques.

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

A correct representation of the short-range contributions such as exchange-repulsion (Erep ) and charge-transfer (Ect ) is essential for the soundness of separable, anisotropic polarizable molecular mechanics potentials. Within the context of the SIBFA procedure, this is aimed at by explicit representations of lone pairs in their expressions. It is necessary to account for their anisotropic behaviors upon performing not only in-plane, but also out-of-plane, variations of a probe molecule or cation interacting with a target molecule or molecular fragment. Thus, Erep and Ect have to reproduce satisfactorily the corresponding anisotropies of their quantum chemical (QC) counterparts. A significant improvement of the out-of-plane dependencies was enabled when the sp2 and sp localized lone-pairs are, even though to a limited extent, delocalized on both sides of the plane, above and below the atom bearer but at the closely similar angles as the in-plane lone pair. We report calibration and validation tests on a series of monoligated complexes of a probe Zn(II) cation with several biochemically relevant ligands. Validations are then performed on several polyligated Zn(II) complexes found in the recognition sites of Zn-metalloproteins. Such calibrations and validations are extended to representative monoligated and polyligated complexes of Mg(II) and Ca(II). It is emphasized that the calibration of all three cations was for each ΔE contribution done on a small training set bearing on a limited number of representative N, O, and S monoligated complexes. Owing to the separable nature of ΔE, a secure transferability is enabled to a diversity of polyligated complexes. For these the relative errors with respect to the target ΔE(QC) values are generally < 3%. Overall, the article proposes a full set of benchmarks that could be useful for force field developers. © 2017 Wiley Periodicals, Inc.

Channeling through Two Stacked Guanine Quartets of One and Two Alkali Cations in the Li+, Na+, K+, and Rb+ Series. Assessment of the Accuracy of the SIBFA Anisotropic Polarizable Molecular Mechanics Potential.

Stacking of guanine quartets (GQs) can trigger the formation of DNA or RNA quadruple helices, which play numerous biochemical roles. The GQs are stabilized by alkali cations, mainly K+ and Na+, which can reside in, or channel through, the central axis of the GQ stems. Further, ion conduction through GQ wires can be leveraged for nanochemistry applications. G-quadruplex systems have been extensively studied by classical molecular dynamics (MD) simulations using pair-additive force fields or by quantum-chemical (QC) calculations. However, the non-polarizable force fields are very approximate, while QC calculations lack the necessary sampling. Thus, ultimate description of GQ systems would require long-enough simulations using advanced polarizable molecular mechanics (MM). However, to perform such calculations, it is first mandatory to evaluate the method's accuracy using benchmark QC. We report such an evaluation for SIBFA polarizable MM, bearing on the channeling (movement) of an alkali cation (Li+, Na+, K+, or Rb+) along the axis of two stacked G quartets interacting with either one or two ions. The QC energy profiles display markedly different features depending upon the cation but can be retrieved in the majority of cases by the SIBFA profiles. An appropriate balance of first-order (electrostatic and short-range repulsion) and second-order (polarization, charge-transfer, and dispersion) contributions within ΔE is mandatory. With two cations in the channel, the relative weights of the second-order contributions increase steadily upon increasing the ion size. In the G8 complexes with two K+ or two Rb+ cations, the sum of polarization and charge-transfer exceeds the first order terms for all ion positions.