“All in All It's Just Another Brick in the Wall” Cooperative Effects of

Jul 31, 2019 - The behavior of the structural parameters of DNA considering different level of methylation in CpG islands is studied by means of full-...
0 downloads 0 Views 5MB Size
Subscriber access provided by BUFFALO STATE

B: Biophysics; Physical Chemistry of Biological Systems and Biomolecules

“All in All It's Just Another Brick in the Wall” Cooperative Effects of Cytosine Methylation on DNA Structure and Dynamics Cécilia Hognon, Vanessa Besancenot, Arnaud Gruez, Stephanie Grandemange, and Antonio Monari J. Phys. Chem. B, Just Accepted Manuscript • DOI: 10.1021/acs.jpcb.9b05835 • Publication Date (Web): 31 Jul 2019 Downloaded from pubs.acs.org on August 2, 2019

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 19 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

“All in All It’s Just Another Brick in the Wall” Cooperative Effects of Cytosine Methylation on DNA Structure and Dynamics Cécilia Hognon,a Vanessa Besancenot,b Arnaud Gruez,c Stephanie Grandemange,b* Antonio Monaria* Université de Lorraine and CNRS, UMR 7019 LPCT, F-54000 Nancy, France. b Université de Lorraine and CNRS, UMR 7039 CRAN, F-54000 Nancy, France. C Université de Lorraine and CNRS, UMR 7356 IMOPA, F-54000 Nancy, France a

AUTHOR INFORMATION Corresponding Author [email protected] [email protected]

ABSTRACT The behavior of the structural parameters of DNA considering different level of methylation in CpG islands is studied by means of full-atom molecular dynamics simulations and electronic circular dichroism, both in an artificial model system and in a gene promoter sequence. It is demonstrated that methylation although intrinsically brings quite local perturbations may, if its level is high enough, induce cooperative effects that strongly modify the DNA backbone torsional parameters altering the helicity as compared to the non-methylated case. Since methylation of CpG island is correlated with the regulation of gene expression, understanding the structural modifications induced in DNA is crucial to characterize all the fine equilibria into play in epigenetics phenomena.

ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 2 of 19

INTRODUCTION The role of nucleic acids in storing the genetic information while allowing its faithful replication during cells mitosis, as well as its decoding and transcription culminating in protein synthesis, is widely recognized. The molecular bases underlying this phenomenon have been firmly established since the seminal discover of the DNA double helix structure in the fifties, setting the stage for the molecular genetic sciences.1 However, and in addition to the former, understanding the fine regulations taking place at cellular level also requires the comprehension of the molecular factors leading to the regulation of gene expressions2–7 and its variability as a result of the cellular cycle or as a response to stress or external stimuli.8,9 The latter represents the domain of epigenetic, i.e. the study of the different reversible modification of DNA composition or structure leading to the regulation of the gene expression.10–12 As such, epigenetics is related to one of the most fundamental processes in biology, while the deregulation in gene expression has been linked to the development of different debilitating diseases including cancers13–16 and neurodegeneration.17–20 Indeed, in a vast amount of cancer cell lines the global transcriptome21,22 and proteome23–25 are highly deregulated as compared to healthy ones, with the overexpression of oncogenes and the inhibition of tumor suppressor genes. Thus, the precise understanding of epigenetics mechanisms could allow the possibility to design specifically tailored epigenetics therapy, targeting in a personalized way the factors related to tumor resistance or progression. However, even if in principle promising, epigenetics therapies are plagued with a plethora of harmful side-effects, mainly due to a global deregulation of gene expressions, that seriously limit their applicability.26–28 From a molecular point of view epigenetics modulations are performed via a series of reversible chemical and structural modification of the DNA structure, particularly in correspondence of the gene regulating sequences.29 By far the most common and celebrated of such modification, as 2 ACS Paragon Plus Environment

Page 3 of 19 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

presented in Figure 1, is the methylation of cytosine (dC) nucleobases to produce 5-methyl-cytosine (5mC).30–34 The former chemical modification takes place mostly in guanine and cytosine rich regions of the genome known as CpG islands. Cytosine methylation is normally associated to gene silencing, i.e. to the repression of its expression, that can be related either to a modified interaction with the corresponding transcription factors or to a subsequent reorganization of the DNA spatial arrangement limiting the accessibility of the gene. Cytosine methylation level is maintained in cell via a rather complex interplay between specific enzymes such as the DNA methyltransferase (DNMT) class, and the DNA lesion repair machinery. Even if fundamental, DNA methylation is not the only epigenetics regulatory factor, other factors involving both nucleic acids and proteins such as chromatine remodeling, histone acetylation, or post-trasductional modifications have to be considered. The complexity of epigenetics phenomena and the subtle network of cross-talks and regulations between different factors necessary to assure its maintenance is one of the reasons behind the modest success of epigenetics therapies so far, since the latter have mostly targeted the DNMT enzymes without proper selectivity towards specific genes.35

3 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

A

Page 4 of 19



B



C

C*



C

