Observing an Induced-fit Mechanism during Sequence-specific DNA Methylation

Base Dynamics in the HhaI Protein Binding Site

Protein-DNA interactions play an important role in numerous biological functions within the living cell. In many of these interactions, the DNA helix is significantly distorted upon protein-DNA complex formation. The HhaI restriction-modification system is one such system, where the methylation target is flipped out of the helix when bound to the methyltransferase. However, the base flipping mechanism is not well understood. The dynamics of the binding site of the HhaI methyltransferase and endonuclease (underlined) within the DNA oligomer [d(G1A2T3A4G5C6G7C8T9A10T11C12)]2 are studied using deuterium solid-state NMR (SSNMR). SSNMR spectra obtained from DNAs deuterated on the base of nucleotides within and flanking the [5'-GCGC-3']2 sequence indicate that all of these positions are structurally flexible. Previously, conformational flexibility within the phosphodiester backbone and furanose ring within the target sequence has been observed and hypothesized to play a role in the distortion mechanism. However, whether that distortion was occurring through an active or passive mechanism remained unclear. These NMR data demonstrate that although the [5'-GCGC-3']2 sequence is dynamic, the target cytosine is not passively flipping out of the double-helix on the millisecond-picosecond time scale. Additionally, although previous studies have shown that both the furanose ring and phosphodiester backbone experience a change in dynamics upon methylation, which may play a role in recognition and cleavage by the endonuclease, our observations here indicate that methylation has no effect on the dynamics of the base itself.

The Arabidopsis (ASHH2) CW domain binds monomethylated K4 of the histone H3 tail through conformational selection

Chromatin post-translational modifications are thought to be important for epigenetic effects on gene expression. Methylation of histone N-terminal tail lysine residues constitutes one of many such modifications, executed by families of histone lysine methyltransferase (HKMTase). One such protein is ASHH2 from the flowering plant Arabidopsis thaliana, equipped with the interaction domain, CW, and the HKMTase domain, SET. The CW domain of ASHH2 is a selective binder of monomethylation at lysine 4 on histone H3 (H3K4me1) and likely helps the enzyme dock correctly onto chromatin sites. The study of CW and related interaction domains has so far been emphasizing lock-key models, missing important aspects of histone-tail CW interactions. We here present an analysis of the ASHH2 CW-H3K4me1 complex using NMR and molecular dynamics, as well as mutation and affinity studies of flexible coils. β-augmentation and rearrangement of coils coincide with changes in the flexibility of the complex, in particular the η1, η3 and C-terminal coils, but also in the β1 and β2 strands and the C-terminal part of the ligand. Furthermore, we show that mutating residues with outlier dynamic behaviour affect the complex binding affinity despite these not being in direct contact with the ligand. Overall, the binding process is consistent with conformational selection. We propose that this binding mechanism presents an advantage when searching for the correct post-translational modification state among the highly modified and flexible histone tails, and also that the binding shifts the catalytic SET domain towards the nucleosome. DATABASES: Structural data are available in the PDB database under the accession code 6QXZ. Resonance assignments for CW42 in its apo- and holo-forms are available in the BMRB database under the accession code 27251.

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.

Structure-guided residence time optimization of a dabigatran reversal agent

Novel oral anticoagulants are effective and safe alternatives to vitamin-K antagonists for anticoagulation therapy. However, anticoagulation therapy in general is associated with an elevated risk of bleeding. Idarucizumab is a reversal agent for the direct thrombin inhibitor, dabigatran etexilate (Pradaxa®) and is currently in Phase 3 studies. Here, we report data on the antibody fragment aDabi-Fab2, a putative backup molecule for idarucizumab. Although aDabi-Fab2 completely reversed effects of dabigatran in a rat model in vivo, we observed significantly reduced duration of action compared to idarucizumab. Rational protein engineering, based on the X-ray structure of aDabi-Fab2, led to the identification of mutant Y103W. The mutant had optimized shape complementarity to dabigatran while maintaining an energetically favored hydrogen bond. It displayed increased affinity for dabigatran, mainly driven by a slower off-rate. Interestingly, the increased residence time translated into longer duration of action in vivo. It was thus possible to further enhance the efficacy of aDabi-Fab2 based on rational design, giving it the potential to serve as a back-up candidate for idarucizumab.

