Near Attack Conformers Research Articles

Extein residues regulate the catalytic function of Spl DnaX intein enzyme by restricting the near-attack conformations of the active-site residues.

Intein enzymes catalyze the splicing of their flanking polypeptide chains and have found tremendous biotechnological applications. Their terminal residues form the catalytic core and participate in the splicing reaction. Hence, the neighboring N- and C-terminal extein residues influence the catalytic rate. As these extein residues vary depending on the substrate identity, we tested the influence of 20 amino acids at these sites in the Spl DnaX intein and observed significant variation of spliced product as well as N- and C-terminus cleavage product formation. We investigated the dependence of these reactions on the extein residues by molecular dynamics (MD) simulations on eight extein variants, and found that the conformational sampling of the active-site residues of the intein enzyme differed among these extein variants. We found that the extein variants that sample higher population of near-attack conformers (NACs) of the active-site residues undergo higher product formation in our activity assays. Ground state conformers that closely resemble the transition state are referred to as NACs. Very good correlation was observed between the NAC populations from the MD simulations of eight extein variants and the corresponding product formation from our activity assays. Furthermore, this molecular detail enabled us to elucidate the mechanistic roles of several conserved active-site residues in the splicing reaction. Overall, this study shows that the catalytic power of Spl DnaX intein enzyme, and most likely other inteins, depends on the efficiency of formation of NACs in the ground state, which is further modulated by the extein residues.

Broken-Symmetry Density Functional Theory Analysis of the Ω Intermediate in Radical S-Adenosyl-l-methionine Enzymes: Evidence for a Near-Attack Conformer over an Organometallic Species.

Radical S-adenosyl-l-methionine (SAM) enzymes are found in all domains of life and catalyze a wide range of biochemical reactions. Recently, an organometallic intermediate, Ω, has been experimentally implicated in the 5'-deoxyadenosyl radical generation mechanism of the radical SAM superfamily. In this work, we employ broken-symmetry density functional theory to evaluate several structural models of Ω. The results show that the calculated hyperfine coupling constants (HFCCs) for the proposed organometallic structure of Ω are inconsistent with the experiment. In contrast, a near-attack conformer of SAM bound to the catalytic [4Fe-4S] cluster, in which the distance between the unique iron and SAM sulfur is ∼3 Å, yields HFCCs that are all within 1 MHz of the experimental values. These results clarify the structure of the ubiquitous Ω intermediate and suggest a paradigm shift reversal regarding the mechanism of SAM cleavage by members of the radical SAM superfamily.

Rational modulation of the enzymatic intermediates for tuning the phosphatase activity of histidine kinase HK853

Histidine kinase (HK) of two-component signal transduction system (TCS) is a potential drug target for treating bacterial infections, and most HKs are bifunctional. We have previously identified the HXXXT motif of HK in HisKA subfamily to perform the phosphatase activity, but the specific working mechanism of the threonine is not well understood. In this paper, we use the phosphate group analog BeF3− to capture the enzymatic intermediates between HK853 and RR468 from Thermotoga maritima during dephosphorylation, and demonstrate that the T264 site is essential for populating capable near attack conformers (NAC) between enzyme and substrate to facilitate catalysis. Importantly, mutations at this site can modulate the phosphatase activity of HK. Our results help to understand the TCS signal transduction mechanisms and provide a reference for drug design.

Simple Design of an Enzyme-Inspired Supported Catalyst Based on a Catalytic Triad

Summary Enzyme active sites afford an intricate interplay of functional groups to mediate complex organic and inorganic reactions. Many hydrolytic enzymes use a catalytic triad comprising three different functional residues—(Ser(-OH), Hist(-imidazole), Asp(-CO 2 H))—that catalyze the hydrolysis of numerous unique substrates. Inspired by this design, we have developed a simple one-step synthesis for preparing a new supported catalytic system in which the three reactive groups of the catalytic triad (alcohol, imidazole, and carboxylate) are incorporated into a single functional unit. These artificial active sites can be coupled to a solid-phase support (Merrifield resin) by copper(I)-catalyzed azide-alkyne cycloaddition "click chemistry," and their effectiveness as esterolysis catalysts was demonstrated. Furthermore, tuning the local hydrophobicity of the resin particles with an approach analogous to the native enzyme hydrophobic pocket increased the catalytic efficiency. Quantum mechanics and molecular dynamics computational modeling were used to probe the catalytic effect and suggested a concerted two-step mechanism and hydrophobic nanoenvironment similar to that of hydrolytic enzymes.

