Ensemble perspective for catalytic promiscuity: calorimetric analysis of the active site conformational landscape of a detoxification enzyme

Enzymological paradigms have shifted recently to acknowledge the biological importance of catalytic promiscuity. However, catalytic promiscuity is a poorly understood property, and no thermodynamic treatment has described the conformational landscape of promiscuous versus substrate-specific enzymes. Here, two structurally similar glutathione transferase (GST, glutathione S-transferase) isoforms with high specificity or high promiscuity are compared. Differential scanning calorimetry (DSC) indicates a reversible low temperature transition for the promiscuous GSTA1-1 that is not observed with substrate-specific GSTA4-4. This transition is assigned to rearrangement of the C terminus at the active site of GSTA1-1 based on the effects of ligands and mutations. Near-UV and far-UV circular dichroism indicate that this transition is due to repacking of tertiary contacts with the remainder of the subunit, rather than "unfolding" of the C terminus per se. Analysis of the DSC data using a modified Landau theory indicates that the local conformational landscape of the active site of GSTA1-1 is smooth, with barrierless transitions between states. The partition function of the C-terminal states is a broad unimodal distribution at all temperatures within this DSC transition. In contrast, the remainder of the GSTA1-1 subunit and the GSTA4-4 protein exhibit folded and unfolded macrostates with a significant energy barrier separating them. Their partition function includes a sharp unimodal distribution of states only at temperatures that yield either folded or unfolded macrostates. At intermediate temperatures the partition function includes a bimodal distribution. The barrierless rearrangement of the GSTA1-1 active site within a local smooth energy landscape suggests a thermodynamic basis for catalytic promiscuity.     

Conformational Plasticity of the Essential Membrane-associated Mannosyltransferase PimA from Mycobacteria

Phosphatidyl-myo-inositol mannosyltransferase A (PimA) is an essential glycosyltransferase (GT) that initiates the biosynthetic pathway of phosphatidyl-myo-inositol mannosides, lipomannan, and lipoarabinomannan, which are key glycolipids/lipoglycans of the mycobacterial cell envelope. PimA belongs to a large family of peripheral membrane-associated GTs for which the understanding of the molecular mechanism and conformational changes that govern substrate/membrane recognition and catalysis remains a major challenge. Here we used single molecule force spectroscopy techniques to study the mechanical and conformational properties of PimA. In our studies, we engineered a polyprotein containing PimA flanked by four copies of the well characterized I27 protein, which provides an unambiguous mechanical fingerprint. We found that PimA exhibits weak mechanical stability albeit displaying β-sheet topology expected to unfold at much higher forces. Notably, PimA unfolds following heterogeneous multiple step mechanical unfolding pathways at low force akin to molten globule states. Interestingly, the ab initio low resolution envelopes obtained from small angle x-ray scattering of the unliganded PimA and the PimA·GDP complexed forms clearly demonstrate that not only the “open” and “closed” conformations of the GT-B enzyme are largely present in solution, but in addition, PimA experiences remarkable flexibility that undoubtedly corresponds to the N-terminal “Rossmann fold” domain, which has been proved to participate in protein-membrane interactions. Based on these results and on our previous experimental data, we propose a model wherein the conformational transitions are important for the mannosyltransferase to interact with the donor and acceptor substrates/membrane.     

Conformational states of syntaxin-1 govern the necessity of N-peptide binding in exocytosis of PC12 cells and Caenorhabditis elegans

Syntaxin-1 is the central SNARE protein for neuronal exocytosis. It interacts with Munc18-1 through its cytoplasmic domains, including the N-terminal peptide (N-peptide). Here we examine the role of the N-peptide binding in two conformational states (“closed” vs. “open”) of syntaxin-1 using PC12 cells and Caenorhabditis elegans. We show that expression of “closed” syntaxin-1A carrying N-terminal single point mutations (D3R, L8A) that perturb interaction with the hydrophobic pocket of Munc18-1 rescues impaired secretion in syntaxin-1–depleted PC12 cells and the lethality and lethargy of unc-64 (C. elegans orthologue of syntaxin-1)-null mutants. Conversely, expression of the “open” syntaxin-1A harboring the same mutations fails to rescue the impairments. Biochemically, the L8A mutation alone slightly weakens the binding between “closed” syntaxin-1A and Munc18-1, whereas the same mutation in the “open” syntaxin-1A disrupts it. Our results reveal a striking interplay between the syntaxin-1 N-peptide and the conformational state of the protein. We propose that the N-peptide plays a critical role in intracellular trafficking of syntaxin-1, which is dependent on the conformational state of this protein. Surprisingly, however, the N-peptide binding mode seems dispensable for SNARE-mediated exocytosis per se, as long as the protein is trafficked to the plasma membrane.     

