Convergent Mechanistic Features between the Structurally Diverse N- and O-Methyltransferases: Glycine N-Methyltransferase and Catechol O-Methyltransferase

Electrostatic Preorganization in Three Distinct Heterogeneous Proteasome β-Subunits

The origin of the enzyme's powerful role in accelerating chemical reactions is one of the most critical and still widely discussed questions. It is already accepted that enzymes impose an electrostatic field onto their substrates by adopting complex three-dimensional structures; therefore, the preorganization of electric fields inside protein active sites has been proposed as a crucial contributor to catalytic mechanisms and rate constant enhancement. In this work, we focus on three catalytically active β-subunits of 20S proteasomes with low sequence identity (∼30%) whose active sites, although situated in an electrostatically miscellaneous environment, catalyze the same chemical reaction with similar catalytic efficiency. Our in silico experiments reproduce the experimentally observed equivalent reactivity of the three sites and show that obliteration of the electrostatic potential in all active sites would deprive the enzymes of their catalytic power by slowing down the chemical process by a factor of 1035. To regain enzymatic efficiency, besides catalytic Thr1 and Lys33 residues, the presence of aspartic acid in position 17 and an aqueous solvent is required, proving that the electrostatic potential generated by the remaining residues is insignificant for catalysis. Moreover, it was found that the gradual decay of atomic charges on Asp17 strongly correlates with the enzyme's catalytic rate deterioration as well as with a change in the charge distributions due to introduced mutations. The computational procedure used and described here may help identify key residues for catalysis in other biomolecular systems and consequently may contribute to the process of designing enzyme-like synthetic catalysts.

Methyl transfer in psilocybin biosynthesis

Psilocybin, the natural hallucinogen produced by Psilocybe ("magic") mushrooms, holds great promise for the treatment of depression and several other mental health conditions. The final step in the psilocybin biosynthetic pathway, dimethylation of the tryptophan-derived intermediate norbaeocystin, is catalysed by PsiM. Here we present atomic resolution (0.9 Å) crystal structures of PsiM trapped at various stages of its reaction cycle, providing detailed insight into the SAM-dependent methylation mechanism. Structural and phylogenetic analyses suggest that PsiM derives from epitranscriptomic N6-methyladenosine writers of the METTL16 family, which is further supported by the observation that bound substrates physicochemically mimic RNA. Inherent limitations of the ancestral monomethyltransferase scaffold hamper the efficiency of psilocybin assembly and leave PsiM incapable of catalysing trimethylation to aeruginascin. The results of our study will support bioengineering efforts aiming to create novel variants of psilocybin with improved therapeutic properties.

Electrostatics as a Guiding Principle in Understanding and Designing Enzymes

Enzyme design faces challenges related to the implementation of the basic principles that govern the catalytic activity in natural enzymes. In this work, we revisit basic electrostatic concepts that have been shown to explain the origin of enzymatic efficiency like preorganization and reorganization. Using magnitudes such as the electrostatic potential and the electric field generated by the protein, we explain how these concepts work in different enzymes and how they can be used to rationalize the consequences of point mutations. We also discuss examples of protein design in which electrostatic effects have been implemented. For the near future, molecular simulations, coupled with the use of machine learning methods, can be used to implement electrostatics as a guiding principle for enzyme designs.

The Glycine N-Methyltransferase Case Study: Another Challenge for QM-Cluster Models?

