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.

Tinker-HP: a massively parallel molecular dynamics package for multiscale simulations of large complex systems with advanced point dipole polarizable force fields.

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.

Polarizable molecular dynamics in a polarizable continuum solvent.

We present, for the first time, scalable polarizable molecular dynamics (MD) simulations within a polarizable continuum solvent with molecular shape cavities and exact solution of the mutual polarization. The key ingredients are a very efficient algorithm for solving the equations associated with the polarizable continuum, in particular, the domain decomposition Conductor-like Screening Model (ddCOSMO), which involves a rigorous coupling of the continuum with the polarizable force field achieved through a robust variational formulation and an effective strategy to solve the coupled equations. The coupling of ddCOSMO with nonvariational force fields, including AMOEBA, is also addressed. The MD simulations are feasible, for real-life systems, on standard cluster nodes; a scalable parallel implementation allows for further acceleration in the context of a newly developed module in Tinker, named Tinker-HP. NVE simulations are stable, and long-term energy conservation can be achieved. This paper is focused on the methodological developments, the analysis of the algorithm, and the stability of the simulations; a proof-of-concept application is also presented to attest to the possibilities of this newly developed technique.