Stochastic ensembles, conformationally adaptive teamwork, and enzymatic detoxification

It has been appreciated for a long time that enzymes exist as conformational ensembles throughout multiple stages of the reactions they catalyze, but there is renewed interest in the functional implications. The energy landscape that results from conformationlly diverse poteins is a complex surface with an energetic topography in multiple dimensions, even at the transition state(s) leading to product formation, and this represents a new paradigm. At the same time there has been renewed interest in conformational ensembles, a new paradigm concerning enzyme function has emerged, wherein catalytic promiscuity has clear biological advantages in some cases. "Useful", or biologically functional, promiscuity or the related behavior of "multifunctionality" can be found in the immune system, enzymatic detoxification, signal transduction, and the evolution of new function from an existing pool of folded protein scaffolds. Experimental evidence supports the widely held assumption that conformational heterogeneity promotes functional promiscuity. The common link between these coevolving paradigms is the inherent structural plasticity and conformational dynamics of proteins that, on one hand, lead to complex but evolutionarily selected energy landscapes and, on the other hand, promote functional promiscuity. Here we consider a logical extension of the overlap between these two nascent paradigms: functionally promiscuous and multifunctional enzymes such as detoxification enzymes are expected to have an ensemble landscape with more states accessible on multiple time scales than substrate specific enzymes. Two attributes of detoxification enzymes become important in the context of conformational ensembles: these enzymes metabolize multiple substrates, often in substrate mixtures, and they can form multiple products from a single substrate. These properties, combined with complex conformational landscapes, lead to the possibility of interesting time-dependent, or emergent, properties. Here we demonstrate these properties with kinetic simulations of nonequilibrium steady state (NESS) behavior resulting from energy landscapes expected for detoxification enzymes. Analogous scenarios with other promiscuous enzymes may be worthy of consideration.     

Active Site Loop Conformation Regulates Promiscuous Activity in a Lactonase from Geobacillus kaustophilus HTA426

Enzyme promiscuity is a prerequisite for fast divergent evolution of biocatalysts. A phosphotriesterase-like lactonase (PLL) from Geobacillus kaustophilus HTA426 (GkaP) exhibits main lactonase and promiscuous phosphotriesterase activities. To understand its catalytic and evolutionary mechanisms, we investigated a “hot spot” in the active site by saturation mutagenesis as well as X-ray crystallographic analyses. We found that position 99 in the active site was involved in substrate discrimination. One mutant, Y99L, exhibited 11-fold improvement over wild-type in reactivity (k_cat/K_m) toward the phosphotriesterase substrate ethyl-paraoxon, but showed 15-fold decrease toward the lactonase substrate δ-decanolactone, resulting in a 157-fold inversion of the substrate specificity. Structural analysis of Y99L revealed that the mutation causes a ∼6.6 Å outward shift of adjacent loop 7, which may cause increased flexibility of the active site and facilitate accommodation and/or catalysis of organophosphate substrate. This study provides for the PLL family an example of how the evolutionary route from promiscuity to specificity can derive from very few mutations, which promotes alteration in the conformational adjustment of the active site loops, in turn draws the capacity of substrate binding and activity.     

Ordered Regions of Channel Nucleoporins Nup62, Nup54, and Nup58 Form Dynamic Complexes in Solution