The methyl transfer reaction between SAM and glycine catalyzed by glycine N-methyltransferase (GNMT) was examined using QM-cluster models generated by Residue Interaction Network ResidUe Selector (RINRUS). RINRUS is a Python-based tool that can build QM-cluster models with rules-based processing of the active site residue interaction network. This way of enzyme model-building allows quantitative analysis of residue and fragment contributions to kinetic and thermodynamic properties of the enzyme. Many residue fragments are important for the GNMT catalytic reaction, such as Gly137, Asn138, and Arg175, which interact with the glycine substrate, and Trp30, Asp85, and Tyr242, which interact with the SAM cofactor. Our study shows that active site fragments that interact with the glycine substrate and the SAM cofactor must both be included in the QM-cluster models. Even though the proposed mechanism is a simple one-step reaction, GNMT may be a rather challenging case study for QM-cluster models because convergence in energetics requires models with >350 atoms. "Maximal" QM-cluster models built with either qualitative contact count ranking or quantitative interaction energies from functional group symmetry adapted perturbation theory provide acceptable results. Hence, important residue fragments that contribute to the energetics of the methyl-transfer reaction in GNMT are correctly identified in the RIN. Observations from this work suggest new directions to better establish an effective approach for constructing atomic-level enzyme models.

High-throughput computational investigation of protein electrostatics and cavity for SAM-dependent methyltransferases

S-adenosyl methionine (SAM)-dependent methyl transferases (MTases) are a ubiquitous class of enzymes catalyzing dozens of essential life processes. Despite targeting a large space of substrates with diverse intrinsic reactivity, SAM MTases have similar catalytic efficiency. While understanding of MTase mechanism has grown tremendously through the integration of structural characterization, kinetic assays, and multiscale simulations, it remains elusive how these enzymes have evolved to fit the diverse chemical needs of their respective substrates. In this work, we performed a high-throughput molecular modeling analysis of 91 SAM MTases to better understand how their properties (i.e., electric field [EF] strength and active site volumes) help achieve similar catalytic efficiency toward substrates of different reactivity. We found that EF strengths have largely adjusted to make the target atom a better methyl acceptor. For MTases that target RNA/DNA and histone proteins, our results suggest that EF strength accommodates formal hybridization state and variation in cavity volume trends with diversity of substrate classes. Metal ions in SAM MTases contribute negatively to EF strength for methyl donation and enzyme scaffolds tend to offset these contributions.

Target Class Profiling of Small-Molecule Methyltransferases

Target class profiling (TCP) is a chemical biology approach to investigate understudied biological target classes. TCP is achieved by developing a generalizable assay platform and screening curated compound libraries to interrogate the chemical biological space of members of an enzyme family. In this work, we took a TCP approach to investigate inhibitory activity across a set of small-molecule methyltransferases (SMMTases), a subclass of methyltransferase enzymes, with the goal of creating a launchpad to explore this largely understudied target class. Using the representative enzymes nicotinamide N-methyltransferase (NNMT), phenylethanolamine N-methyltransferase (PNMT), histamine N-methyltransferase (HNMT), glycine N-methyltransferase (GNMT), catechol O-methyltransferase (COMT), and guanidinoacetate N-methyltransferase (GAMT), we optimized high-throughput screening (HTS)-amenable assays to screen 27,574 unique small molecules against all targets. From this data set, we identified a novel inhibitor which selectively inhibits the SMMTase HNMT and demonstrated how this platform approach can be leveraged for a targeted drug discovery campaign using the example of HNMT.

Biocatalytic C3-Indole Methylation-A Useful Tool for the Natural-Product-Inspired Stereoselective Synthesis of Pyrroloindoles

Enantioselective synthesis of bioactive compounds bearing a pyrroloindole framework is often laborious. In contrast, there are several S-adenosyl methionine (SAM)-dependent methyl transferases known for stereo- and regioselective methylation at the C3 position of various indoles, directly leading to the formation of the desired pyrroloindole moiety. Herein, the SAM-dependent methyl transferase PsmD from Streptomyces griseofuscus, a key enzyme in the biosynthesis of physostigmine, is characterized in detail. The biochemical properties of PsmD and its substrate scope were demonstrated. Preparative scale enzymatic methylation including SAM regeneration was achieved for three selected substrates after a design-of-experiment optimization.

Structural diversity, bioactivities, and biosynthesis of natural diterpenoid alkaloids

