Enforcing Local DNA Kinks by Sequence-Selective Trisintercalating Oligopeptides of a Tricationic Porphyrin: A Polarizable Molecular Dynamics Study.

Bisacridinyl-bisarginyl porphyrin (BABAP) is a trisintercalating derivative of a tricationic porphyrin, formerly designed and synthesized in order to selectively target and photosensitize the ten-base pair palindromic sequence d(CGGGCGCCCG)2 . We resorted to the previously derived (Far et al., 2004) lowest energy-minimized (EM) structure of the BABAP complex with this sequence as a starting point. We performed polarizable molecular dynamics (MD) on this complex. It showed, over a 150 ns duration, the persistent binding of the Arg side-chain on each BABAP arm to the two G bases upstream from the central porphyrin intercalation site. We subsequently performed progressive shortenings of the connector chain linking the Arg-Gly backbone to the acridine, from n=6 methylenes to 4, followed by removal of the Gly backbone and further connector shortenings, from n=4 to n=1. These resulted into progressive deformations ('kinks') of the DNA backbone. In its most accented kinked structure, the DNA backbone was found to have a close overlap with that of DNA bound to Cre recombinase, with, at the level of one acridine intercalation site, negative roll and positive tilt values consistent with those experimentally found for this DNA at its own kinked dinucleotide sequence. Thus, in addition to their photosensitizing properties, some BABAP derivatives could induce sequence-selective, controlled DNA deformations, which are targets for cleavage by endonucleases or for repair enzymes.

Targeting the Major Groove of the Palindromic d(GGCGCC)2 Sequence by Oligopeptide Derivatives of Anthraquinone Intercalators.

GC-rich sequences are recurring motifs in oncogenes and retroviruses and could be targeted by noncovalent major-groove therapeutic ligands. We considered the palindromic sequence d(G1G2C3G4C5C6)2, and designed several oligopeptide derivatives of the anticancer intercalator mitoxantrone. The stability of their complexes with an 18-mer oligonucleotide encompassing this sequence in its center was validated using polarizable molecular dynamics. We report the most salient structural features of two novel compounds, having a dialkylammonium group as a side chain on both arms. The anthraquinone ring is intercalated in the central d(CpG)2 sequence with its long axis perpendicular to that of the two base pairs. On each strand, this enables each ammonium group to bind in-register to O6/N7 of the two facing G bases upstream. We subsequently designed tris-intercalating derivatives, each dialkylammonium substituted with a connector to an N9-aminoacridine intercalator extending our target range from a six- to a ten-base-pair palindromic sequence, d(C1G2G3G4C5G6C7C8C9G10)2. The structural features of the complex of the most promising derivative are reported. The present design strategy paves the way for designing intercalator-oligopeptide derivatives with even higher selectivity, targeting an increased number of DNA bases, going beyond ten.

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.

Intermolecular interactions of the extended recognition site of VIM-2 metallo-ß-lactamase with 1,2,4-triazole-3-thione inhibitors. Validations of a polarizable molecular mechanics potential by ab initio QC.

Molecular dynamics on the complexes of inhibitors with Zn-metalloproteins are a privileged area of applications of polarizable molecular mechanics potentials. With which accuracy could these reproduce the QC intermolecular interaction energies in the two mono-zinc cores and in the dizinc core, toward full-fledged MD simulations on the entire protein complexes? We considered the complexes of the extended recognition site of a Zn-dependent metallo-ß-lactamase, VIM-2, produced by bacteria responsible for nosocomial infections, with five newly synthesized inhibitors sharing an original dizinc binding group, 1,2,4-triazole-3-thione (TZT). We considered the energy-minimized structures of each of the five VIM-2 complexes obtained with the SIBFA potential. Energy decomposition analyses (EDA) at the HF level enabled to compare the QC and the SIBFA ΔE values and their contributions in the zinc cores, with and without TZT, totaling 30 complexes. With one exception, the ΔE(QC) values were reproduced with relative errors <1.5%. We next considered the complex of the entire inhibitors with an extended model of VIM-2 recognition site, totaling up to 280 atoms. ΔE(SIBFA) could closely reproduce ΔE(QC). EDA analyses were resumed on the complexes of each inhibitor arm with its interacting VIM-2 residues. As a last step, EDA results at correlated levels were analyzed for the mono- and dizinc sites enabling comparisons with dispersion-augmented ΔE(SIBFA) and correlated multipoles and polarizabilities. Closely reproducing ΔE(QC) and the contrasting trends of its individual contributions should enable for dependable free energy perturbation studies and comparisons to recent experimental ΔG values, limiting as much as possible the reliance on error compensations.