Dynamic conformational change regulates the protein-DNA recognition: an investigation on binding of a Y-family polymerase to its target DNA

Protein-DNA recognition is a central biological process that governs the life of cells. A protein will often undergo a conformational transition to form the functional complex with its target DNA. The protein conformational dynamics are expected to contribute to the stability and specificity of DNA recognition and therefore may control the functional activity of the protein-DNA complex. Understanding how the conformational dynamics influences the protein-DNA recognition is still challenging. Here, we developed a two-basin structure-based model to explore functional dynamics in Sulfolobus solfataricus DNA Y-family polymerase IV (DPO4) during its binding to DNA. With explicit consideration of non-specific and specific interactions between DPO4 and DNA, we found that DPO4-DNA recognition is comprised of first 3D diffusion, then a short-range adjustment sliding on DNA and finally specific binding. Interestingly, we found that DPO4 is under a conformational equilibrium between multiple states during the binding process and the distributions of the conformations vary at different binding stages. By modulating the strength of the electrostatic interactions, the flexibility of the linker, and the conformational dynamics in DPO4, we drew a clear picture on how DPO4 dynamically regulates the DNA recognition. We argue that the unique features of flexibility and conformational dynamics in DPO4-DNA recognition have direct implications for low-fidelity translesion DNA synthesis, most of which is found to be accomplished by the Y-family DNA polymerases. Our results help complete the description of the DNA synthesis process for the Y-family polymerases. Furthermore, the methods developed here can be widely applied for future investigations on how various proteins recognize and bind specific DNA substrates.

Protein-DNA recognition is crucial for many key biological processes in cells. Protein often undergoes large-scale conformational change during DNA recognition. However, the physical and global understanding of flexible protein-DNA binding is still challenging. Here, we developed a theoretical approach to investigate binding of a Y-family DNA polymerase to its target DNA during the DNA synthesis process. The results of electrostatic-controlled multi-step DNA binding process accompanied with multi-state conformational transition of protein occurring throughout are in remarkable agreement with experiments. During the process of protein-DNA recognition, the flexibility is found to facilitate both the conformational transition of protein (intra-chain dynamics) and DNA binding (inter-chain dynamics) simultaneously. Therefore, we provided a quantitative description of protein-DNA binding mechanism that flexibility or conformational change regulates DNA recognition dynamically, leading to high efficiency and specificity of function for protein-DNA recognition.

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.

Distal structural elements coordinate a conserved base flipping network

One of the most dramatic illustrations of enzymatic promotion of a high-energy intermediate is observed in DNA modification and repair enzymes where an individual base is rotated (flipped) 180° around the deoxyribose-phosphate backbone and into the active site. While the end states have been extensively characterized, experimental techniques have yet to yield a full description of the base flipping process and the role played by the enzyme. The C5 cytosine methyltransferase M.HhaI coordinates an ensemble of reciprocal DNA and enzyme rearrangements to efficiently flip the target cytosine from the DNA helix. We sought to understand the role of individual amino acids during base flipping. Our results demonstrate that M.HhaI initiates base flipping before closure of the catalytic loop and utilizes the conserved serine 85 in the catalytic loop to accelerate flipping and maintain distortion of the DNA backbone. Serine 87, which forms specific contacts within the DNA helix after base flipping, is not involved in the flipping process or in maintaining the catalytically competent complex. At the base of the catalytic loop, glycine 98 acts as a hinge to allow conformational dynamism of the loop and mutation to alanine inhibits stabilization of the closed loop. Our results illustrate how an enzyme utilizes numerous, distal residues in concert to transform substrate recognition into catalysis.

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.

Direct observation of cytosine flipping and covalent catalysis in a DNA methyltransferase

