Journal of Chongqing University-Eng. Ed.
Vol. 2 No. 1 June 2003
Molecular modeling study of the effect of base methylation on internal
interactions and motions in DNA and implication to B-Z conformation change
*
CAI Congzhong
1,2
, WANG Wanlu
1
,
CHEN Yuzong
1,2
1
College of Mathematics and Sciences, Chongqing University, Chongqing 400044, P. R. China
2
Department of Computational Science, National University of Singapore, Singapore 117543
Received 18 November 2002; revised 24 February 2003
Abstract: Methylation in the bases of DNA is known to induce B-Z conformation change. In this work, molecular mechanics and
normal mode analysis are used to probe how certain methylation affects the internal interactions and thermodynamic motions in
the DNA double helixes in both B and Z conformations, and its implication to B-Z conformation change. By molecular modeling
with Insight II, two cases involving cytosine C5 and guanine C8 methylation on both B and Z-form DNA duplex d(CGCGCG)
2
are
studied in comparison with the corresponding unmethylated duplexes. The internal interaction energies computed based on a
molecular mechanics force field and the entropies due to internal motions computed according to a normal mode analysis are in
fare agreement with respective observed thermodynamic quantities. The analysis on the computed individual energy terms
suggests that the observed B-Z conformation change induced by methylation is primarily driven by enthalpic factors. A
combination of changes in Van der Waals interaction, electrostatic interaction and hydrogen bonding likely contributes to the
change of enthalpy that favors Z-conformation in the methylated states.
Keywords: molecular modeling; DNA base methylation; conformation change; entropy; enthalpy
CAI Congzhong (
蔡从中
): Male; Born 1966; PhD candidate; Asst. Prof.; Research interests: bioinformatics, biophysics, material
physics and scanning probe microscope.
*
Funded by the International Joint Research Project of Chongqing University and National University of Singapore (ARF-151-000-014-
112), and the Basic and Applied Research Foundation of Chongqing University.
1. Introduction
Experimental studies have shown that a DNA
molecule with alternating pyrimidine-purine sequence
can adopt a left-handed, double-helical Z-DNA
conformation [1,2]. Such structural changes of DNA
occur as a consequence of environmental conditions
such as at salt concentration above 4 mol/L NaCl [1,2]
or through certain chemical modification such as
methylation [1] or bromination [2,3] of bases at
physiological conditions. These structural changes and
chemical modifications are of importance in certain
biological processes and thus are subjects of a number
of investigations [1-8].
Molecular modeling has been separately used to
construct methylated DNA structures
[4,5], to analyze
the contribution of internal interactions to the stability
of DNA duplexes
[6,7], and to examine the
contribution of internal motions to B-Z transition [8].
The success of molecular modeling in these separate
studies have raised an interesting question as to
whether it can be used to probe the mechanism of the
observed effect of methylation on DNA B-Z confor-
mation change [1,2].
In the present work, in an attempt to probe possible
mechanism of the effect of cytosine C5 and guanine C8
methylation on B-Z transition in DNA respectively,
molecular modeling techniques are used to model the
3-dimensional structures of both normal and methylated
DNAs. Effects of the methylation of every previously
mentioned bases on the internal interactions and
internal motions in both B and Z forms of DNA are
analyzed to determine how they might affect the
internal energy and other thermodynamic quantities
which are important to determination of equilibrium
state between different conformations. The results are
used to assess the effect on B-Z conformation change.
2. Methods
2.1 Structures
The structures of both B- and Z-form d(CGCGCG)
2
were generated from Insight II [9] in which the
standard parameters for helical rise and twist angle
were used. The structures of both B and Z forms
d(m
5
CGCGm
5
CG)
2
and d(CGCm
8
GCG)
2
were
generated by adding methyl groups to the relevant
positions of the corresponding B and Z form