Generalized Poisson-Boltzmann calculation of the distribution of electrolyte ions around the B- and Z-conformers of DNA

A solution to the three-dimensional Poisson-Boltzmann equation, generalized to include the finite size of the ions, is presented for the environment of DNA in the B- and Z-conformations. The results clearly indicate that despite the lower linear charge density of the left-handed Z-conformer, there is a higher concentration of Na⁺ at the immediate surface of the Z-from than at the surface of the B-form. The average concentration of counterions within a 12-Å radius of the DNA is, nonetheless, higher for the B- than for the Z-form. Calculations of the electrostatic interactions of these conformers with an environment of 0.01M monovalent salt show that the salt exerts a greater stabilizing effect on the left-handed conformer than on the right-handed form.     

Differential electrostatic stabilization of A-, B-, and Z-forms of DNA

The static accessibility discrete charge algorithm for protein charge interactions is extended to the case of linear polyelectrolytes. In this model, the effective dielectric value between surface charge sites depends predominantly on the solvent ionic strength and the solvent accessibilities of the charge sites. This treatment accounts for the phenomena of specific ion binding in the context of a general electrostatic effect [Matthew and Richards (1982) Biochemistry 21 , 4989]. Specific ion sites are determined by locating areas of high electrostatic potential at the solvent interface of the macromolecule. At a given ionic strength the calculated potential at a site is taken to describe a binding constant and therefore the ion site occupancy. For a 20-base-pair fragment of B-DNA, net charge of ?40, 16 ion sites are indicated in the minor groove. The partial occupancy of each site increases from 0.2 to 0.5 as the ionic strength is increased from 0.01 to 0.50. Over the same range of ionic strength, the electrostatic free energy of this charge array is calculated to change from +0.6 to ?0.05 kcal/bp. Parallel behavior is predicted for A- and Z-DNA charge geometries. The most stable configuration, based on electrostatic criteria, at high ionic strength (I = 0.1–0.5) is that of Z-DNA. In this range, the ratio of “bound” sodium to phosphate is predicted to be less than 0.4.     

The effect of methylation of the 6 oxygen of guanine on the structure and stability of double helical DNA

