S-adenosylmethionine conformations in solution and in protein complexes: conformational influences of the sulfonium group

Engineering of Methionine Adenosyltransferase toward Mitigated Product Inhibition for Efficient Production of S-Adenosylmethionine

S-Adenosylmethionine (SAM) is a crucial metabolic intermediate playing irreplaceable roles in organismal activities. However, the synthesis of SAM by methionine adenosyltransferase (MAT) is hindered by low conversion due to severe product inhibition. Herein structure-guided semirational engineering was conducted on MAT from Escherichia coli (EcMAT) to mitigate the product inhibitory effect. Compared with the wild-type EcMAT, the best variant E56Q/Q105R exhibited an 8.13-fold increase in half maximal inhibitory concentration and a 4.46-fold increase in conversion (150 mM ATP and l-methionine), leading to a SAM titer of 47.02 g/L. Another variant, E56N/Q105R, showed superior thermostability with an impressive 85.30-fold increase in half-life (50 °C) value. Furthermore, molecular dynamics (MD) simulation results demonstrate that the alleviation in product inhibitory effect could be attributed to facilitated product release. This study offers molecular insights into the mitigated product inhibition, and provides valuable guidance for engineering MAT toward enhanced catalytic performance.

An optimized purification protocol for enzymatically synthesized S-adenosyl-L-methionine (SAM) for applications in solution state infrared spectroscopic studies

S-adenosyl-L-methionine (SAM) is an abundant biomolecule used by methyltransferases to regulate a wide range of essential cellular processes such as gene expression, cell signaling, protein functions, and metabolism. Despite considerable effort, there remain many specificity challenges associated with designing small molecule inhibitors for methyltransferases, most of which exhibit off-target effects. Interestingly, NMR evidence suggests that SAM undergoes conformeric exchange between several states when free in solution. Infrared spectroscopy can detect different conformers of molecules if present in appreciable populations. When SAM is noncovalently bound within enzyme active sites, the nature and the number of different conformations of the molecule are likely to be altered from when it is free in solution. If there are unique structures or different numbers of conformers between different methyltransferase active sites, solution-state information may provide promising structural leads to increase inhibitor specificity for a particular methyltransferase. Toward this goal, frequencies measured in SAM's infrared spectra must be assigned to the motions of specific atoms via isotope incorporation at discrete positions. The incorporation of isotopes into SAM's structure can be accomplished via an established enzymatic synthesis using isotopically labeled precursors. However, published protocols produced an intense and highly variable IR signal which overlapped with many of the signals from SAM rendering comparison between isotopes challenging. We observed this intense absorption to be from co-purifying salts and the SAM counterion, producing a strong, broad signal at 1100 cm-1. Here, we report a revised SAM purification protocol that mitigates the contaminating salts and present the first IR spectra of isotopically labeled CD3-SAM. These results provide a foundation for isotopic labeling experiments of SAM that will define which atoms participate in individual molecular vibrations, as a means to detect specific molecular conformations.

Comparative Analysis of Sulfonium-π, Ammonium-π, and Sulfur-π Interactions and Relevance to SAM-Dependent Methyltransferases

We report the measurement and analysis of sulfonium-π, thioether-π, and ammonium-π interactions in a β-hairpin peptide model system, coupled with computational investigation and PDB analysis. These studies indicated that the sulfonium-π interaction is the strongest and that polarizability contributes to the stronger interaction with sulfonium relative to ammonium. Computational studies demonstrate that differences in solvation of the trimethylsulfonium versus the trimethylammonium group also contribute to the stronger sulfonium-π interaction. In comparing sulfonium-π versus sulfur-π interactions in proteins, analysis of SAM- and SAH-bound enzymes in the PDB suggests that aromatic residues are enriched in close proximity to the sulfur of both SAM and SAH, but the populations of aromatic interactions of the two cofactors are not significantly different, with the exception of the Me-π interactions in SAM, which are the most prevalent interaction in SAM but are not possible for SAH. This suggests that the weaker interaction energies due to loss of the cation-π interaction in going from SAM to SAH may contribute to turnover of the cofactor.

Restriction of S-adenosylmethionine conformational freedom by knotted protein binding sites

S-adenosylmethionine (SAM) is one of the most important enzyme substrates. It is vital for the function of various proteins, including large group of methyltransferases (MTs). Intriguingly, some bacterial and eukaryotic MTs, while catalysing the same reaction, possess significantly different topologies, with the former being a knotted one. Here, we conducted a comprehensive analysis of SAM conformational space and factors that affect its vastness. We investigated SAM in two forms: free in water (via NMR studies and explicit solvent simulations) and bound to proteins (based on all data available in the PDB and on all-atom molecular dynamics simulations in water). We identified structural descriptors—angles which show the major differences in SAM conformation between unknotted and knotted methyltransferases. Moreover, we report that this is caused mainly by a characteristic for knotted MTs compact binding site formed by the knot and the presence of adenine-binding loop. Additionally, we elucidate conformational restrictions imposed on SAM molecules by other protein groups in comparison to conformational space in water.