5’ - CGCGCGCGCGCGCGCGCGCG - 3’ 3’ - GCGCGCGCGCGCGCGCGCGC - 5’



D

5’ - CGCGCGCGCGC*GCGCGCGCG - 3’ 3’ - GCGCGCGCGCGCGCGCGCGC - 5’



E 5’ - CGC*GC*GC*GC*GC*GC*GC*GC*GC*G - 3’



3’ - GCGCGCGCGCGCGCGCGCGC - 5’

F 5’ - CGC*GC*GC*GC*GC*GC*GC*GC*GC*G - 3’ 3’ - GC*GC*GC*GC*GC*GC*GC*GC*GC*GC - 5’



Figure 1. Molecular sketch of cytosine (A) and 5mC (B). The model strands used for the subsequent molecular simulations are also reported: un-methylated (C), singly-methylated (D), hemi-methylated (E), and fully methylated (F). Note that the position of 5mC is indicated in red as C*.



As already underlined the biological role of the DNA methylation is exerted also via the modification of the intrinsic structural properties of DNA double strands that may lead to chromatin rearrangements. As such it is fundamental to characterize precisely at a molecular and atomistic level the subtle modifications induced by methylation in DNA intrinsic dynamics and stability. Indeed, the passage from dC to 5mC may seam a rather innocent modification that should have only marginal effects on the DNA properties since it should not impact neither the p-stacking properties of the nucleobases nor the Watson and Crick hydrogen bond pairing.36,37 However, DNA inherent flexibility and its complex macromolecular structure may lead to the appearance of complicated conformational landscapes as the result of the presence of damages in its sequence.38 Furthermore, even slightly different lesions may give raise to totally different structural properties and conformational flexibility.39 Finally one has also to consider that in biological relevant situations DNA methylation takes place in or intergenic or intrageneic with CpG islands, leading to the accumulation of the chemical modifications in spatially close regions of the DNA sequence, that in turn may induce the 4 ACS Paragon Plus Environment

Page 5 of 19 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

coupling of the individual deformations leading to a cooperative and cumulative effect as it was evidenced in case of oxidative DNA cluster lesions.40,41 The role of DNA methylation, and the related change in the DNA mechanical properties, is indeed fundamental in assuring epigenetics regulation.42–47 Some important computational studies have been performed concerning from the one side the interaction between DNA and histones, and hence nucleosomal compaction and the role of methylation,48 and on the other side the interaction with specific enzymes and the sensitivity to methylation experienced by transcription factors.49 In addition, combined molecular force assay, single-molecule spectroscopy, and molecular dynamic simulations have pinpointed the role of methylation in affecting DNA strand separation.50 As a matter of fact methylation has also been shown to significantly alter the interaction of DNA with different materials such as graphene quantum dots.51 In this contribution we rely on extend full atom molecular dynamics (MD) to unravel at an atomistic level the effects of DNA methylation on both model poly(dCdG) strands (Figure 1) and a methylated consensus region extracted from the adenomatous polyposis coli (APC) gene promoter, whose methylation status has been largely described while its aberration is correlated to cancer development affecting different organs and including breast, lungs, and digestive tract.52–55 In addition electronic circular dichroism (ECD) is used to confirm the presence of non-negligible structural modification of DNA.56–62 This will allow to point out the fundamental effects of methylation on the native B-DNA structure and the alteration of its native structure, hence filling a gap in the understanding of the molecular, and in a certain instance mechanical, epigenetics factors. Our results clearly underline the fact that DNA methylation induces important conformational effects that are particularly captured by significant differences observed in the backbone structure and flexibility. Indeed, even though the effect of one single methylation is rather marginal the cooperativity in CpG islands induces modifications that are easily confirmed by the observed change of 5 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 6 of 19

the ECD profile. The differences in the intrinsic properties of methylated and non-methylated DNA sequences may help setting the base to understand on the one side the differential interactions with transcription factors and on the other one the perturbation of the coiling leading to chromatin remodeling.44,47 EXPERIMENTAL AND COMPUTATIONAL DETAILS To characterize the impact of 5-mC on the global structural and dynamical properties of DNA we performed classical molecular dynamics on a 20 base pairs model B-DNA double strand of general sequence poly-(dC)-(dG) and on the 17 base pairs consensus sequence of the APC promoting gene (Figure 1), following different methylation patterns, namely fully methylated, hemimethylated or singly-methylated. All the systems have been built and set-up using the NAB utility distributed in the AmberTools package. DNA was described using the amberf99 force field 63 including the bsc1 corrections.64,65 The force field for 5mC was parameterized coherently with the amber force field, following the usual Amber antechamber procedure. The Restricted Electrostatic Potential (RESP)66 procedure was used consistently, the ground state geometry of the model 5mC nucleotide was optimized at density functional theory (DFT) level using the standard 6-31G basis set and the B3LYP functional.67 (See SI) The different DNA strands have been solvated in a cubic water box described by TIP3P force field68,69 and K+ cations were added to ensure electroneutrality of the simulation box. The H Mass Repartition (HMR) algorithm has been used consistently to speed up the simulation and allow a more important sampling.70 HMR allows to increase the time step used to integrate Newton’s equations of motion to 4 fs by artificially scaling the mass of all non water hydrogen atoms from 1.008 Da to 3.024 Da. 6 ACS Paragon Plus Environment

