Mycolactone is the exotoxin produced by Mycobacterium ulcerans and is the virulence factor behind the neglected tropical disease Buruli ulcer. The toxin has a broad spectrum of biological effects ...within the host organism, stemming from its interaction with at least two molecular targets and the inhibition of protein uptake into the endoplasmic reticulum. Although it has been shown that the toxin can passively permeate into host cells, it is clearly lipophilic. Association with lipid carriers would have substantial implications for the toxin's distribution within a host organism, delivery to cellular targets, diagnostic susceptibility, and mechanisms of pathogenicity. Yet the toxin's interactions with, and distribution in, lipids are unknown. Herein we have used coarse-grained molecular dynamics simulations, guided by all-atom simulations, to study the interaction of mycolactone with pure and mixed lipid membranes. Using established techniques, we calculated the toxin's preferential localization, membrane translocation, and impact on membrane physical and dynamical properties. The computed water-octanol partition coefficient indicates that mycolactone prefers to be in an organic phase rather than in an aqueous environment. Our results show that in a solvated membrane environment the exotoxin mainly localizes in the water-membrane interface, with a preference for the glycerol moiety of lipids, consistent with the reported studies that found it in lipid extracts of the cell. The calculated association constant to the model membrane is similar to the reported association constant for Wiskott-Aldrich syndrome protein. Mycolactone is shown to modify the physical properties of membranes, lowering the transition temperature, compressibility modulus, and critical line tension at which pores can be stabilized. It also shows a tendency to behave as a linactant, a molecule that localizes at the boundary between different fluid lipid domains in membranes and promotes inter-mixing of domains. This property has implications for the toxin's cellular access, T-cell immunosuppression, and therapeutic potential.
Significance Proton transfer is a fundamental mechanism at the core of many enzyme-catalyzed reactions. It is also exquisitely sensitive to a number of factors, including pH, electrostatics, proper ...active-site geometry, and chemistry. Carbonic anhydrase has evolved a fast and efficient way to conduct protons through a combination of hydrophilic amino acid side chains that coordinate a highly ordered H-bonded water network. This study uses a powerful approach, combining NMR solution studies with neutron protein crystallography, to determine the effect of pH and divalent cations on key residues involved in proton transfer in human carbonic anhydrase. The results have broad implications for our understanding of proton transfer and how subtle changes in ionization and H-bonding interactions can modulate enzyme catalysis.
Human carbonic anhydrase II (HCA II) uses a Zn-bound OH ⁻/H ₂O mechanism to catalyze the reversible hydration of CO ₂. This catalysis also involves a separate proton transfer step, mediated by an ordered solvent network coordinated by hydrophilic residues. One of these residues, Tyr7, was previously shown to be deprotonated in the neutron crystal structure at pH 10. This observation indicated that Tyr7 has a perturbed pK ₐ compared with free tyrosine. To further probe the pK ₐ of this residue, NMR spectroscopic measurements of ¹³CTyr-labeled holo HCA II (with active-site Zn present) were preformed to titrate all Tyr residues between pH 5.4–11.0. In addition, neutron studies of apo HCA II (with Zn removed from the active site) at pH 7.5 and holo HCA II at pH 6 were conducted. This detailed interrogation of tyrosines in HCA II by NMR and neutron crystallography revealed a significantly lowered pK ₐ of Tyr7 and how pH and Tyr proximity to Zn affect hydrogen-bonding interactions.
A 1.1 Å resolution, room-temperature X-ray structure and a 2.1 Å resolution neutron structure of a chitin-degrading lytic polysaccharide monooxygenase domain from the bacterium Jonesia denitrificans ...(JdLPMO10A) show a putative dioxygen species equatorially bound to the active site copper. Both structures show an elongated density for the dioxygen, most consistent with a Cu(II)-bound peroxide. The coordination environment is consistent with Cu(II). In the neutron and X-ray structures, difference maps reveal the N-terminal amino group, involved in copper coordination, is present as a mixed ND2 and ND–, suggesting a role for the copper ion in shifting the pK a of the amino terminus.
The Protein Crystallography Station (PCS), located at the Los Alamos Neutron Scattering Center (LANSCE), was the first macromolecular crystallography beamline to be built at a spallation neutron ...source. Following testing and commissioning, the PCS user program was funded by the Biology and Environmental Research program of the Department of Energy Office of Science (DOE-OBER) for 13 years (2002-2014). The PCS remained the only dedicated macromolecular neutron crystallography station in North America until the construction and commissioning of the MaNDi and IMAGINE instruments at Oak Ridge National Laboratory, which started in 2012. The instrument produced a number of research and technical outcomes that have contributed to the field, clearly demonstrating the power of neutron crystallo-graphy in helping scientists to understand enzyme reaction mechanisms, hydrogen bonding and visualization of H-atom positions, which are critical to nearly all chemical reactions. During this period, neutron crystallography became a technique that increasingly gained traction, and became more integrated into macromolecular crystallography through software developments led by investigators at the PCS. This review highlights the contributions of the PCS to macromolecular neutron crystallography, and gives an overview of the history of neutron crystallography and the development of macromolecular neutron crystallography from the 1960s to the 1990s and onwards through the 2000s.
Our knowledge of the metabolites in cells and their reactions is far from complete as revealed by metabolomic measurements that detect many more small molecules than are documented in metabolic ...databases. Here, we develop an approach for predicting the reactivity of small-molecule metabolites in enzyme-catalyzed reactions that combines expert knowledge, computational chemistry and machine learning.
We classified 4843 reactions documented in the KEGG database, from all six Enzyme Commission classes (EC 1-6), into 80 reaction classes, each of which is marked by a characteristic functional group transformation. Reaction centers and surrounding local structures in substrates and products of these reactions were represented using SMARTS. We found that each of the SMARTS-defined chemical substructures is widely distributed among metabolites, but only a fraction of the functional groups in these substructures are reactive. Using atomic properties of atoms in a putative reaction center and molecular properties as features, we trained support vector machine (SVM) classifiers to discriminate between functional groups that are reactive and non-reactive. Classifier accuracy was assessed by cross-validation analysis. A typical sensitivity TP/(TP+FN) or specificity TN/(TN+FP) is ≈0.8. Our results suggest that metabolic reactivity of small-molecule compounds can be predicted with reasonable accuracy based on the presence of a potentially reactive functional group and the chemical features of its local environment.
The classifiers presented here can be used to predict reactions via a web site (http://cellsignaling.lanl.gov/Reactivity/). The web site is freely available.
The clustering of genes in a pathway and the co-location of functionally related genes is widely recognized in prokaryotes. We used these characteristics to predict the metabolic involvement for a ...Transcriptional Regulator (TR) of unknown function, identified and confirmed its biological activity.
A software tool that identifies the genes encoded within a defined genomic neighborhood for the subject TR and its homologs was developed. The output lists of genes in the genetic neighborhoods, their annotated functions, the reactants/products, and identifies the metabolic pathway in which the encoded-proteins function. When a set of TRs of known function was analyzed, we observed that their homologs frequently had conserved genomic neighborhoods that co-located the metabolically related genes regulated by the subject TR. We postulate that TR effectors are metabolites in the identified pathways; indeed the known effectors were present. We analyzed Bxe_B3018 from Burkholderia xenovorans, a TR of unknown function and predicted that this TR was related to the glycine, threonine and serine degradation. We tested the binding of metabolites in these pathways and for those that bound, their ability to modulate TR binding to its specific DNA operator sequence. Using rtPCR, we confirmed that methylglyoxal was an effector of Bxe_3018.
These studies provide the proof of concept and validation of a systematic approach to the discovery of the biological activity for proteins of unknown function, in this case a TR. Bxe_B3018 is a methylglyoxal responsive TR that controls the expression of an operon composed of a putative efflux system.
Three high‐resolution X‐ray crystal structures of malate dehydrogenase (MDH; EC 1.1.1.37) from the methylotroph Methylobacterium extorquens AM1 are presented. By comparing the structures of apo MDH, ...a binary complex of MDH and NAD+, and a ternary complex of MDH and oxaloacetate with ADP‐ribose occupying the pyridine nucleotide‐binding site, conformational changes associated with the formation of the catalytic complex were characterized. While the substrate‐binding site is accessible in the enzyme resting state or NAD+‐bound forms, the substrate‐bound form exhibits a closed conformation. This conformational change involves the transition of an α‐helix to a 310‐helix, which causes the adjacent loop to close the active site following coenzyme and substrate binding. In the ternary complex, His284 forms a hydrogen bond to the C2 carbonyl of oxaloacetate, placing it in a position to donate a proton in the formation of (2S)‐malate.
Crystal structures of apo malate dehydrogenase (MDH) from Methylobacterium extorquens, MDH bound to NAD+, and MDH with oxaloacetate and ADP‐ribose revealed conformational changes, closing the active site upon coenzyme and substrate binding. In the ternary complex, His284 is in position to donate a proton in the formation of (2S)‐malate.
The most abundant carbohydrate product of cellulosic biomass pyrolysis is the anhydrosugar levoglucosan (1,6-anhydro-β-d-glucopyranose), which can be converted to glucose 6-phosphate by levoglucosan ...kinase (LGK). In addition to the canonical kinase phosphotransfer reaction, the conversion requires cleavage of the 1,6-anhydro ring to allow ATP-dependent phosphorylation of the sugar O6 atom. Using x-ray crystallography, we show that LGK binds two magnesium ions in the active site that are additionally coordinated with the nucleotide and water molecules to result in ideal octahedral coordination. To further verify the metal binding sites, we co-crystallized LGK in the presence of manganese instead of magnesium and solved the structure de novo using the anomalous signal from four manganese atoms in the dimeric structure. The first metal is required for catalysis, whereas our work suggests that the second is either required or significantly promotes the catalytic rate. Although the enzyme binds its sugar substrate in a similar orientation to the structurally related 1,6-anhydro-N-acetylmuramic acid kinase (AnmK), it forms markedly fewer bonding interactions with the substrate. In this orientation, the sugar is in an optimal position to couple phosphorylation with ring cleavage. We also observed a second alternate binding orientation for levoglucosan, and in these structures, ADP was found to bind with lower affinity. These combined observations provide an explanation for the high Km of LGK for levoglucosan. Greater knowledge of the factors that contribute to the catalytic efficiency of LGK can be used to improve applications of this enzyme for levoglucosan-derived biofuel production.
Engineered microbes can be used for producing value-added chemicals from renewable feedstocks, relieving the dependency on nonrenewable resources such as petroleum. These microbes often are composed ...of synthetic metabolic pathways; however, one major problem in establishing a synthetic pathway is the challenge of precisely controlling competing metabolic routes, some of which could be crucial for fitness and survival. While traditional gene deletion and/or coarse overexpression approaches do not provide precise regulation, cis-repressors (CRs) are RNA-based regulatory elements that can control the production levels of a particular protein in a tunable manner. Here, we describe a protocol for a generally applicable fluorescence-activated cell sorting technique used to isolate eight subpopulations of CRs from a semidegenerate library in Escherichia coli, followed by deep sequencing that permitted the identification of 15 individual CRs with a broad range of protein production profiles. Using these new CRs, we demonstrated a change in production levels of a fluorescent reporter by over two orders of magnitude and further showed that these CRs are easily ported from E. coli to Pseudomonas putida. We next used four CRs to tune the production of the enzyme PpsA, involved in pyruvate to phosphoenolpyruvate (PEP) conversion, to alter the pool of PEP that feeds into the shikimate pathway. In an engineered P. putida strain, where carbon flux in the shikimate pathway is diverted to the synthesis of the commodity chemical cis,cis-muconate, we found that tuning PpsA translation levels increased the overall titer of muconate. Therefore, CRs provide an approach to precisely tune protein levels in metabolic pathways and will be an important tool for other metabolic engineering efforts.
Until recently, engineering strategies for altering gene expression have focused on transcription control using strong inducible promoters or one of several methods to knock down wasteful genes. ...Recently, synthetic riboregulators have been developed for translational regulation of gene expression. Here, we report a new modular synthetic riboregulator class that has the potential to finely tune protein expression and independently control the concentration of each enzyme in an engineered metabolic pathway. This development is important because the most straightforward approach to altering the flux through a particular metabolic step is to increase or decrease the concentration of the enzyme. Our design includes a cis-repressor at the 5′ end of the mRNA that forms a stem-loop helix, occluding the ribosomal binding sequence and blocking translation. A trans-expressed activating-RNA frees the ribosomal-binding sequence, which turns on translation. The overall architecture of the riboregulators is designed using Watson–Crick base-pairing stability. We describe here a cis-repressor that can completely shut off translation of antibiotic-resistance reporters and a trans-activator that restores translation. We have established that it is possible to use these riboregulators to achieve translational control of gene expression over a wide dynamic range. We have also found that a targeting sequence can be modified to develop riboregulators that can, in principle, independently regulate translation of many genes. In a selection experiment, we demonstrated that by subtly altering the sequence of the trans-activator it is possible to alter the ratio of the repressed and activated states and to achieve intermediate translational control.