The topology of a folded polypeptide chain has great impact on the resulting protein function and its interaction with ligands. Interestingly, topological constraints appear to affect binding of one of the most ubiquitous substrates in the cell, S-adenosylmethionine (SAM), to its target proteins. Here, we demonstrate how binding sites of specific proteins restrict SAM conformational freedom in comparison to its unbound state, with a special interest in proteins with non-trivial topology, including an exciting group of knotted methyltransferases. Using a vast array of computational methods combined with NMR experiments, we identify key structural features of knotted methyltransferases that impose unorthodox SAM conformations. We compare them with the characteristics of standard, unknotted SAM binding proteins. These results are significant for understanding differences between analogous, yet topologically different enzymes, as well as for future rational drug design.

Nonstationary Two-Dimensional Nuclear Magnetic Resonance: A Method for Studying Reaction Mechanisms in Situ

Nuclear magnetic resonance spectroscopy (NMR) is a versatile tool of chemical analysis allowing one to determine structures of molecules with atomic resolution. Particularly informative are two-dimensional (2D) experiments that directly identify atoms coupled by chemical bonds or a through-space interaction. Thus, NMR could potentially be powerful tool to study reactions in situ and explain their mechanisms. Unfortunately, 2D NMR is very time-consuming and thus often cannot serve as a "snapshot" technique for in situ reaction monitoring. Particularly difficult is the case of spectra, in which resonance frequencies vary in the course of reaction. This leads to resolution and sensitivity loss, often hindering the detection of transient products. In this paper we introduce a novel approach to correct such nonstationary 2D NMR signals and raise the detection limits over 10 times. We demonstrate success of its application for studying the mechanism of the reaction of AgSO₄-induced synthesis of diphenylmethane-type compounds. Several reactions occur in the studied mixture of benzene and toluene, all with rather low yield and leading to compounds with similar chemical shifts. Nevertheless, with the use of a proposed 2D NMR approach we were able to describe complex mechanisms of diphenylmethane formation involving AgSO₄-induced toluene deprotonation and formation of benzyl carbocation, followed by nucleophilic attacks.

Smyd2 conformational changes in response to p53 binding: role of the C-terminal domain

Smyd2 lysine methyltransferase regulates monomethylation of histone and nonhistone lysine residues using S‐adenosylmethionine cofactor as the methyl donor. The nonhistone interactors include several tumorigenic targets, including p53. Understanding this interaction would allow the structural principles that underpin Smyd2‐mediated p53 methylation to be elucidated. Here, we performed μ‐second molecular dynamics (MD) simulations on binary Smyd2‐cofactor and ternary Smyd2‐cofactor‐p53 peptide complexes. We considered both unmethylated and monomethylated p53 peptides (at Lys370 and Lys372). The results indicate that (a) the degree of conformational freedom of the C‐terminal domain of Smyd2 is restricted by the presence of the p53 peptide substrate, (b) the Smyd2 C‐terminal domain shows distinct dynamic properties when interacting with unmethylated and methylated p53 peptides, and (c) Lys372 methylation confines the p53 peptide conformation, with detectable influence on Lys370 accessibility to the cofactor. These MD results are therefore of relevance for studying the biology of p53 in cancer progression.

Structural and Functional Characterization of Sulfonium Carbon-Oxygen Hydrogen Bonding in the Deoxyamino Sugar Methyltransferase TylM1

The N-methyltransferase TylM1 from Streptomyces fradiae catalyzes the final step in the biosynthesis of the deoxyamino sugar mycaminose, a substituent of the antibiotic tylosin. The high-resolution crystal structure of TylM1 bound to the methyl donor S-adenosylmethionine (AdoMet) illustrates a network of carbon-oxygen (CH···O) hydrogen bonds between the substrate's sulfonium cation and residues within the active site. These interactions include hydrogen bonds between the methyl and methylene groups of the AdoMet sulfonium cation and the hydroxyl groups of Tyr14 and Ser120 in the enzyme. To examine the functions of these interactions, we generated Tyr14 to phenylalanine (Y14F) and Ser120 to alanine (S120A) mutations to selectively ablate the CH···O hydrogen bonding to AdoMet. The TylM1 S120A mutant exhibited a modest decrease in its catalytic efficiency relative to that of the wild type (WT) enzyme, whereas the Y14F mutation resulted in an approximately 30-fold decrease in catalytic efficiency. In contrast, site-specific substitution of Tyr14 by the noncanonical amino acid p-aminophenylalanine partially restored activity comparable to that of the WT enzyme. Correlatively, quantum mechanical calculations of the activation barrier energies of WT TylM1 and the Tyr14 mutants suggest that substitutions that abrogate hydrogen bonding with the AdoMet methyl group impair methyl transfer. Together, these results offer insights into roles of CH···O hydrogen bonding in modulating the catalytic efficiency of TylM1.

Improved NOE fitting for flexible molecules based on molecular mechanics data – a case study with S-adenosylmethionine

Using the important biomolecule S-adenosyl methionine as an exemplar, we provide a new, enhanced approach for fitting MD data to high-accuracy NOE data, providing improvements in structure determination.

Computational Study of the Degradation of S-Adenosyl Methionine in Water

The degradation of S-adenosyl methionine (SAM) to homoserine-γ-lactone (HSL) and methyltioadenine (MTA) in water is studied with MD simulations. The AM1 Hamiltonian is used for the quantum part and the flexible AMBER force field for the H₂O molecules. The MD simulations predict the free energy barrier for the degradation reaction to be between 109 and 112 kJ mol^-1 and an overall gain in free energy of -26 kJ mol^-1. The high barrier and the low energy gain of this reaction can be linked to interactions among the carboxylate group of the SAM molecule and solvent H₂O molecules, which are not observed on the product side. Hence, the H₂O molecules effectively slow down the reaction that otherwise would be much faster.