Page 7 of 19 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

After an equilibration and thermalization of 36 ns, 500 ns for the model systems bearing different methylation levels and 200 ns for the APC promoting gene have been run in the NPT thermodynamic ensemble at the temperature of 300 K and the pressure of 1 atm. All MD simulations have been performed using the NAMD 2.12 code,71 and results have been analyzed and visualized with VMD72 and Curves+ suites.73 Circular dichroism (CD) was performed out using a Chirascan CD from Applied Photophysics. All the data have been collected at 20 °C with 0.5 nm intervals in the wavelength range of 180−320 nm, using a temperature-controlled chamber. A 0.01 cm cuvette containing 30 μL of DNA sample at 100 M was used for all the measurements. All the measurement were replicated at least three times, and sample spectra were corrected for buffer background by subtracting the average spectrum of buffer alone.

Figure 2. Time series of the RMSD for the different model DNA strands.



7 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 8 of 19

RESULTS AND DISCUSSION In order to asses the effects of methylation on the DNA structural parameters in a systematic and rather unbiased way we first compare the behavior of a native alternating poly-(dC)-(dG) strand with different degrees of methylation (Figure 1). In the case of the non-methylated DNA we can observe a dynamic totally coherent with the one exhibited by B-DNA double strands with a rather small root mean square deviation (RMSD) as reported in Figure 2. The inclusion of only one 5mC positioned at the center of the strand, as expected, does not alter the RMSD, whose time series is almost superposed with the previous one, since the modification induced are of low magnitude and more importantly extremely local. On the contrary some differences appear when one consider the hemi-methylated model, i.e. the system in which one strand is totally methylated while the other one is left unchanged, is considered. As expected the discrepancy is also increased when a fully methylated DNA is considered. From the analysis of a global index such as the RMSD one may already infer the cumulative effects induced by methylation leading to an accumulation of structural deformation. However, more specific and local parameters are necessary to proper characterize the intrinsic behavior of poly-methylated DNA and more importantly its difference compared to the native B-DNA. As confirmed by our simulations the methyl group in 5mC is directed directly into the major groove of DNA. Hence, we expect to observe important modifications of the major groove parameters depth upon methylation.

8 ACS Paragon Plus Environment

Page 9 of 19 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Figure 3. Time series of the major groove width in correspondence of the different bases for the 4 DNA model strands as defined in Figure 1. A) native DNA, b) singly-methylated, C) hemi-methylated, D) fully-methylated.



A



C O5’

P



O3’ C3’



B



G G C 5mC



O5’ P O3’ C3’



O5’ P

Unmethylated Methylated



Figure 4. A) Definition of the torsional angle z. B) Superposition of the backbone of a DNA strand in correspondence of methylated and non-methylated sites showing the different orientation of z. C) Superposition of methylated and nonmethylated DNA strand highlighting the global effects on the DNA backbone.

9 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 10 of 19

The time series along the dynamics of this parameter for the different model strands are reported in Figure 3. Indeed, an increase of the major groove depth can be observed, coherently with the presence of a more sterical hindered group and of backbone reorganization. Interestingly enough, once again this modification is quite local as observed for the single methylated strand in which the corresponding increase in the depth value is observed only in correspondence of 5mC. While the other structural parameters describing the base coupling are almost unaltered upon methylation, as can be seen from the data reported in SI, the backbone torsional parameters are much more affected and globally induce an important deviation of the backbone helicity from the one of the native B-DNA as illustrated pictorially in Figure 4. In particular it appears that the most affected parameters is the one phosphate torsional angle z defined in Figure 4A, that leads to a different orientation of the sugar and phosphate moieties (Figure 4B) and reflects on a global distortion of the helicity (Figure 4C). A

B

C

D

Figure 5. Distribution of the average z angle over all the bases of a single strand for the un-methylated DNA (A), singlemethylated DNA (B), hemi-methylated DNA (C), and fully methylated DNA (D). In case of hemi-methylated and singlymethylated DNA the strand bearing 5mC is represented in red, while the ones having dC in black.



10 ACS Paragon Plus Environment

Page 11 of 19 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

In Figure 5 we report the distribution for the average value of the torsional angle z over all the sequence. We may observe that while in the case of the native B-DNA the value is centered around -100° while a large distribution peaking around 70° is observed for the methylated case. Interestingly in the case of the hemy-methylated system a very different behavior for the two strands is observed with the non-methylated one experiences a very sharp peak around -100° while the distribution for the methylated strand is wide and again centered around 70°. Interestingly in the case of a singly-methylated strand a double population emerges with two peaks at -100° and 100°. On the other hand no significant modification for the other parameters such as the torsional angle e is observed.

C

A

5’ - CGGCGCACGTGACCGCG - 3’ 3’ - GCCGCGTGCACTGGCGC- 5’

B

