The Recognition Pathway for the DNA Cytosine Methyltransferase M.HhaI,

A perspective on the molecular simulation of DNA from structural and functional aspects

As genetic material, DNA not only carries genetic information by sequence, but also affects biological functions ranging from base modification to replication, transcription and gene regulation through its structural and dynamic properties and variations. The motion and structural properties of DNA involved in related biological processes are also multi-scale, ranging from single base flipping to local DNA deformation, TF binding, G-quadruplex and i-motif formation, TAD establishment, compartmentalization and even chromosome territory formation, just to name a few. The sequence-dependent physical properties of DNA play vital role in all these events, and thus it is interesting to examine how simple sequence information affects DNA and the formation of the chromatin structure in these different hierarchical orders. Accordingly, molecular simulations can provide atomistic details of interactions and conformational dynamics involved in different biological processes of DNA, including those inaccessible by current experimental methods. In this perspective, which is mainly based on our recent studies, we provide a brief overview of the atomistic simulations on how the hierarchical structure and dynamics of DNA can be influenced by its sequences, base modifications, environmental factors and protein binding in the context of the protein-DNA interactions, gene regulation and structural organization of chromatin. We try to connect the DNA sequence, the hierarchical structures of DNA and gene regulation.

Single-molecule study on conformational dynamics of M.HhaI

We found that apo DNA methyltransferase M.HhaI under the physiological salt concentration does not possess the structure characterized by X-ray crystallography; instead, it interchanges between prefolded and unfolded states. Only after binding to the substrate, it transforms into a crystal-structure-like state. Flipping rates of its catalytic loop were directly measured.

Huge conformational rearrangements in M.HhaI were observed by a single-molecule study.

Role of Conformational Fluctuations of Protein toward Methylation in DNA by Cytosine-5-methyltransferase

Methylation of cytosine is the common epigenetic modification in genomes ranging from bacteria to mammals, and aberrant methylation leads to human diseases including cancer. Recognition of a cognate DNA sequence by DNA methyltransferases and flipping of a target base into the enzyme active site pocket are the key steps in DNA methylation. Using molecular dynamics simulations and enhanced sampling techniques here we elucidate the role of conformational fluctuations of protein and active or passive involvement of protein elements that mediate base flipping and formation of the closed catalytic complex. The free energy profiles for the flipping of target cytosine into the enzyme active site support the major groove base eversion pathway; and the results show that the closed state of enzyme increases the free energy barrier, whereas the open state reduces it. We found that the interactions of the key loop residues of protein with cognate DNA altered the protein motions, and modulation of protein fluctuations relates to the closed catalytic complex formation. Methylation of cytosine in the active site of the closed complex destabilizes the interactions of catalytic loop residues with cognate DNA and reduces the stability of the closed state. Our study provides microscopic insights on the base flipping mechanism coupled with enzyme's loop motions and provides evidence for the role of conformational fluctuations of protein in the enzyme-catalyzed DNA processing mechanism.

NMR-based investigations into target DNA search processes of proteins

To perform their function, transcription factors and DNA-repair/modifying enzymes must first locate their targets in the vast presence of nonspecific, but structurally similar sites on genomic DNA. Before reaching their targets, these proteins stochastically scan DNA and dynamically move from one site to another on DNA. Solution NMR spectroscopy provides unique atomic-level insights into the dynamic DNA-scanning processes, which are difficult to gain by any other experimental means. In this review, we provide an introductory overview on the NMR methods for the structural, dynamic, and kinetic investigations of target DNA search by proteins. We also discuss advantages and disadvantages of these NMR methods over other methods such as single-molecule techniques and biochemical approaches.

Metadynamics simulation study on the conformational transformation of HhaI methyltransferase: an induced-fit base-flipping hypothesis

DNA methyltransferases play crucial roles in establishing and maintenance of DNA methylation, which is an important epigenetic mark. Flipping the target cytosine out of the DNA helical stack and into the active site of protein provides DNA methyltransferases with an opportunity to access and modify the genetic information hidden in DNA. To investigate the conversion process of base flipping in the HhaI methyltransferase (M.HhaI), we performed different molecular simulation approaches on M.HhaI-DNA-S-adenosylhomocysteine ternary complex. The results demonstrate that the nonspecific binding of DNA to M.HhaI is initially induced by electrostatic interactions. Differences in chemical environment between the major and minor grooves determine the orientation of DNA. Gln237 at the target recognition loop recognizes the GCGC base pair from the major groove side by hydrogen bonds. In addition, catalytic loop motion is a key factor during this process. Our study indicates that base flipping is likely to be an “induced-fit” process. This study provides a solid foundation for future studies on the discovery and development of mechanism-based DNA methyltransferases regulators.