The effect of methylation of the O-6 position of guanine in short segments of double helical DNA has been investigated by molecular mechanical simulations on the sequences d(CGCGCG)2, d(CGC[OMG]CG)2, d(CGT[OMG]CG)2, d(CGC[OMC]CG/(CGCGCG), d(CGC[OMG]CG/d(CGTGCG), d(CGCGAATTCGCG)2 and d(CGCGAATTC[OMG]CG)2. Guanines methylated at the O-6 position are found to form hydrogen bonds of roughly equal strength to cytosine and thymine. The optimum structure of these modified base pairs are not dramatically different from normal GC pairs, but both involve some bifurcation of the proton donors of cytosine (4NH2) or thymine (3NH) between the guanine N3 and O6 groups.     

Influence of sequence-dependent cytosine protonation and methylation on DNA triplex stability

To investigate cytosine protonation and its influence on the sequence-dependent thermal stability of DNA triplexes in detail, we have employed homo- and heteronuclear NMR experiments on specifically (15)N-labeled oligodeoxynucleotides that were designed to fold into intramolecular triple helices of the pyrimidine motif under appropriate conditions. These experiments reveal that cytosines in central positions of the triplex are significantly protonated even at neutral pH. However, semiprotonation points for individual cytosine bases as determined from pH-dependent measurements show considerable differences depending on their position. Thus, protonation is disfavored for adjacent cytosines or for cytosines at the triplex termini, resulting in a smaller contribution to the overall free energy of the triple helical system. In contrast, protonation of the base upon substitution of 5-methylcytosine for cytosine in the triplex third strand is only affected to a minor extent, and triplex stabilization by the methyl substituent is shown to primarily arise from stacking energies and/or hydrophobic effects.     

The effect of gamma radiation on DNA methylation

The effect of 60Co gamma radiation on DNA methylation was studied in four cultured cell lines. In all cases a dose-dependent decrease in 5-methylcytosine was observed at 24, 48, and 72 h postexposure to 0.5-10 Gy. Nuclear DNA methyltransferase activity decreased while cytoplasmic activity increased in irradiated (10 Gy) V79A03 cells as compared to controls. No DNA demethylase activity was detected in the nuclei of control or irradiated V79A03 cells. Additionally, gamma radiation resulted in the differentiation of C-1300 N1E-115 cells, a mouse neuroblastoma line, in a dose- and time-dependent manner. These results are consistent with the hypothesis that (1) genes may be turned on following radiation via a mechanism involving hypomethylation of cytosine and (2) radiation-induced hypomethylation results from decreased intranuclear levels of DNA methyltransferase.     

DNA methylation cues in nucleosome geometry,stability and unwrapping

Cytosine methylation at the 5-carbon position is an essential DNA epigenetic mark in many eukaryotic organisms. Although countless structural and functional studies of cytosine methylation have been reported, our understanding of how it influences the nucleosome assembly, structure, and dynamics remains obscure. Here, we investigate the effects of cytosine methylation at CpG sites on nucleosome dynamics and stability. By applying long molecular dynamics simulations on several microsecond time scale, we generate extensive atomistic conformational ensembles of full nucleosomes. Our results reveal that methylation induces pronounced changes in geometry for both linker and nucleosomal DNA, leading to a more curved, under-twisted DNA, narrowing the adjacent minor grooves, and shifting the population equilibrium of sugar-phosphate backbone geometry. These DNA conformational changes are associated with a considerable enhancement of interactions between methylated DNA and the histone octamer, doubling the number of contacts at some key arginines. H2A and H3 tails play important roles in these interactions, especially for DNA methylated nucleosomes. This, in turn, prevents a spontaneous DNA unwrapping of 3–4 helical turns for the methylated nucleosome with truncated histone tails, otherwise observed in the unmethylated system on several microseconds time scale.     

DNA methylation polymorphism and stability in Chinese indica hybrid rice

Conventional rice breeding has long focused on exploiting the DNA sequence diversity. However, epigenetic diversity, reflected particularly in DNA methylation, can also contribute to phenotypic variation and should not be overlooked in rice breeding. In this study, 20 parental lines of indica rice, which are widely used in hybrid rice breeding in China, were analyzed to investigate variations of DNA methylation and its inheritance. The results revealed a wide diversity in DNA methylation among these breeding lines. A positive correlation was seen between DNA methylation and genetic diversity. Furthermore, some of the methylated DNA was inherited in the subsequent generation, regardless of whether they were produced by selfing or hybrid-crossing. This study provides insight into the methylation patterns in rice, and suggests the importance of epigenetic diversity in rice breeding.     

The effect of 5-azadeoxycytidine on cell growth and DNA methylation

By growing cells in the presence of 3 mM thymidine and 5-azadeoxycytidine up to 20% of DNA cytosines have been substituted with azacytosine. No substitution was obtained on incubating with 5-methyldeoxycytidine. Azacytosine-substituted DNA has a very low level of 5-methylcytosine and cells, which survive azadeoxycytidine treatments maintain this low level of methylation in the absence of the drug. The DNA of such cells is undermethylated fairly evenly in all classes of DNA e.g., satellite and unique DNA. Incubation of cells in azadeoxycytidine leads to high cell mortality which is not related to the lack of DNA methylation but may be linked to the altered interactions of proteins with the substituted DNA. This effect, rather than reduced DNA methylation, may be the cause of differentiative changes observed on treatment of cells with 5-azacytidine.     

In silico study on the effect of surface lysines and arginines on the electrostatic interactions and protein stability

Charged amino acids are mostly exposed on a protein surface, thereby forming a network of interactions with the surrounding amino acids as well as with water. In particular, positively charged arginine and lysine have different side chain geometries and provide a different number of potential electrostatic interactions. This study reports a comparative analysis of the difference in the number of two representative electrostatic interactions, such as salt-bridges and hydrogen bonds, contributed by surface arginine and lysine, as well as their effect on protein stability using molecular modeling and dynamics simulation techniques. Two in silico variants, the R variant with all arginines and the K variant with all lysines on the protein surface, were modeled by mutating all the surface lysines to arginines and the surface arginines to lysines, respectively, for each of the 10 model proteins. A structural comparison of the respective two variants showed that the majority of R variants possessed more salt-bridges and hydrogen bond interactions than the K variants, indicating that arginine provides a higher probability of electrostatic interactions than lysine owing to its side chain geometry. Molecular dynamics simulations of these variants revealed the R variants to be more stable than the K variants at room temperature but this effect was not prominent under protein denaturating conditions, such as 353 and 333 K with 8 M urea. These results suggest that the arginine residues on a protein surface contribute to the protein stability slightly more than lysine by enhancing the electrostatic interactions.     

The effect of HhaI methylation on DNA local structure.

The effect of methylating the central cytosine in the sequence GCGC on the sensitivity to digestion by deoxyribonuclease I was investigated. Examination of five such sites shows that the extra methyl group renders the bond on its 5' side more susceptible to cleavage, and this is interpreted by suggesting that it causes small changes in the local DNA phosphate orientation.     

Electrostatic interactions in proteins: Calculations of the electrostatic term of free energy and the electrostatic potential field

The pH-dependence of the electrostatic energy of interactions between titratable groups is calculated for some well studied globular proteins: basic pancreatic trypsin inhibitor, sperm whale myoglobin and tuna cytochrome c. The calculations are carried out using a semi-empirical appraach in terms of the macroscopic model based on the Kirkwood-Tanford theory. The results are discussed in the light of their physicochemical and biological properties. It was found that the pH-dependence of the electrostatic energy correlates with the III–IV transition of cytochrome c. The electrostatic field of the cysteine proteinase inhibitor, cystatin, was calculated in two ways. In the first one, the electrostatic field created by the pH dependent charges of the ionizable groups and peptide dipoles was calculated using the approach proposed. In the second one, the finite-difference method was used. The results obtained by the two methods are in overall agreement. The calculated field was discussed in terms of the binding of cystatin to papain.     

DNA methylation analysis reveals the effect of arsenic on gestational diabetes mellitus

BackgroundArsenic (As) exposure is one of the risk factors for gestational diabetes mellitus (GDM). This study aimed to explore the effect of As-exposure on DNA methylation in GDM and to establish a risk assessment model of GDM in As exposed pregnant women.MethodWe collected elbow vein blood of pregnant women before delivery to measure As concentration and DNA methylation data. Then compared the DNA methylation data and established a nomogram.ResultWe identified a total of 10 key differentially methylated CpGs (DMCs) and found 6 corresponding genes. Functions were enriched in Hippo signaling pathway, cell tight junction, prophetic acid metabolism, ketone body metabolic process, and antigen processing and presentation. A nomogram was established that can predict GDM risks (c-index = 0.595, s:p = 0.973).ConclusionWe found 6 genes associated with GDM with high As exposure. The prediction of the nomograms has been proven to be effective.