5’ - CGGC*GCAC*GTGACC*GCG - 3’ 3’ - GCC*GC*GTGCACTGGCGC- 5’

D

Figure 6. Sequence of the un-methylated (A) and methylated (B) APC promoting region. Distribution of the average torsional z angle for the un-methylated (C) and methylated (D) sequence, respectively.



The MD simulations on the model DNA strand have clearly confirmed that methylation induces cooperative structural modifications that strongly alter the helical parameters. However, those conclusions have been reported based on a model in which the methylation level is artificially high and hence could suffer from an over interpretation. In order to analyze a more biological relevant model we considered a fragment of the promoting sequence of the APC gene (Figure 6) that is composed of seventeen base pairs and harbors two CpG islands per strand. The same protocol as for the model system has been used, i.e.

11 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 12 of 19

considering non-methylated, hemi-methylated, and fully-methylated strands. Importantly, the same strands have also been experimentally produced and considered for spectroscopical characterization From the analysis of the MD trajectories (Figure 6) we can infer that the behavior underlined for the model system is indeed confirmed and important modifications of the backbone torsional angle are present upon methylation or hemi-methylation involving the relevant CpG islands. In particular, and coherently with what was observed for the poly-(dC)-poly(dG) strands methylation is accompanied by the inversion of the torsions z and the increase of its distribution width and hence by the corresponding cumulative effects on the backbone properties. One technique of choice to highlight possible structural deformations in biological macromolecules is indeed ECD, due to its sensitivity able to monitor even slight changes in global arrangement of organized chiral structures. Hence, we performed ECD for the promoting gene strands having different methylation level and results are collected in Figure 7. The helical structure of DNA is evidenced in the ECD spectra by the negative peaks at around 210 nm and 250 nm, and the positive peak around 280 nm that constitute a very well-known signature of the excitonic coupling between the B-DNA bases. Methylation, on the other hand, has only a very limited effect on the values of the absorption maxima as also confirmed by the absorption spectrum reported in SI. This result is expected since the modifications on the DNA nucleobases induced by methylation are too small to allow for a significant changes of the excitonic coupling strength. On the other hand, the intensity of the ECD positive peak at around 280 nm clearly increases with the methylation level, while the absolute value of the intensity of the negative peak at 250 nm follows a less clear trend. Indeed, the continuum of wavelength shifts between maxima (210 nm, 250 nm and 280 nm) of non-methylated, hemi-methylated, and fully-methylated strands also reflects a DNA structural change. Such spectral signatures are coherent with a partial modification of the helical parameters of the B-DNA structure upon methylation and with a slightly more important rigidity of methylated DNA, as evidenced by the RMSD. Indeed, ECD spectra is particularly sensitive to modifications happening in macromolecules secondary structure and hence to all the structural parameters of DNA affecting its 12 ACS Paragon Plus Environment

Page 13 of 19

helicity. In the present case the inversion of the backbone angles provoked by methylation and evidenced in the MD simulations may correlate well with the observed spectral signature.

50

Unmethylated Hemimethylated (5’−3’) Hemimethylated (3’−5’) Methylated

40

Δε/n x 104 (M−1.cm−1.n−1)

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

30 20 10 0

−10

200

220

240 260 280 Wavelength λ (nm)

300

320

Figure 7. Experimental ECD for the un-methylated, hemi-methylated, and fully-methylated APC promoting region.

CONCLUSIONS By resorting to extended full atom MD simulations we have clearly demonstrated that even if the inclusion of one single 5mC has little to no effect on the global structure, extended methylation of CpG islands clearly correlates with cooperative effects and hence significantly alters the DNA backbone torsional parameters, especially the z angle, in addition to modifying the major groove accessibility. The results have been confirmed both in model systems consisting of poly-(dC)-poly(dG) strands that have also allowed to underlined the differential effects related to either full- or hemy-methylation. Importantly, the same results have been confirmed in the case of a gene consensus sequence containing CpG islands whose methylation level has been varied accordingly. ECD on the gene sequence has also experimentally confirmed the alteration of the helicity compared to the ideal B-DNA upon methylation. Our work highlights, for the first time, the crucial and cooperative effects induced by methylation on the intrinsic structural parameters and DNA and hence can be fundamental in understanding the role played 13 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 14 of 19

by methylation in regulating gene expression by modifying the intrinsic DNA flexibility and hence, at a larger scale, allowing not only preferential interactions with gene promoting factors, but also in possibly driving chromatin reorganization that is known to be a crucial epigenetic marker. In addition, by highlighting the differential dynamics of methylated DNA sequences of biological relevance, our work may also be fundamental to allow for a better rationalization or development of novel drugs for epigenetic treatment having DNA as a target. In the future we plan to extend the present study considering DNA methylation in coiled DNA, either by histones in nucleosomes or by histones-like (HU). Furthermore, and to allow a better one-to-one mapping between simulation and experience, from the one hand ECD will be simulated resorting to hybrid quantum mechanical molecular mechanics (QM/MM) approaches and on the other hand the structural of methylated DNA will be explored using small-angle neutron scattering (SANS) techniques. ASSOCIATED CONTENT Time series of the torsional parameters for the different model systems and the APC promoting gene, absorption spectrum of the methylated, hemi-methylated and non-methylated APC gene, force field parameters for the 5mC nucleotide. The Supporting Information is available free of charge on the ACS Publications website. ACKNOWLEDGMENT The authors acknowledge the LUE IMPACT “Biomolécules” for financial supports. C.H. Ph.D. scholarship is funded by the Université de Lorraine. Molecular Modeling and simulation was performed on the LPCT local computing resources.