Understanding protein arginine methyltransferase 1 (PRMT1) product specificity from molecular dynamics

Protein arginine methyltransferases (PRMTs) catalyze the post-translational methylation of specific arginyl groups within targeted proteins to regulate fundamental biological responses in eukaryotic cells. The major Type I PRMT enzyme, PRMT1, strictly generates monomethyl arginine (MMA) and asymmetric dimethylarginine (ADMA), but not symmetric dimethylarginine (SDMA). Multiple diseases can arise from the dysregulation of PRMT1, including heart disease and cancer, which underscores the need to elucidate the origin of product specificity. Molecular dynamics (MD) simulations were carried out for WT PRMT1 and its M48F, H293A, H293S, and H293S-M48F mutants bound with S-adenosylmethionine (AdoMet) and the arginine substrate in an unmethylated or methylated form. Experimental site-directed mutagenesis and analysis of the resultant products were also performed. Two specific PRMT1 active site residues, Met48 and His293, have been determined to play a key role in dictating product specificity, as: (1) the single mutation of Met48 to Phe enabled PRMT1 to generate MMA, ADMA, and a limited amount of SDMA; (2) the single mutation of His293 to Ser formed the expected MMA and ADMA products only; whereas (3) the double mutant H293S-M48F-PRMT1 produced SMDA as the major product with limited amounts of MMA and ADMA. Calculating the formation of near-attack conformers resembling SN2 transition states leading to either the ADMA or SDMA products finds that Met48 and His293 may enable WT PRMT1 to yield ADMA exclusively by precluding MMA from binding in an orientation more conducive to SDMA formation, i.e., the methyl group bound at the arginine Nη2 position.

Computational design of a lipase for catalysis of the Diels-Alder reaction

Combined molecular docking, molecular dynamics (MD) and density functional theory (DFT) studies have been employed to study catalysis of the Diels-Alder reaction by a modified lipase. Six variants of the versatile enzyme Candida Antarctica lipase B (CALB) have been rationally engineered in silico based on the specific characteristics of the pericyclic addition. A kinetic analysis reveals that hydrogen bond stabilization of the transition state and substrate binding are key components of the catalytic process. In the case of substrate binding, which has the greater potential for optimization, both binding strength and positioning of the substrates are important for catalytic efficiency. The binding strength is determined by hydrophobic interactions and can be tuned by careful selection of solvent and substrates. The MD simulations show that substrate positioning is sensitive to cavity shape and size, and can be controlled by a few rational mutations. The well-documented S105A mutation is essential to enable sufficient space in the vicinity of the oxyanion hole. Moreover, bulky residues on the edge of the active site hinders the formation of a sandwich-like nearattack conformer (NAC), and the I189A mutation is needed to obtain enough space above the face of the α,β-double bond on the dienophile. The double mutant S105A/I189A performs quite well for two of three dienophiles. Based on binding constants and NAC energies obtained from MD simulations combined with activation energies from DFT computations, relative catalytic rates (v(cat)/v(uncat)) of up to 103 are predicted.

Molecular basis for the high-affinity binding and stabilization of firefly luciferase by PTC124