Enzyme-promoted base flipping controls DNA methylation fidelity

A quantitative understanding of how conformational transitions contribute to enzyme catalysis and specificity remains a fundamental challenge. A suite of biophysical approaches was used to reveal several transient states of the enzyme-substrate complexes of the model DNA cytosine methyltransferase M.HhaI. Multidimensional, transverse relaxation-optimized nuclear magnetic resonance (NMR) experiments show that M.HhaI has the same conformation with noncognate and cognate DNA sequences. The high-affinity cognatelike mode requires the formation of a subset of protein-DNA interactions that drive the flipping of the target base from the helix to the active site. Noncognate substrates lacking these interactions undergo slow base flipping, and fluorescence tracking of the catalytic loop corroborates the NMR evidence of a loose, nonspecific binding mode prior to base flipping and subsequent closure of the catalytic loop. This slow flipping transition defines the rate-limiting step for the methylation of noncognate sequences. Additionally, we present spectroscopic evidence of an intermediate along the base flipping pathway that has been predicted but never previously observed. These findings provide important details of how conformational rearrangements are used to balance specificity with catalytic efficiency.

Evolutionary transitions to new DNA methyltransferases through target site expansion and shrinkage

DNA-binding and modifying proteins show high specificity but also exhibit a certain level of promiscuity. Such latent promiscuous activities comprise the starting points for new protein functions, but this hypothesis presents a paradox: a new activity can only evolve if it already exists. How then, do novel activities evolve? DNA methyltransferases, for example, are highly divergent in their target sites, but how transitions toward novel sites occur remains unknown. We performed laboratory evolution of the DNA methyltransferase M.HaeIII. We found that new target sites emerged primarily through expansion of the original site, GGCC, and the subsequent shrinkage of evolved expanded sites. Variants evolved for sites that are promiscuously methylated by M.HaeIII [GG(^A/_T)CC and GGCGCC] carried mutations in ‘gate-keeper’ residues. They could thereby methylate novel target sites such as GCGC and GGATCC that were neither selected for nor present in M.HaeIII. These ‘generalist’ intermediates were further evolved to obtain variants with novel target specificities. Our results demonstrate the ease by which new DNA-binding and modifying specificities evolve and the mechanism by which they occur at both the protein and DNA levels.

Differential expression and cellular localization of novel isoforms of the tendon biomarker tenomodulin

Tenomodulin (Tnmd, also called Tendin) is classified as a type II transmembrane glycoprotein and is highly expressed in developing as well as in mature tendons. Along with scleraxis (scx), Tnmd is a candidate marker gene for tenocytes. Its function is unknown, but it has been reported to have anti-angiogenic properties. Results in a knockout mouse model did not substantiate that claim. It has homology to chondromodulin-I. Single nucleotide polymorphisms of TNMD have been associated with obesity, macular degeneration, and Alzheimer's disease in patients. In the present study, three Tnmd isoforms with deduced molecular weights of 20.3 (isoform II), 25.4 (isoform III), and 37.1 (isoform I) kDa were proposed and verified by Western blot from cells with green fluorescent protein-linked, overexpressed constructs, tissue, and by qPCR of isoforms from human tissues and cultured cells. Overexpression of each Tnmd isoform followed by immunofluorescence imaging showed that isoforms I and II had perinuclear localization while isoform III was cytoplasmic. Results of qPCR demonstrated differential expression of each Tnmd isoform in patient's specimens taken from flexor carpi radialis, biceps brachii, and flexor digitorum profundus tendons. Knockdown of Tnmd increased the expression of both scleraxis (scx) and myostatin, indicating a potential negative feedback loop between Tnmd and its regulators. Knockdown of all Tnmd isoforms simultaneously also reduced tenocyte proliferation. I-TASSER protein three-dimensional conformation modeling predictions indicated each Tnmd isoform had different structures and potential functions: isoform 1, modeled as a cytosine methyltransferase; isoform 2, a SUMO-1-like SENP-1 protease; and isoform 3, an α-syntrophin, plextrin homology domain scaffolding protein. Further functional studies with each Tnmd isoform may help us to better understand regulation of tenocyte proliferation, tendon development, response to injury and strain, as well as mechanisms in tendinoses. These results may indicate novel therapeutic targets in specific tenomodulin isoforms as well as treatments for tendon diseases.