REFERENCES (1)

Watson, J. D.; Crick, F. H. C. Molecular Structure of Nucleic Acids: A Structure for Deoxyribose Nucleic Acid. Nature 1953, 171, 737–738. 14 ACS Paragon Plus Environment

Page 15 of 19 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

(2) (3) (4) (5) (6) (7) (8) (9) (10) (11) (12) (13) (14) (15) (16) (17) (18) (19) (20) (21) (22) (23) (24)

The Journal of Physical Chemistry

Le Novere, N. Quantitative and Logic Modelling of Molecular and Gene Networks. Nat. Rev. Genet. 2015, 16 (3), 146–158. Holoch, D.; Moazed, D. RNA-Mediated Epigenetic Regulation of Gene Expression. Nat. Rev. Genet. 2015, 16 (2), 71–84. Jonas, S.; Izaurralde, E. Towards a Molecular Understanding of MicroRNA-Mediated Gene Silencing. Nat. Rev. Genet. 2015, 16 (7), 421–433. Gibcus, J.; Dekker, J. The Context of Gene Expression Regulation. F1000 Biol. Rep. 2012, 4, 8. Aranda, S.; Mas, G.; Di Croce, L. Regulation of Gene Transcription by Polycomb Proteins. Sci. Adv. 2015, 1 (11), e1500737. Ennen, M.; Klotz, R.; Touche, N.; Pinel, S.; Barbieux, C.; Besancenot, V.; Brunner, E.; Thiebaut, D.; Jung, A. C.; Ledrappier, S.; et al. DDB2: A Novel Regulator of Nf-Kb and Breast Tumor Invasion. Cancer Res. 2013, 73 (16), 5040–5052. De Nadal, E.; Ammerer, G.; Posas, F. Controlling Gene Expression in Response to Stress. Nat. Rev. Genet. 2011, 12 (12), 833–845. Murray, J. I. Diverse and Specific Gene Expression Responses to Stresses in Cultured Human Cells. Mol. Biol. Cell 2004, 15 (5), 2361–2374. Weinhold, B. Epigenetics: The Science of Change. Environ. Health Perspect. 2006, 114 (3), A1607. Deans, C.; Maggert, K. A. What Do You Mean, “Epigenetic”? Genetics 2015, 199 (4), 887–896. Rozek, L. S.; Dolinoy, D. C.; Sartor, M. A.; Omenn, G. S. Epigenetics: Relevance and Implications for Public Health; Annual Reviews, 2014; Vol. 35. Dawson, M. A.; Kouzarides, T. Cancer Epigenetics: From Mechanism to Therapy. Cell 2012, 150 (1), 12–27. Fardi, M.; Solali, S.; Farshdousti Hagh, M. Epigenetic Mechanisms as a New Approach in Cancer Treatment: An Updated Review. Genes Dis. 2018, 5 (4), 304–311. Jones, P. A.; Ohtani, H.; Chakravarthy, A.; De Carvalho, D. D. Epigenetic Therapy in ImmuneOncology. Nat. Rev. Cancer 2019, 19 (3), 151–161. Nebbioso, A.; Tambaro, F. P.; Dell’Aversana, C.; Altucci, L. Cancer Epigenetics: Moving Forward. PLoS Genetics. Public Library of Science 2018, p e1007362. Hwang, J. Y.; Aromolaran, K. A.; Zukin, R. S. The Emerging Field of Epigenetics in Neurodegeneration and Neuroprotection. Nat. Rev. Neurosci. 2017, 18 (6), 347–361. Lardenoije, R.; Iatrou, A.; Kenis, G.; Kompotis, K.; Steinbusch, H. W. M.; Mastroeni, D.; Coleman, P.; Lemere, C. A.; Hof, P. R.; van den Hove, D. L. A.; et al. The Epigenetics of Aging and Neurodegeneration. Prog. Neurobiol. 2015, 131, 21–64. Berson, A.; Nativio, R.; Berger, S. L.; Bonini, N. M. Epigenetic Regulation in Neurodegenerative Diseases. Trends Neurosci. 2018, 41 (9), 587–598. Hadizadeh Esfahani, A.; Sverchkova, A.; Saez-Rodriguez, J.; Schuppert, A. A.; Brehme, M. A Systematic Atlas of Chaperome Deregulation Topologies across the Human Cancer Landscape. PLoS Comput. Biol. 2018, 14 (1), e1005890. Wang, W.; Zhan, L.; Guo, D.; Xiang, Y.; Zhang, Y. U.; Tian, M.; Han, Z. Transcriptome Analysis of Pancreatic Cancer Cell Response to Treatment with Grape Seed Proanthocyanidins. Oncol. Lett. 2019, 17 (2), 1741–1749. Dominguez, D.; Tsai, Y. H.; Gomez, N.; Jha, D. K.; Davis, I.; Wang, Z. A High-Resolution Transcriptome Map of Cell Cycle Reveals Novel Connections between Periodic Genes and Cancer. Cell Res. 2016, 26 (8), 946–962. Ananthi, S.; Lakshmi, C. N. P.; Atmika, P.; Anbarasu, K.; Mahalingam, S. Global Quantitative Proteomics Reveal Deregulation of Cytoskeletal and Apoptotic Signalling Proteins in Oral Tongue Squamous Cell Carcinoma. Sci. Rep. 2018, 8 (1), 1567. Pappá, K. I.; Lygirou, V.; Kontostathi, G.; Zoidakis, J.; Makridakis, M.; Vougas, K.; Daskalakis, G.; Polyzos, A.; Anagnou, N. P. Proteomic Analysis of Normal and Cancer Cervical Cell Lines 15 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