Firefly luciferase (FLuc), an ATP-dependent bioluminescent reporter enzyme, is broadly used in chemical biology and drug discovery assays. PTC124 (Ataluren; (3-[5-(2-fluorophenyl)-1,2,4-oxadiazol-3-yl]benzoic acid) discovered in an FLuc-based assay targeting nonsense codon suppression, is an unusually potent FLuc-inhibitor. Paradoxically, PTC124 and related analogs increase cellular FLuc activity levels by posttranslational stabilization. In this study, we show that FLuc inhibition and stabilization is the result of an inhibitory product formed during the FLuc-catalyzed reaction between its natural substrate, ATP, and PTC124. A 2.0 A cocrystal structure revealed the inhibitor to be the acyl-AMP mixed-anhydride adduct PTC124-AMP, which was subsequently synthesized and shown to be a high-affinity multisubstrate adduct inhibitor (MAI; K(D) = 120 pM) of FLuc. Biochemical assays, liquid chromatography/mass spectrometry, and near-attack conformer modeling demonstrate that formation of this novel MAI is absolutely dependent upon the precise positioning and reactivity of a key meta-carboxylate of PTC124 within the FLuc active site. We also demonstrate that the inhibitory activity of PTC124-AMP is relieved by free coenzyme A, a component present at high concentrations in luciferase detection reagents used for cell-based assays. This explains why PTC124 can appear to increase, instead of inhibit, FLuc activity in cell-based reporter gene assays. To our knowledge, this is an unusual example in which the "off-target" effect of a small molecule is mediated by an MAI mechanism.

Molecular dynamics studies of ground state and intermediate of the hyperthermophilic indole-3-glycerol phosphate synthase.

Indole-3-glycerol phosphate synthase catalyzes the terminal ring closure step in tryptophan biosynthesis. In this paper, we compare the results from molecular dynamics (MD) simulations of enzyme-bound substrate at 298, 333, 363, and 385 K and the enzyme-bound intermediate at 385 K, solvated in TIP3P water box with a CHARMM force field. Results from MD simulations agree with experimental studies supporting the observation that Lys-110 is the general acid. Based on its location in the active site during the MD simulations, Glu-210 warrants classification as the general base instead of the previously proposed Glu-159. We find that the relative population of the reactive enzyme-substrate Michaelis conformers [near attack conformers (NACs)] with temperature correlates well (correlation coefficient of 0.96) with the relative activity of this thermophilic enzyme. At higher temperature, the enzyme-substrate electrostatic interaction favors the binding of the substrate in NAC conformation, whereas, at lower temperature, the substrate is distorted and bound in a nonreactive conformation. This change is reflected in the approximately 1,100-fold increase in population of NACs at 385 K relative to 298 K. The easily determined population of NACs at given temperature tells much about the thermophilic property of the enzyme. Thus, the hyperthermophilic enzyme has evolved to have optimum activity at high temperatures, and, with lowering of the temperature, the electrostatic interaction at the active site is enhanced and the structure is deformed. This model can be regarded as a general explanation for the activity of hyperthermophilic enzymes.

Anticorrelated motions as a driving force in enzyme catalysis: the dehydrogenase reaction.

Molecular dynamics and cross-correlation analysis of the horse liver alcohol dehydrogenase HLADH.NAD(+).PhCH(2)O(-) complex has established anticorrelated motions between the NAD(+)-binding domain and other portions of the enzyme. Four pairs of anticorrelated interactions are (i and ii) cofactor-binding domain: C(alpha) of V292 and the CG1 of V203 with C7 of PhCH(2)O(-); (iii) cofactor-binding domain: amide carbonyl oxygen of I318 with amide N of H67; and (iv) cofactor domain: C(alpha) of T178 with carbonyl oxygen of L141. The average distances between pairs are 9.2 A for i, 8.2 A for ii, 14.7 A for iii, and 18.2 A for iv. The motions of i and ii are most important in the approximately 0.5 A pushing of C4 of NAD(+) toward C7 of PhCH(2)O(-) to form push near-attack conformer (NACs). The motions of iv are less so, and those of iii are not important. Seventy-five quantum mechanics/molecular mechanics calculations of the energies of reaction were carried out without structural restrictions from different stages of the molecular dynamics trajectory. Of the 71 conformations, the 29 fulfilling NAC criteria were associated with the lowest energies of activation. Thus, anticorrelated motions from the NAD(+)-binding domain by way of the neighboring V292 and V203 have a pushing motion, which moves the C4 of NAD(+) toward the H-C7 of the substrate. Longer-range anticorrelated motions involving the cofactor-binding domain have no or very little influence on NAC formation.

Molecular recognition of aldehydes by aldehyde dehydrogenase and mechanism of nucleophile activation.

Experimental structural data on the state of substrates bound to class 3 Aldehyde Dehydrogenases (ALDH3A1) is currently unknown. We have utilized molecular mechanics (MM) simulations, in conjunction with new force field parameters for aldehydes, to study the atomic details of benzaldehyde binding to ALDH3A1. Our results indicate that while the nucleophilic Cys243 must be in the neutral state to form what are commonly called near-attack conformers (NACs), these structures do not correlate with increased complexation energy calculated with the MM-Generalized Born Molecular Volume (GBMV) method. The negatively charged Cys243 (thiolate form) of ALDH3A1 also binds benzaldehyde in a stable conformation but in this complex the sulfur of Cys243 is oriented away from benzaldehyde yet yields the most favorable MM-GBMV complexation energy. The identity of the general base, Glu209 or Glu333, in ALDHs remains uncertain. The MM simulations reveal structural and possible functional roles for both Glu209 and Glu333. Structures from the MM simulations that would support either glutamate residue as the general base were further examined with Hybrid Quantum Mechanical (QM)/MM simulations. These simulations show that, with the PM3/OPLS potential, Glu209 must go through a step-wise mechanism to activate Cys243 through an intervening water molecule while Glu333 can go through a more favorable concerted mechanism for the same activation process.

Transition state stabilization by general acid catalysis, water expulsion, and enzyme reorganization in Medicago savita chalcone isomerase.

In aqueous solution, Medicago savita chalcone isomerase (CHI) enhances the reaction rate for the unimolecular rearrangement of chalcone (CHN) into flavanone by seven orders of magnitude. Conformations of CHN and their relative free energies in water and CHI were investigated by the thermodynamic perturbation method. In water, CHN adopts two conformations (I and II) with conformation I being higher in energy than conformation II by 3 kcal/mol. Only I can give rise to a near attack conformer (NAC) where the nucleophile O2' and the electrophile C9 are placed in proximity. In CHI, I binds less tightly than II by approximately 2 kcal/mol, resulting in the free energy for NAC formation being approximately 2 kcal/mol higher in the enzyme than in water. This unfavorable feature in the ground state of the CHI reaction requires the predominant catalytic advantage to be taken in the step of NAC --> transition state (TS). From the molecular dynamics simulations of apo-CHI, CHI complexed with CHN (CHI.CHN) and CHI.TS, we found: (i) Lys-97-general-acid catalysis of the O2'(-) nucleophilic addition; (ii) expulsion of three water molecules in the process of TS formation; (iii) release of enzyme structural distortion on TS formation. In the conclusion, CHI's remarkable efficiency of stabilizing the TS and its relatively poor ability in organizing the ground state is compared with chorismate mutase whose catalytic prowess, when compared with water, originates predominantly from the enhanced NAC population at the active site.

A substrate variant as a high-affinity, reversible inhibitor: insight from the X-ray structure of cilastatin bound to membrane dipeptidase

An analysis of the X-ray structure of cilastatin bound to membrane dipeptidase, together with docking studies, is presented here to reveal how a simple amide may act as a high-affinity, reversible, amidase inhibitor. Cilastatin binds as a normal substrate and is orientated in a perfect near-attack conformer for formation of a tetrahedral intermediate with the zinc-bound water/hydroxide. This intermediate is fated, however, only to revert to its starting components as scission of the amide bond is prevented by the precise fit of cilastatin within the active site. The cilastatin alkyl end groups that are tightly buttressed against amino acid residues on opposite sides of the active site, are aligned along the C–N reaction coordinate axis thereby preventing collapse of the intermediate via rupture of the C–N bond. Such a feature could have more general applicability in the explicit design of substrate variants as selective, tight-binding, and reversible inhibitors.

Reaction mechanism of soluble epoxide hydrolase: insights from molecular dynamics simulations.

Molecular dynamics simulations have been performed to gain insights into the catalytic mechanism of the hydrolysis of epoxides to vicinal diols by soluble epoxide hydrolase (sEH). The binding of a substrate, 1S,2S-trans-methylstyrene oxide, was studied in two conformations in the active site of the enzyme. It was found that only one is likely to be found in the active enzyme. In the preferred conformation the phenyl group of the substrate is pi-sandwiched between two aromatic residues, Tyr381 and His523, whereas the other conformation is pi-stacked with only one aromatic residue, Trp334. Two simulations were carried out to 1 ns for each conformation to evaluate the protonation state of active site residue His523. It was found that a protonated histidine is essential for keeping the active site from being disrupted. Long time scale, 4 ns, molecular dynamics simulation was done for the structure with the most likely combination of binding conformation and protonation state of His523. Near Attack Conformers (NACs) are present 5.3% of the time and nucleophilic attack on either epoxide carbon atom, approximately 75% on C(1) and approximately 25% on C(2), is found. A maximum of one hydrogen bond between the epoxide oxygen and either of the active site tyrosines, Tyr465 and Tyr381, is present, in agreement with experimental mutagenesis results that reveal a slight loss in activity if one tyrosine is mutated and essential loss of all activity upon double mutation of the two tyrosines in question. It was found that a hydrogen bond from Tyr465 to the substrate oxygen is essential for controlling the regioselectivity of the reaction. Furthermore, a relationship between the presence of this hydrogen bond and the separation of reactants was found. Two groups of amino acid segments were identified each as moving collectively. Furthermore, an overall anti-correlation was found between the movements of these two individually collectively moving groups, made up by parts of the cap-region, including the two tyrosines, and the site of the catalytic triad, respectively. This overall anti-correlated collective domain motion is, perhaps, involved in the conversion of E.NAC to E.TS.

Protein Engineering of Nitrile Hydratase Activity of Papain: Molecular Dynamics Study of a Mutant and Wild-Type Enzyme

The mechanism of hydrolysis of the nitrile (N-acetyl-phenylalanyl-2-amino-propionitrile, I) catalyzed by Gln19Glu mutant of papain has been studied by nanosecond molecular dynamics (MD) simulations. MD simulations of the complex of mutant enzyme with I and of mutant enzyme covalently attached to both neutral (II) and protonated (III) thioimidate intermediates were performed. An MD simulation with the wild-type enzyme.I complex was undertaken as a reference. The ion pair between protonated His159 and thiolate of Cys25 is coplanar, and the hydrogen bonding interaction S(-)(25).HD1-ND1(159) is observed throughout MD simulation of the mutant enzyme.I complex. Such a sustained hydrogen bond is absent in nitrile-bound wild-type papain due to the flexibility of the imidazole ring of His159. The nature of the residue at position 19 plays a critical role in the hydrolysis of the covalent thioimidate intermediate. When position 19 represents Glu, the imidazolium ion of His159-ND1(+).Cys25-S(-) ion pair is distant, on average, from the nitrile nitrogen of substrate I. Near attack conformers (NACs) have been identified in which His159-ImH(+) is positioned to initiate a general acid-catalyzed addition of Cys-S(-) to nitrile. Though Glu19-CO(2)H is distant from nitrile nitrogen in the mutant.I structure, MD simulations of the mutant.II covalent adduct finds Glu19-CO(2)H hydrogen bonded to the thioimide nitrogen of II. This hydrogen bonded species is much less stable than the hydrogen bonded Glu19-CO(2)(-) with mutant-bound protonated thioimidate (III). This observation supports Glu19-CO(2)H general acid catalysis of the formation of mutant.III. This is the commitment step in the Gln19Glu mutant catalysis of nitrile hydrolysis.

The mechanism of cis-trans isomerization of prolyl peptides by cyclophilin.

The mechanism of cis-trans isomerization of prolyl peptides catalyzed by cyclophilin (CyP) was studied computationally via molecular dynamics (MD) simulations of the transition state (TS) and the cis and trans forms of the ground state (GS), when bound to CyP and when free in aqueous solution. The MD simulations include four enzyme-bound species of tetrapeptide (Suc-Ala-XC([double bond]O)-NPro-Phe-pNA; X = Gly, Trp, Ala, and Leu). In water, the prolyl amide bond is favorably planar with the presence of conformers exhibiting +/-20 degrees twist of the C-N dihedral. In the active site a hydrogen bond between the cis-prolyl amide carbonyl O and the backbone amide N-H of Asn102 retains the 20 degrees twist of the C-N dihedral. The TS structure is characterized by a 90 degrees twist of the amide C-N bond and a more favorable interaction with Asn102 due to the shorter distance between Asn102(HN) and the amide carbonyl O. The conformational change of cis --> TS also involves pyramidalization of the amide N, which results in the formation of a hydrogen bond between the amide N and the guanidino group of Arg55. Both Asn102 and Arg55 are held in the same position in CyP.cis-isomer as in CyP.TS. In the ligand-free CyP the Arg55 guanidino group is highly disorganized and Asn102 is displaced 1 A from the position in the ligand-bound CyP. Thus, the organization of Arg55 and Asn102 occurs upon substrate binding. The geometrical complimentarity of the organized enzyme structure to the TS structure is a result of preferential binding of the proline N and the amide carbonyl of the TS compared to that of GS. However, the N-terminal part (Suc-Ala) becomes repositioned in the TS such that two hydrogen bonds disappear, one hydrogen bond appears and two other hydrogen bonds becomes weaker on the conversion of CyP.cis to CyP.TS. During this conversion, total hydrophobic contact between enzyme and the peptide is preserved. Thus, the interaction energies of GS and TS with enzyme are, as a whole, much alike. This does not support the contention that TS is bound more tightly than GS by K(m)/K(TS) = 10(6) in the cis --> trans reaction. Repositioning of the N-terminal part of the peptide on CyP.TS formation becomes more pronounced when the substrate X residue is changed from Gly < Trp < Ala < Leu. We propose that the larger turning of the N-terminus is responsible for the larger value of the experimentally observed Delta S(++) and Delta H(++), which sum up to little change in Delta G(++). The positioning of the Arg55 and the degree of 20 degrees twist of the amide C-N bond are considered as criteria for Near Attack Conformers (NACs) in cis-trans isomerization. NACs account for approximately 30% of the total GS populations of the cis-isomer. Similar NAC populations were observed with four different substrates. This is consistent with the insensitivity of enzymatic activity to the nature of the X residue. Also, the NAC population in CyP.trans-AAPF was comparable to that in CyP.cis-AAPF, in accord with similar experimentally measured rates of the cis --> trans and trans --> cis reaction in CyP. These NACs, found in CyP.cis and CyP.trans, resemble only one of the four possible TS configurations in the water reaction. The identity of this TS structure (syn/exo) is in accord with experimentally determined KIE values in the enzymatic reaction. However, the geometry of the active site was also complementary to another TS structure (anti/exo) that was not detected in the active site by the same KIE measurements, implying that the geometrical fitness of the TS cannot be a single determining factor for enzymatic reactions.

Dynamic structures of horse liver alcohol dehydrogenase (HLADH): results of molecular dynamics simulations of HLADH-NAD(+)-PhCH(2)OH, HLADH-NAD(+)-PhCH(2)O(-), and HLADH-NADH-PhCHO.

Molecular dynamics simulations of the oxidation of benzyl alcohol by horse liver alcohol dehydrogenase (HLADH) have been carried out. The following three states have been studied: HLADH.PhCH(2)OH.NAD(+) (MD1), HLADH.PhCH(2)O(-).NAD(+) (MD2), and HLADH.PhCHO.NADH (MD3). MD1, MD2, and MD3 simulations were carried out on one of the subunits of the dimeric enzyme covered in a 32-A-radius sphere of TIP3P water centered on the active site. The proton produced on ionization of the alcohol when HLADH.PhCH(2)OH.NAD(+) --> HLADH.PhCH(2)O(-).NAD(+) is transferred from the active site to solvent water via a hydrogen bonding network consisting of serine48 hydroxyl, ribose 2'- and 3'-hydroxyl groups, and Hist51. Hydrogen bonding of the 3'OH of ribose to Ile269 carbonyl maintains this proton in position to be transferred to water. Molecular dynamic simulations have been employed to track water1287 from the TIP3 water pool to the active site, thus exhibiting the mode of entrance of water to the active site. With time the water1287 accumulates in two different positions in order to accept the proton from the ribose 3'-OH and from His51. There can be identified two structural substates for proton passage. In the first substate the imidazole Ne2 of His51 is adjacent to the nicotinamide ribose C2'-OH and hydrogen bonding distances for proton transfer through the hydrogen bonded relay series PhCH(2)OH...Ser48-OH...Ribose2'-OH...His51...OH(2) (path 1) average 2.0, 2.0, and 2.1 A and (for His51...OH(2)) minimal distances less or equal to 2.5 A. The structure for path 1 is present 20% of the time span. And in the second substate, there are two possible proton passages: path 1 as before and path 2. Path 2 involves the hydrogen-bonded relay series PhCH(2)OH...Ser48-OH...Ribose2'-OH...Ribose3'-OH...His51.OH(2) with the average bonding distances being 2.0, 2.0, 2.1, and 2.0 A and (for His51...OH(2)) minimal distances less or equal to 2.5 A (20% probability of the time span), respectively. During the molecular dynamics simulation the NAD(+) ribose conformations have stabilized at the C2'-endo-C3'-exo or the C2'-endo conformations. With the C2'-endo conformation the first and second substates are able to persist for different time spans, while with the C2'-endo-C3'-exo conformation the only possible pathway involves the first substate. For both first and second substates the fluctuation of the distances between the ribose-OH protons and N epsilon 2 of His51 imidazole ring is partially contributed by the "windshield wiper" motion of the His51 imidazole ring. Since the imidazole of His-51 contributes only about 10-fold to activity, as estimated from the decrease in activity upon substitution with a Gln, there must be an alternate route for the proton to pass to solvent without going through this histidine. A third pathway involves ribose C3'-OH and Ile-269. In MD2, near attack conformers (NACs) for hydride transfer from PhCH(2)O(-) to NAD(+) represent approximately 60% of E.S conformers. The molecular dynamic study of MD3 at mildly basic pH reveals that reactive ground state conformers (NACs) for hydride transfer from NADH to PhCHO amount to 12 mol % of conformers. In MD3, anisotropic bending of the dihydronicotinamide ring of NADH (average value of alpha(c) = 4.0 degrees and alpha(n) = 0.5 degrees, respectively) is observed.

Comparison of the Dynamics for Ground-State and Transition-State Structures in the Active Site of Catechol O-Methyltransferase

Three 1-ns molecular dynamics (MD) simulations have been used to study the active-site structure of catechol O-methyltransferase when containing catechol, catecholate, or the transition state. It was found that the physical properties of the enzyme active sites are very similar when containing either catecholate or the transition state. The calculated root-mean-squared deviation and positional fluctuations of the active-site residues within 10 Å of the methyl group of S-adenosyl-l-methionine (AdoMet) for the catecholate and transition-state simulations are similar, and the average solvent accessible surface areas for these residues are almost identical. It was found in the ground-state simulation with catecholate that interactions between AdoMet and the enzyme are critical in positioning AdoMet into near attack conformers (NACs) which resemble the transition-state geometry. The CB of Tyr68 influences the distance between the methyl carbon of AdoMet and the ionized oxygen of catecholate. Interactions between the backbone amide carbonyls of Met40 and Asp141 control the angular positioning of the AdoMet relative to catecholate. When the interactions dissipated, the angle formed by the sulfur and methyl carbon of AdoMet and ionized oxygen of catecholate decreased and no longer fulfilled the NAC criteria. Comparisons of these 3 MD simulations in combination with results from previous quantum mechanical calculations from this lab suggest that the catalytic power of this enzyme in going from E·S to E·TS is not due to transition-state stabilization but from the ability of the active site of catechol O-methyltransferase to arrange the reactants into conformers that closely resemble the transition state. Overall, in E·S going to E·TS one would also consider as a driving force catecholate and AdoMet desolvation.