Enzymatic mechanism and product specificity of SET-domain protein lysine methyltransferases.

Molecular dynamics and hybrid quantum mechanics/molecular mechanics have been used to investigate the mechanisms of (+)AdoMet methylation of protein-Lys-NH(2) catalyzed by the lysine methyltransferase enzymes: histone lysine monomethyltransferase SET7/9, Rubisco large-subunit dimethyltransferase, viral histone lysine trimethyltransferase, and the Tyr245Phe mutation of SET7/9. At neutrality in aqueous solution, primary amines are protonated. The enzyme reacts with Lys-NH(3)(+) and (+)AdoMet species to provide an Enz.Lys-NH(3)(+).(+)AdoMet complex. The close positioning of two positive charges lowers the pK(a) of the Lys-NH(3)(+) entity, a water channel appears, and the proton escapes to the aqueous solvent; then the reaction Enz.Lys-NH(2).(+)AdoMet --> Enz.Lys-N(Me)H(2)(+).AdoHcy occurs. Repeat of the sequence provides dimethylated lysine, and another repeat yields a trimethylated lysine. The sequence is halted at monomethylation when the conformation of the Enz.Lys-N(Me)H(2)(+).(+)AdoMet has the methyl positioned to block formation of a water channel. The sequence of reactions stops at dimethylation if the conformation of Enz.Lys-N(Me)(2)H(+).(+)AdoMet has a methyl in position, which forbids the formation of the water channel.

Development of potential anticancer agents that target the telomere sequence.

The immortality of cancer cells is due to the relatively high concentration of telomerase enzyme that maintains the telomere sequence during cell division. Human telomeric DNA consists of repeats of the sequence d(5'-TTAGGG-3'). Deoxyribonucleic guanidine (DNG) is a DNA analog in which positively charged guanidine [-NH-C(NH2+)-NH-] replaces the negatively charged phosphodiester of DNA. The synthesized DNG hexamer AgAgTgCgCpC and dodecamer AgAgTgCgCgCAgAgTgCgCpC are complementary to the non-coding telomere sequence d(5'-TTAGGG-3'). We have found that binding of the complementary DNG hexamer to the telomere is favored over that of DNA telomere by 10(2.5)-fold and binding the dodecamer with 2-mismatched DNA is favored by 10(5)-fold. We have shown that DNG binding to RNA is favored over binding to DNA. A complementary complex of DNG with RNA at the active site of telomerase enzyme would be very stable.

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.

Complexation of single strand telomere and telomerase RNA template polyanions by deoxyribonucleic guanidine (DNG) polycations: plausible anticancer agents.

Cancerous cell immortality is due to relatively high concentrations of telomerase enzyme which maintains telomere sequence during cell division. Deoxyribonucleic guanidine (DNG) is a positively charged DNA analog in which guanidine replaces the phosphordiester linkage of DNA. Mixed sequences of DNG and DNA oligonucleotides are referred to as chimera. Complexation of DNG and chimeric polycations with the complementary negatively charged non-coding telomere single strand d(5'-TTAGGG-3')(n) and the 11-base telomeric RNA template (5'-CUAACCCUAAC-3') in the active site of telomerase has been studied. Calculated by ensemble sampling simulations in GBMV solvent model, we found that binding of complementary DNG hexamer with telomere is favored over that of DNA-telomere by approximately 10(6)-fold and binding of chimera hexamer is favored by approximately 10(4)-fold. Binding of complementary DNG with telomeric RNA is favored by 43 kcal/mol over telomere substrate binding with telomeric RNA.

Incorporation of positively charged ribonucleic guanidine linkages into oligodeoxyribonucleotides: Development of potent antisense agents.

Oligodeoxynucleic acid (21-mer) containing both negatively charged phosphate and positively charged ribonucleic guanidine linkages (RNG/DNA chimera) have been synthesized. DNA binding characteristics and nuclease resistance of RNG/DNA chimeras have been evaluated. Using the bcr-abl oncogene (cause of chronic myeloid leukemia) as a target, the binding of a 21-mer RNG/DNA chimera that includes six RNG's is more than 103.5 stronger than the binding of 21-mer composed solely of DNA.

Binding properties of positively charged deoxynucleic guanidine (DNG), AgTgAgTgAgT and DNG/DNA chimeras to DNA.

