How a single 5-methylation of cytosine regulates the recognition of C/EBPβ transcription factor: a molecular dynamic simulation study

Comprehensive analysis unveils altered binding kinetics of 5-/6-methylCytosine/adenine modifications in R2R3-DNA system

Recent studies have shown that DNA methylation is an important epigenetic marker. Two prominent forms are methylation of the C5 position of cytosine and methylation of the C6 position of adenine. Given the vital significance of DNA methylation, investigating the mechanisms that influence protein binding remains a compelling pursuit. This study used molecular dynamics simulations to investigate the binding patterns of R2R3 protein and four differentially methylated DNAs. The alanine scanning combined with interaction entropy method was used to identify key residues that respond to different methylation patterns. The order of protein binding ability to DNA is as follows: unmethylated DNA > A11 methylation (5'-A6mAC-3') (6m2A system) > A10 methylation (5'-6mAAC-3') (6m1A system) > both A10 and A11 methylation (5'-6mA6mAC-3') (6mAA system) > C12 methylation (5'-AA5mC-3') (5mC system). All methylation systems lead to the sixth α helix (H6) (residues D105 to L116) moving away from the binding interface, and in the 5mC and 6m1A systems, the third α helix (H3) (residues G54 to L65) exhibits a similar trend. When the positively charged amino acids in H3 and H6 move away from the binding interface, their electrostatic and van der Waals interactions with the negatively charged DNA are weakened. Structural changes induced by methylation contributed to the destabilization of the hydrogen bond network near the original binding site, except for the 6m2A system. Moreover, there is a positive correlation between the number of methylated sites and the probability of distorting the DNA structure. Our study explores how different methylation patterns affect binding and structural adaptability, and have implications for drug discovery and understanding diseases related to abnormal methylation.

Interaction preferences between protein side chains and key epigenetic modifications 5-methylcytosine, 5-hydroxymethycytosine and N6-methyladenine

Covalent modifications of standard DNA/RNA nucleobases affect epigenetic regulation of gene expression by modulating interactions between nucleic acids and protein readers. We derive here the absolute binding free energies and analyze the binding modalities between key modified nucleobases 5-methylcytosine (5mC), 5-hydroxymethylcytosine (5hmC) and N⁶-methyladenine (m⁶A) and all non-prolyl/non-glycyl protein side chains using molecular dynamics simulations and umbrella sampling in both water and methanol, the latter mimicking the low dielectric environment at the dehydrated nucleic-acid/protein interfaces. We verify the derived affinities by comparing against a comprehensive set of high-resolution structures of nucleic-protein complexes involving 5mC. Our analysis identifies protein side chains that are highly tuned for detecting cytosine methylation as a function of the environment and can thus serve as microscopic readers of epigenetic marks. Conversely, we show that the relative ordering of sidechain affinities for 5hmC and m⁶A does not differ significantly from those for their precursor bases, cytosine and adenine, respectively, especially in the low dielectric environment. For those two modified bases, the effect is more nuanced and manifests itself primarily at the level of absolute changes in the binding free energy. Our results contribute towards establishing a quantitative foundation for understanding, predicting and modulating the interactions between modified nucleic acids and proteins at the atomistic level.

Methylation of ERβ 5'-untranslated region attenuates its inhibitory effect on ERα gene transcription and promotes the initiation and progression of papillary thyroid cancer

Normal thyroid tissue displays a prevalent expression of ERβ than ERα, which drastically turns upside down in the initiation and progression of papillary thyroid cancer (PTC). The underlying molecular mechanism of this phenomenon remains unclear. Here, we demonstrated that ERα and ERβ were coexpressed in human thyroid tissues and cells. ERα mRNA (A-1) and ERβ mRNA (0N-1), transcribed from Promoter A of ERα gene and Promoter 0N of ERβ gene, respectively, were the major mRNA isoforms which mainly contributed to total ERα mRNA and total ERβ mRNA in human thyroid-derived cell lines and tissues. The expression levels of ERα mRNA (A-1) and total ERα mRNA were gradually increased, and those of ERβ mRNA (0N-1) and total ERβ mRNA were decreased by degree in the initiation and progression of PTC. No aberrant DNA methylation of ERα 5'-untranslated region was involved in its up-regulation; however, aberrant DNA methylation in Promoter 0N and Exon 0N of ERβ gene was found to be involved in its down-regulation in the initiation and progression of PTC. ERβ can repress ERα gene transcription via recruitment of NCoR and displacement of RNA polymerase II at the Sp1 site in ERα Promoter A-specific region in thyroid-derived cells. It is suggested that DNA methylation of CpG islands in Promoter 0N and Exon 0N of ERβ gene leads to a decreased ERβ gene expression, which attenuates its inhibitory effect on ERα gene transcription and results in an increased ERα gene expression, cell proliferation, initiation, and progression of PTC.

