Molecular Mechanisms of DNA Replication and Repair Machinery: Insights from Microscopic Simulations

Modelling DNA damage-repair and beyond

The paper presents a review of mechanistic modelling studies of DNA damage and DNA repair, and consequences to follow in mammalian cell nucleus. We hypothesize DNA deletions are consequences of repair of double strand breaks leading to the modifications of genome that play crucial role in long term development of genetic inheritance and diseases. The aim of the paper is to review formation mechanisms underlying naturally occurring DNA deletions in the human genome and their potential relevance for bridging the gap between induced DNA double strand breaks and deletions in damaged human genome from endogenous and exogenous events. The model of the cell nucleus presented enables simulation of DNA damage at molecular level identifying the spectrum of damage induced in all chromosomal territories and loops. Our mechanistic modelling of DNA repair for double stand breaks (DSB), single strand breaks (SSB) and base damage (BD), shows the complexity of DNA damage is responsible for the longer repair times and the reason for the biphasic feature of mammalian cells repair curves. In the absence of experimentally determined data, the mechanistic model of repair predicts the in vivo rate constants for the proteins involved in the repair of DSB, SSB, and of BD.

Protein-Nucleic Acid Interactions for RNA Polymerase II Elongation Factors by Molecular Dynamics Simulations

RNA polymerase II (Pol II) forms a complex with elongation factors to proceed to the elongation stage of the transcription process. In this work, we studied the elongation factor SPT5 and explored the protein-nucleic acid interactions for the isolated systems of KOW1 and KOW4 domains of SPT5 with DNA and RNA, respectively. We performed molecular dynamics (MD) simulations using three commonly used force fields that are CHARMM c36m, AMBER ff14sb, and ff19sb. Simulations showed strong protein-nucleic acid interactions and low electrostatic binding free energies for all force fields used. RNA was found to be highly dynamic with all force fields, while DNA had relatively more stable conformations with the AMBER force fields compared to that with CHARMM. Furthermore, we performed MD simulations of the complete elongation complex using CHARMM c36m and AMBER ff19sb force fields to compare the dynamics and interactions with the isolated systems. Similarly, strong KOW1 and DNA interactions were observed in the complete elongation complex simulations and DNA was further stabilized by a network of interactions involving SPT5-KOW1, SPT4, and rpb2 of Pol II. Overall, our study showed that the differences between CHARMM and AMBER force fields strongly affect the dynamics of the nucleic acids. CHARMM provides highly flexible DNA, while AMBER largely stabilizes the DNA structure. Although the presence of the entire interaction network stabilized the DNA and decreased the differences in the results from the two force fields, the discrepancies of the force fields for smaller systems may reflect their problems in generating accurate dynamics of nucleic acids.

Single-molecule biophysics experiments in silico: Toward a physical model of a replisome

The interpretation of single-molecule experiments is frequently aided by computational modeling of biomolecular dynamics. The growth of computing power and ongoing validation of computational models suggest that it soon may be possible to replace some experiments outright with computational mimics. Here, we offer a blueprint for performing single-molecule studies in silico using a DNA-binding protein as a test bed. We demonstrate how atomistic simulations, typically limited to sub-millisecond durations and zeptoliter volumes, can guide development of a coarse-grained model for use in simulations that mimic single-molecule experiments. We apply the model to recapitulate, in silico, force-extension characterization of protein binding to single-stranded DNA and protein and DNA replacement assays, providing a detailed portrait of the underlying mechanics. Finally, we use the model to simulate the trombone loop of a replication fork, a large complex of proteins and DNA.

•

Coarse-grained model derived from all-atom simulation recapitulates experiments

•

Model reproduces the elastic response to force and exchange dynamics

•

Model reveals structure of intermediate states usually inaccessible to experiment

•

Model applied to viral replisome with trombone loop containing tens of SSB proteins

Antitumor activity of novel POLA1-HDAC11 dual inhibitors

Hybrid molecules targeting simultaneously DNA polymerase α (POLA1) and histone deacetylases (HDACs) were designed and synthesized to exploit a potential synergy of action. Among a library of screened molecules, MIR002 and GEM144 showed antiproliferative activity at nanomolar concentrations on a panel of human solid and haematological cancer cell lines. In vitro functional assays confirmed that these molecules inhibited POLA1 primer extension activity, as well as HDAC11. Molecular docking studies also supported these findings. Mechanistically, MIR002 and GEM144 induced acetylation of p53, activation of p21, G1/S cell cycle arrest, and apoptosis. Oral administration of these inhibitors confirmed their antitumor activity in in vivo models. In human non-small cancer cell (H460) xenografted in nude mice MIR002 at 50 mg/kg, Bid (qd × 5 × 3w) inhibited tumor growth (TGI = 61%). More interestingly, in POLA1 inhibitor resistant cells (H460-R9A), the in vivo combination of MIR002 with cisplatin showed an additive antitumor effect with complete disappearance of tumor masses in two animals at the end of the treatment. Moreover, in two human orthotopic malignant pleural mesothelioma xenografts (MM473 and MM487), oral treatments with MIR002 and GEM144 confirmed their significant antitumor activity (TGI = 72-77%). Consistently with recent results that have shown an inverse correlation between POLA1 expression and type I interferon levels, MIR002 significantly upregulated interferon-α in immunocompetent mice.

Novel adamantyl retinoid-related molecules with POLA1 inhibitory activity

Atypical retinoids (AR) or retinoid-related molecules (RRMs) represent a promising class of antitumor compounds. Among AR, E-3-(3'-adamantan-1-yl-4'-hydroxybiphenyl-4-yl)acrylic acid (adarotene), has been extensively investigated. In the present work we report the results of our efforts to develop new adarotene-related atypical retinoids endowed also with POLA1 inhibitory activity. The effects of the synthesized compounds on cell growth were determined on a panel of human and hematological cancer cell lines. The most promising compounds showed antitumor activity against several tumor histotypes and increased cytotoxic activity against an adarotene-resistant cell line, compared to the parent molecule. The antitumor activity of a selected compound was evaluated on HT-29 human colon carcinoma and human mesothelioma (MM487) xenografts. Particularly significant was the in vivo activity of the compound as a single agent compared to adarotene and cisplatin, against pleural mesothelioma MM487. No reduction of mice body weight was observed, thus suggesting a higher tolerability with respect to the parent compound adarotene.

Molecular dynamics simulations of DNA-DNA and DNA-protein interactions

The all-atom molecular dynamics method can characterize the molecular-level interactions in DNA and DNA-protein systems with unprecedented resolution. Recent advances in computational technologies have allowed the method to reveal the unbiased behavior of such systems at the microseconds time scale, whereas enhanced sampling approaches have matured enough to characterize the interaction free energy with quantitative precision. Here, we describe recent progress toward increasing the realism of such simulations by refining the accuracy of the molecular dynamics force field, and we highlight recent application of the method to systems of outstanding biological interest.