Covering: 2009 to 2018. Diterpenoid alkaloids, originating from the amination of natural tetracyclic diterpenes, are a diverse class of compounds having complex structural features with many stereocenters. The important pharmacological activities and structural complexity of the diterpenoid alkaloids have long interested scientists due to their medicinal uses, infamous toxicity, and unique biosynthesis. Since 2009, 373 diterpenoid alkaloids, assigned to 46 skeletons, have been isolated and identified from plants mostly in the Ranunculaceae family. The names, classes, molecular weight, molecular formula, NMR data, and plant sources of these diterpene alkaloids are collated here. This review will be a detailed update of the naturally occurring diterpene alkaloids reported from the plant kingdom from 2009-2018, providing an in-depth discussion of their diversity, biological activities, pharmacokinetics, toxicity, application, evolution, and biosynthesis.

Mechanistic Studies on CysS - A Vitamin B12-Dependent Radical SAM Methyltransferase Involved in the Biosynthesis of the tert-Butyl Group of Cystobactamid

Cobalamin (Cbl)-dependent radical S-adenosylmethionine (SAM) methyltransferases catalyze methylation reactions at non-nucleophilic centers in a wide range of substrates. CysS is a Cbl-dependent radical SAM methyltransferase involved in cystobactamid biosynthesis. This enzyme catalyzes the sequential methylation of a methoxy group to form ethoxy, i-propoxy, s-butoxy, and t-butoxy groups on a p-aminobenzoate peptidyl carrier protein thioester intermediate. This biosynthetic strategy enables the host myxobacterium to biosynthesize a combinatorial antibiotic library of 25 cystobactamid analogues. In this Article, we describe three experiments to elucidate how CysS uses Cbl, SAM, and a [4Fe-4S] cluster to catalyze iterative methylation reactions: a cyclopropylcarbinyl rearrangement was used to trap the substrate radical and to estimate the rate of the radical substitution reaction involved in the methyl transfer; a bromoethoxy analogue was used to explore the active site topography; and deuterium isotope effects on the hydrogen atom abstraction by the adenosyl radical were used to investigate the kinetic significance of the hydrogen atom abstraction. On the basis of these experiments, a revised mechanism for CysS is proposed.

Hydrogen deuterium exchange defines catalytically linked regions of protein flexibility in the catechol O-methyltransferase reaction

Human catechol O-methyltransferase (COMT) has emerged as a model for understanding enzyme-catalyzed methyl transfer from S-adenosylmethionine (AdoMet) to small-molecule catecholate acceptors. Mutation of a single residue (tyrosine 68) behind the methyl-bearing sulfonium of AdoMet was previously shown to impair COMT activity by interfering with methyl donor-acceptor compaction within the activated ground state of the wild type enzyme [J. Zhang, H. J. Kulik, T. J. Martinez, J. P. Klinman, Proc. Natl. Acad. Sci. U.S.A. 112, 7954-7959 (2015)]. This predicts the involvement of spatially defined protein dynamical effects that further tune the donor/acceptor distance and geometry as well as the electrostatics of the reactants. Here, we present a hydrogen/deuterium exchange (HDX)-mass spectrometric study of wild type and mutant COMT, comparing temperature dependences of HDX against corresponding kinetic and cofactor binding parameters. The data show that the impaired Tyr68Ala mutant displays similar breaks in Arrhenius plots of both kinetic and HDX properties that are absent in the wild type enzyme. The spatial resolution of HDX below a break point of 15-20 °C indicates changes in flexibility across ∼40% of the protein structure that is confined primarily to the periphery of the AdoMet binding site. Above 20 °C, Tyr68Ala behaves more like WT in HDX, but its rate and enthalpic barrier remain significantly altered. The impairment of catalysis by Tyr68Ala can be understood in the context of a mutationally induced alteration in protein motions that becomes manifest along and perpendicular to the primary group transfer coordinate.

Moving Through Barriers in Science and Life