Three out of ∼30 nucleoporins, Nup62, Nup54, and Nup58, line the nuclear pore channel. These “channel” nucleoporins each contain an ordered region of ∼150–200 residues, which is predicted to be segmented into 3–4 α-helical regions of ∼40–80 residues. Notably, these segmentations are evolutionarily conserved between uni- and multicellular eukaryotes. Strikingly, the boundaries of these segments match our previously reported mapping and crystal data, which collectively identified two “cognate” segments of Nup54, each interacting with cognate segments, one in Nup58 and the other one in Nup62. Because Nup54 and Nup58 cognate segments form crystallographic hetero- or homo-oligomers, we proposed that these oligomers associate into inter-convertible “mid-plane” rings: a single large ring (40–50 nm diameter, consisting of eight hetero-dodecamers) or three small rings (10–20 nm diameter, each comprising eight homo-tetramers). Each “ring cycle” would recapitulate “dilation” and “constriction” of the nuclear pore complex''s central transport channel. As for the Nup54·Nup62 interactome, it forms a 1:2 triple helix (“finger”), multiples of which project alternately up and down from mid-plane ring(s). Collectively, our previous crystal data suggested a copy number of 128, 64, and 32 for Nup62, Nup54, and Nup58, respectively, that is, a 4:2:1 stoichiometry. Here, we carried out solution analysis utilizing the entire ordered regions of Nup62, Nup54, and Nup58, and demonstrate that they form a dynamic “triple complex” that is heterogeneously formed from our previously characterized Nup54·Nup58 and Nup54·Nup62 interactomes. These data are consistent both with our crystal structure-deduced copy numbers and stoichiometries and also with our ring cycle model for structure and dynamics of the nuclear pore channel.     

Structural Insights into the Catalytic Mechanism of Synechocystis Magnesium Protoporphyrin IX O-Methyltransferase (ChlM)

Magnesium protoporphyrin IX O-methyltransferase (ChlM) catalyzes transfer of the methyl group from S-adenosylmethionine to the carboxyl group of the C13 propionate side chain of magnesium protoporphyrin IX. This reaction is the second committed step in chlorophyll biosynthesis from protoporphyrin IX. Here we report the crystal structures of ChlM from the cyanobacterium Synechocystis sp. PCC 6803 in complex with S-adenosylmethionine and S-adenosylhomocysteine at resolutions of 1.6 and 1.7 Å, respectively. The structures illustrate the molecular basis for cofactor and substrate binding and suggest that conformational changes of the two “arm” regions may modulate binding and release of substrates/products to and from the active site. Tyr-28 and His-139 were identified to play essential roles for methyl transfer reaction but are not indispensable for cofactor/substrate binding. Based on these structural and functional findings, a catalytic model is proposed.     

The Novel New Jersey Eyewitness Instruction Induces Skepticism but Not Sensitivity

In recent decades, social scientists have shown that the reliability of eyewitness identifications is much worse than laypersons tend to believe. Although courts have only recently begun to react to this evidence, the New Jersey judiciary has reformed its jury instructions to notify jurors about the frailties of human memory, the potential for lineup administrators to nudge witnesses towards suspects that they police have already identified, and the advantages of alternative lineup procedures, including blinding of the administrator. This experiment tested the efficacy of New Jersey’s jury instruction. In a 2×2 between-subjects design, mock jurors (N = 335) watched a 35-minute murder trial, wherein identification quality was either “weak” or “strong” and either the New Jersey or a “standard” instruction was delivered. Jurors were more than twice as likely to convict when the standard instruction was used (OR = 2.55; 95% CI = 1.37–4.89, p < 0.001). The New Jersey instruction, however, did not improve juror''s ability to discern quality; rather, jurors receiving those instructions indiscriminatingly discounted “weak” and “strong” testimony in equal measure.     

Modularity of the hydrophobic core and evolution of functional diversity in fold A glycosyltransferases

Hydrophobic cores are fundamental structural properties of proteins typically associated with protein folding and stability; however, how the hydrophobic core shapes protein evolution and function is poorly understood. Here, we investigated the role of conserved hydrophobic cores in fold-A glycosyltransferases (GT-As), a large superfamily of enzymes that catalyze formation of glycosidic linkages between diverse donor and acceptor substrates through distinct catalytic mechanisms (inverting versus retaining). Using hidden Markov models and protein structural alignments, we identify similarities in the phosphate-binding cassette (PBC) of GT-As and unrelated nucleotide-binding proteins, such as UDP-sugar pyrophosphorylases. We demonstrate that GT-As have diverged from other nucleotide-binding proteins through structural elaboration of the PBC and its unique hydrophobic tethering to the F-helix, which harbors the catalytic base (xED-Asp). While the hydrophobic tethering is conserved across diverse GT-A fold enzymes, some families, such as B3GNT2, display variations in tethering interactions and core packing. We evaluated the structural and functional impact of these core variations through experimental mutational analysis and molecular dynamics simulations and find that some of the core mutations (T336I in B3GNT2) increase catalytic efficiency by modulating the conformational occupancy of the catalytic base between “D-in” and acceptor-accessible “D-out” conformation. Taken together, our studies support a model of evolution in which the GT-A core evolved progressively through elaboration upon an ancient PBC found in diverse nucleotide-binding proteins, and malleability of this core provided the structural framework for evolving new catalytic and substrate-binding functions in extant GT-A fold enzymes.     

