A quantum mechanics/molecular mechanics study of the catalytic mechanism and product specificity of viral histone lysine methyltransferase

There are three reaction steps in the S-adenosylmethionine (AdoMet) methylation of lysine-NH2 catalyzed by a methyltransferase. They are (i) combination of enzyme.Lys-NH3+ with AdoMet, (ii) substrate ionization to provide enzyme.AdoMet.Lys-NH2, and (iii) methyl transfer providing enzyme.AdoHcy.Lys-N(Me)H2+ and the dissociation of AdoHcy. In this study of the viral histone methyltransferase (vSET), we find that substrate ionization of vSET.Lys27-NH3+, vSET.Lys27-N(Me)H2+, and vSET.Lys27-N(Me)2H+ takes place upon combination with AdoMet. The presence of a water channel allows dissociation of a proton to the solvent. There is no water channel in the absence of AdoMet. That the formation of a water channel is combined with AdoMet binding was first discovered in our investigation of Rubisco large subunit methyltransferase. Via a quantum mechanics/molecular mechanics (QM/MM) approach, the calculated free energy barrier (DeltaG++) of the first methyl transfer reaction catalyzed by vSET [Lys27-NH2 + AdoMet --> Lys27-N(Me)H2+ + AdoHcy] equals 22.5 +/- 4.3 kcal/mol, which is in excellent agreement with the free energy barrier (21.7 kcal/mol) calculated from the experimental rate constant (0.047 min-1). The calculated DeltaG++ of the second methyl transfer reaction [AdoMet + Lys27-N(Me)H --> AdoHcy + Lys27-N(Me)2H+] at the QM/MM level is 22.6 +/- 3.6 kcal/mol, which is in agreement with the value of 22.4 kcal/mol determined from the experimental rate constant (0.015 min-1). The third methylation [Lys27-N(Me)2 + AdoMet --> Lys27-N(Me)3+ + AdoHcy] is associated with a DeltaG++ of 23.1 +/- 4.0 kcal/mol, which is in agreement with the value of 23.0 kcal/mol determined from the experimental rate constant (0.005 min-1). Our computations establish that the first, second, and third methyl transfer steps catalyzed by vSET are linear SN2 reactions with the bond making being approximately 50% associative.     

Product specificity and mechanism of protein lysine methyltransferases: insights from the histone lysine methyltransferase SET8

Molecular dynamics simulations employing a molecular mechanics (MM) force field and hybrid quantum mechanics (QM) and MM (QM/MM) have been carried out to investigate the product specificity and mechanism of the histone H4 lysine 20 (H4-K20) methylation by human histone lysine methyltransferase SET8. At neutral pH, the target lysine is available to only the enzyme in the protonated state. The first step in the methylation reaction must be deprotonation of the lysine target which is followed by the (+)AdoMet methylation of the neutral lysine [Enz.Lys-CH(2)-NH(3)(+).(+)AdoMet --> H(+) + Enz.Lys-CH(2)-NH(2).(+)AdoMet -->--> Enz.Lys-CH(2)-N(Me)H(2)(+).AdoHcy]. The electrostatic interactions between two positive charges on (+)AdoMet and Lys20-NH(3)(+) decrease the pK(a) of Lys20-NH(3)(+). Upon formation of Enz.Lys-NH(3)(+).(+)AdoMet, a water channel by which the proton escapes to the outer solvent phase is formed. The formation of a water channel for the escape of a proton from Lys20-N(Me)H(2)(+) in Enz.Lys20-N(Me)H(2)(+).(+)AdoMet is not formed because the methyl substituent blocks the starting of the water channel. Thus, a second methylation does not take place. The dependence of the occurrence of methyl transfer on the formation of a water channel in SET8 is in accord with our previous reports on product specificity by histone lysine monomethyltransferase SET7/9, large subunit lysine dimethyltransferase (LSMT), and viral histone lysine trimethyltransferase (vSET). The average value of the experimental DeltaG(E)() for the six lysine methyl transfer reactions catalyzed by vSET, LSMT, and SET7/9 with p53 as a substrate is 22.1 +/- 1.0 kcal/mol, and the computed average (DeltaG(C)()) is 22.2 +/- 0.8 kcal/mol. In this study, the computed free energy barrier of the methyl transfer reaction [Lys20-NH(2) + (+)AdoMet --> Lys20-N(Me)H(2)(+) + AdoHcy] catalyzed by SET8 is 20.8 kcal/mol. This is in agreement with the value of 20.6 kcal/mol calculated from the experimental rate constant (0.43 +/- 0.02 min(-1)). Our bond-order computations establish that the H4-K20 monomethylation in SET8 is a concerted linear S(N)2 displacement reaction.     