This first serious attempt at an autobiographical accounting has forced me to sit still long enough to compile my thoughts about a long personal and scientific journey. I especially hope that my trajectory will be of interest and perhaps beneficial to much younger women who are just getting started in their careers. To paraphrase from Virginia Woolf's writings in A Room of One's Own at the beginning of the 20th century, "for most of history Anonymous was a Woman." However, Ms. Woolf is also quoted as saying "nothing has really happened until it has been described," a harbinger of the enormous historical changes that were about to be enacted and recorded by women in the sciences and other disciplines. The progress in my chosen field of study-the chemical basis of enzyme action-has also been remarkable, from the first description of an enzyme's 3D structure to a growing and deep understanding of the origins of enzyme catalysis.

A single residue determines substrate preference in benzylisoquinoline alkaloid N-methyltransferases

N-methylation is a recurring feature in the biosynthesis of many plant specialized metabolites, including alkaloids. A crucial step in the conserved central pathway that provides intermediates for the biosynthesis of benzylisoquinoline alkaloids (BIAs) involves conversion of the secondary amine (S)-coclaurine into the tertiary amine (S)-N-methylcoclaurine by coclaurine N-methyltransferase (CNMT). Subsequent enzymatic steps yield the core intermediate (S)-reticuline, from which various branch pathways for the biosynthesis of major BIAs such as morphine, noscapine and sanguinarine diverge. An additional N-methylation yielding quaternary BIAs is catalyzed by reticuline N-methyltransferase (RNMT), such as in the branch pathway leading to the taxonomically widespread and ecologically significant alkaloid magnoflorine. Despite their functional differences, analysis of primary sequence information has been unable to accurately distinguish between CNMT-like and RNMT-like enzymes, necessitating laborious in vitro screening. Furthermore, despite a recent emphasis on structural characterization of BIA NMTs, the features and mechanisms underlying the CNMT-RNMT functional dichotomy were unknown. We report the identification of structural variants tightly correlated with function in known BIA NMTs and show through reciprocal mutagenesis that a single residue acts as a switch between CNMT- and RNMT-like functions. We use yeast in vivo screening to show that this discovery allows for accurate prediction of activity strictly from primary sequence information and, on this basis, improve the annotation of previously reported putative BIA NMTs. Our results highlight the unusually short mutational distance separating ancestral CNMT-like enzymes from more evolutionarily advanced RNMT-like enzymes, and thus help explain the widespread yet sporadic occurrence of quaternary BIAs in plants. While this is the first report of structural variants controlling mono-versus di-methylation activity among plant NMT enzymes, comparison with bacterial MT enzymes also suggests possible convergent evolution.

Replication of the Enzymatic Temperature Dependency of the Primary Hydride Kinetic Isotope Effects in Solution: Caused by the Protein-Controlled Rigidity of the Donor-Acceptor Centers?

The change from the temperature independence of the primary (1°) H/D kinetic isotope effects (KIEs) in wild-type enzyme-catalyzed H-transfer reactions (ΔE_a = E_aD - E_aH ∼ 0) to a strong temperature dependence with the mutated enzymes (ΔE_a ≫ 0) has recently been frequently observed. This has prompted some enzymologists to develop new H-tunneling models to correlate ΔE_a with the donor-acceptor distance (DAD) at the tunneling-ready state (TRS) as well as the protein thermal motions/dynamics that sample the short DAD_TRS's for H-tunneling to occur. While extensive evidence supporting or disproving the thermally activated DAD sampling concept has emerged, a comparable study of the simpler bimolecular H-tunneling reactions in solution has not been carried out. In particular, small ΔE_a's (∼0) have not been found. In this paper, we report a study of the hydride-transfer reactions from four NADH models to the same hydride acceptor in acetonitrile. The ΔE_a's were determined: 0.37 (small), 0.60, 0.99, and 1.53 kcal/mol (large). The α-secondary (2°) KIEs on the acceptor that serve as a ruler for the rigidity of reaction centers were previously reported or determined. All possible productive reactant complex (PRC) configurations were computed to provide insight into the structures of the TRS's. Relationships among structures, 2° KIEs, DAD_PRC's, and ΔE_a's were discussed. The more rigid system with more suppressed 2° C-H vibrations at the TRS and more narrowly distributed DAD_PRC's in PRCs gave a smaller ΔE_a. The results replicated the trend observed in enzymes versus mutated enzymes and appeared to support the concepts of different thermally activated DAD_TRS sampling processes in response to the rigid versus flexible donor-acceptor centers.