The melting properties of hexameric oligonucleotide AgTgAgTgAgT, in which the phosphodiester linkages of the DNA have been replaced by guanidium linkages, have been evaluated. Using the juvenile esterase gene as a target, the binding of a 20-mer DNG/DNA chimera that includes AgTgAgTgAgT is more than 10(5.7) stronger than the binding of 20-mer composed solely of DNA.

Low-frequency normal modes in horse liver alcohol dehydrogenase and motions of residues involved in the enzymatic reaction.

Normal mode analysis using the elastic network model has provided characteristics and directions of the low-frequency large domain motions of horse liver alcohol dehydrogenase. Three normal modes (mode 1, mode 7, and mode 8) were identified as representative domain motions that may promote the onset of Near Attack Conformers or facilitate the product to be released. The pattern of the atomic displacement for some key residues (such as Val292 and Val203) revealed in this study is in line with experimental structural and kinetic studies and theoretical simulations.

An investigation of the ionic and solvation patterns of dsDNG versus dsDNA by use of molecular dynamics simulations.

Molecular dynamics has been employed to analyze the counterion and water atmospheres around the deoxynucleic guanidine (DNG) duplex G(12)-C(12). These features are compared to those of DNA G(12)-C(12). The chloride counterions of cationic DNG demonstrate fewer penetrations and only fleeting residence times in the minor groove, as opposed to the multi-nanosecond visits seen by sodium ions in DNA minor grooves. The 10 ns simulations also show the differences in hydration patterns between the DNG and DNA duplexes.

Inhibiting transcription factor/DNA complexes using fluorescent microgonotropens (FMGTs).

Fluorescent microgonotropens (FMGTs) are A/T selective, minor groove-binding bisbenzimidazole ligands. Basic side chains extending from these agents electrostatically contact the major groove side of the phosphodiester backbone of DNA, endowing them with high binding affinity. Here, we evaluate the potential of these agents as inhibitors of transcription factor (TF) binding to DNA and explore whether their ability to contact both grooves enhances their inhibitory activity. A series of FMGTs (L2-L5), with polyamine tails of varying lengths and degrees of branching, were compared to an analog lacking these basic side chains (L1), and the classical bisbenzimidazole Hoechst 33342 for effects on TF complex formation on the c-fos serum response element (SRE). Although L1 could not inhibit TF/SRE interactions, L2-L5 did so at submicromolar concentrations. Moreover, the FMGTs were up to 50 times more potent than Hoechst 33342 in inhibiting TF complex formation in electrophoretic mobility shift assays. The FMGTs also inhibited c-fos promoter-driven cell-free transcription and topoisomerase II activity in nuclei. These studies establish the potential of FMGTs as TF/DNA complex inhibitors in cell-free systems, provide insight into the relationship between their structure and biological activities, and demonstrate the benefits of functionalizing minor groove binding-agents with major groove-contacting groups.

Sequence specific DNA binding of Ets-1 transcription factor: molecular dynamics study on the Ets domain--DNA complexes.