Emergence,Development, and Maturity of the Gonad of Two Species of Chitons “Sea Cockroach” (Mollusca: Polyplacophora) through the Early Life Stages

This study describes and recognises, using histological and microscopical examinations on a morphometrical basis, several gonad traits through the early life stages of Chiton articulatus and C. albolineatus. Gonadal ontogenesis, gonad development stages, sexual differentiation, onset of the first sexual maturity, and growth sequences or “early life stages” were determined. In addition, allometry between lengths and body weight pooled for both sexes per each chiton were calculated using equation Y = aX^b. A total of 125 chitons (4≤TL≤40 mm, in total length “TL”) were used. All allometric relations showed a strong positive correlation (r), close to 1, with b-values above three, indicating an isometric growth. Gonadal ontogenesis and gonad development stages were categorised into three periods (“Pw” without gonad, “Pe” gonad emergence, and “Pf” gonadal sac formed) and four stages (“S0” gametocytogenesis, “S1” gametogenesis, “S2” mature, and “S3” spawning), respectively. Compound digital images were attained for each process. Periods and stages are overlapped among them and between species, with the following overall confidence intervals in TL: Pw 6.13–14.32 mm, Pe 10.32–16.93 mm, Pf 12.99–25.01 mm, S0 16.08–24.34 mm (females) and 19.51–26.60 mm (males), S1 27.15–35.63 mm (females) and 23.45–32.27 mm (males), S2 24.48–40.24 mm (females) and 25.45–32.87 mm (males). Sexual differentiation (in S0) of both chitons occurs first as a female then as a male; although, males reach the onset of the first sexual maturity earlier than females, thus for C. articulatus males at 17 mm and females at 32 mm, and for C. albolineatus males at 23.5 mm and females at 28 mm, all in TL. Four early life stages (i.e., subjuvenile, juvenile, subadult, and adult) are described and proposed to distinguish growth sequences. Our results may be useful to diverse disciplines, from developmental biology to fisheries management.     

Structure of Oxalacetate Acetylhydrolase,a Virulence Factor of the Chestnut Blight Fungus

Oxalacetate acetylhydrolase (OAH), a member of the phosphoenolpyruvate mutase/isocitrate lyase superfamily, catalyzes the hydrolysis of oxalacetate to oxalic acid and acetate. This study shows that knock-out of the oah gene in Cryphonectria parasitica, the chestnut blight fungus, reduces the ability of the fungus to form cankers on chestnut trees, suggesting that OAH plays a key role in virulence. OAH was produced in Escherichia coli and purified, and its catalytic rates were determined. Oxalacetate is the main OAH substrate, but the enzyme also acts as a lyase of (2R,3S)-dimethyl malate with ∼1000-fold lower efficacy. The crystal structure of OAH was determined alone, in complex with a mechanism-based inhibitor, 3,3-difluorooxalacetate (DFOA), and in complex with the reaction product, oxalate, to a resolution limit of 1.30, 1.55, and 1.65 Å, respectively. OAH assembles into a dimer of dimers with each subunit exhibiting an (α/β)₈ barrel fold and each pair swapping the 8th α-helix. An active site “gating loop” exhibits conformational disorder in the ligand-free structure. To obtain the structures of the OAH·ligand complexes, the ligand-free OAH crystals were soaked briefly with DFOA or oxalacetate. DFOA binding leads to ordering of the gating loop in a conformation that sequesters the ligand from the solvent. DFOA binds in a gem-diol form analogous to the oxalacetate intermediate/transition state. Oxalate binds in a planar conformation, but the gating loop is largely disordered. Comparison between the OAH structure and that of the closely related enzyme, 2,3-dimethylmalate lyase, suggests potential determinants of substrate preference.     