Structure-function studies of tetrahydroprotoberberine N-methyltransferase reveal the molecular basis of stereoselective substrate recognition

Benzylisoquinoline alkaloids (BIAs) are a structurally diverse class of plant-specialized metabolites that have been particularly well-studied in the order Ranunculales. The N-methyltransferases (NMTs) in BIA biosynthesis can be divided into three groups according to substrate specificity and amino acid sequence. Here, we report the first crystal structures of enzyme complexes from the tetrahydroprotoberberine NMT (TNMT) subclass, specifically for GfTNMT from the yellow horned poppy (Glaucium flavum). GfTNMT was co-crystallized with the cofactor S-adenosyl-l-methionine (d _min = 1.6 Å), the product S-adenosyl-l-homocysteine (d _min = 1.8 Å), or in complex with S-adenosyl-l-homocysteine and (S)-cis-N-methylstylopine (d _min = 1.8 Å). These structures reveal for the first time how a mostly hydrophobic L-shaped substrate recognition pocket selects for the (S)-cis configuration of the two central six-membered rings in protoberberine BIA compounds. Mutagenesis studies confirm and functionally define the roles of several highly-conserved residues within and near the GfTNMT-active site. The substrate specificity of TNMT enzymes appears to arise from the arrangement of subgroup-specific stereospecific recognition elements relative to catalytic elements that are more widely-conserved among all BIA NMTs. The binding mode of protoberberine compounds to GfTNMT appears to be similar to coclaurine NMT, with the isoquinoline rings buried deepest in the binding pocket. This binding mode differs from that of pavine NMT, in which the benzyl ring is bound more deeply than the isoquinoline rings. The insights into substrate recognition and catalysis provided here form a sound basis for the rational engineering of NMT enzymes for chemoenzymatic synthesis and metabolic engineering.

The Structure of the Bifunctional Everninomicin Biosynthetic Enzyme EvdMO1 Suggests Independent Activity of the Fused Methyltransferase-Oxidase Domains

Members of the orthosomycin family of natural products are decorated polysaccharides with potent antibiotic activity and complex biosynthetic pathways. The defining feature of the orthosomycins is an orthoester linkage between carbohydrate moieties that is necessary for antibiotic activity and is likely formed by a family of conserved oxygenases. Everninomicins are octasaccharide orthosomycins produced by Micromonospora carbonacea that have two orthoester linkages and a methylenedioxy bridge, three features whose formation logically requires oxidative chemistry. Correspondingly, the evd gene cluster encoding everninomicin D encodes two monofunctional nonheme iron, α-ketoglutarate-dependent oxygenases and one bifunctional enzyme with an N-terminal methyltransferase domain and a C-terminal oxygenase domain. To investigate whether the activities of these domains are linked in the bifunctional enzyme EvdMO1, we determined the structure of the N-terminal methyltransferase domain to 1.1 Å and that of the full-length protein to 3.35 Å resolution. Both domains of EvdMO1 adopt the canonical folds of their respective superfamilies and are connected by a short linker. Each domain's active site is oriented such that it faces away from the other domain, and there is no evidence of a channel connecting the two. Our results support EvdMO1 working as a bifunctional enzyme with independent catalytic activities.

Probing the transcriptome of Aconitum carmichaelii reveals the candidate genes associated with the biosynthesis of the toxic aconitine-type C19-diterpenoid alkaloids