(25) (26) (27) (28) (29) (30) (31) (32) (33) (34) (35) (36) (37) (38) (39) (40) (41) (42) (43) (44)

Page 16 of 19

Reveals Deregulation of Cytoskeleton-Associated Proteins. Cancer Genomics and Proteomics 2017, 14 (4), 253–266. Pappa, K. I.; Kontostathi, G.; Makridakis, M.; Lygirou, V.; Zoidakis, J.; Daskalakis, G.; Anagnou, N. P. High Resolution Proteomic Analysis of the Cervical Cancer Cell Lines Secretome Documents Deregulation of Multiple Proteases. Cancer Genomics Proteomics 2017, 14 (6), 507–521. Sharma, S.; Kelly, T. K.; Jones, P. A. Epigenetics in Cancer. Carcinogenesis. March 13, 2009, pp 27–36. Yoo, C. B.; Jones, P. A. Epigenetic Therapy of Cancer: Past, Present and Future. Nat. Rev. Drug Discov. 2006, 5 (1), 37–50. Riley, L. B.; Anderson, D. W. Cancer Epigenetics. In Handbook of Epigenetics; American Cancer Society, 2011; Vol. 60, pp 521–534. Becker, P. B.; Workman, J. L. Nucleosome Remodeling and Epigenetics. Cold Spring Harb. Perspect. Biol. 2013, 5 (9), a017905–a017905. Breiling, A.; Lyko, F. Epigenetic Regulatory Functions of DNA Modifications: 5-Methylcytosine and Beyond. Epigenetics and Chromatin. BioMed Central 2015, p 24. Kumar, S.; Chinnusamy, V.; Mohapatra, T. Epigenetics of Modified DNA Bases: 5Methylcytosine and Beyond. Front. Genet. 2018, 9, 640. Oakeley, E. J. DNA Methylation Analysisa Review of Current Methodologies. Pharmacol. Ther. 1999, 84 (3), 389–400. Wyatt, G. R. Occurrence of 5-Methyl-Cytosine in Nucleic Acids. Nature 1950, 166 (4214), 237– 238. Hergersberg, M. Biological Aspects of Cytosine Methylation in Eukaryotic Cells. Experientia 1991, 47 (11–12), 1171–1185. Cheray, M.; Pacaud, R.; Hervouet, E.; Vallette, F.; Cartron, P.-F. DNMT Inhibitors in Cancer, Current Treatments and Future Promising Approach: Inhibition of Specific DNMT-Including Complexes. Epigenetic Diagnosis Ther. 2016, 1 (1), 37–48. Fonseca Guerra, C.; Bickelhaupt, F. M.; Snijders, J. G.; Baerends, E. J. The Nature of the Hydrogen Bond in DNA Base Pairs: The Role of Charge Transfer and Resonance Assistance. Chem. - A Eur. J. 1999, 5 (12), 3581–3594. Poater, J.; Swart, M.; Guerra, C. F.; Bickelhaupt, F. M. Selectivity in DNA Replication. Interplay of Steric Shape, Hydrogen Bonds, π-Stacking and Solvent Effects. Chem. Commun. 2011, 47 (26), 7326. Dumont, E.; Monari, A. Understanding DNA under Oxidative Stress and Sensitization: The Role of Molecular Modeling. Front. Chem. 2015, 3, 43. Dehez, F.; Gattuso, H.; Bignon, E.; Morell, C.; Dumont, E.; Monari, A. Conformational Polymorphism or Structural Invariance in DNA Photoinduced Lesions: Implications for Repair Rates. Nucleic Acids Res. 2017, 45 (7), 3654–3662. Bignon, E.; Gattuso, H.; Morell, C.; Dehez, F.; Georgakilas, A. G.; Monari, A.; Dumont, E. Correlation of Bistranded Clustered Abasic DNA Lesion Processing with Structural and Dynamic DNA Helix Distortion. Nucleic Acids Res. 2016, 44 (18), 8588–8599. Gattuso, H.; Durand, E.; Bignon, E.; Morell, C.; Georgakilas, A. G.; Dumont, E.; Chipot, C.; Dehez, F.; Monari, A. Repair Rate of Clustered Abasic DNA Lesions by Human Endonuclease: Molecular Bases of Sequence Specificity. J. Phys. Chem. Lett. 2016, 7 (19), 3760–3765. Thliveris, A.; Samowitz, W.; Matsunami, N.; Groden, J.; White, R. Demonstration of Promoter Activity and Alternative Splicing in the Region 5′ to Exon 1 of the APC Gene. Cancer Res. 1994, 54 (11), 2991–2995. Deng, G.; Song, G. A.; Pong, E.; Sleisenger, M.; Kim, Y. S. Promoter Methylation Inhibits APC Gene Expression by Causing Changes in Chromatin Conformation and Interfering with the Binding of Transcription Factor CCAAT-Binding Factor. Cancer Res. 2004, 64 (8), 2692–2698. Geiman, T. M.; Robertson, K. D. Chromatin Remodeling, Histone Modifications, and DNA 16 ACS Paragon Plus Environment