Methylation of the five position of cytosine in DNA plays important roles in epigenetic regulation in diverse organisms including humans. The transfer of methyl groups from the cofactor S-adenosyl-l-methionine is carried out by methyltransferase enzymes. Using the paradigm bacterial methyltransferase M.HhaI we demonstrate, in a chemically unperturbed system, the first direct real-time analysis of the key mechanistic events—the flipping of the target cytosine base and its covalent activation; these changes were followed by monitoring the hyperchromicity in the DNA and the loss of the cytosine chromophore in the target nucleotide, respectively. Combined with studies of M.HhaI variants containing redesigned tryptophan fluorophores, we find that the target base flipping and the closure of the mobile catalytic loop occur simultaneously, and the rate of this concerted motion inversely correlates with the stability of the target base pair. Subsequently, the covalent activation of the target cytosine is closely followed by but is not coincident with the methyl group transfer from the bound cofactor. These findings provide new insights into the temporal mechanism of this physiologically important reaction and pave the way to in-depth studies of other base-flipping systems.

Molecular and enzymatic profiles of mammalian DNA methyltransferases: structures and targets for drugs

DNA methylation is an epigenetic event involved in a variety array of processes that may be the foundation of genetic phenomena and diseases. DNA methyltransferase is a key enzyme for cytosine methylation in DNA, and can be divided into two functional families (Dnmt1 and Dnmt3) in mammals. All mammalian DNA methyltransferases are encoded by their own single gene, and consisted of catalytic and regulatory regions (except Dnmt2). Via interactions between functional domains in the regulatory or catalytic regions and other adaptors or cofactors, DNA methyltransferases can be localized at selective areas (specific DNA/nucleotide sequence) and linked to specific chromosome status (euchromatin/heterochromatin, various histone modification status). With assistance from UHRF1 and Dnmt3L or other factors in Dnmt1 and Dnmt3a/Dnmt3b, mammalian DNA methyltransferases can be recruited, and then specifically bind to hemimethylated and unmethylated double-stranded DNA sequence to maintain and de novo setup patterns for DNA methylation. Complicated enzymatic steps catalyzed by DNA methyltransferases include methyl group transferred from cofactor Ado-Met to C5 position of the flipped-out cytosine in targeted DNA duplex. In the light of the fact that different DNA methyltransferases are divergent in both structures and functions, and use unique reprogrammed or distorted routines in development of diseases, design of new drugs targeting specific mammalian DNA methyltransferases or their adaptors in the control of key steps in either maintenance or de novo DNA methylation processes will contribute to individually treating diseases related to DNA methyltransferases.

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.

Folding kinetics of recognition loop peptides from a photolyase and cryptochrome-DASH

Cryptochromes (CRY) and photolyases (PL) use a common flavin adenine dinucleotide cofactor and homologous protein scaffold to accomplish numerous, seemingly dissimilar functions. PL repairs UV-damaged DNA in a mechanism requiring light and DNA base flipping. CRY cannot repair DNA, and instead function in core biological processes including plant photomorphogenesis, circadian rhythm, and magnetoreception. One subclass, CRY-DASH, does catalyze repair of single-stranded DNA; compromised base flipping may deactivate its tight binding to duplex DNA substrates. We recently demonstrated that the a "recognition loop" involved in DNA binding by both PL and CRY-DASH is among the most flexible regions in the two proteins, and exhibits especially heightened dynamics in CRY-DASH. Here, we establish that these distinct dynamics are encoded by the loop sequences: we quantify the flexibility of the isolated loop peptides through the kinetics and activation parameters for their folding. Mirroring the dynamics within the proteins, the CRY-DASH recognition loop peptide folds 2.5-fold faster than its counterpart in PL, predominantly due to a lower enthalpy of activation. We propose that these distinct dynamics are functionally significant in DNA recognition. Binding duplex DNA in the catalytically-active base-flipped conformation imposes significant order on the recognition loop, and a corresponding entropic penalty. This may be surmounted by the more preorganized PL recognition loop, but may impose too large a barrier for the more dynamic loop in CRY-DASH. These results suggest that evolution of protein dynamics, through local sequence tuning in the recognition loop, may be an important mechanism for functional diversification in PL and CRY.