Calibration of the dianionic phosphate group: Validation on the recognition site of the homodimeric enzyme phosphoglucose isomerase.

We calibrate and validate the parameters necessary to represent the dianionic phosphate group (DPG) in molecular mechanics. DPG is an essential fragment of signaling biological molecules and protein-binding ligands. It is a constitutive fragment of biosensors, which bind to the dimer interface of phosphoglucose isomerase (PGI), an intracellular enzyme involved in sugar metabolism, as well as an extracellular protein known as autocrine motility factor (AMF) closely related to metastasis formation. Our long-term objective is to design DPG-based biosensors with enhanced affinities for AMF/PGI cancer biomarker in blood. Molecular dynamics with polarizable potentials could be used toward this aim. This requires to first evaluate the accuracy of such potentials upon representing the interactions of DPG with its PGI ligands and tightly bound water molecules. Such evaluations are done by comparisons with high-level ab initio quantum chemistry (QC) calculations. We focus on the Sum of Interactions Between Fragments Ab initio computed (SIBFA) polarizable molecular mechanics procedure. We present first the results of the DPG calibration. This is followed by comparisons between ΔE(SIBFA) and ΔE(QC) regarding bi-molecular complexes of DPG with the main-chain and side-chain PGI residues, which bind to it in the recognition site. We then consider DPG complexes with an increasing number of PGI residues. The largest QC complexes encompass the entirety of the recognition site, with six structural water molecules totaling up to 211 atoms. A persistent and satisfactory agreement could be shown between ΔE(SIBFA) and ΔE(QC). These validations constitute an essential first step toward large-scale molecular dynamics simulations of DPG-based biosensors bound at the PGI dimer interface. © 2020 Wiley Periodicals, Inc.

The inhibition process of HIV-1 integrase by diketoacids molecules: Understanding the factors governing the better efficiency of dolutegravir.

The Human Immunodeficiency Virus-1 integrase is responsible for the covalent insertion of a newly synthesized double-stranded viral DNA into the host cells, and is an emerging target for antivirus drug design. Raltegravir (RAL) and elvitegravir (EVG) are the first two integrase strand transfer inhibitors used in therapy. However, treated patients eventually develop detrimental resistance mutations. By contrast, a recently approved drug, dolutegravir (DTG), presents a high barrier to resistance. This study aims to understand the increased efficiency of DTG upon focusing on its interaction properties with viral DNA. The results showed DTG to be involved in more extended interactions with viral DNA than EVG. Such interactions involve the halobenzene and scaffold of DTG and EVG and bases 5'G-43', 3'A35'and 3'C45'.

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.

Complexes of a Zn-metalloenzyme binding site with hydroxamate-containing ligands. A case for detailed benchmarkings of polarizable molecular mechanics/dynamics potentials when the experimental binding structure is unknown.