Histone lysine methyltransferase SET7/9: formation of a water channel precedes each methyl transfer

Molecular dynamics (MD) simulations and hybrid quantum mechanics/molecular mechanics (QM/MM) calculations have been carried out in an investigation of histone lysine methyltransferase (SET7/9). Proton dissociation (SET7/9.Lys4-NH3+.AdoMet --> SET7/9.Lys4-NH2.AdoMet + H+) must be prior to the methylation by S-adenosylmethionine (AdoMet). We find that a water channel is formed to allow escape of the proton to solvent. The water channel appears in the presence of AdoMet, but is not present in the species SET7/9.Lys4-NH3+ or SET7/9.Lys4-N(Me)H2+.AdoHcy. A water channel is not formed in the ground state of SET7/9.Lys4-N(Me)H2+.AdoMet, and the second methyl transfer does not occur. The structure of SET7/9.Lys4-N(Me)H2+.AdoMet includes a greater distance (6.1 +/- 0.3 A) between Cgamma(AdoMet) and N(MeLys4) than is present in SET7/9.Lys4-NH3+.AdoMet (5.7 +/- 0.2 A). The electrostatic interactions between the positive charges on AdoMet and SET7/9.Lys4-NH3+ decrease the pKa of the latter from 10.9 +/- 0.4 to 8.2 +/- 0.6, and this is not seen in the SET7/9.Lys4-N(Me)H2+.AdoMet species. The formation, or not, of a water channel, the distance between Sdelta(AdoMet) and N(Lys4), and the angle Sdelta(AdoMet)-Cgamma(AdoMet)-N(Lys4) determine whether methyl transfer can occur. By QM/MM, the calculated free energy barrier of the methyl transfer reaction in the SET7/9 [Lys4-NH2 + AdoMet --> Lys4-N(Me)H2+ + AdoHcy] complex is DeltaG++ = 19.0 +/- 1.6 kcal/mol. This DeltaG++ is in agreement with the value of 20.9 kcal/mol calculated from the experimental rate constant (0.24 min(-1)).     

Mechanism of multiple lysine methylation by the SET domain enzyme Rubisco LSMT

SET domain protein methyltransferases catalyze the transfer of methyl groups from the cofactor S-adenosylmethionine (AdoMet) to specific lysine residues of protein substrates, such as the N-terminal tails of histones H3 and H4 and the large subunit of the Rubisco holoenzyme complex. The crystal structures of pea Rubisco large subunit methyltransferase (LSMT) in ternary complexes with either lysine or epsilon-N-methyllysine (MeLys) and the product S-adenosylhomocysteine (AdoHcy) were determined to resolutions of 2.65 and 2.55 A, respectively. The zeta-methyl group of MeLys is bound to the enzyme via carbon-oxygen hydrogen bonds that play a key role in catalysis. The methyl donor and acceptor are aligned in a linear geometry for S(N)2 nucleophilic transfer of the methyl group during catalysis. Differences in hydrogen bonding between the MeLys epsilon-amino group and Rubisco LSMT and SET7/9 explain why Rubisco LSMT generates multiply methylated Lys, wheras SET7/9 generates only MeLys.     

Carbinolamine formation and dehydration in a DNA repair enzyme active site