Escherichia coli DNA adenine methyltransferase: intrasite processivity and substrate-induced dimerization and activation

Methylation of GATC sites in Escherichia coli by DNA adenine methyltransferase (EcoDam) is essential for proper DNA replication timing, gene regulation, and mismatch repair. The low cellular concentration of EcoDam and the high number of GATC sites in the genome (approximately 20000) support the reliance on methylation efficiency-enhancing strategies such as extensive intersite processivity. Here, we present evidence that EcoDam has evolved other unique mechanisms of activation not commonly observed with restriction-modification methyltransferases. EcoDam dimerizes on short, synthetic DNA, resulting in enhanced catalysis; however, dimerization is not observed on large genomic DNA where the potential for intersite processive methylation precludes any dimerization-dependent activation. An activated form of the enzyme is apparent on large genomic DNA and can also be achieved with high concentrations of short, synthetic substrates. We suggest that this activation is inherent on polymeric DNA where either multiple GATC sites are available for methylation or the partitioning of the enzyme onto nonspecific DNA is favored. Unlike other restriction-modification methyltransferases, EcoDam carries out intrasite processive catalysis whereby the enzyme-DNA complex methylates both strands of an unmethylated GATC site prior to dissociation from the DNA. This occurs with short 21 bp oligonucleotides and is highly dependent upon salt concentrations. Kinetic modeling which invokes enzyme activation by both dimerization and excess substrate provides mechanistic insights into key regulatory checkpoints for an enzyme involved in multiple, diverse biological pathways.

Coupling sequence-specific recognition to DNA modification

Enzymes that modify DNA are faced with significant challenges in specificity for both substrate binding and catalysis. We describe how single hydrogen bonds between M.HhaI, a DNA cytosine methyltransferase, and its DNA substrate regulate the positioning of a peptide loop which is approximately 28 A away. Stopped-flow fluorescence measurements of a tryptophan inserted into the loop provide real-time observations of conformational rearrangements. These long-range interactions that correlate with substrate binding and critically, enzyme turnover, will have broad application to enzyme specificity and drug design for this medically relevant class of enzymes.

Dynamic energy landscape view of coupled binding and protein conformational change: induced-fit versus population-shift mechanisms

Allostery, the coupling between ligand binding and protein conformational change, is the heart of biological network and it has often been explained by two representative models, the induced-fit and the population-shift models. Here, we clarified for what systems one model fits better than the other by performing molecular simulations of coupled binding and conformational change. Based on the dynamic energy landscape view, we developed an implicit ligand-binding model combined with the double-basin Hamiltonian that describes conformational change. From model simulations performed for a broad range of parameters, we uncovered that each of the two models has its own range of applicability, stronger and longer-ranged interaction between ligand and protein favors the induced-fit model, and weaker and shorter-ranged interaction leads to the population-shift model. We further postulate that the protein binding to small ligand tends to proceed via the population-shift model, whereas the protein docking to macromolecules such as DNA tends to fit the induced-fit model.

Backbone dynamics in the DNA HhaI protein binding site

The dynamics of the phosphodiester backbone in the [5'-GCGC-3'] 2 moiety of the DNA oligomer [d(G 1A 2T 3A 4 G 5 C 6 G 7 C 8T 9A 10T 11C 12)] 2 are studied using deuterium solid-state NMR (SSNMR). SSNMR spectra obtained from DNAs nonstereospecifically deuterated on the 5' methylene group of nucleotides within the [5'-GCGC-3'] 2 moiety indicated that all of these positions are structurally flexible. Previous work has shown that methylation reduces the amplitude of motion in the phosphodiester backbone and furanose ring of the same DNA, and our observations indicate that methylation perturbs backbone dynamics through not only a loss of mobility but also a change of direction of motion. These NMR data indicate that the [5'-GCGC-3'] 2 moiety is dynamic, with the largest amplitude motions occurring nearest the methylation site. The change of orientation of this moiety in DNA upon methylation may make the molecule less amenable to binding to the HhaI endonuclease.