Page 17 of 19 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

(45) (46) (47) (48) (49) (50) (51) (52) (53) (54) (55)

(56) (57) (58) (59) (60)

(61) (62)

The Journal of Physical Chemistry

Methylation - How Does It All Fit Together? J. Cell. Biochem. 2002, 87 (2), 117–125. Kim, S.; Kaang, B. K. Epigenetic Regulation and Chromatin Remodeling in Learning and Memory. Exp. Mol. Med. 2017, 49 (1), e281–e281. Laird-Offringa, I. A.; Sanchez-Cespedes, M. 12 – Epigenetic Events in Lung Cancer: Chromatin Remodeling and DNA Methylation. In IASLC Thoracic Oncology; Content Repository Only!, 2018; pp 104-116.e5. Meehan, R. R.; Pennings, S.; Stancheva, I. Lashings of DNA Methylation, Forkfuls of Chromatin Remodeling. Genes Dev. 2001, 15 (24), 3231–3236. Portella, G.; Battistini, F.; Orozco, M. Understanding the Connection between Epigenetic DNA Methylation and Nucleosome Positioning from Computer Simulations. PLoS Comput. Biol. 2013, 9 (11), e1003354. Machado, A. C. D.; Zhou, T.; Rao, S.; Goel, P.; Rastogi, C.; Lazarovici, A.; Bussemaker, H. J.; Rohs, R. Evolving Insights on How Cytosine Methylation Affects Protein-DNA Binding. Brief. Funct. Genomics 2015, 14 (1), 61–73. Severin, P. M. D.; Zou, X.; Gaub, H. E.; Schulten, K. Cytosine Methylation Alters DNA Mechanical Properties. Nucleic Acids Res. 2011, 39 (20), 8740–8751. Rafiei, S.; Dadmehr, M.; Hosseini, M.; Kermani, H. A.; Ganjali, M. R. A Fluorometric Study on the Effect of DNA Methylation on DNA Interaction with Graphene Quantum Dots. Methods Appl. Fluoresc. 2019, 7 (2), 025001. Tsuchiya, T.; Tamura, G.; Sato, K.; Endoh, Y.; Sakata, K.; Jin, Z.; Motoyama, T.; Usuba, O.; Kimura, W.; Nishizuka, S.; et al. Distinct Methylation Patterns of Two APC Gene Promoters in Normal and Cancerous Gastric Epithelia. Oncogene 2000, 19 (32), 3642–3646. Jin, Z.; Tamura, G.; Tsuchiya, T.; Sakata, K.; Kashiwaba, M.; Osakabe, M.; Motoyama, T. Adenomatous Polyposis Coli (APC) Gene Promoter Hypermethylation in Primary Breast Cancers. Br. J. Cancer 2001, 85 (1), 69–73. Dulaimi, E.; Hillinck, J.; De Caceres, I. I.; Al-Saleem, T.; Cairns, P. Tumor Suppressor Gene Promoter Hypermethylation in Serum of Breast Cancer Patients. Clin. Cancer Res. 2004, 10 (18 I), 6189–6193. Brabender, J.; Usadel, H.; Danenberg, K. D.; Metzger, R.; Schneider, P. M.; Lord, R. V; Wickramasinghe, K.; Lum, C. E.; Park, J. M.; Salonga, D.; et al. Adenomatous Polyposis Coli Gene Promoter Hypermethylation in Non-Small Cell Lung Cancer Is Associated with Survival. Oncogene 2001, 20 (27), 3528–3532. Vorlíčková, M.; Kejnovská, I.; Bednářová, K.; Renčiuk, D.; Kypr, J. Circular Dichroism Spectroscopy of DNA: From Duplexes to Quadruplexes. Chirality 2012, 24 (9), 691–698. Pérez-Arnaiz, C.; Busto, N.; Leal, J. M.; García, B. New Insights into the Mechanism of the DNA/Doxorubicin Interaction. J. Phys. Chem. B 2014, 118 (5), 1288–1295. Vorlíčkova, M.; Kejnovská, I.; Sagi, J.; Renčiuk, D.; Bednářová, K.; Motlová, J.; Kypr, J. Circular Dichroism and Guanine Quadruplexes. Methods 2012, 57 (1), 64–75. Gattuso, H.; Spinello, A.; Terenzi, A.; Assfeld, X.; Barone, G.; Monari, A. Circular Dichroism of DNA G-Quadruplexes: Combining Modeling and Spectroscopy to Unravel Complex Structures. J. Phys. Chem. B 2016, 120 (12), 3113–3121. Lauria, A.; Bonsignore, R.; Terenzi, A.; Spinello, A.; Giannici, F.; Longo, A.; Almerico, A. M.; Barone, G. Nickel(II), Copper(II) and Zinc(II) Metallo-Intercalators: Structural Details of the DNA-Binding by a Combined Experimental and Computational Investigation. Dalton Trans. 2014, 43 (16), 6108–6119. Gattuso, H.; Assfeld, X.; Monari, A. Modeling DNA Electronic Circular Dichroism by QM/MM Methods and Frenkel Hamiltonian. Theor. Chem. Acc. 2015, 134 (3), 36. Norman, P.; Parello, J.; Polavarapu, P. L.; Linares, M. Predicting Near-UV Electronic Circular Dichroism in Nucleosomal DNA by Means of DFT Response Theory. Phys. Chem. Chem. Phys. 2015, 17 (34), 21866–21879. 17 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