In order to suggest detailed mechanistic hypotheses for the formation and dehydration of a key carbinolamine intermediate in the T4 pyrimidine dimer glycosylase (T4PDG) reaction, we have investigated these reactions using steered molecular dynamics with a coupled quantum mechanics-molecular mechanics potential (QM/MM). We carried out simulations of DNA abasic site carbinolamine formation with and without a water molecule restrained to remain within the active site quantum region. We recovered potentials of mean force (PMF) from thirty replicate reaction trajectories using Jarzynski averaging. We demonstrated feasible pathways involving water, as well as those independent of water participation. The water-independent enzyme-catalyzed reaction had a bias-corrected Jarzynski-average barrier height of approximately (6.5 kcal mol(-1) (27.2 kJ mol(-1)) for the carbinolamine formation reaction and 44.5 kcal mol(-1) (186 kJ mol(-1)) for the reverse reaction at this level of representation. When the proton transfer was facilitated with an intrinsic quantum water, the barrier height was approximately 15 kcal mol(-1) (62.8 kJ mol(-1)) in the forward (formation) reaction and 19 kcal mol(-1) (79.5 kJ mol(-1)) for the reverse. In addition, two modes of unsteered (free dynamics) carbinolamine dehydration were observed: in one, the quantum water participated as an intermediate proton transfer species, and in the other, the active site protonated glutamate hydrogen was directly transferred to the carbinolamine oxygen. Water-independent unforced proton transfer from the protonated active site glutamate carboxyl to the unprotonated N-terminal amine was also observed. In summary, complex proton transfer events, some involving water intermediates, were studied in QM/MM simulations of T4PDG bound to a DNA abasic site. Imine carbinolamine formation was characterized using steered QM/MM molecular dynamics. Dehydration of the carbinolamine intermediate to form the final imine product was observed in free, unsteered, QM/MM dynamics simulations, as was unforced acid-base transfer between the active site carboxylate and the N-terminal amine.     

Bacteriophage T4 Dam DNA-[N6-adenine]methyltransferase. Kinetic evidence for a catalytically essential conformational change in the ternary complex.

We carried out a steady state kinetic analysis of the bacteriophage T4 DNA-[N6-adenine]methyltransferase (T4 Dam) mediated methyl group transfer reaction from S-adenosyl-l-methionine (AdoMet) to Ade in the palindromic recognition sequence, GATC, of a 20-mer oligonucleotide duplex. Product inhibition patterns were consistent with a steady state-ordered bi-bi mechanism in which the order of substrate binding and product (methylated DNA, DNA(Me) and S-adenosyl-l-homocysteine, AdoHcy) release was AdoMet downward arrow DNA downward arrow DNA(Me) upward arrow AdoHcy upward arrow. A strong reduction in the rate of methylation was observed at high concentrations of the substrate 20-mer DNA duplex. In contrast, increasing substrate AdoMet concentration led to stimulation in the reaction rate with no evidence of saturation. We propose the following model. Free T4 Dam (initially in conformational form E) randomly interacts with substrates AdoMet and DNA to form a ternary T4 Dam-AdoMet-DNA complex in which T4 Dam has isomerized to conformational state F, which is specifically adapted for catalysis. After the chemical step of methyl group transfer from AdoMet to DNA, product DNA(Me) dissociates relatively rapidly (k(off) = 1.7 x s(-1)) from the complex. In contrast, dissociation of product AdoHcy proceeds relatively slowly (k(off) = 0.018 x s(-1)), indicating that its release is the rate-limiting step, consistent with kcat = 0.015 x s(-1). After AdoHcy release, the enzyme remains in the F conformational form and is able to preferentially bind AdoMet (unlike form E, which randomly binds AdoMet and DNA), and the AdoMet-F binary complex then binds DNA to start another methylation cycle. We also propose an alternative pathway in which the release of AdoHcy is coordinated with the binding of AdoMet in a single concerted event, while T4 Dam remains in the isomerized form F. The resulting AdoMet-F binary complex then binds DNA, and another methylation reaction ensues. This route is preferred at high AdoMet concentrations.     

Ionic selectivity, saturation, and block in a K+-selective channel from sarcoplasmic reticulum

The open-channel conductance properties of a voltage-gated channel from sarcoplasmic reticulum were studied in planar phospholipid membranes. The channel is ideally selective for K+ over Cl- and for K+ over Ca++. In symmetrical 1 M solutions, the single-channel conductance (in pmho) falls in the order: K+ (214) > NH4+ (157) > Rb+ (125) > Na+ (72) > La+ (8.1) > Cs+ (< 3). In neutral bilayers, the channel conductance saturates with ion activity according to a rectangular hyperbolic relation, with half-saturation activities of 54 mM for K+ and 34 mM for Na+. Under symmetrical salt conditions, the K+:Na+ channel conductance ratio increases with salt activity, but the permeability ratio, measured by single-channel bi-ionic potentials, is constant between 20 mM and 2.5 M salt; the permeability ratio is equal to the conductance ratio in the limit of low-salt concentration. The channel conductance varies < 5% in the voltage range -100 to +70 mV. The maximum conductance varies K+ and Na+ is only weakly temperature dependent (delta H++ = 4.6 and 5.3 kcal/mol, respectively), but that of Li+ varies strongly with temperature (delta H++ = 13 kcal/mol). The channel's K+ conductance is blocked asymmetrically by Cs+, and this block is competitive with K+. The results are consistent with an Eyring-type barriers as it permeates the channel. The data conform to Lüger's (1973. Biochem. Biophys. Acta. 311:423-441) predictions for a "pure" single-ion channel.     

Bacteriophage T4 Dam DNA-(N6-adenine)-methyltransferase. Processivity and orientation to the methylation target

We carried out steady state and pre-steady state (burst) kinetic analyses of the bacteriophage T4 Dam DNA-(N(6)-adenine)-methyltransferase (MTase)-mediated methyl group transfer from S-adenosyl-l-methionine (AdoMet) to Ade in oligonucleotide duplexes containing one or two specific GATC sites with different combinations of methylated and unmodified targets. We compared the results for ligated 40-mer duplexes with those of the mixtures of the two unligated duplexes used to generate the 40-mers. The salient results are as follows: (i) T4 Dam MTase modifies 40-mer duplexes in a processive fashion. (ii) During processive movement, T4 Dam rapidly exchanges product S-adenosyl-l-homocysteine (AdoHcy) for substrate AdoMet without dissociating from the DNA duplex. (iii) T4 Dam processivity is consistent with an ordered bi-bi mechanism AdoMet downward arrow DNA downward arrow DNA(Me) upward arrow AdoHcy upward arrow. However, in contrast to the steady state, here DNA(Me) upward arrow signifies departure from a methylated site GMTC upward arrow without physically dissociating from the DNA. (iv) Following methyl transfer at one site and linear diffusion to a hemimethylated site, a reconstituted T4 Dam-AdoMet complex rapidly reorients itself to the (productive) unmethylated strand. T4 Dam-AdoHcy cannot reorient at an enzymatically created GMTC site. (v) The inhibition potential of fully methylated sites 5'-GMTC/5'-GMTC is much lower for a long DNA molecule compared with short single-site duplexes.     

Molecular dynamics simulations of the intramolecular proton transfer and carbanion stabilization in the pyridoxal 5'-phosphate dependent enzymes L-dopa decarboxylase and alanine racemase

Molecular dynamics simulations using a combined quantum mechanical and molecular mechanical (QM/MM) potential have been carried out to investigate the internal proton transfer equilibrium of the external aldimine species in l-dopa decarboxylase, and carbanion stabilization by the enzyme cofactor in the active site of alanine racemase. Solvent effects lower the free energy of the O-protonated PLP tautomer both in aqueous solution and in the active site, resulting a free energy difference of about -1 kcal/mol relative to the N-protonated Schiff base in the enzyme. The external aldimine provides the dominant contribution to lowering the free energy barrier for the spontaneous decarboxylation of l-dopa in water, by a remarkable 16 kcal/mol, while the enzyme l-dopa decarboxylase further lowers the barrier by 8 kcal/mol. Kinetic isotope effects were also determined using a path integral free energy perturbation theory on the primary (13)C and the secondary (2)H substitutions. In the case of alanine racemase, if the pyridine ring is unprotonated as that in the active site, there is destabilizing contribution to the formation of the α-carbanion in the gas phase, although when the pyridine ring is protonated the contribution is stabilizing. In aqueous solution and in alanine racemase, the α-carbanion is stabilized both when the pyridine ring is protonated and unprotonated. The computational studies illustrated in this article show that combined QM/MM simulations can help provide a deeper understanding of the mechanisms of PLP-dependent enzymes. This article is part of a Special Issue entitled: Pyridoxal Phosphate Enzymology.     

H nuclear magnetic resonance study of restricted internal rotation of N-6,N-6-dimethyladenine in aqueous solution.

Kinetics of internal rotation about the C(6)-N(6) bond of N-6,N-6-dimethyladenine (M2-6A) was investigated by -1H nuclear magnetic resonance line-shape analysis of the methyl resonances (220 MHz). Rates of rotation were determined for M2-6A deuterated at N(1) and for neutral M2-6A. Activation parameters for monodeuterated M2-6A at 22 degrees are Ea = 13.8kcal/mol, log A = 12.6, incrementG++=14.9 kcal/mol, incrementH++ = 13.1 kcal/mol, incrementS++ = minus 5.8 eu; for neutral M2-6A: Ea = 15.5 kcal/mol, log A = 14.9, incrementG++ = 12.6 kcal/mol, incrementH++ = 14.9 kcal/mol, incrementS++ =7.8 eu. Vertical stacking of bases interferes with internal rotation of the dimethylamino group.     

The enthalpy change for the reduction of nicotinamide–adenine dinucleotide

The heat of the reaction NAD(+)+propan-2-ol=NADH+acetone+H(+) was determined to be 42.5+/-0.6kJ/mol (10.17+/-0.15kcal/mol) from equilibrium measurements at 9-42 degrees C catalysed by yeast alcohol dehydrogenase. With the aid of thermochemical data for acetone and propan-2-ol the values of DeltaH=-29.2kJ/mol (-6.99kcal/mol) and DeltaG(0)=22.1kJ/mol (5.28kcal/mol) are derived for the reduction of NAD (NAD(+)+H(2)=NADH+H(+)). These values are consistent with analogous but less accurate data for the ethanol-acetaldehyde reaction. Thermodynamic data for the reduction of NAD and NADP are summarized.     

DNA (cytosine-N4-)- and -(adenine-N6-)-methyltransferases have different kinetic mechanisms but the same reaction route. A comparison of M.BamHI and T4 Dam

We studied the kinetics of methyl group transfer by the BamHI DNA-(cytosine-N(4)-)-methyltransferase (MTase) from Bacillus amyloliquefaciens to a 20-mer oligodeoxynucleotide duplex containing the palindromic recognition site GGATCC. Under steady state conditions the BamHI MTase displayed a simple kinetic behavior toward the 20-mer duplex. There was no apparent substrate inhibition at concentrations much higher than the K(m) for either DNA (100-fold higher) or S-adenosyl-l-methionine (AdoMet) (20-fold higher); this indicates that dead-end complexes did not form in the course of the methylation reaction. The DNA methylation rate was analyzed as a function of both substrate and product concentrations. It was found to exhibit product inhibition patterns consistent with a steady state random bi-bi mechanism in which the dominant order of substrate binding and product release (methylated DNA, DNA(Me), and S-adenosyl-l-homocysteine, AdoHcy) was Ado-Met DNA DNA(Me) AdoHcy. The M.BamHI kinetic scheme was compared with that for the T4 Dam (adenine-N(6)-)-MTase. The two differed with respect to an effector action of substrates and in the rate-limiting step of the reaction (product inhibition patterns are the same for the both MTases). From this we conclude that the common chemical step in the methylation reaction, methyl transfer from AdoMet to a free exocyclic amino group, is not sufficient to dictate a common kinetic scheme even though both MTases follow the same reaction route.     

Structure and mechanism of mRNA cap (guanine-N7) methyltransferase

A suite of crystal structures is reported for a cellular mRNA cap (guanine-N7) methyltransferase in complex with AdoMet, AdoHcy, and the cap guanylate. Superposition of ligand complexes suggests an in-line mechanism of methyl transfer, albeit without direct contacts between the enzyme and either the N7 atom of guanine (the attacking nucleophile), the methyl carbon of AdoMet, or the sulfur of AdoMet/AdoHcy (the leaving group). The structures indicate that catalysis of cap N7 methylation is accomplished by optimizing proximity and orientation of the substrates, assisted by a favorable electrostatic environment. The enzyme-ligand structures, together with new mutational data, fully account for the biochemical specificity of the cap guanine-N7 methylation reaction, an essential and defining step of eukaryotic mRNA synthesis.     

QM/MM study of the active site of free papain and of the NMA-papain complex.

Hybrid quantum mechanical/molecular mechanical (QM/MM) calculations using restricted and unrestricted Hartree-Fock and B3LYP ab initio (QM) and Amber force field (MM), respectively, have been applied to study the catalytic site of papain in both free and substrate bonded forms. Ab initio geometry optimizations have been performed for the active site of papain and the N-methyl-acetamide (NMA)-papain complex within the molecular mechanical treatment of the protein environment. A covalent tetrahedral intermediate structure could be obtained only when the amide N atom of the substrate molecule was protonated through a proton transfer from the His-159 in the catalytic site. Our results support the previous assumption that a proton transfer from His-159 to the amide N atom of the substrate occurs prior to or concerted with the nucleophilic attack of the Cys-25 sulfur atom to the carbonyl group of the substrate. The electron correlation effect will reduce the proton transfer barrier. Therefore, this proton transfer can be easily observed in the B3LYP/6-31G* calculations. The HF/6-31G* method overestimates the reaction barrier against this proton transfer. The sulfur atom of Cys-25 and the imidazole ring of His-159 are found to be coplanar in the free form of the enzyme. However, the rotation of the imidazole ring of His-159 was observed during the formation of the tetrahedral intermediate. Without the papain environment, the coplanar thiolate-imidazolium ion pair RS-...ImH+ is much less stable than the neutral form of RSH....Im. Within the protein environment, however, the thiolate-imidazolium ion pair becomes more stable than its neutral form by 4.1 and 0.4 kcal/mol in HF/6-31G* and B3LYP/6-31G* calculations, respectively. The barrier of proton transfer from S-H group of Cys-25 to the imidazole ring of His-159 was reduced from 22.0 kcal/mol to 15.2 kcal/mol by the protein environment in HF/6-31G* calculations. This barrier is found to be much smaller (2.5 kcal/mol) in B3LYP/6-31G* calculations.     

Hydration of carbon dioxide by carbonic anhydrase: internal proton transfer of Zn2+-bound HCO3-

Proton transfer within HCO3- has been examined under various conditions through molecular orbital methods: partial retention of diatomic differential overlap and 4-31G self-consistent field programs. These conditions include the absence or presence of Zn2+, Zn2+(NH3)3, or a water ligand on Zn2+. In addition, 4-31G+ and some MP2/4-31G results are obtained. The use of Be2+ to simulate Zn2+ reproduces reaction pathways and energy barriers, except for marginal cases. The barrier of 35.6 kcal/mol for direct internal proton transfer is reduced to 3.5 kcal/mol when one water molecule, not bound to Zn2+, is included for proton relay and to 1.4 kcal/mol when two such water molecules are included. In the enzyme, either Thr-199 or solvent molecules could perform this relay function. Our results favor this facilitated proton transfer over a mechanism in which Zn2+-bound OH- attacks CO2, a bidentate intermediate forms, and the OH moiety of the resulting HCO3- dissociates from Zn2+, thus leaving one of the oxygens of the original CO2 as a ligand to Zn2+.     

Exploring the RING-Catalyzed Ubiquitin Transfer Mechanism by MD and QM/MM Calculations

Ubiquitylation is a universal mechanism for controlling cellular functions. A large family of ubiquitin E3 ligases (E3) mediates Ubiquitin (Ub) modification. To facilitate Ub transfer, RING E3 ligases bind both the substrate and ubiquitin E2 conjugating enzyme (E2) linked to Ub via a thioester bond to form a catalytic complex. The mechanism of Ub transfer catalyzed by RING E3 remains elusive. By employing a combined computational approach including molecular modeling, molecular dynamics (MD) simulations, and quantum mechanics/molecular mechanics (QM/MM) calculations, we characterized this catalytic mechanism in detail. The three-dimensional model of dimeric RING E3 ligase RNF4 RING, E2 ligase UbcH5A, Ub and the substrate SUMO2 shows close contact between the substrate and Ub transfer catalytic center. Deprotonation of the substrate lysine by D117 on UbcH5A occurs with almost no energy barrier as calculated by MD and QM/MM calculations. Then, the side chain of the activated lysine gets close to the thioester bond via a conformation change. The Ub transfer pathway begins with a nucleophilic addition that forms an oxyanion intermediate of a 4.23 kcal/mol energy barrier followed by nucleophilic elimination, resulting in a Ub modified substrate by a 5.65 kcal/mol energy barrier. These results provide insight into the mechanism of RING-catalyzed Ub transfer guiding the discovery of Ub system inhibitors.     

Synthesis of S-adenosyl-L-methionine analogs and their use for sequence-specific transalkylation of DNA by methyltransferases

Here we describe a one-step synthetic procedure for the preparation of S-adenosyl-L-methionine (AdoMet) analogs with extended carbon chains replacing the methyl group. These AdoMet analogs function as efficient cofactors for DNA methyltransferases (MTases), and we provide a protocol for sequence-specific transfer of extended side chains from these AdoMet analogs to DNA by DNA MTases. Direct chemoselective allylation or propargylation of S-adenosyl-L-homocysteine (AdoHcy) at sulfur is achieved under the acidic conditions needed to protect other nucleophilic positions in AdoHcy. The unsaturated bonds in beta position to the sulfonium center of the resulting AdoMet analogs are designed to stabilize the transition state formed upon DNA MTase-catalyzed nucleophilic attack at the carbon next to the sulfonium center and lead to efficient transfer of the extended side chains to DNA. Using these protocols, sequence-specific functionalized DNA can be obtained within one to two weeks.     

S-adenosyl-L-methionine-dependent methyl transfer: observable precatalytic intermediates during DNA cytosine methylation

S-adenosyl-L-methionine- (AdoMet-) dependent methyltransferases are widespread, play critical roles in diverse biological pathways, and are antibiotic and cancer drug targets. Presently missing from our understanding of any AdoMet-dependent methyl-transfer reaction is a high-resolution structure of a precatalytic enzyme/AdoMet/DNA complex. The catalytic mechanism of DNA cytosine methylation was studied by structurally and functionally characterizing several active site mutants of the bacterial enzyme M.HhaI. The 2.64 A resolution protein/DNA/AdoMet structure of the inactive C81A M.HhaI mutant suggests that active site water, an approximately 13 degree tilt of the target base toward the active site nucleophile, and the presence or absence of the cofactor methylsulfonium are coupled via a hydrogen-bonding network involving Tyr167. The active site in the mutant complex is assembled to optimally align the pyrimidine for nucleophilic attack and subsequent methyl transfer, consistent with previous molecular dynamics ab initio and quantum mechanics/molecular mechanics calculations. The mutant/DNA/AdoHcy structure (2.88 A resolution) provides a direct comparison to the postcatalytic complex. A third C81A ternary structure (2.22 A resolution) reveals hydrolysis of AdoMet to adenosine in the active site, further validating the coupling between the methionine portion of AdoMet and ultimately validating the structural observation of a prechemistry/postchemistry water network. Disruption of this hydrogen-bonding network by a Tyr167 to Phe167 mutation does not alter the kinetics of nucleophilic attack or methyl transfer. However, the Y167F mutant shows detectable changes in kcat, caused by the perturbed kinetics of AdoHcy release. These results provide a basis for including an extensive hydrogen-bonding network in controlling the rate-limiting product release steps during cytosine methylation.