Aconitum carmichaelii has long been used as a traditional Chinese medicine, and its processed lateral roots are known commonly as fuzi. Aconitine-type C₁₉-diterpenoid alkaloids accumulating in the lateral roots are some of the main toxicants of this species, yet their biosynthesis remains largely unresolved. As a first step towards understanding the biosynthesis of aconitine-type C₁₉-diterpenoid alkaloids, we performed de novo transcriptome assembly and analysis of rootstocks and leaf tissues of Aconitum carmichaelii by next-generation sequencing. A total of 525 unigene candidates were identified as involved in the formation of C₁₉-diterpenoid alkaloids, including those encoding enzymes in the early steps of diterpenoid alkaloids scaffold biosynthetic pathway, such as ent-copalyl diphosphate synthases, ent-kaurene synthases, kaurene oxidases, cyclases, and key aminotransferases. Furthermore, candidates responsible for decorating of diterpenoid alkaloid skeletons were discovered from transcriptome sequencing of fuzi, such as monooxygenases, methyltransferase, and BAHD acyltransferases. In addition, 645 differentially expressed genes encoding transcription factors potentially related to diterpenoid alkaloids accumulation underground were documented. Subsequent modular domain structure phylogenetics and differential expression analysis led to the identification of BAHD acyltransferases possibly involved in the formation of acetyl and benzoyl esters of diterpenoid alkaloids, associated with the acute toxicity of fuzi. The transcriptome data provide the foundation for future research into the molecular basis for aconitine-type C₁₉-diterpenoid alkaloids biosynthesis in A. carmichaelii.

Chemical and Biochemical Perspectives of Protein Lysine Methylation

Protein lysine methylation is a distinct posttranslational modification that causes minimal changes in the size and electrostatic status of lysine residues. Lysine methylation plays essential roles in regulating fates and functions of target proteins in an epigenetic manner. As a result, substrates and degrees (free versus mono/di/tri) of protein lysine methylation are orchestrated within cells by balanced activities of protein lysine methyltransferases (PKMTs) and demethylases (KDMs). Their dysregulation is often associated with neurological disorders, developmental abnormalities, or cancer. Methyllysine-containing proteins can be recognized by downstream effector proteins, which contain methyllysine reader domains, to relay their biological functions. While numerous efforts have been made to annotate biological roles of protein lysine methylation, limited work has been done to uncover mechanisms associated with this modification at a molecular or atomic level. Given distinct biophysical and biochemical properties of methyllysine, this review will focus on chemical and biochemical aspects in addition, recognition, and removal of this posttranslational mark. Chemical and biophysical methods to profile PKMT substrates will be discussed along with classification of PKMT inhibitors for accurate perturbation of methyltransferase activities. Semisynthesis of methyllysine-containing proteins will also be covered given the critical need for these reagents to unambiguously define functional roles of protein lysine methylation.

Insights on the Origin of Catalysis on Glycine N-Methyltransferase from Computational Modeling

The origin of enzyme catalysis remains a question of debate despite much intense study. We report a QM/MM theoretical study of the SN2 methyl transfer reaction catalyzed by a glycine N-methyltransferase (GNMT) and three mutants to test whether recent experimental observations of rate-constant reductions and variations in inverse secondary α-3H kinetic isotope effects (KIEs) should be attributed to changes in the methyl donor-acceptor distance (DAD): Is catalysis due to a compression effect? Semiempirical (AM1) and DFT (M06-2X) methods were used to describe the QM subset of atoms, while OPLS-AA and TIP3P classical force fields were used for the protein and water molecules, respectively. The computed activation free energies and KIEs are in good agreement with experimental data, but the mutations do not meaningfully affect the DAD: Compression cannot explain the experimental variations on KIEs. On the contrary, electrostatic properties in the active site correlate with the catalytic activity of wild type and mutants. The plasticity of the enzyme moderates the effects of the mutations, explaining the rather small degree of variation in KIEs and reactivities.