Structure and Function of the DUF2233 Domain in Bacteria and in the Human Mannose 6-Phosphate Uncovering Enzyme

DUF2233, a domain of unknown function (DUF), is present in many bacterial and several viral proteins and was also identified in the mammalian transmembrane glycoprotein N-acetylglucosamine-1-phosphodiester α-N-acetylglucosaminidase (“uncovering enzyme” (UCE)). We report the crystal structure of BACOVA_00430, a 315-residue protein from the human gut bacterium Bacteroides ovatus that is the first structural representative of the DUF2233 protein family. A notable feature of this structure is the presence of a surface cavity that is populated by residues that are highly conserved across the entire family. The crystal structure was used to model the luminal portion of human UCE (hUCE), which is involved in targeting of lysosomal enzymes. Mutational analysis of several residues in a highly conserved surface cavity of hUCE revealed that they are essential for function. The bacterial enzyme (BACOVA_00430) has ∼1% of the catalytic activity of hUCE toward the substrate GlcNAc-P-mannose, the precursor of the Man-6-P lysosomal targeting signal. GlcNAc-1-P is a poor substrate for both enzymes. We conclude that, for at least a subset of proteins in this family, DUF2233 functions as a phosphodiester glycosidase.     

Insights into diterpene cyclization from structure of bifunctional abietadiene synthase from Abies grandis

Abietadiene synthase from Abies grandis (AgAS) is a model system for diterpene synthase activity, catalyzing class I (ionization-initiated) and class II (protonation-initiated) cyclization reactions. Reported here is the crystal structure of AgAS at 2.3 Å resolution and molecular dynamics simulations of that structure with and without active site ligands. AgAS has three domains (α, β, and γ). The class I active site is within the C-terminal α domain, and the class II active site is between the N-terminal γ and β domains. The domain organization resembles that of monofunctional diterpene synthases and is consistent with proposed evolutionary origins of terpene synthases. Molecular dynamics simulations were carried out to determine the effect of substrate binding on enzymatic structure. Although such studies of the class I active site do lead to an enclosed substrate-Mg²⁺ complex similar to that observed in crystal structures of related plant enzymes, it does not enforce a single substrate conformation consistent with the known product stereochemistry. Simulations of the class II active site were more informative, with observation of a well ordered external loop migration. This “loop-in” conformation not only limits solvent access but also greatly increases the number of conformational states accessible to the substrate while destabilizing the nonproductive substrate conformation present in the “loop-out” conformation. Moreover, these conformational changes at the class II active site drive the substrate toward the proposed transition state. Docked substrate complexes were further assessed with regard to the effects of site-directed mutations on class I and II activities.     

Substrate-induced Changes in Protease Active Site Conformation Impact on Subsequent Reactions with Substrates

Enzymatic catalysis of biochemical reactions is essential to all living systems. The “lock and key” and “induced fit” models were early contributions to our understanding of the mechanisms involved in the reaction between an enzyme and its substrate. However, whether a given substrate-induced conformation is rigid or remains flexible has not yet been determined. By measuring the enzyme activity and intrinsic fluorescence of a nonspecific Eisenia fetida protease-I with different chromogenic substrates, we show that in subsequent reactions of protease with substrates, both the “lock and key” and “induced fit” mechanisms are used depending on the degree of conformational change required. Chromozym-Th- or chromosym-Ch-induced protease conformations were unable to bind chromozym-U. The chromosym-U-induced protease conformation remained flexible and could be further induced by chromozym-Th and chromozym-Ch. When low concentrations of guanidine HCl were used to disturb the conformation of the enzyme, only small changes in intrinsic fluorescence of the chromozym-Th-induced protease were detected, in contrast to the native enzyme whose intrinsic fluorescence markedly increased. This indicates that the substrate-induced enzyme was relatively rigid compared with the native protease. Utilizing a lock and key mechanism for secondary substrate reactions may have adaptive value in that it facilitates high efficiency in enzymatic reactions.     

Rapid Identification of Sequences for Orphan Enzymes to Power Accurate Protein Annotation