Conformational Dynamics of Lysine Methyltransferase Smyd2. Insights into the Different Substrate Crevice Characteristics of Smyd2 and Smyd3

Smyd2, the SET and MYND domain containing protein lysine methyltransferase, targets histone and nonhistone substrates. Methylation of nonhistone substrates has direct implications in cancer development and progression. Dynamic regulation of Smyd2 activity and the structural basis of broad substrate specificity still remain elusive. Herein, we report on extensive molecular dynamics simulations on a full length Smyd2 in the presence and absence of AdoMet cofactor (covering together 1.3 μs of sampling), and the accompanying conformational transitions. Additionally, dynamics of the C-terminal domain (CTD) and structural features of substrate crevices of Smyd2 and Smyd3 are compared. The CTD of Smyd2 exhibits conformational flexibility in both states. In the holo form, however, it undergoes larger hinge motions resulting in more opened configurations than the apo form, which is confined around the partially open starting X-ray configuration. AdoMet binding triggers increased elasticity of the CTD leading Smyd2 to adopt fully opened configurations, which completely exposes the substrate binding crevice. These long-range concerted motions highlight Smyd2's ability to target substrates of varying sizes. Substrate crevices of Smyd2 and Smyd3 show distinct features in terms of spatial, hydration, and electrostatic properties that emphasize their characteristic modes of substrates interaction and entry pathways for inhibitor binding. On the whole, our study shows how the elasticity and hinge motion of the CTD regulate its functional role and underpin the basis of broad substrate specificity of Smyd2. We also highlight the specific structural principles that guide substrate and inhibitor binding to Smyd2 and Smyd3.

Nucleosome Binding Alters the Substrate Bonding Environment of Histone H3 Lysine 36 Methyltransferase NSD2

Nuclear receptor-binding SET domain protein 2 (NSD2) is a histone H3 lysine 36 (H3K36)-specific methyltransferase enzyme that is overexpressed in a number of cancers, including multiple myeloma. NSD2 binds to S-adenosyl-l-methionine (SAM) and nucleosome substrates to catalyze the transfer of a methyl group from SAM to the ε-amino group of histone H3K36. Equilibrium binding isotope effects and density functional theory calculations indicate that the SAM methyl group is sterically constrained in complex with NSD2, and that this steric constraint is released upon nucleosome binding. Together, these results show that nucleosome binding to NSD2 induces a significant change in the chemical environment of enzyme-bound SAM.

Smyd3 open & closed lock mechanism for substrate recruitment: The hinge motion of C-terminal domain inferred from μ-second molecular dynamics simulations

BACKGROUND

The human lysine methyltransferase Smyd3, a member of the SET and MYND domain containing protein family, harbors methylation activity on both histone and non-histone targets in a tightly regulated manner. The mechanism of how Smyd3 dynamically regulates substrate recognition is still not fully unveiled.

METHODS

Here, we employed molecular dynamics simulations on full length human Smyd3, performed to a total of 1.2 μ-second, in the presence (holo) and absence (apo) of the S-Adenosyl methionine (AdoMet) cofactor. The dynamical features of Smyd3 in apo and holo states have been examined and compared via examining geometrical and electrostatic properties.

RESULTS

The results show a distinct dynamics of the C-terminal domain (CTD) in the two states. In the apo state, the CTD undergoes a large hinge like motion and samples more opened configurations, thus acting like a loosened clamp and resulting in expanded substrate binding crevice. In the holo state, the CTD exhibits a restricted motion while the overall structure remains compact, mimicking a closed clamp. This leads to a localized increase in the negative potential at the substrate binding cleft. Further, solvent accessibility of critical residues at the target lysine access channel, important for methylation activity, is increased.

CONCLUSIONS

We postulate that AdoMet cofactor acts like a key and locks Smyd3 in a closed conformation. In effect, the cofactor binding restricts the elasticity of the CTD, presenting a compact substrate binding cleft with high negative potential, which may have implications on substrate recruitment via long range electrostatics.

GENERAL SIGNIFICANCE

The deletion of the CTD from Smyd3 has been shown to abolish the basal histone methylation activity. Our study highlights the importance of the CTD elasticity in shaping the substrate binding site for recognition and supports the previously proposed role of the CTD in stabilizing the active site for methylation activity.

Label-free colorimetric detection of biothiols utilizing SAM and unmodified Au nanoparticles

Herein, a sensitive and selective sensor for biothiols based on colorimetric assay is reported. S-adenosyl-L-methionine (SAM) could induce the selective aggregation of unmodified gold nanoparticles (AuNPs) by electrostatic interaction. In the presence of biothiols, such as glutathione (GSH), homocysteine (Hcy), and cysteine (Cys), AuNPs prefer to react with thiols of biothiols rather than SAM due to the formation of Au-S bond. Thus, the AuNPs turn from the aggregation to the dispersion state, and the corresponding color variation in the process of anti-aggregation of AuNPs can be used for the quantitative screening of biothiols through UV-vis spectroscopy or by the naked eye. Under optimized conditions, a good linear relationship in the range of 0.4-1.2 µM is obtained for Cys, 0.2-0.9 µM for GSH, and 0.6-3.0 µM for Hcys. The detection limits of this assay for GSH, Cys and Hcys are 35.8 nM, 21.7 nM, and 62.4 nM, respectively. This colorimetric assay exhibits rapid operation (within 5 min), high selectivity and sensitivity towards biothiols with tunable dynamic ranges.