Zn-metalloproteins are a major class of targets for drug design. They constitute a demanding testing ground for polarizable molecular mechanics/dynamics aimed at extending the realm of quantum chemistry (QC) to very long-duration molecular dynamics (MD). The reliability of such procedures needs to be demonstrated upon comparing the relative stabilities of competing candidate complexes of inhibitors with the recognition site stabilized in the course of MD. This could be necessary when no information is available regarding the experimental structure of the inhibitor-protein complex. Thus, this study bears on the phosphomannose isomerase (PMI) enzyme, considered as a potential therapeutic target for the treatment of several bacterial and parasitic diseases. We consider its complexes with 5-phospho-d-arabinonohydroxamate and three analog ligands differing by the number and location of their hydroxyl groups. We evaluate the energy accuracy expectable from a polarizable molecular mechanics procedure, SIBFA. This is done by comparisons with ab initio quantum-chemistry (QC) calculations in the following cases: (a) the complexes of the four ligands in three distinct structures extracted from the entire PMI-ligand energy-minimized structures, and totaling up to 264 atoms; (b) the solvation energies of several energy-minimized complexes of each ligand with a shell of 64 water molecules; (c) the conformational energy differences of each ligand in different conformations characterized in the course of energy-minimizations; and (d) the continuum solvation energies of the ligands in different conformations. The agreements with the QC results appear convincing. On these bases, we discuss the prospects of applying the procedure to ligand-macromolecule recognition problems. © 2016 Wiley Periodicals, Inc.

Scalable improvement of SPME multipolar electrostatics in anisotropic polarizable molecular mechanics using a general short-range penetration correction up to quadrupoles.

We propose a general coupling of the Smooth Particle Mesh Ewald SPME approach for distributed multipoles to a short-range charge penetration correction modifying the charge-charge, charge-dipole and charge-quadrupole energies. Such an approach significantly improves electrostatics when compared to ab initio values and has been calibrated on Symmetry-Adapted Perturbation Theory reference data. Various neutral molecular dimers have been tested and results on the complexes of mono- and divalent cations with a water ligand are also provided. Transferability of the correction is adressed in the context of the implementation of the AMOEBA and SIBFA polarizable force fields in the TINKER-HP software. As the choices of the multipolar distribution are discussed, conclusions are drawn for the future penetration-corrected polarizable force fields highlighting the mandatory need of non-spurious procedures for the obtention of well balanced and physically meaningful distributed moments. Finally, scalability and parallelism of the short-range corrected SPME approach are addressed, demonstrating that the damping function is computationally affordable and accurate for molecular dynamics simulations of complex bio- or bioinorganic systems in periodic boundary conditions.

Could the "Janus-like" properties of the halobenzene CX bond (X=Cl, Br) be leveraged to enhance molecular recognition?

The CX bond in halobenzenes (XCl, Br) exhibits a dual character, being electron-deficient along the CX direction, and electron-rich on its flanks. We sought to amplify both features by resorting to electron-withdrawing and electron-donating substituents, respectively. This was done by quantum chemistry (QC) computations in the recognition sites of three protein targets: farnesyl transferase, coagulation factor Xa, and the HIV-1 integrase. In this context, some substituents, notably fluorine, CF3 , and NHCH3 , afforded significant overall gains in the binding energies as compared to the parent halobenzene, in the 2-5 kcal/mol range. In fact, we found that some di- and up to tetra-substitutions enabled even larger gains than those they contribute separately owing to many-body effects. Moreover, desolvation was also found to be a key contributor to the energy balances. As a consequence, some particular substituents, contributing to reduce the halobenzene dipole moment, accordingly reduced solvation: this factor acted in synergy with their enhancement of the intermolecular interaction energies along and around the CX bond. We could thus leverage the "Janus-like" properties of such a bond and the fact that it can be tuned and possibly amplified by well-chosen substituents. We propose a simple yet rigorous computational strategy resorting to QC to prescreen novel substituted halobenzenes. The QC results on the recognition sites then set benchmarks to validate polarizable molecular mechanics/dynamics approaches used to handle the entirety of the inhibitor-protein complex.

Quantum-chemistry based calibration of the alkali metal cation series (Li(+)-Cs(+)) for large-scale polarizable molecular mechanics/dynamics simulations.

The alkali metal cations in the series Li(+)-Cs(+) act as major partners in a diversity of biological processes and in bioinorganic chemistry. In this article, we present the results of their calibration in the context of the SIBFA polarizable molecular mechanics/dynamics procedure. It relies on quantum-chemistry (QC) energy-decomposition analyses of their monoligated complexes with representative O-, N-, S-, and Se- ligands, performed with the aug-cc-pVTZ(-f) basis set at the Hartree-Fock level. Close agreement with QC is obtained for each individual contribution, even though the calibration involves only a limited set of cation-specific parameters. This agreement is preserved in tests on polyligated complexes with four and six O- ligands, water and formamide, indicating the transferability of the procedure. Preliminary extensions to density functional theory calculations are reported.