Molecular mechanism of methyl-dependent and spatial-specific DNA recognition of c-Jun homodimer

DNA methylation is important in regulation of gene expression and normal development because it alters the interplay between protein and DNA. Experiments have shown that a single 5-methylcytosine at different CpG sites (mCpG) might have different effects on specific recognition, but the atomistic origin and dynamic details are largely unclear. In this work, we investigated the mechanism of monomethylation at different CpG sites in the cognate motif and the cooperativity of full methylation. By constructing four models of c-Jun/Jun protein binding to the 5[Formula: see text]-XGAGTCA-3[Formula: see text] (X represents C or methylated C) motif, we characterized the dynamics of the contact interface using the all-atom molecular dynamics method. Free energy analysis of MM/GBSA suggests that regardless of whether the C12pG13 site of the bottom strand is methylated, the effects from mC25 of the top strand are dominant and can moderately enhance the binding by [Formula: see text] 31 kcal/mol, whereas mC12 showed a relatively small contribution, in agreement with the experimental data. Remarkably, we found that this spatial-specific influence was induced by different regulatory rules. The influence of the mC25 site is mainly mediated by steric hindrance. The additional methyl group leads to the conformational changes in nearby residues and triggers an obvious structural bending in the protein, which results in the formation of a new T-Asn-C triad that enhances the specific recognition of TCA half-sites. The substitution of the methyl group at the mC12 site of the bottom strand breaks the original H-bonds directly. Such changes in electrostatic interactions also lead to the remote allosteric effects of protein by multifaceted interactions but have negligible contributions to binding. Although these two influence modes are different, they can both fine-tune the local environment, which might produce remote allosteric effects through protein-protein interactions. Further analysis reveals that the discrepancies in these two modes are primarily due to their location. Moreover, when both sites are methylated, the major determinant of binding specificity depends on the context and the location of the methylation site, which is the result of crosstalk and cooperativity.

DNA methylation changes involved in the tumor increase in F2 males born to gestationally arsenite-exposed F1 male mice

Multigenerational adverse effects from the environment such as nutrition and chemicals are among important concerns in environmental health issues. Previously, we have found that arsenite exposure of only F0 females during their pregnancy increases hepatic tumors in the F2 males in C3H mice. In the current study, we investigated the association of DNA methylation with the hepatic tumor increase in the F2 males of the arsenite group. Reduced-representation bisulfite sequencing analysis newly identified that DNA methylation levels of regions around the transcriptional start sites of Tmem54 and Cd74 were decreased and the expression of these genes were significantly increased in the hepatic tumors of F2 males of the arsenite group. The associations between DNA methylation in these regions and gene expression changes were confirmed by treatment of murine hepatoma cell lines and hepatic stellate cell line with 5-aza-2'-deoxycytidine. Overexpression of Cd74 in Hepa1c1c7 cells increased Trib3 expression and suppressed the expression of tumor suppressor genes Id3 and Atoh8. Human database analysis using the Cancer Genome Atlas indicated that TMEM54, CD74, and TRIB3 were significantly increased and that ATOH8 was decreased in hepatocellular carcinoma. The data also showed that high expression of TMEM54 and TRIB3 and low expression of ATOH8 were associated with poor survival. These results suggested that an increase in Tmem54 and Cd74 expression via DNA methylation reduction was involved in the tumor increase in the F2 male offspring by gestational arsenite exposure of F0 females. This study also suggested that genes downstream of Cd74 were involved in tumorigenesis.

Molecular mechanisms of the protein-protein interaction-regulated binding specificity of basic-region leucine zipper transcription factors

It is well known that the DNA-binding specificity of transcription factors (TFs) is influenced by protein-protein interactions (PPIs). However, the underlying molecular mechanisms remain largely unknown. In this work, we adopted the cAMP-response element-binding protein (CREB) of the basic leucine zipper (bZIP) TF family as a model system, and a workflow of combined bioinformatics and molecular modeling analysis of protein-DNA interaction was tested. First, the multiple sequence alignment and SDPsite method were used to find potential bZIP family binding specificity determining positions (SDPs) within the protein-protein interaction region. Second, the mutation system was analyzed using molecular dynamics simulation. Molecular mechanics Poisson-Boltzmann surface area (MM/PBSA) free energy calculations confirmed the enhancement of the binding affinity of the mutation, which was in agreement with experimental results. The root mean square fluctuation (RMSF) and hydrogen bonding changes suggested an open and close protein dimerization process after the system was mutated, which resulted in the change of the hydrogen bonding of the protein-DNA interface and a slight conformational change. We believe that this work will contribute to understanding the protein-protein interaction-regulated binding specificity of bZIP transcription factors.