Molecular dynamics (MD) simulations for Ets-1 ETS domain-DNA complexes were performed to investigate the mechanism of sequence-specific recognition of the GGAA DNA core by the ETS domain. Employing the crystal structure of the Ets-1 ETS domain-DNA complex as a starting structure we carried out MD simulations of: (i). the complex between Ets-1 ETS domain and a 14 base-pair DNA containing GGAA core sequence (ETS-GGAA); (ii). the complex between the ETS domain and a DNA having single base-pair mutation, GGAG sequence (ETS-GGAG); and (iii). the 14 base-pair DNA alone (GGAA). Comparative analyses of the MD structures of ETS-GGAA and ETS-GGAG reveal that the DNA bending angles and the ETS domain-DNA phosphate interactions are similar in these complexes. These results support that the GGAA core sequence is distinguished from the mutated GGAG sequence by a direct readout mechanism in the Ets-1 ETS domain-DNA complex. Further analyses of the direct contacts in the interface between the helix-3 region of Ets-1 and the major groove of the core DNA sequence clearly show that the highly conserved arginine residues, Arg391 and Arg394, play a critical role in binding to the GGAA core sequence. These arginine residues make bidentate contacts with the nucleobases of GG dinucleotides in GGAA core sequence. In ETS-GGAA, the hydroxyl group of Tyr395 is hydrogen bonded to N7 nitrogen of A(3) (the third adenosine in the GGAA core), while the hydroxyl group makes a contact with N4 nitrogen of C(4') (the complementary nucleotide of the fourth guanosine G(4) in the GGAG sequence) in the ETS-GGAG complex. We have found that this difference in behavior of Tyr395 results in the relatively large motion of helix-3 in the ETS-GGAG complex, causing the collapse of bidentate contacts between Arg391/Arg394 and the GG dinucleotides in the GGAG sequence.

Determination of enzyme mechanisms by molecular dynamics: studies on quinoproteins, methanol dehydrogenase, and soluble glucose dehydrogenase.

Molecular dynamics (MD) simulations have been carried out to study the enzymatic mechanisms of quinoproteins, methanol dehydrogenase (MDH), and soluble glucose dehydrogenase (sGDH). The mechanisms of reduction of the orthoquinone cofactor (PQQ) of MDH and sGDH involve concerted base-catalyzed proton abstraction from the hydroxyl moiety of methanol or from the 1-hydroxyl of glucose, and hydride equivalent transfer from the substrate to the quinone carbonyl carbon C5 of PQQ. The products of methanol and glucose oxidation are formaldehyde and glucolactone, respectively. The immediate product of PQQ reduction, PQQH- [-HC5(O-)-C4(=O)-] and PQQH [-HC5(OH)-C4(=O)-] converts to the hydroquinone PQQH2 [-C5(OH)=C4(OH)-]. The main focus is on MD structures of MDH * PQQ * methanol, MDH * PQQH-, MDH * PQQH, sGDH * PQQ * glucose, sGDH * PQQH- (glucolactone, and sGDH * PQQH. The reaction PQQ-->PQQH- occurs with Glu 171-CO2- and His 144-Im as the base species in MDH and sGDH, respectively. The general-base-catalyzed hydroxyl proton abstraction from substrate concerted with hydride transfer to the C5 of PQQ is assisted by hydrogen-bonding to the C5=O by Wat1 and Arg 324 in MDH and by Wat89 and Arg 228 in sGDH. Asp 297-COOH would act as a proton donor for the reaction PQQH(-)-->PQQH, if formed by transfer of the proton from Glu 171-COOH to Asp 297-CO2- in MDH. For PQQH-->PQQH2, migration of H5 to the C4 oxygen may be assisted by a weak base like water (either by crystal water Wat97 or bulk solvent, hydrogen-bonded to Glu 171-CO2- in MDH and by Wat89 in sGDH).

Envisioning the loop movements and rotation of the two subdomains of dihydrofolate reductase by elastic normal mode analysis.

Normal mode analysis, using the elastic network model, has been employed to envision the low frequency normal mode motion trends in the structures of five intermediates and a transition state in the kinetic pathway of E. coli dihydrofolate reductase (DHFR). Five of the reaction pathway analog structures and a crystal structure resembling the transition state, using X-ray analyses determined by Kraut et al., have been adapted as structural models. The motions that poise pathways of the M20 loop transitions from closed to occluded conformations and sub domain rotation to close the substrate cleft, have been predicted and envisioned for the first time by this study. Pathway entries to the movement of the substrate binding cleft helices are also envisioned. These motions play roles in transition structure stabilization and in regulating the release of the product tetrahydrofolate (THF). The motions observed push the ground state conformation of each intermediate towards a higher energy sub state conformation. A set of conserved residues involved in the catalytic reactions and conformational changes, previously studied by kinetic, theoretical and NMR, have been analyzed. The importance of these motions in terms of protein dynamics are revealed and envisioned by the normal mode analysis. Additional residues are proposed as candidates for further study of their potential promotional function.

Mechanism of product specificity of AdoMet methylation catalyzed by lysine methyltransferases: transcriptional factor p53 methylation by histone lysine methyltransferase SET7/9.

The catalysis by SET7/9 histone lysine methyltransferase of AdoMet N-methylation of the transcriptional factor p53-Lys4-NH 2 has been investigated with particular attention paid to the means of product specificity. After formation of the SET7/9.p53-Lys4-NH 3 (+).AdoMet complex, the following events occur: (i) the appearance of a water channel, (ii) a depronation of p53-Lys4-NH 3 (+) via this water channel into the aqueous solvent, and (iii) AdoMet methylation of p53-Lys4-NH 2 to form p53-Lys4-N(Me)H 2 (+). The formation of a water channel does not occur on formation of the SET7/9.p53-Lys4-NH 3 (+), SET7/9.p53-Lys4-N(Me)H 2 (+).AdoHcy, or SET7/9.p53-Lys4-N(Me)H 2 (+).AdoMet complex. Without a water channel, the substrate p53-Lys4-N(Me)H is not available because the proton dissociation p53-Lys4-N(Me)H 2 (+) --> p53-Lys4-N(Me)H + H (+) does not occur. The lack of formation of a water channel is due to the positioning of the methyl substituent of the SET7/9.p53-Lys4-N(Me)H 2 (+).AdoMet complex. By quantum mechanics/molecular mechanics, the computed free energy barrier of the methyl transfer reaction [p53-Lys4-NH 2 + AdoMet --> p53-Lys4-N(Me)H 2 (+) + AdoHcy] in the SET7/9 complex is Delta G (++) = 20.1 +/- 2.9 kcal/mol. This Delta G (++) is in agreement with the value of 20.9 kcal/mol calculated from the experimental rate constant (1.2 +/- 0.1 min (-1)). Our bond-order computations establish that the methyl transfer reaction in protein lysine methyltransferases occurs via a linear S N2 associative reaction with bond making of approximately 50%.

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)).