Estimating the binding ability of onium ions with CO2 and π systems: a computational investigation

The impact of increasing methyl substitution on onium ions in their complexes with CO₂ and aromatic systems has been analyzed using DFT calculations.

DNA cytosine methylation: structural and thermodynamic characterization of the epigenetic marking mechanism

DNA cytosine methyltransferases regulate the expression of the genome through the precise epigenetic marking of certain cytosines with a methyl group, and aberrant methylation is a hallmark of human diseases including cancer. Targeting these enzymes for drug design is currently a high priority. We have utilized ab initio quantum mechanical/molecular mechanical (QM/MM) molecular dynamics (MD) simulations to investigate extensively the reaction mechanism of the representative DNA methyltransferase HhaI (M.HhaI) from prokaryotes, whose overall mechanism is shared with the mammalian enzymes. We obtain for the first time full free energy profiles for the complete reaction, together with reaction dynamics in atomistic detail. Our results show an energetically preferred mechanism in which nucleophilic attack of cytosine C5 on the S-adenosyl-L-methionine (AdoMet) methyl group is concerted with formation of the Michael adduct between a conserved Cys in the active site with cytosine C6. Spontaneous and reversible proton transfer between a conserved Glu in the active site and cytosine N3 at the transition state was observed in our simulations, revealing the chemical participation of this Glu residue in the catalytic mechanism. Subsequently, the β-elimination of the C5 proton utilizes as base an OH(-) derived from a conserved crystal water that is part of a proton wire water channel, and this syn β-elimination reaction is the rate-limiting step. Design of novel cytosine methylation inhibitors would be advanced by our structural and thermodynamic characterization of the reaction mechanism.

Direct evidence for methyl group coordination by carbon-oxygen hydrogen bonds in the lysine methyltransferase SET7/9

SET domain lysine methyltransferases (KMTs) are S-adenosylmethionine (AdoMet)-dependent enzymes that catalyze the site-specific methylation of lysyl residues in histone and non-histone proteins. Based on crystallographic and cofactor binding studies, carbon-oxygen (CH · · · O) hydrogen bonds have been proposed to coordinate the methyl groups of AdoMet and methyllysine within the SET domain active site. However, the presence of these hydrogen bonds has only been inferred due to the uncertainty of hydrogen atom positions in x-ray crystal structures. To experimentally resolve the positions of the methyl hydrogen atoms, we used NMR (1)H chemical shift coupled with quantum mechanics calculations to examine the interactions of the AdoMet methyl group in the active site of the human KMT SET7/9. Our results indicated that at least two of the three hydrogens in the AdoMet methyl group engage in CH · · · O hydrogen bonding. These findings represent direct, quantitative evidence of CH · · · O hydrogen bond formation in the SET domain active site and suggest a role for these interactions in catalysis. Furthermore, thermodynamic analysis of AdoMet binding indicated that these interactions are important for cofactor binding across SET domain enzymes.

Crystal structures of NodS N-methyltransferase from Bradyrhizobium japonicum in ligand-free form and as SAH complex

NodS is an S-adenosyl-L-methionine (SAM)-dependent N-methyltransferase that is involved in the biosynthesis of Nod factor (NF) in rhizobia, which are bacterial symbionts of legume plants. NF is a modified chitooligosaccharide (COS) signal molecule that is recognized by the legume host, where it initiates symbiotic processes leading to atmospheric nitrogen fixation. We report the crystal structure of recombinant NodS protein from Bradyrhizobium japonicum, which infects lupine and serradella legumes. Two crystal forms--ligand-free NodS and NodS in complex with S-adenosyl-L-homocysteine, which is a by-product of the methylation reaction--were obtained, and their structures were refined to resolutions of 2.43 Å and 1.85 Å, respectively. Although the overall fold (consisting of a seven-stranded β-sheet flanked by layers of helices) is similar to those of other SAM-dependent methyltransferases, NodS has specific features reflecting the unique character of its oligosaccharide substrate. In particular, the N-terminal helix and its connecting loop get ordered upon SAM binding, thereby closing the methyl donor cavity and shaping a long surface canyon that is clearly the binding site for the acceptor molecule. Comparison of the two structural forms of NodS suggests that there are also other conformational changes taking place upon the binding of the donor substrate. As an enzyme that methylates a COS substrate, NodS is the first example among all SAM-dependent methyltransferases to have its three-dimensional structure elucidated. Gaining insight about how NodS binds its donor and acceptor substrates helps to better understand the mechanism of NodS activity and the basis of its functional difference in various rhizobia.

The major product ion of S-adenosyl-L-methionine arises from a neighbouring group reaction