A supervised fitting approach to force field parametrization with application to the SIBFA polarizable force field.

A supervised, semiautomated approach to force field parameter fitting is described and applied to the SIBFA polarizable force field. The I-NoLLS interactive, nonlinear least squares fitting program is used as an engine for parameter refinement while keeping parameter values within a physical range. Interactive fitting is shown to avoid many of the stability problems that frequently afflict highly correlated, nonlinear fitting problems occurring in force field parametrizations. The method is used to obtain parameters for the H2O, formamide, and imidazole molecular fragments and their complexes with the Mg(2+) cation. Reference data obtained from ab initio calculations using an auc-cc-pVTZ basis set exploit advances in modern computer hardware to provide a more accurate parametrization of SIBFA than has previously been available.

Polarizable molecular mechanics studies of Cu(I)/Zn(II) superoxide dismutase: bimetallic binding site and structured waters.

The existence of a network of structured waters in the vicinity of the bimetallic site of Cu/Zn-superoxide dismutase (SOD) has been inferred from high-resolution X-ray crystallography. Long-duration molecular dynamics (MD) simulations could enable to quantify the lifetimes and possible interchanges of these waters between themselves as well as with a ligand diffusing toward the bimetallic site. The presence of several charged or polar ligands makes it necessary to resort to second-generation polarizable potentials. As a first step toward such simulations, we benchmark in this article the accuracy of one such potential, sum of interactions between fragments Ab initio computed (SIBFA), by comparisons with quantum mechanics (QM) computations. We first consider the bimetallic binding site of a Cu/Zn-SOD, in which three histidines and a water molecule are bound to Cu(I) and three histidines and one aspartate are bound to Zn(II). The comparisons are made for different His6 complexes with either one or both cations, and either with or without Asp and water. The total net charges vary from zero to three. We subsequently perform preliminary short-duration MD simulations of 296 waters solvating Cu/Zn-SOD. Six representative geometries are selected and energy-minimized. Single-point SIBFA and QM computations are then performed in parallel on model binding sites extracted from these six structures, each of which totals 301 atoms including the closest 28 waters from the Cu metal site. The ranking of their relative stabilities as given by SIBFA is identical to the QM one, and the relative energy differences by both approaches are fully consistent. In addition, the lowest-energy structure, from SIBFA and QM, has a close overlap with the crystallographic one. The SIBFA calculations enable to quantify the impact of polarization and charge transfer in the ranking of the six structures. Five structural waters, which connect Arg141 and Glu131, are endowed with very high dipole moments (2.7-3.0 Debye), equal and larger than the one computed by SIBFA in ice-like arrangements (2.7 D).

Synthesis and structure-activity relationship of non-peptidic antagonists of neuropilin-1 receptor.

Neuropilins (NRPs) are VEGF-A165 co-receptors over-expressed in tumor cells, and considered as targets in angiogenic-related pathologies. We previously identified compound 1, the first non-peptidic antagonist of the VEGF-A165/NRP binding, which exhibits in vivo anti-angiogenic and anti-tumor activities. We report here the synthesis and biological evaluations of new antagonists structurally-related to compound 1. Among these molecules, 4a, 4c and 4d show cytotoxic effects on HUVEC and MDA-MB-31 cells, and antagonize VEGF-A165/NRP-1 binding. This study confirmed our key structure-activity relationships hypothesis and paved the way to compound 1 'hit to lead' optimization.

Substituent-modulated affinities of halobenzene derivatives to the HIV-1 integrase recognition site. Analyses of the interaction energies by parallel quantum chemical and polarizable molecular mechanics.