Catalytic mechanism and product specificity of rubisco large subunit methyltransferase: QM/MM and MD investigations.

Molecular dynamics (MD) simulations and hybrid quantum mechanics/molecular mechanics (QM/MM) calculations have been carried out in an investigation of Rubisco large subunit methyltransferase (LSMT). It was found that the appearance of a water channel is required for the stepwise methylation by S-adenosylmethionine (AdoMet). The water channel appears in the presence of AdoMet (LSMT.Lys-NH3+.AdoMet), but is not present immediately after methyl transfer (LSMT.Lys-N(Me)H2+.AdoHcy). The water channel allows proton dissociation from both LSMT.AdoMet.Lys-NH3+ and LSMT.AdoMet.Lys-N(Me)H2+. The water channel does not appear for proton dissociation from LSMT.AdoMet.Lys-N(Me)2H+, and a third methyl transfer does not occur. By QM/MM, the calculated free energy barrier of the first methyl transfer reaction catalyzed by LSMT (Lys-NH2 + AdoMet --> Lys-N(Me)H2+ + AdoHcy) is DeltaG++ = 22.8 +/- 3.3 kcal/mol. This DeltaG++ is in remarkable agreement with the value 23.0 kcal/mol calculated from the experimental rate constant (6.2 x 10-5 s-1). The calculated DeltaG++ of the second methyl transfer reaction (AdoMet + Lys-N(Me)H --> AdoHcy + Lys-N(Me)2H+) at the QM/MM level is 20.5 +/- 3.6 kcal/mol, which is in agreement with the value 22.0 kcal/mol calculated from the experimental rate constant (2.5 x 10-4 s-1). The third methyl transfer (Lys-N(Me)2 + AdoMet --> Lys-N(Me)3+ + AdoHcy) is associated with an allowed DeltaG++ of 25.9 +/- 3.2 kcal/mol. However, this reaction does not occur because a water channel does not form to allow the proton dissociation of Lys-N(Me)2H+. Future studies will determine whether the product specificity of lysine (mono, di, and tri) methyltransferases is determined by the formation of water channels.

Temperature-dependent structure of the E x S complex of Bacillus stearothermophilus alcohol dehydrogenase.

