ABSTRACT
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.
Subject(s)
Molecular Dynamics Simulation , Porphyrins , Porphyrins/chemistry , DNA/chemistry , Oligopeptides , AcridinesABSTRACT
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.
Subject(s)
Intercalating Agents , Oligopeptides , Intercalating Agents/pharmacology , Intercalating Agents/chemistry , Mitoxantrone/pharmacology , DNA/chemistry , Molecular Dynamics Simulation , Nucleic Acid ConformationABSTRACT
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.
Subject(s)
Molecular Dynamics Simulation , Water , Static ElectricityABSTRACT
Following our previous work ( Chem. Sci. 2021, 12, 4889-4907), we study the structural dynamics of the SARS-CoV-2 Main Protease dimerization interface (apo dimer) by means of microsecond adaptive sampling molecular dynamics simulations (50 µs) using the AMOEBA polarizable force field (PFF). This interface is structured by a complex H-bond network that is stable only at physiological pH. Structural correlations analysis between its residues and the catalytic site confirms the presence of a buried allosteric site. However, noticeable differences in allosteric connectivity are observed between PFFs and non-PFFs. Interfacial polarizable water molecules are shown to appear at the heart of this discrepancy because they are connected to the global interface H-bond network and able to adapt their dipole moment (and dynamics) to their diverse local physicochemical microenvironments. The water-interface many-body interactions appear to drive the interface volume fluctuations and to therefore mediate the allosteric interactions with the catalytic cavity.
Subject(s)
Molecular Dynamics Simulation , SARS-CoV-2/metabolism , Viral Matrix Proteins/chemistry , Water/chemistry , Allosteric Site , COVID-19/pathology , COVID-19/virology , Catalytic Domain , Dimerization , Humans , Hydrogen Bonding , Hydrogen-Ion Concentration , SARS-CoV-2/isolation & purification , Viral Matrix Proteins/metabolismABSTRACT
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.
Subject(s)
Density Functional Theory , Magnesium/chemistry , Organophosphates/chemistry , Water/chemistry , Zinc/chemistry , Molecular Dynamics SimulationABSTRACT
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.
ABSTRACT
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.
Subject(s)
Enzyme Inhibitors/chemistry , Enzyme Inhibitors/pharmacology , Molecular Dynamics Simulation , Thiones/chemistry , Thiones/pharmacology , beta-Lactamases/metabolism , Binding Sites , Models, Molecular , Molecular Structure , Protein Conformation , beta-Lactamases/chemistryABSTRACT
Resistance to ß-lactam antibiotics in Gram-negatives producing metallo-ß-lactamases (MBLs) represents a major medical threat and there is an extremely urgent need to develop clinically useful inhibitors. We previously reported the original binding mode of 5-substituted-4-amino/H-1,2,4-triazole-3-thione compounds in the catalytic site of an MBL. Moreover, we showed that, although moderately potent, they represented a promising basis for the development of broad-spectrum MBL inhibitors. Here, we synthesized and characterized a large number of 4-amino-1,2,4-triazole-3-thione-derived Schiff bases. Compared to the previous series, the presence of an aryl moiety at position 4 afforded an average 10-fold increase in potency. Among 90 synthetic compounds, more than half inhibited at least one of the six tested MBLs (L1, VIM-4, VIM-2, NDM-1, IMP-1, CphA) with Ki values in the µM to sub-µM range. Several were broad-spectrum inhibitors, also inhibiting the most clinically relevant VIM-2 and NDM-1. Active compounds generally contained halogenated, bicyclic aryl or phenolic moieties at position 5, and one substituent among o-benzoic, 2,4-dihydroxyphenyl, p-benzyloxyphenyl or 3-(m-benzoyl)-phenyl at position 4. The crystallographic structure of VIM-2 in complex with an inhibitor showed the expected binding between the triazole-thione moiety and the dinuclear centre and also revealed a network of interactions involving Phe61, Tyr67, Trp87 and the conserved Asn233. Microbiological analysis suggested that the potentiation activity of the compounds was limited by poor outer membrane penetration or efflux. This was supported by the ability of one compound to restore the susceptibility of an NDM-1-producing E. coli clinical strain toward several ß-lactams in the presence only of a sub-inhibitory concentration of colistin, a permeabilizing agent. Finally, some compounds were tested against the structurally similar di-zinc human glyoxalase II and found weaker inhibitors of the latter enzyme, thus showing a promising selectivity towards MBLs.
Subject(s)
Schiff Bases/pharmacology , Thiones/pharmacology , Triazoles/pharmacology , beta-Lactamase Inhibitors/pharmacology , beta-Lactamases/metabolism , Crystallography, X-Ray , Escherichia coli/drug effects , Escherichia coli Proteins/metabolism , Human Umbilical Vein Endothelial Cells , Humans , Microbial Sensitivity Tests , Protein Binding , Pseudomonas aeruginosa/chemistry , Schiff Bases/chemical synthesis , Schiff Bases/metabolism , Thiones/chemical synthesis , Thiones/metabolism , Triazoles/chemical synthesis , Triazoles/metabolism , beta-Lactamase Inhibitors/chemical synthesis , beta-Lactamase Inhibitors/metabolismABSTRACT
The HIV-1 integrase (IN) is a major target for the design of novel anti-HIV inhibitors. Among these, three inhibitors which embody a halobenzene ring derivative (HR) in their structures are presently used in clinics. High-resolution X-ray crystallography of the complexes of the IN-viral DNA transient complex bound to each of the three inhibitors showed in all cases the HR ring to interact within a confined zone of the viral DNA, limited to the highly conserved 5'CpA 3'/5'TpG 3' step. The extension of its extracyclic CX bond is electron-depleted, owing to the existence of the "sigma-hole." It interacts favorably with the electron-rich rings of base G4. We have sought to increase the affinity of HR derivatives for the G4/C16 base pair. We thus designed thirteen novel derivatives and computed their Quantum Chemistry (QC) intermolecular interaction energies (ΔE) with this base-pair. Most compounds had ΔE values significantly more favorable than those of the HR of the most potent halobenzene drug presently used in clinics, Dolutegravir. This should enable the improvement in a modular piece-wise fashion, the affinities of halogenated inhibitors for viral DNA (vDNA). In view of large scale polarizable molecular dynamics simulations on the entirety of the IN-vDNA-inhibitor complexes, validations of the SIBFA polarizable method are also reported, in which the evolution of each ΔE(SIBFA) contribution is compared to its QC counterpart along this series of derivatives.
ABSTRACT
To fight the increasingly worrying bacterial resistance to antibiotics, the discovery and development of new therapeutics is urgently needed. Here, we report on a new series of 1,2,4-triazole-3-thione compounds as inhibitors of metallo-ß-lactamases (MBLs), which represent major resistance determinants to ß-lactams, and especially carbapenems, in Gram-negative bacteria. These molecules are stable analogs of 4-amino-1,2,4-triazole-derived Schiff bases, where the hydrazone-like bond has been reduced (hydrazine series) or the 4-amino group has been acylated (hydrazide series); the synthesis and physicochemical properties thereof are described. The inhibitory potency was determined on the most clinically relevant acquired MBLs (IMP-, VIM-, and NDM-types subclass B1 MBLs). When compared with the previously reported hydrazone series, hydrazine but not hydrazide analogs showed similarly potent inhibitory activity on VIM-type enzymes, especially VIM-2 and VIM-4, with Ki values in the micromolar to submicromolar range. One of these showed broad-spectrum inhibition as it also significantly inhibited VIM-1 and NDM-1. Restoration of ß-lactam activity in microbiological assays was observed for one selected compound. Finally, the binding to the VIM-2 active site was evaluated by isothermal titration calorimetry and a modeling study explored the effect of the linker structure on the mode of binding with this MBL.
Subject(s)
Drug Resistance, Bacterial/drug effects , Gram-Negative Bacteria/drug effects , Thiones/chemistry , Triazoles/chemistry , beta-Lactamase Inhibitors/pharmacology , beta-Lactamases/metabolism , Anti-Bacterial Agents/chemistry , Anti-Bacterial Agents/pharmacology , Bacterial Infections/microbiology , Bacterial Infections/prevention & control , Biocatalysis/drug effects , Carbapenems/chemistry , Carbapenems/pharmacology , Gram-Negative Bacteria/metabolism , Humans , Microbial Sensitivity Tests , Molecular Structure , beta-Lactamase Inhibitors/chemistry , beta-Lactams/chemistry , beta-Lactams/pharmacologyABSTRACT
Using polarizable (AMOEBA) and nonpolarizable (CHARMM) force fields, we compare the relative free energy stability of two extreme conformations of the HIV-1 nucleocapsid protein NCp7 that had been previously experimentally advocated to prevail in solution. Using accelerated sampling techniques, we show that they differ in stability by no more than 0.75-1.9 kcal/mol depending on the reference protein sequence. While the extended form appears to be the most probable structure, both forms should thus coexist in water explaining the differing NMR findings.
Subject(s)
gag Gene Products, Human Immunodeficiency Virus/chemistry , Entropy , Magnetic Resonance Spectroscopy , Molecular Dynamics SimulationABSTRACT
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.
Subject(s)
Density Functional Theory , Glucose-6-Phosphate Isomerase/chemistry , Phosphates/chemistry , Anions/chemistry , Calibration , Glucose-6-Phosphate Isomerase/metabolismABSTRACT
We present Tinker-HP, a massively MPI parallel package dedicated to classical molecular dynamics (MD) and to multiscale simulations, using advanced polarizable force fields (PFF) encompassing distributed multipoles electrostatics. Tinker-HP is an evolution of the popular Tinker package code that conserves its simplicity of use and its reference double precision implementation for CPUs. Grounded on interdisciplinary efforts with applied mathematics, Tinker-HP allows for long polarizable MD simulations on large systems up to millions of atoms. We detail in the paper the newly developed extension of massively parallel 3D spatial decomposition to point dipole polarizable models as well as their coupling to efficient Krylov iterative and non-iterative polarization solvers. The design of the code allows the use of various computer systems ranging from laboratory workstations to modern petascale supercomputers with thousands of cores. Tinker-HP proposes therefore the first high-performance scalable CPU computing environment for the development of next generation point dipole PFFs and for production simulations. Strategies linking Tinker-HP to Quantum Mechanics (QM) in the framework of multiscale polarizable self-consistent QM/MD simulations are also provided. The possibilities, performances and scalability of the software are demonstrated via benchmarks calculations using the polarizable AMOEBA force field on systems ranging from large water boxes of increasing size and ionic liquids to (very) large biosystems encompassing several proteins as well as the complete satellite tobacco mosaic virus and ribosome structures. For small systems, Tinker-HP appears to be competitive with the Tinker-OpenMM GPU implementation of Tinker. As the system size grows, Tinker-HP remains operational thanks to its access to distributed memory and takes advantage of its new algorithmic enabling for stable long timescale polarizable simulations. Overall, a several thousand-fold acceleration over a single-core computation is observed for the largest systems. The extension of the present CPU implementation of Tinker-HP to other computational platforms is discussed.
ABSTRACT
In the context of the SIBFA polarizable molecular mechanics/dynamics (PMM/PMD) procedure, we report the calibration and a series of validation tests for the 1,2,4-triazole-3-thione (TZT) heterocycle. TZT acts as the chelating group of inhibitors of dizinc metallo-ß-lactamases (MBL), an emerging class of Zn-dependent bacterial enzymes, which by cleaving the ß-lactam bond of most ß-lactam antibiotics are responsible for the acquired resistance of bacteria to these drugs. Such a study is indispensable prior to performing PMD simulations of complexes of TZT-based inhibitors with MBL's, on account of the anchoring role of TZT in the dizinc MBL recognition site. Calibration was done by comparisons to energy decomposition analyses (EDA) of high-level ab initio QC computations of the TZT complexes with two probes: Zn(II), representative of "soft" dications, and water, representative of dipolar molecules. We performed distance variations of the approach of each probe to each of the two TZT atoms involved in Zn ligation, the S atom and the N atom ortho to it, so that each SIBFA contribution matches its QC counterpart. Validations were obtained by performing in- and out-of-plane angular variations of Zn(II) binding in monoligated Zn(II)-TZT complexes. The most demanding part of this study was then addressed. How well does ΔE(SIBFA) and its individual contributions compare to their QC counterparts in the dizinc binding site of one MBL, L1, whose structure is known from high-resolution X-ray crystallography? Six distinct complexes were considered, namely each separate monozinc site, and the dizinc site, whether ligated or unligated by TZT. Despite the large magnitude of the interaction energies, in all six complexes ΔE(SIBFA) can match ΔE(QC) with relative errors <2% and the proper balance of individual energy contributions. The computations were extended to the dizinc site of another MBL, VIM-2, and its complexes with two other TZT analogues. ΔE(SIBFA) faithfully reproduced ΔE(QC) in terms of magnitude, ranking of the three ligands, and trends of the separate energy contributions. A preliminary extension to correlated calculations is finally presented. All these validations should enable a secure design of a diversity of TZT-containing MBL inhibitors: a structurally and energetically correct anchoring of TZT should enable all other inhibitor groups to in turn optimize their interactions with the other target MBL residues.
Subject(s)
Quantum Theory , Triazoles/chemistry , Zinc/chemistry , beta-Lactamase Inhibitors/chemistry , beta-Lactamases/chemistry , Calibration , Crystallography, X-Ray , Reproducibility of Results , Triazoles/pharmacology , Zinc/metabolism , beta-Lactamase Inhibitors/pharmacology , beta-Lactamases/metabolismABSTRACT
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.
ABSTRACT
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'.
Subject(s)
HIV Integrase Inhibitors/pharmacology , HIV Integrase/metabolism , Heterocyclic Compounds, 3-Ring/pharmacology , Keto Acids/pharmacology , DNA, Viral/drug effects , Dose-Response Relationship, Drug , Fluorescence Polarization , HIV Integrase Inhibitors/chemistry , Heterocyclic Compounds, 3-Ring/chemistry , Keto Acids/chemistry , Models, Molecular , Molecular Conformation , Oxazines , Piperazines , Pyridones , Structure-Activity RelationshipABSTRACT
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.
Subject(s)
G-Quadruplexes , Guanine/chemistry , Lithium/chemistry , Potassium/chemistry , Rubidium/chemistry , Sodium/chemistry , Alkalies/chemistry , Anisotropy , Cations, Monovalent/chemistry , Models, Molecular , ThermodynamicsABSTRACT
Roscovitine is an antitumor purine inhibitor of cyclin-dependent kinase CDK5, for which it displays submicromolar affinity. It reached phase IIb clinical trials in 2007. The search for analogues with improved kinase affinities led recently to an isomer, finisterine, having a nearly 10-fold greater affinity for both CDK5 and CDK9. It solely differs by the displacement of one nitrogen atom in the purine ring, from position 6 to position 9. This has no incidence on the intermolecular interaction of either drug with the neighboring sites that anchor the ring in the recognition site. Quantum chemistry calculations combined with conformational and topological analyses of the impact of the purine ring isomerization of roscovitine and finisterine on its conformational stability show that the modified electronic conjugation, on the other hand, results in a stiffening of the rotational barrier around the extracyclic C-NH bond of finisterine, vicinal to N9, and to which an aryl ring is connected, along with a loosening of the barrier around an extracyclic C6-C bond connecting to a shorter, hydrophobic arm. The first effect is proposed to lead to a lesser hydration entropy of solvation in the case of finisterine, thus to a facilitated desolvation term in the overall energy balances.
ABSTRACT
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.
Subject(s)
Hydroxamic Acids/chemistry , Metalloproteins/chemistry , Molecular Dynamics Simulation , Quantum Theory , Sugar Phosphates/chemistry , Zinc/chemistry , Binding Sites , Ligands , Metalloproteins/metabolism , Zinc/metabolismABSTRACT
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.