The C-X bond of halobenzenes (X = Cl, Br) has a dual character, its electron density being depleted in its prolongation and built-up on its sides. We have recently considered three protein or nucleic acid recognition sites of halobenzenes and quantified the energy gains that either electron-attracting substituents or electron-donating ones contribute due to such a character (El Hage et al., paper in revision). Nonadditivity was found to impact the total interaction energies. We focus here on one recognition site, that of the HIV-1 integrase, in which the halobenzene ring of the drug elvitegravir is sandwiched between a guanine and a cytosine base. We perform energy-decomposition analyses of the ab initio quantum-chemistry (QC) binding energies of the parent halobenzene ring and its derivatives with this G-C base pair. In these complexes, the nonadditivity of ΔE could be traced back mostly to the polarization contribution Epol. In view of large-scale applications to the entirety of the complex formed between the integrase, the viral DNA, and the whole drug, the analyses were performed in parallel with a polarizable molecular mechanics method, SIBFA. This method could faithfully reproduce most features of the QC energies. This is due to its use of QC-derived distributed multipoles and polarizabilities, which enable us to account for both nonisotropy and nonadditivity.

S/G-1: an ab initio force-field blending frozen Hermite Gaussian densities and distributed multipoles. Proof of concept and first applications to metal cations.

We demonstrate as a proof of principle the capabilities of a novel hybrid MM'/MM polarizable force field to integrate short-range quantum effects in molecular mechanics (MM) through the use of Gaussian electrostatics. This lead to a further gain in accuracy in the representation of the first coordination shell of metal ions. It uses advanced electrostatics and couples two point dipole polarizable force fields, namely, the Gaussian electrostatic model (GEM), a model based on density fitting, which uses fitted electronic densities to evaluate nonbonded interactions, and SIBFA (sum of interactions between fragments ab initio computed), which resorts to distributed multipoles. To understand the benefits of the use of Gaussian electrostatics, we evaluate first the accuracy of GEM, which is a pure density-based Gaussian electrostatics model on a test Ca(II)-H2O complex. GEM is shown to further improve the agreement of MM polarization with ab initio reference results. Indeed, GEM introduces nonclassical effects by modeling the short-range quantum behavior of electric fields and therefore enables a straightforward (and selective) inclusion of the sole overlap-dependent exchange-polarization repulsive contribution by means of a Gaussian damping function acting on the GEM fields. The S/G-1 scheme is then introduced. Upon limiting the use of Gaussian electrostatics to metal centers only, it is shown to be able to capture the dominant quantum effects at play on the metal coordination sphere. S/G-1 is able to accurately reproduce ab initio total interaction energies within closed-shell metal complexes regarding each individual contribution including the separate contributions of induction, polarization, and charge-transfer. Applications of the method are provided for various systems including the HIV-1 NCp7-Zn(II) metalloprotein. S/G-1 is then extended to heavy metal complexes. Tested on Hg(II) water complexes, S/G-1 is shown to accurately model polarization up to quadrupolar response level. This opens up the possibility of embodying explicit scalar relativistic effects in molecular mechanics thanks to the direct transferability of ab initio pseudopotentials. Therefore, incorporating GEM-like electron density for a metal cation enable the introduction of nonambiguous short-range quantum effects within any point-dipole based polarizable force field without the need of an extensive parametrization.

Could an anisotropic molecular mechanics/dynamics potential account for sigma hole effects in the complexes of halogenated compounds?

Halogenated compounds are gaining an increasing importance in medicinal chemistry and materials science. Ab initio quantum chemistry (QC) has unraveled the existence of a "sigma hole" along the C-X (X = F, Cl, Br, I) bond, namely, a depletion of electronic density prolonging the bond, concomitant with a build-up on its sides, both of which are enhanced along the F < Cl < Br < I series. We have evaluated whether these features were intrinsically built-in in an anisotropic, polarizable molecular mechanics (APMM) procedure such as SIBFA (sum of interactions between fragments ab initio computed). For that purpose, we have computed the interaction energies of fluoro-, chloro-, and bromobenzene with two probes: a divalent cation, Mg(II), and water approaching X through either one H or its O atom. This was done by parallel QC energy-decomposition analyses (EDA) and SIBFA computations. With both probes, the leading QC contribution responsible for the existence of the sigma hole is the Coulomb contribution E(c). For all three halogenated compounds, and with both probes, the in- and out-of-plane angular features of E(c) were closely mirrored by the SIBFA electrostatic multipolar contribution (E(MTP)). Resorting to such a contribution thus dispenses with empirically-fitted "extra", off-centered partial atomic charges as in classical molecular mechanics/dynamics.