Previous studies have shown that low-energy collision-induced dissociation (CID) of the important sulfonium ion metabolite S-adenosyl-L-methionine (AdoMet, m/z 399) yields five main product ions: an ion at m/z 250 arising from methionine loss; ions at m/z 102 and 298, which arise via cleavage of the gamma C-S bond of methionine; and ions at m/z 136 and 264, which arise via loss of protonated and neutral adenine, respectively. These metabolomics studies have, however, either totally ignored the mechanisms that govern the formation of the major product ion at m/z 250 (Gellekink H, van Oppenraaij-Emmerzaal D, van Rooij A, Struys EA, den Heijer M, Blom HJ. Clin. Chem. 2005; 51: 1487), or have proposed an oxonium ion structure that must arise via a rearrangement involving a 1,2 hydride shift (Cataldi TRI, Bianco G, Abate S, Mattia D. Rapid Commun. Mass Spectrom. 2009; 23: 3465). Here DFT calculations on a model system are used to examine potential mechanisms for the formation of the major product ion of AdoMet. These calculations suggest that a neighbouring group mechanism is preferred over a 1,2 hydride shift mechanism.

Structural biology of S-adenosylmethionine decarboxylase

S-adenosylmethionine decarboxylase (AdoMetDC) is a critical enzyme in the polyamine biosynthetic pathway and a subject of many structural and biochemical investigations for anti-cancer and anti-parasitic therapy. The enzyme undergoes an internal serinolysis reaction as a post-translational modification to generate the active site pyruvoyl group for the decarboxylation process. The crystal structures of AdoMetDC from Homo sapiens, Solanum tuberosum, Thermotoga maritima, and Aquifex aeolicus have been determined. Numerous crystal structures of human AdoMetDC and mutants have provided insights into the mechanism of autoprocessing, putrescine activation, substrate specificity, and inhibitor design to the enzyme. The comparison of the human and potato enzyme with the T. maritima and A. aeolicus enzymes supports the hypothesis that the eukaryotic enzymes evolved by gene duplication and fusion. The residues implicated in processing and activity are structurally conserved in all forms of the enzyme, suggesting a divergent evolution of AdoMetDC.

An investigation of the catalytic mechanism of S-adenosylmethionine synthetase by QM/MM calculations

Catalysis by S-adenosylmethionine synthetase has been investigated by quantum mechanical/molecular mechanical calculations, exploiting structures of the active crystalline enzyme. The transition state energy of +19.1 kcal/mol computed for a nucleophilic attack of the methionyl sulfur on carbon-5' of the nucleotide was indistinguishable from the experimental (solution) value when the QM residues were an uncharged histidine that hydrogen bonds to the leaving oxygen-5' and an aspartate that chelates a Mg2+ ion, and was similar (+18.8 kcal/mol) when the QM region also included the active site arginine and lysines. The computed energy difference between reactant and product was also consistent with their equimolar abundance in co-crystals. The calculated geometrical changes support catalysis of a S(N)2 reaction through hydrogen bonding of the liberated oxygen-5' to the histidine, charge neutralization by the two Mg2+ ions, and stabilization of the product sulfonium cation through a close, non-bonded, contact between the sulfur and the ribose oxygen-4'.

Role of the sulfonium center in determining the ligand specificity of human s-adenosylmethionine decarboxylase

S-Adenosylmethionine decarboxylase (AdoMetDC) is a key enzyme in the polyamine biosynthetic pathway. Inhibition of this pathway and subsequent depletion of polyamine levels is a viable strategy for cancer chemotherapy and for the treatment of parasitic diseases. Substrate analogue inhibitors display an absolute requirement for a positive charge at the position equivalent to the sulfonium of S-adenosylmethionine. We investigated the ligand specificity of AdoMetDC through crystallography, quantum chemical calculations, and stopped-flow experiments. We determined crystal structures of the enzyme cocrystallized with 5'-deoxy-5'-dimethylthioadenosine and 5'-deoxy-5'-(N-dimethyl)amino-8-methyladenosine. The crystal structures revealed a favorable cation-pi interaction between the ligand and the aromatic side chains of Phe7 and Phe223. The estimated stabilization from this interaction is 4.5 kcal/mol as determined by quantum chemical calculations. Stopped-flow kinetic experiments showed that the rate of the substrate binding to the enzyme greatly depends on Phe7 and Phe223, thus supporting the importance of the cation-pi interaction.

New insights into the design of inhibitors of human S-adenosylmethionine decarboxylase: studies of adenine C8 substitution in structural analogues of S-adenosylmethionine

S-adenosylmethionine decarboxylase (AdoMetDC) is a critical enzyme in the polyamine biosynthetic pathway and depends on a pyruvoyl group for the decarboxylation process. The crystal structures of the enzyme with various inhibitors at the active site have shown that the adenine base of the ligands adopts an unusual syn conformation when bound to the enzyme. To determine whether compounds that favor the syn conformation in solution would be more potent AdoMetDC inhibitors, several series of AdoMet substrate analogues with a variety of substituents at the 8-position of adenine were synthesized and analyzed for their ability to inhibit hAdoMetDC. The biochemical analysis indicated that an 8-methyl substituent resulted in more potent inhibitors, yet most other 8-substitutions provided no benefit over the parent compound. To understand these results, we used computational modeling and X-ray crystallography to study C(8)-substituted adenine analogues bound in the active site.

Structure of Staphylococcus aureus 5'-methylthioadenosine/S-adenosylhomocysteine nucleosidase