Origins of Enzyme Catalysis: Experimental Findings for C-H Activation, New Models, and Their Relevance to Prevailing Theoretical Constructs

The physical basis for enzymatic rate accelerations is a subject of great fundamental interest and of direct relevance to areas that include the de novo design of green catalysts and the pursuit of new drug regimens. Extensive investigations of C-H activating systems have provided considerable insight into the relationship between an enzyme's overall structure and the catalytic chemistry at its active site. This Perspective highlights recent experimental data for two members of distinct, yet iconic C-H activation enzyme classes, lipoxygenases and prokaryotic alcohol dehydrogenases. The data necessitate a reformulation of the dominant textbook definition of biological catalysis. A multidimensional model emerges that incorporates a range of protein motions that can be parsed into a combination of global stochastic conformational thermal fluctuations and local donor-acceptor distance sampling. These motions are needed to achieve a high degree of precision with regard to internuclear distances, geometries, and charges within the active site. The available model also suggests a physical framework for understanding the empirical enthalpic barrier in enzyme-catalyzed processes. We conclude by addressing the often conflicting interface between computational and experimental chemists, emphasizing the need for computation to predict experimental results in advance of their measurement.

Kinetic Isotope Effects and Transition State Structure for Human Phenylethanolamine N-Methyltransferase

Phenylethanolamine N-methyltransferase (PNMT) catalyzes the S-adenosyl-l-methionine (SAM)-dependent conversion of norepinephrine to epinephrine. Epinephrine has been associated with critical processes in humans including the control of respiration and blood pressure. Additionally, PNMT activity has been suggested to play a role in hypertension and Alzheimer's disease. In the current study, labeled SAM substrates were used to measure primary methyl-14C and 36S and secondary methyl-3H, 5'-3H, and 5'-14C intrinsic kinetic isotope effects for human PNMT. The transition state of human PNMT was modeled by matching kinetic isotope effects predicted via quantum chemical calculations to intrinsic values. The model provides information on the geometry and electrostatics of the human PNMT transition state structure and indicates that human PNMT catalyzes the formation of epinephrine through an early SN2 transition state in which methyl transfer is rate-limiting.

Structural and Functional Studies of Pavine N-Methyltransferase from Thalictrum flavum Reveal Novel Insights into Substrate Recognition and Catalytic Mechanism

Benzylisoquinoline alkaloids (BIAs) are produced in a wide variety of plants and include many common analgesic, antitussive, and anticancer compounds. Several members of a distinct family of S-adenosylmethionine (SAM)-dependent N-methyltransferases (NMTs) play critical roles in BIA biosynthesis, but the molecular basis of substrate recognition and catalysis is not known for NMTs involved in BIA metabolism. To address this issue, the crystal structure of pavine NMT from Thalictrum flavum was solved using selenomethionine-substituted protein (d_min = 2.8 Å). Additional structures were determined for the native protein (d_min = 2.0 Å) as well as binary complexes with SAM (d_min = 2.3 Å) or the reaction product S-adenosylhomocysteine (d_min = 1.6 Å). The structure of a complex with S-adenosylhomocysteine and two molecules of tetrahydropapaverine (THP; one as the S conformer and a second in the R configuration) (d_min = 1.8 Å) revealed key features of substrate recognition. Pavine NMT converted racemic THP to laudanosine, but the enzyme showed a preference for (±)-pavine and (S)-reticuline as substrates. These structures suggest the involvement of highly conserved residues at the active site. Mutagenesis of three residues near the methyl group of SAM and the nitrogen atom of the alkaloid acceptor decreased enzyme activity without disrupting the structure of the protein. The binding site for THP provides a framework for understanding substrate specificity among numerous NMTs involved in the biosynthesis of BIAs and other specialized metabolites. This information will facilitate metabolic engineering efforts aimed at producing medicinally important compounds in heterologous systems, such as yeast.