The power of genome sequencing depends on the ability to understand what those genes and their proteins products actually do. The automated methods used to assign functions to putative proteins in newly sequenced organisms are limited by the size of our library of proteins with both known function and sequence. Unfortunately this library grows slowly, lagging well behind the rapid increase in novel protein sequences produced by modern genome sequencing methods. One potential source for rapidly expanding this functional library is the “back catalog” of enzymology – “orphan enzymes,” those enzymes that have been characterized and yet lack any associated sequence. There are hundreds of orphan enzymes in the Enzyme Commission (EC) database alone. In this study, we demonstrate how this orphan enzyme “back catalog” is a fertile source for rapidly advancing the state of protein annotation. Starting from three orphan enzyme samples, we applied mass-spectrometry based analysis and computational methods (including sequence similarity networks, sequence and structural alignments, and operon context analysis) to rapidly identify the specific sequence for each orphan while avoiding the most time- and labor-intensive aspects of typical sequence identifications. We then used these three new sequences to more accurately predict the catalytic function of 385 previously uncharacterized or misannotated proteins. We expect that this kind of rapid sequence identification could be efficiently applied on a larger scale to make enzymology’s “back catalog” another powerful tool to drive accurate genome annotation.     

Higher Landing Accuracy in Expert Pilots is Associated with Lower Activity in the Caudate Nucleus

The most common lethal accidents in General Aviation are caused by improperly executed landing approaches in which a pilot descends below the minimum safe altitude without proper visual references. To understand how expertise might reduce such erroneous decision-making, we examined relevant neural processes in pilots performing a simulated landing approach inside a functional MRI scanner. Pilots (aged 20–66) were asked to “fly” a series of simulated “cockpit view” instrument landing scenarios in an MRI scanner. The scenarios were either high risk (heavy fog–legally unsafe to land) or low risk (medium fog–legally safe to land). Pilots with one of two levels of expertise participated: Moderate Expertise (Instrument Flight Rules pilots, n = 8) or High Expertise (Certified Instrument Flight Instructors or Air-Transport Pilots, n = 12). High Expertise pilots were more accurate than Moderate Expertise pilots in making a “land” versus “do not land” decision (CFII: d′ = 3.62±2.52; IFR: d′ = 0.98±1.04; p<.01). Brain activity in bilateral caudate nucleus was examined for main effects of expertise during a “land” versus “do not land” decision with the no-decision control condition modeled as baseline. In making landing decisions, High Expertise pilots showed lower activation in the bilateral caudate nucleus (0.97±0.80) compared to Moderate Expertise pilots (1.91±1.16) (p<.05). These findings provide evidence for increased “neural efficiency” in High Expertise pilots relative to Moderate Expertise pilots. During an instrument approach the pilot is engaged in detailed examination of flight instruments while monitoring certain visual references for making landing decisions. The caudate nucleus regulates saccade eye control of gaze, the brain area where the “expertise” effect was observed. These data provide evidence that performing “real world” aviation tasks in an fMRI provide objective data regarding the relative expertise of pilots and brain regions involved in it.     

Potentiation of Ligand Binding through Cooperative Effects in Monoamine Oxidase B

Crystallographic and biochemical studies have been employed to identify the binding site and mechanism for potentiation of imidazoline binding in human monoamine oxidase B (MAO B). 2-(2-Benzofuranyl)-2-imidazoline (2-BFI) inhibits recombinant human MAO B with a K_i of 8.3 ± 0.6 μm, whereas tranylcypromine-inhibited MAO B binds 2-BFI with a K_d of 9 ± 2 nm, representing an increase in binding energy Δ(ΔG) of −3.9 kcal/mol. Crystal structures show the imidazoline ligand bound in a site that is distinct from the substrate-binding cavity. Contributions to account for the increase in binding affinity upon tranylcypromine inhibition include a conformational change in the side chain of Gln²⁰⁶ and a “closed conformation” of the side chain of Ile¹⁹⁹, forming a hydrophobic “sandwich” with the side chain of Ile³¹⁶ on each face of the benzofuran ring of 2-BFI. Data with the I199A mutant of human MAO B and failure to observe a similar binding potentiation with rat MAO B, where Ile³¹⁶ is replaced with a Val residue, support an allosteric mechanism where the increased binding affinity of 2-BFI results from a cooperative increase in H-bond strength through formation of a more hydrophobic milieu. These insights should prove valuable in the design of high affinity and specific reversible MAO B inhibitors.