5'-Methylthioadenosine/S-adenosylhomocysteine nucleosidase (MTAN) catalyzes the irreversible cleavage of the glycosidic bond in 5'-methylthioadenosine (MTA) and S-adenosylhomocysteine (SAH) and plays a key role in four metabolic processes: biological methylation, polyamine biosynthesis, methionine recycling and bacterial quorum sensing. The absence of the nucleosidase in mammalian species has implicated this enzyme as a target for antimicrobial drug design. MTAN from the pathogenic bacterium Staphylococcus aureus (SaMTAN) has been kinetically characterized and its structure has been determined in complex with the transition-state analogue formycin A (FMA) at 1.7 A resolution. A comparison of the SaMTAN-FMA complex with available Escherichia coli MTAN structures shows strong conservation of the overall structure and in particular of the active site. The presence of an extra water molecule, which forms a hydrogen bond to the O4' atom of formycin A in the active site of SaMTAN, produces electron withdrawal from the ribosyl group and may explain the lower catalytic efficiency that SaMTAN exhibits when metabolizing MTA and SAH relative to the E. coli enzyme. The implications of this structure for broad-based antibiotic design are discussed.

How do SET-domain protein lysine methyltransferases achieve the methylation state specificity? Revisited by Ab initio QM/MM molecular dynamics simulations

A distinct protein lysine methyltransferase (PKMT) only transfers a certain number of methyl group(s) to its target lysine residue in spite of the fact that a lysine residue can be either mono-, di-, or tri-methylated. In order to elucidate how such a remarkable product specificity is achieved, we have carried out ab initio quantum mechanical/molecular mechanical (QM/MM) molecular dynamics simulations on two SET-domain PKMTs: SET7/9 and Rubisco large subunit methyltransferase (LSMT). The results indicate that the methylation state specificity is mainly controlled by the methyl-transfer reaction step, and confirm that SET7/9 is a mono-methyltransferase while LSMT has both mono-and di-methylation activities. It is found that the binding of the methylated lysine substrate in the active site of SET7/ 9 opens up the cofactor AdoMet binding channel so that solvent water molecules get access to the active site. This disrupts the catalytic machinery of SET7/9 for the di-methylation reaction, which leads to a higher activation barrier, whereas for the LSMT, its active site is more spacious than that of SET7/9, so that the methylated lysine substrate can be accommodated without interfering with its catalytic power. These detailed insights take account of protein dynamics and are consistent with available experimental results as well as recent theoretical findings regarding the catalytic power of SET7/9.

Reactivity and reaction order in acylhomoserine lactone formation by Pseudomonas aeruginosa RhlI

The formation of N-butyrylhomoserine lactone catalyzed by RhlI has been investigated by transient-state kinetic methods. A single intermediate, assigned to N-butyryl- S-adenosylmethionine, was observed. Under single-turnover conditions, the intermediate formed with a rate constant of 4.0 +/- 0.2 s (-1) and decayed with a rate constant of 3.7 +/- 0.2 s (-1). No other intermediates were detected, demonstrating that the RhlI reaction proceeds via acylation of S-adenosylmethionine, followed by lactonization. S-Adenosylhomocysteine acted as a pseudosubstrate, in that it did not undergo either acylation or lactonization but did induce the deacylation of butyryl-acyl carrier protein. The K m for S-adenosylhomocysteine was approximately 15-fold higher than the K m for S-adenosylmethionine. The reactivities of acylated and unacylated sulfonium ions that were analogues of S-adenosylmethionine were investigated by computational methods. The calculations indicated that acylation of the substrate amino group activated the substrate for lactonization, by allowing the carboxyl group oxygen to approach more closely the methylene carbon to which it adds. This observation provides a satisfying chemical rationale for the order of the individual reactions in the catalytic cycle.

Ab initio quantum mechanical/molecular mechanical molecular dynamics simulation of enzyme catalysis: the case of histone lysine methyltransferase SET7/9

To elucidate enzyme catalysis through computer simulation, a prerequisite is to reliably compute free energy barriers for both enzyme and solution reactions. By employing on-the-fly Born-Oppenheimer molecular dynamics simulations with the ab initio quantum mechanical/molecular mechanical approach and the umbrella sampling method, we have determined free energy profiles for the methyl-transfer reaction catalyzed by the histone lysine methyltransferase SET7/9 and its corresponding uncatalyzed reaction in aqueous solution, respectively. Our calculated activation free energy barrier for the enzyme catalyzed reaction is 22.5 kcal/mol, which agrees very well with the experimental value of 20.9 kcal/mol. The difference in potential of mean force between a corresponding prereaction state and the transition state for the solution reaction is computed to be 30.9 kcal/mol. Thus, our simulations indicate that the enzyme SET7/9 plays an essential catalytic role in significantly lowering the barrier for the methyl-transfer reaction step. For the reaction in solution, it is found that the hydrogen bond network near the reaction center undergoes a significant change, and there is a strong shift in electrostatic field from the prereaction state to the transition state, whereas for the enzyme reaction, such an effect is much smaller and the enzyme SET7/9 is found to provide a preorganized electrostatic environment to facilitate the methyl-transfer reaction. Meanwhile, we find that the transition state in the enzyme reaction is a little more dissociative than that in solution.

Catalytic mechanism and product specificity of the histone lysine methyltransferase SET7/9: an ab initio QM/MM-FE study with multiple initial structures