(63) (64) (65) (66) (67) (68) (69) (70) (71) (72) (73)

Page 18 of 19

Cornell, W. D.; Cieplak, P.; Bayly, C. I.; Gould, I. R.; Merz, K. M.; Ferguson, D. M.; Spellmeyer, D. C.; Fox, T.; Caldwell, J. W.; Kollman, P. A. A Second Generation Force Field for the Simulation of Proteins, Nucleic Acids, and Organic Molecules. J. Am. Chem. Soc. 1995, 117 (19), 5179–5197. Ivani, I.; Dans, P. D.; Noy, A.; Pérez, A.; Faustino, I.; Hospital, A.; Walther, J.; Andrio, P.; Goñi, R.; Balaceanu, A.; et al. Parmbsc1: A Refined Force Field for DNA Simulations. Nat. Methods 2015, 13 (1), 55–58. Dans, P. D.; Ivani, I.; Hospital, A.; Portella, G.; González, C.; Orozco, M. How Accurate Are Accurate Force-Fields for B-DNA? Nucleic Acids Res. 2017, 45 (7), 4217–4230. Wang, J.; Cieplak, P.; Kollman, P. A. How Well Does a Restrained Electrostatic Potential (RESP) Model Perform in Calculating Conformational Energies of Organic and Biological Molecules? J. Comput. Chem. 2000, 21 (12), 1049–1074. Becke, A. B3LYP. J. Chem. Phys. 1993, 98, 5648–5652. Jorgensen, W. L.; Chandrasekhar, J.; Madura, J. D.; Impey, R. W.; Klein, M. L. Comparison of Simple Potential Functions for Simulating Liquid Water. J. Chem. Phys. 1983, 79 (2), 926–935. Mark, P.; Nilsson, L. Structure and Dynamics of the TIP3P, SPC, and SPC/E Water Models at 298 K. J. Phys. Chem. A 2001, 105 (43), 9954–9960. Hopkins, C. W.; Le Grand, S.; Walker, R. C.; Roitberg, A. E. Long-Time-Step Molecular Dynamics through Hydrogen Mass Repartitioning. J. Chem. Theory Comput. 2015, 11 (4), 1864– 1874. Phillips, J. C.; Braun, R.; Wang, W.; Gumbart, J.; Tajkhorshid, E.; Villa, E.; Chipot, C.; Skeel, R. D.; Kalé, L.; Schulten, K. Scalable Molecular Dynamics with NAMD. J. Comput. Chem. 2005, 26 (16), 1781–1802. Humphrey, W.; Dalke, A.; Schulten, K. VMD: Visual Molecular Dynamics. J. Mol. Graph. 1996, 14 (1), 33–38, 27–28. Lavery, R.; Moakher, M.; Maddocks, J. H.; Petkeviciute, D.; Zakrzewska, K. Conformational Analysis of Nucleic Acids Revisited: Curves+. Nucleic Acids Res. 2009, 37 (17), 5917–5929.

18 ACS Paragon Plus Environment

Page 19 of 19

TOC GRAPHICS

fects ive Ef ion Cooperat e Methylat osin By e r u of Cyt t c D A stru ON DN and EC lations u im S MD 50 40

Unmethylated Hemimethylated (5’−3’) Hemimethylated (3’−5’) Methylated

30 20

4

−1

−1 −1

De/n x 10 (M .cm .n )

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

10 0

−10

200

220

240 260 280 Wavelength l (nm)

300

320

19 ACS Paragon Plus Environment