Molecular drivers of base flipping during sequence-specific DNA methylation

One step at a time: Substrates containing nucleotide analogues lacking sequence-specific contacts to the C5 methyltransferase M.HhaI were used to probe the role of individual interactions in effecting conformational transitions during base flipping. A segregation of duties, that is, specific recognition and chemomechanical force for base flipping and active site assembly, within the enzyme is confirmed.

Dnmt1 structure and function

Dnmt1, the principal DNA methyltransferase in mammalian cells, is a large and a highly dynamic enzyme with multiple regulatory features that can control DNA methylation in cells. This chapter highlights how insights into Dnmt1 structure and function can advance our understanding of DNA methylation in cells. The allosteric site(s) on Dnmt1 can regulate processes of de novo and maintenance DNA methylation in cells. Remaining open questions include which molecules, by what mechanism, bind at the allosteric site(s) in cells? Different phosphorylation sites on Dnmt1 can change its activity or ability to bind DNA target sites. Thirty-one different molecules are currently known to have physical and/or functional interaction with Dnmt1 in cells. The Dnmt1 structure and enzymatic mechanism offer unique insights into those interactions. The interacting molecules are involved in chromatin organization, DNA repair, cell cycle regulation, and apoptosis and also include RNA polymerase II, some RNA-binding proteins, and some specific Dnmt1-inhibitory RNA molecules. Combined insights from studies of different enzymatic features of Dnmt1 offer novel ideas for development of drug candidates, and can be used in selection of promising drug candidates from more than 15 different compounds that have been identified as possible inhibitors of DNA methylation in cells.

Determinants of precatalytic conformational transitions in the DNA cytosine methyltransferase M.HhaI

The DNA methyltransferase M.HhaI is an excellent model for understanding how recognition of a nucleic acid substrate is translated into site-specific modification. In this study, we utilize direct, real-time monitoring of the catalytic loop position via engineered tryptophan fluorescence reporters to dissect the conformational transitions that occur in both enzyme and DNA substrate prior to methylation of the target cytosine. Using nucleobase analogues in place of the target and orphan bases, the kinetics of the base flipping and catalytic loop closure rates were determined, revealing that base flipping precedes loop closure as the rate-determining step prior to methyl transfer. To determine the mechanism by which individual specific hydrogen bond contacts at the enzyme-DNA interface mediate these conformational transitions, nucleobase analogues lacking hydrogen bonding groups were incorporated into the recognition sequence to disrupt the major groove recognition elements. The consequences of binding, loop closure, and catalysis were determined for four contacts, revealing large differences in the contribution of individual hydrogen bonds to DNA recognition and conformational transitions on the path to catalysis. Our results describe how M.HhaI utilizes direct readout contacts to accelerate extrication of the target base that offer new insights into the evolutionary history of this important class of enzymes.

Homology modeling and molecular dynamics simulations of HgiDII methyltransferase in complex with DNA and S-adenosyl-methionine: catalytic mechanism and interactions with DNA

M.HgiDII is a methyltransferase (MTase) from Herpetosiphon giganteus that recognizes the sequence GTCGAC. This enzyme belongs to a group of MTases that share a high degree of amino acid similarity, albeit none of them has been thoroughly characterized. To study the catalytic mechanism of M.HgiDII and its interactions with DNA, we performed molecular dynamics simulations with a homology model of M.HgiDII complexed with DNA and S-adenosyl-methionine. Our results indicate that M.HgiDII may not rely only on Glu119 to activate the cytosine ring, which is an early step in the catalysis of cytosine methylation; apparently, Arg160 and Arg162 may also participate in the activation by interacting with cytosine O2. Another residue from the catalytic site, Val118, also played a relevant role in the catalysis of M.HgiDII. Val118 interacted with the target cytosine and kept water molecules from accessing the region of the catalytic pocket where Cys79 interacts with cytosine, thus preventing water-mediated disruption of interactions in the catalytic site. Specific recognition of DNA was mediated mainly by amino acids of the target recognition domain, although some amino acids (loop 80-88) of the catalytic domain may also contribute to DNA recognition. These interactions involved direct contacts between M.HgiDII and DNA, as well as indirect contacts through water bridges. Additionally, analysis of sequence alignments with closely related MTases helped us to identify a motif in the TRD of M.HgiDII that may be relevant to specific DNA recognition.