Histone lysine methylation is emerging as an important mechanism to regulate chromatin structure and gene activity. To provide theoretical understanding of its reaction mechanism and product specificity, ab initio quantum mechanical/molecular mechanical free energy (QM/MM-FE) calculations and molecular dynamics simulations have been carried out to investigate the histone lysine methyltransferase SET7/9. It is found that the methyl-transfer reaction catalyzed by SET7/9 is a typical in-line S(N)2 nucleophilic substitution reaction with a transition state of 70% dissociative character. The calculated average free energy barrier at the MP2(6-31+G) QM/MM level is 20.4 +/- 1.1 kcal/mol, consistent with the activation barrier of 20.9 kcal/mol estimated from the experimental reaction rate. The barrier fluctuation has a strong correlation with the nucleophilic attack distance and angle in the reactant complex. The calculation results show that the product specificity of SET7/9 as a monomethyltransferase is achieved by disrupting the formation of near-attack conformations for the dimethylation reaction.

Insight into enzymatic C-F bond formation from QM and QM/MM calculations

The C-F bond-forming step in the fluorinase, the only native fluorination enzyme characterized to date, has been studied. The enzyme catalyzes the reaction between S-adenosyl-L-methionine (SAM) and fluoride ions to form 5'-fluoro-5'-deoxyadenosine (5'-FDA) and L-methionine. To obtain an insight into the mechanism of this unusual enzymatic reaction and to elucidate the role of the enzyme in catalysis, we have explored the conformational energetics of SAM and the intrinsic reactivity patterns of SAM and fluoride with DFT (BP86) and continuum solvent methods, before investigating the full enzymatic system with combined DFT/CHARMM calculations. We find that the enzymatic reaction follows an S(N)2 reaction mechanism, concurring with the intrinsic reactivity preferences in solution. The formation of sulfur ylides is thermodynamically strongly disfavored, and an alternative elimination-addition mechanism involving the concerted anti-Markovnikov addition of HF to an enol ether is energetically viable, but kinetically prohibitive. The S(N)2 activation energy is 92 (112) kJ mol(-)(1) in solution, but only 53 (63) kJ mol(-1) in the enzyme, and the reaction energy in the enzyme is -25 (-34) kJ mol(-1) (values in parentheses are B3LYP single-point energies). The fluorinase thus lowers the barrier for C-F bond formation by 39 (49) kJ mol(-)(1). A decomposition analysis shows that the major role of the enzyme is in the preparation and positioning of the substrates.

Structural snapshots of MTA/AdoHcy nucleosidase along the reaction coordinate provide insights into enzyme and nucleoside flexibility during catalysis

MTA/AdoHcy nucleosidase (MTAN) irreversibly hydrolyzes the N9-C1' bond in the nucleosides, 5'-methylthioadenosine (MTA) and S-adenosylhomocysteine (AdoHcy) to form adenine and the corresponding thioribose. MTAN plays a vital role in metabolic pathways involving methionine recycling, biological methylation, polyamine biosynthesis, and quorum sensing. Crystal structures of a wild-type (WT) MTAN complexed with glycerol, and mutant-enzyme and mutant-product complexes have been determined at 2.0A, 2.0A, and 2.1A resolution, respectively. The WT MTAN-glycerol structure provides a purine-free model and in combination with the previously solved thioribose-free MTAN-ADE structure, we now have separate apo structures for both MTAN binding subsites. The purine and thioribose-free states reveal an extensive enzyme-immobilized water network in their respective binding subsites. The Asp197Asn MTAN-MTA and Glu12Gln MTAN-MTR.ADE structures are the first enzyme-substrate and enzyme-product complexes reported for MTAN, respectively. These structures provide representative snapshots along the reaction coordinate and allow insight into the conformational changes of the enzyme and the nucleoside substrate. A "catalytic movie" detailing substrate binding, catalysis, and product release is presented.

Characterization of a methionine adenosyltransferase over-expressing strain in the trypanosomatid Leishmania donovani

Methionine adenosyltransferase (MAT: EC 2.5.1.6) catalyzes the synthesis of S-adenosylmethionine (AdoMet) in two sequential steps, AdoMet formation and subsequent tripolyphosphate (PPPi) cleavage, induced by AdoMet. In pursuit of a better understanding of the biological function of the enzyme, the MAT gene was cloned into vector PX63NEO to induce episomal overexpression in leishmania parasites. Neomycin-selected clones originated a strain of such overexpressing parasites that accumulated more than 3-fold AdoMet than mock-transfected cells and showed over ten times the wild type MAT activity, concurring with a significant accumulation of the MAT protein during the early logarithmic phase and MAT transcripts throughout the growth cycle. The rate of AdoMet efflux, practically nil in the control promastigotes, was exceptionally high in the MAT-overexpressing parasites, whilst growth in this strain was comparable to development in control cells, i.e., it was not affected by deleterious hypermethylation. Moreover, the modified strain was 10-fold more resistant to sinefungin, a S-adenosylmethionine-like antibiotic, than control cells. The effects of overexpression on polyamine metabolism and transport were likewise studied.

Structural rationale for the affinity of pico- and femtomolar transition state analogues of Escherichia coli 5'-methylthioadenosine/S-adenosylhomocysteine nucleosidase

Immucillin and DADMe-Immucillin inhibitors are tight binding transition state mimics of purine nucleoside phosphorylases (PNP). 5'-Methylthioadenosine/S-adenosylhomocysteine nucleosidase (MTAN) is proposed to form a similar transition state structure as PNP. The companion paper describes modifications of the Immucillin and DADMe-Immucillin inhibitors to better match transition state features of MTAN and have led to 5'-thio aromatic substitutions that extend the inhibition constants to the femtomolar range (Singh, V., Evans, G. B., Lenz, D. H., Mason, J., Clinch, K., Mee, S., Painter, G. F., Tyler, P. C., Furneaux, R. H., Lee, J. E., Howell, P. L., and Schramm, V. L. (2005) J. Biol. Chem. 280, 18265-18273). 5'-Methylthio-Immucillin A (MT-ImmA) and 5'-methylthio-DADMe-Immucillin A (MT-DADMe-ImmA) exhibit slow-onset inhibition with K(i)(*) of 77 and 2 pm, respectively, and were selected for structural analysis as the parent compounds of each class of transition state analogue. The crystal structures of Escherichia coli MTAN complexed with MT-ImmA and MT-DADMe-ImmA were determined to 2.2 A resolution and compared with the existing MTAN inhibitor complexes. These MTAN-transition state complexes are among the tightest binding enzyme-ligand complexes ever described and analysis of their mode of binding provides extraordinary insight into the structural basis for their affinity. The MTAN-MT-ImmA complex reveals the presence of a new ion pair between the 4'-iminoribitol atom and the nucleophilic water (WAT3) that captures key features of the transition state. Similarly, in the MTAN-MT-DADMe-ImmA complex a favorable hydrogen bond or ion pair interaction between the cationic 1'-pyrrolidine atom and WAT3 is crucial for tight affinity. Distance analysis of the nucleophile and leaving group show that MT-ImmA is a mimic of an early transition state, while MT-DADMe-ImmA is a better mimic of the highly dissociated transition state of E. coli MTAN.

Decarboxylases involved in polyamine biosynthesis and their inactivation by nitric oxide

Polyamines are ubiquitous cellular components that are involved in normal and neoplastic growth. Polyamine biosynthesis is very highly regulated in mammalian cells by the activities of two key decarboxylases acting on ornithine and S-adenosylmethionine. Recent studies, which include crystallographic analysis of the recombinant human proteins, have provided a detailed knowledge of their structure and function. Ornithine decarboxylase is a PLP-requiring decarboxylase, whereas S-adenosylmethionine decarboxylase (AdoMetDC) contains a covalently bound pyruvate prosthetic group. Both enzymes have a key cysteine residue, which is involved in protonation of the Schiff base intermediate C(alpha) to form the product. These residues, Cys360 in ornithine decarboxylase (ODC) and Cys82 in AdoMetDC, react readily with nitric oxide (NO), which is therefore a potent inactivator of polyamine synthesis. The inactivation of these enzymes may mediate some of the antiproliferative actions of NO.

Structure of 1-aminocyclopropane-1-carboxylate synthase in complex with an amino-oxy analogue of the substrate: implications for substrate binding

The crystal structure of 1-aminocyclopropane-1-carboxylate (ACC) synthase in complex with the substrate analogue [2-(amino-oxy)ethyl](5'-deoxyadenosin-5'-yl)(methyl)sulfonium (AMA) was determined at 2.01-A resolution. The crystallographic results show that a covalent adduct (oxime) is formed between AMA (an amino-oxy analogue of the natural substrate S-adenosyl-L-methionine (SAM)) and the pyridoxal 5'-phosphate (PLP) cofactor of ACC synthase. The oxime formation is supported by spectroscopic data. The ACC synthase-AMA structure provides reliable and detailed information on the binding mode of the natural substrate of ACC synthase and complements previous structural and functional work on this enzyme.

Structure of Escherichia coli 5'-methylthioadenosine/ S-adenosylhomocysteine nucleosidase inhibitor complexes provide insight into the conformational changes required for substrate binding and catalysis

5'-Methylthioadenosine/S-adenosylhomocysteine (MTA/AdoHcy) nucleosidase is a key enzyme in a number of critical biological processes in many microbes. This nucleosidase catalyzes the irreversible hydrolysis of the N(9)-C(1') bond of MTA or AdoHcy to form adenine and the corresponding thioribose. The key role of the MTA/AdoHcy nucleosidase in biological methylation, polyamine biosynthesis, methionine recycling, and bacterial quorum sensing has made it an important antimicrobial drug target. The crystal structures of Escherichia coli MTA/AdoHcy nucleosidase complexed with the transition state analog, formycin A (FMA), and the nonhydrolyzable substrate analog, 5'-methylthiotubercidin (MTT) have been solved to 2.2- and 2.0-A resolution, respectively. These are the first MTA/AdoHcy nucleosidase structures to be solved in the presence of inhibitors. These structures clearly identify the residues involved in substrate binding and catalysis in the active site. Comparisons of the inhibitor complexes to the adenine-bound MTA/AdoHcy nucleosidase (Lee, J. E., Cornell, K. A., Riscoe, M. K., and Howell, P. L. (2001) Structure (Camb.) 9, 941-953) structure provide evidence for a ligand-induced conformational change in the active site and the substrate preference of the enzyme. The enzymatic mechanism has been re-examined.