The catalytic chemistry of the thermophilic Bacillus stearothermophilus alcohol dehydrogenase (HtADH) closely resembles that of mesophilic horse liver alcohol dehydrogenase (HLADH). Molecular dynamics (MD) simulations of the htADH x NAD+ x EtO- complex at 298, 323, and 348 K show that the structure of the ligated Zn2+...EtO- complex varies slightly with change in temperature. The MD-created Boltzmann distribution of htADH x NAD+ x EtO- structures establishes the formation of multiple states which increase in number with a decrease in temperature. The motions of the cofactor domain are highly correlated with the motions of NAD+ at the optimal growth temperature (348 K), with NAD+ being pushed toward the substrate by Val260. With a decrease in temperature, the motion together of the cofactor and substrate is reversed, and at 298 K, the nicotinamide ring of the cofactor moves away from the substrate. Both the distance between and the angle of approach of C4 of NAD+ and HD of EtO- become distorted from those of the reactive conformation. The percentages of ground state present as the reactive conformation at different temperatures are approximately correlated with the kcat for the htADH enzymatic reaction. The rate constant for the htADH x NAD+ x EtOH --> htADH x NAD+ x EtO- proton dissociation, which is mediated by Thr40-OH, becomes slower at lower temperatures. The time-dependent distance between EtO- and Thr40-OH reveals that the Thr40 hydroxyl group sways between the substrate and NAD+ ribose 2'-hydroxyl group at the optimal enzyme growth temperature, and this movement is effectively frozen out as the temperature decreases. The temperature dependence of active site conformations is due to the change in both long-range and short-range motions of the E x S complex.

Diels-Alder ribozyme catalysis: a computational approach.

A computational comparison of the Diels-Alder reaction of a maleimide and an anthracene in water and the active site of the ribozyme Diels-Alderase is reported. During the course of the catalyzed reaction, the maleimide is held in the hydrophobic pocket while the anthracene approaches to the maleimide through the back passage of the active site. The active site is so narrow that the anthracene has to adopt a tilted approach angle toward maleimide. The conformation of the active site changes marginally at different states of the reaction. Active site dynamics contribution to catalysis has been ruled out. The active site stabilizes the product more than the transition state (TS). The reaction coordinates of the ribozyme reaction in TS, RC1-CD1 and RC4-CD2, are 2.35 and 2.33 A, respectively, compared to 2.37 and 2.36 A in water. The approach angle of anthracene toward maleimide is twisted by 18 degrees in the TS structure of ribozyme reaction while no twisted angle is found in TS of the reaction in water. The free energy barriers for reactions in both ribozyme and water were obtained by umbrella sampling combined with SCCDFTB/MM. The calculated free energy barriers for the ribozyme and water reactions are in good agreement with the experimental values. As expected, Mulliken charges of the atoms involved in the ribozyme reaction change in a similar manner as that of the reaction in water. The proficiency of the Diels-Alder ribozyme reaction originates from the active site holding the two reactants in reactive conformations, in which the reacting atoms are brought together in van der Waals distances and reactants approach to each other at an appropriate angle.

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.

Mechanism of methanol oxidation by quinoprotein methanol dehydrogenase.

At neutral pH, oxidation of CH(3)OH --> CH(2)O by an o-quinone requires general-base catalysis and the reaction is endothermic. The active-site -CO(2)(-) groups of Glu-171 and Asp-297 (Glu-171-CO(2)(-) and Asp-297-CO(2)(-)) have been considered as the required general base catalysts in the bacterial o-quinoprotein methanol dehydrogenase (MDH) reaction. Based on quantum mechanics/molecular mechanics (QM/MM) calculations, the free energy for MeOH reduction of o-PQQ when MeOH is hydrogen bonded to Glu-171-CO(2)(-) and the crystal water (Wat1) is hydrogen bonded to Asp-297-CO(2)(-) is DeltaG++ = 11.7 kcal/mol, which is comparable with the experimental value of 8.5 kcal/mol. The calculated DeltaG++ when MeOH is hydrogen bonded to Asp-297-CO(2)(-) is >50 kcal/mol. The Asp-297-CO(2)(-)...Wat1 complex is very stable. Molecular dynamics (MD) simulations on MDH.PQQ.Wat1 complex in TIP3P water for 5 ns does not result in interchange of Asp-297-CO(2)(-) bound Wat1 for a solvent water. Starting with Wat1 removed and MeOH hydrogen bonded to Asp-297-CO(2)(-), we find that MeOH returns to be hydrogen bonded to Glu-171-CO(2)(-) and Asp-297-CO(2)(-) coordinates to Ca(2+) during 3 ns simulation. The Asp-297-CO(2)(-)...Wat1 of reactant complex does play a crucial role in catalysis. By QM/MM calculation DeltaG++ = 1.1 kcal/mol for Asp-297-CO(2)(-) general-base catalysis of Wat1 hydration of the immediate CH(2)==O product --> CH(2)(OH)(2). By this means, the endothermic oxidation-reduction reaction is pulled such that the overall conversion of MeOH to CH(2)(OH)(2) is exothermic.