Potent inhibitors of CDK5 derived from roscovitine: synthesis, biological evaluation and molecular modelling.

Cyclin dependent kinase 5 (CDK5) is a serine/threonine kinase belonging to the cyclin dependent kinase (CDK) family. CDK5 is involved in numerous neuronal diseases (including Alzheimer's or Parkinson's diseases, stroke, traumatic brain injury), pain signaling and cell migration. In the present Letter, we describe syntheses and biological evaluations of new 2,6,9-trisubstituted purines, structurally related to roscovitine, a promising CDK inhibitor currently in clinical trials (CDK1/Cyclin B, IC(50)=350 nM; CDK5/p25, IC(50)=200 nM). These new molecules were synthesized using an original Buchwald-Hartwig catalytic procedure; several compounds (3j, 3k, 3l, 3e, 4k, 6b, 6c) displayed potent kinase inhibitory potencies against CDK5 (IC(50) values ranging from 17 to 50 nM) and showed significant cell death inducing activities (IC(50) values ranging from 2 to 9 µM on SH-SY5Y). The docking of the inhibitors into the ATP binding domain of the CDK5 catalytic site highlighted the discriminatory effect of a hydrogen bond involving the CDK5 Lys-89. In addition, the calculated final energy balances for complexation measured for several inhibitors is consistent with the ranking of the IC(50) values. Lastly, we observed that several compounds exhibit submicromolar activities against DYRK1A (dual specificity, tyrosine phosphorylation regulated kinase 1A), a kinase involved in Down syndrome and Alzheimer's disease (3g, 3h, 4m; IC(50) values ranging from 300 to 400 nM).

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.

The reaction mechanism of type I phosphomannose isomerases: new information from inhibition and polarizable molecular mechanics studies.

Type I phosphomannose isomerases (PMIs) are zinc-dependent metalloenzymes involved in the reversible isomerization of D-mannose 6-phosphate (M6P) and D-fructose 6-phosphate (F6P). 5-Phospho-D-arabinonohydroxamic acid (5PAH), an inhibitor endowed with nanomolar affinity for yeast (Type I) and Pseudomonas aeruginosa (Type II) PMIs (Roux et al., Biochemistry 2004; 43:2926-2934), strongly inhibits human (Type I) PMI (for which we report an improved expression and purification procedure), as well as Escherichia coli (Type I) PMI. Its K(i) value of 41 nM for human PMI is the lowest value ever reported for an inhibitor of PMI. 5-Phospho-D-arabinonhydrazide, a neutral analogue of the reaction intermediate 1,2-cis-enediol, is about 15 times less efficient at inhibiting both enzymes, in accord with the anionic nature of the postulated high-energy reaction intermediate. Using the polarizable molecular mechanics, sum of interactions between fragments ab initio computed (SIBFA) procedure, computed structures of the complexes between Candida albicans (Type I) PMI and the cyclic substrate ß-D-mannopyranose 6-phosphate (ß-M6P) and between the enzyme and the high-energy intermediate analogue inhibitor 5PAH are reported. Their analysis allows us to identify clearly the nature of each individual active site amino acid and to formulate a hypothesis for the overall mechanism of the reaction catalyzed by Type I PMIs, that is, the ring-opening and isomerization steps, respectively. Following enzyme-catalyzed ring-opening of ß-M6P by zinc-coordinated water and Gln111 ligands, Lys136 is identified as the probable catalytic base involved in proton transfer between the two carbon atoms C1 and C2 of the substrate D-mannose 6-phosphate.