Putative rhamnogalacturonan-II glycosyltransferase identified through callus gene editing bypasses embryo lethality.

Rhamnogalacturonan II (RG-II) is a structurally complex and conserved domain of the pectin present in the primary cell walls of vascular plants. Borate crosslinking of RG-II is required for plants to grow and develop normally. Mutations that alter RG-II structure also affect crosslinking and are lethal or severely impair growth. Thus, few genes involved in RG-II synthesis have been identified. Here we developed a method to generate viable loss-of-function Arabidopsis (Arabidopsis thaliana) mutants in callus tissue via CRISPR/Cas9-mediated gene editing. We combined this with a candidate gene approach to characterize the male gametophyte defective 2 (MPG2) gene that encodes a putative family GT29 glycosyltransferase. Plants homozygous for this mutation do not survive. We showed that in the callus mutant cell walls, RG-II does not crosslink normally because it lacks 3-deoxy-D-manno-octulosonic acid (Kdo) and thus cannot form the α-L-Rhap-(1→5)-α-D-kdop-(1→ sidechain. We suggest that MGP2 encodes an inverting RG-II CMP-ß-Kdo transferase (RCKT1). Our discovery provides further insight into the role of sidechains in RG-II dimerization. Our method also provides a viable strategy for further identifying proteins involved in the biosynthesis of RG-II.

Complete integration of carbene-transfer chemistry into biosynthesis.

Biosynthesis is an environmentally benign and renewable approach that can be used to produce a broad range of natural and, in some cases, new-to-nature products. However, biology lacks many of the reactions that are available to synthetic chemists, resulting in a narrower scope of accessible products when using biosynthesis rather than synthetic chemistry. A prime example of such chemistry is carbene-transfer reactions1. Although it was recently shown that carbene-transfer reactions can be performed in a cell and used for biosynthesis2,3, carbene donors and unnatural cofactors needed to be added exogenously and transported into cells to effect the desired reactions, precluding cost-effective scale-up of the biosynthesis process with these reactions. Here we report the access to a diazo ester carbene precursor by cellular metabolism and a microbial platform for introducing unnatural carbene-transfer reactions into biosynthesis. The α-diazoester azaserine was produced by expressing a biosynthetic gene cluster in Streptomyces albus. The intracellularly produced azaserine was used as a carbene donor to cyclopropanate another intracellularly produced molecule-styrene. The reaction was catalysed by engineered P450 mutants containing a native cofactor with excellent diastereoselectivity and a moderate yield. Our study establishes a scalable, microbial platform for conducting intracellular abiological carbene-transfer reactions to functionalize a range of natural and new-to-nature products and expands the scope of organic products that can be produced by cellular metabolism.

Structural and biochemical insight into a modular ß-1,4-galactan synthase in plants.

Rhamnogalacturonan I (RGI) is a structurally complex pectic polysaccharide with a backbone of alternating rhamnose and galacturonic acid residues substituted with arabinan and galactan side chains. Galactan synthase 1 (GalS1) transfers galactose and arabinose to either extend or cap the ß-1,4-galactan side chains of RGI, respectively. Here we report the structure of GalS1 from Populus trichocarpa, showing a modular protein consisting of an N-terminal domain that represents the founding member of a new family of carbohydrate-binding module, CBM95, and a C-terminal glycosyltransferase family 92 (GT92) catalytic domain that adopts a GT-A fold. GalS1 exists as a dimer in vitro, with stem domains interacting across the chains in a 'handshake' orientation that is essential for maintaining stability and activity. In addition to understanding the enzymatic mechanism of GalS1, we gained insight into the donor and acceptor substrate binding sites using deep evolutionary analysis, molecular simulations and biochemical studies. Combining all the results, a mechanism for GalS1 catalysis and a new model for pectic galactan side-chain addition are proposed.

Efficient production of oxidized terpenoids via engineering fusion proteins of terpene synthase and cytochrome P450.

The functionalization of terpenes using cytochrome P450 enzymes is a versatile route to the production of useful derivatives that can be further converted to value-added products. Many terpenes are hydrophobic and volatile making their availability as a substrate for P450 enzymes significantly limited during microbial production. In this study, we developed a strategy to improve the accessibility of terpene molecules for the P450 reaction by linking terpene synthase and P450 together. As a model system, fusion proteins of 1,8-cineole synthase (CS) and P450cin were investigated and it showed an improved hydroxylation of the monoterpenoid 1,8-cineole up to 5.4-fold. Structural analysis of the CS-P450cin fusion proteins by SEC-SAXS indicated a dimer formation with preferred orientations of the active sites of the two domains. We also applied the enzyme fusion strategy to the oxidation of a sesquiterpene epi-isozizaene and the fusion enzymes significantly improved albaflavenol production in engineered E. coli. From the analysis of positive and negative examples of the fusion strategy, we proposed key factors in structure-based prediction and evaluation of fusion enzymes. Developing fusion enzymes for terpene synthase and P450 presents an efficient strategy toward oxidation of hydrophobic terpene compounds. This strategy could be widely applicable to improve the biosynthetic titer of the functionalized products from hydrophobic terpene intermediates.

An iron (II) dependent oxygenase performs the last missing step of plant lysine catabolism.

Despite intensive study, plant lysine catabolism beyond the 2-oxoadipate (2OA) intermediate remains unvalidated. Recently we described a missing step in the D-lysine catabolism of Pseudomonas putida in which 2OA is converted to D-2-hydroxyglutarate (2HG) via hydroxyglutarate synthase (HglS), a DUF1338 family protein. Here we solve the structure of HglS to 1.1 Å resolution in substrate-free form and in complex with 2OA. We propose a successive decarboxylation and intramolecular hydroxylation mechanism forming 2HG in a Fe(II)- and O2-dependent manner. Specificity is mediated by a single arginine, highly conserved across most DUF1338 proteins. An Arabidopsis thaliana HglS homolog coexpresses with known lysine catabolism enzymes, and mutants show phenotypes consistent with disrupted lysine catabolism. Structural and biochemical analysis of Oryza sativa homolog FLO7 reveals identical activity to HglS despite low sequence identity. Our results suggest DUF1338-containing enzymes catalyze the same biochemical reaction, exerting the same physiological function across bacteria and eukaryotes.

Structural Mechanism of Regioselectivity in an Unusual Bacterial Acyl-CoA Dehydrogenase.

Terminal alkenes are easily derivatized, making them desirable functional group targets for polyketide synthase (PKS) engineering. However, they are rarely encountered in natural PKS systems. One mechanism for terminal alkene formation in PKSs is through the activity of an acyl-CoA dehydrogenase (ACAD). Herein, we use biochemical and structural analysis to understand the mechanism of terminal alkene formation catalyzed by an γ,δ-ACAD from the biosynthesis of the polyketide natural product FK506, TcsD. While TcsD is homologous to canonical α,ß-ACADs, it acts regioselectively at the γ,δ-position and only on α,ß-unsaturated substrates. Furthermore, this regioselectivity is controlled by a combination of bulky residues in the active site and a lateral shift in the positioning of the FAD cofactor within the enzyme. Substrate modeling suggests that TcsD utilizes a novel set of hydrogen bond donors for substrate activation and positioning, preventing dehydrogenation at the α,ß position of substrates. From the structural and biochemical characterization of TcsD, key residues that contribute to regioselectivity and are unique to the protein family were determined and used to identify other putative γ,δ-ACADs that belong to diverse natural product biosynthetic gene clusters. These predictions are supported by the demonstration that a phylogenetically distant homologue of TcsD also regioselectively oxidizes α,ß-unsaturated substrates. This work exemplifies a powerful approach to understand unique enzymatic reactions and will facilitate future enzyme discovery, inform enzyme engineering, and aid natural product characterization efforts.

Massively Parallel Fitness Profiling Reveals Multiple Novel Enzymes in Pseudomonas putida Lysine Metabolism.

Despite intensive study for 50 years, the biochemical and genetic links between lysine metabolism and central metabolism in Pseudomonas putida remain unresolved. To establish these biochemical links, we leveraged random barcode transposon sequencing (RB-TnSeq), a genome-wide assay measuring the fitness of thousands of genes in parallel, to identify multiple novel enzymes in both l- and d-lysine metabolism. We first describe three pathway enzymes that catabolize l-2-aminoadipate (l-2AA) to 2-ketoglutarate (2KG), connecting d-lysine to the TCA cycle. One of these enzymes, P. putida 5260 (PP_5260), contains a DUF1338 domain, representing a family with no previously described biological function. Our work also identified the recently described coenzyme A (CoA)-independent route of l-lysine degradation that results in metabolization to succinate. We expanded on previous findings by demonstrating that glutarate hydroxylase CsiD is promiscuous in its 2-oxoacid selectivity. Proteomics of selected pathway enzymes revealed that expression of catabolic genes is highly sensitive to the presence of particular pathway metabolites, implying intensive local and global regulation. This work demonstrated the utility of RB-TnSeq for discovering novel metabolic pathways in even well-studied bacteria, as well as its utility a powerful tool for validating previous research.IMPORTANCEP. putida lysine metabolism can produce multiple commodity chemicals, conferring great biotechnological value. Despite much research, the connection of lysine catabolism to central metabolism in P. putida remained undefined. Here, we used random barcode transposon sequencing to fill the gaps of lysine metabolism in P. putida We describe a route of 2-oxoadipate (2OA) catabolism, which utilizes DUF1338-containing protein P. putida 5260 (PP_5260) in bacteria. Despite its prevalence in many domains of life, DUF1338-containing proteins have had no known biochemical function. We demonstrate that PP_5260 is a metalloenzyme which catalyzes an unusual route of decarboxylation of 2OA to d-2-hydroxyglutarate (d-2HG). Our screen also identified a recently described novel glutarate metabolic pathway. We validate previous results and expand the understanding of glutarate hydroxylase CsiD by showing that can it use either 2OA or 2KG as a cosubstrate. Our work demonstrated that biological novelty can be rapidly identified using unbiased experimental genetics and that RB-TnSeq can be used to rapidly validate previous results.

Computational design of self-assembling cyclic protein homo-oligomers.

Self-assembling cyclic protein homo-oligomers play important roles in biology, and the ability to generate custom homo-oligomeric structures could enable new approaches to probe biological function. Here we report a general approach to design cyclic homo-oligomers that employs a new residue-pair-transform method to assess the designability of a protein-protein interface. This method is sufficiently rapid to enable the systematic enumeration of cyclically docked arrangements of a monomer followed by sequence design of the newly formed interfaces. We use this method to design interfaces onto idealized repeat proteins that direct their assembly into complexes that possess cyclic symmetry. Of 96 designs that were characterized experimentally, 21 were found to form stable monodisperse homo-oligomers in solution, and 15 (four homodimers, six homotrimers, six homotetramers and one homopentamer) had solution small-angle X-ray scattering data consistent with the design models. X-ray crystal structures were obtained for five of the designs and each is very close to their corresponding computational model.

Principles for designing proteins with cavities formed by curved ß sheets.

Active sites and ligand-binding cavities in native proteins are often formed by curved ß sheets, and the ability to control ß-sheet curvature would allow design of binding proteins with cavities customized to specific ligands. Toward this end, we investigated the mechanisms controlling ß-sheet curvature by studying the geometry of ß sheets in naturally occurring protein structures and folding simulations. The principles emerging from this analysis were used to design, de novo, a series of proteins with curved ß sheets topped with α helices. Nuclear magnetic resonance and crystal structures of the designs closely match the computational models, showing that ß-sheet curvature can be controlled with atomic-level accuracy. Our approach enables the design of proteins with cavities and provides a route to custom design ligand-binding and catalytic sites.

Computational design of a homotrimeric metalloprotein with a trisbipyridyl core.

Metal-chelating heteroaryl small molecules have found widespread use as building blocks for coordination-driven, self-assembling nanostructures. The metal-chelating noncanonical amino acid (2,2'-bipyridin-5yl)alanine (Bpy-ala) could, in principle, be used to nucleate specific metalloprotein assemblies if introduced into proteins such that one assembly had much lower free energy than all alternatives. Here we describe the use of the Rosetta computational methodology to design a self-assembling homotrimeric protein with [Fe(Bpy-ala)3]2+ complexes at the interface between monomers. X-ray crystallographic analysis of the homotrimer showed that the design process had near-atomic-level accuracy: The all-atom rmsd between the design model and crystal structure for the residues at the protein interface is ∼1.4 Å. These results demonstrate that computational protein design together with genetically encoded noncanonical amino acids can be used to drive formation of precisely specified metal-mediated protein assemblies that could find use in a wide range of photophysical applications.

De novo design of protein homo-oligomers with modular hydrogen-bond network-mediated specificity.

In nature, structural specificity in DNA and proteins is encoded differently: In DNA, specificity arises from modular hydrogen bonds in the core of the double helix, whereas in proteins, specificity arises largely from buried hydrophobic packing complemented by irregular peripheral polar interactions. Here, we describe a general approach for designing a wide range of protein homo-oligomers with specificity determined by modular arrays of central hydrogen-bond networks. We use the approach to design dimers, trimers, and tetramers consisting of two concentric rings of helices, including previously not seen triangular, square, and supercoiled topologies. X-ray crystallography confirms that the structures overall, and the hydrogen-bond networks in particular, are nearly identical to the design models, and the networks confer interaction specificity in vivo. The ability to design extensive hydrogen-bond networks with atomic accuracy enables the programming of protein interaction specificity for a broad range of synthetic biology applications; more generally, our results demonstrate that, even with the tremendous diversity observed in nature, there are fundamentally new modes of interaction to be discovered in proteins.

Structural and Biochemical Characterization of the Early and Late Enzymes in the Lignin ß-Aryl Ether Cleavage Pathway from Sphingobium sp. SYK-6.

There has been great progress in the development of technology for the conversion of lignocellulosic biomass to sugars and subsequent fermentation to fuels. However, plant lignin remains an untapped source of materials for production of fuels or high value chemicals. Biological cleavage of lignin has been well characterized in fungi, in which enzymes that create free radical intermediates are used to degrade this material. In contrast, a catabolic pathway for the stereospecific cleavage of ß-aryl ether units that are found in lignin has been identified in Sphingobium sp. SYK-6 bacteria. ß-Aryl ether units are typically abundant in lignin, corresponding to 50-70% of all of the intermonomer linkages. Consequently, a comprehensive understanding of enzymatic ß-aryl ether (ß-ether) cleavage is important for future efforts to biologically process lignin and its breakdown products. The crystal structures and biochemical characterization of the NAD-dependent dehydrogenases (LigD, LigO, and LigL) and the glutathione-dependent lyase LigG provide new insights into the early and late enzymes in the ß-ether degradation pathway. We present detailed information on the cofactor and substrate binding sites and on the catalytic mechanisms of these enzymes, comparing them with other known members of their respective families. Information on the Lig enzymes provides new insight into their catalysis mechanisms and can inform future strategies for using aromatic oligomers derived from plant lignin as a source of valuable aromatic compounds for biofuels and other bioproducts.

Structural Basis of Stereospecificity in the Bacterial Enzymatic Cleavage of ß-Aryl Ether Bonds in Lignin.

Lignin is a combinatorial polymer comprising monoaromatic units that are linked via covalent bonds. Although lignin is a potential source of valuable aromatic chemicals, its recalcitrance to chemical or biological digestion presents major obstacles to both the production of second-generation biofuels and the generation of valuable coproducts from lignin's monoaromatic units. Degradation of lignin has been relatively well characterized in fungi, but it is less well understood in bacteria. A catabolic pathway for the enzymatic breakdown of aromatic oligomers linked via ß-aryl ether bonds typically found in lignin has been reported in the bacterium Sphingobium sp. SYK-6. Here, we present x-ray crystal structures and biochemical characterization of the glutathione-dependent ß-etherases, LigE and LigF, from this pathway. The crystal structures show that both enzymes belong to the canonical two-domain fold and glutathione binding site architecture of the glutathione S-transferase family. Mutagenesis of the conserved active site serine in both LigE and LigF shows that, whereas the enzymatic activity is reduced, this amino acid side chain is not absolutely essential for catalysis. The results include descriptions of cofactor binding sites, substrate binding sites, and catalytic mechanisms. Because ß-aryl ether bonds account for 50-70% of all interunit linkages in lignin, understanding the mechanism of enzymatic ß-aryl ether cleavage has significant potential for informing ongoing studies on the valorization of lignin.

Structure of the OsSERK2 leucine-rich repeat extracellular domain.

Somatic embryogenesis receptor kinases (SERKs) are leucine-rich repeat (LRR)-containing integral membrane receptors that are involved in the regulation of development and immune responses in plants. It has recently been shown that rice SERK2 (OsSERK2) is essential for XA21-mediated resistance to the pathogen Xanthomonas oryzae pv. oryzae. OsSERK2 is also required for the BRI1-mediated, FLS2-mediated and EFR-mediated responses to brassinosteroids, flagellin and elongation factor Tu (EF-Tu), respectively. Here, crystal structures of the LRR domains of OsSERK2 and a D128N OsSERK2 mutant, expressed as hagfish variable lymphocyte receptor (VLR) fusions, are reported. These structures suggest that the aspartate mutation does not generate any significant conformational change in the protein, but instead leads to an altered interaction with partner receptors.

Molecular simulations provide new insights into the role of the accessory immunoglobulin-like domain of Cel9A.

Cel9A from the thermoacidophilic bacterium Alicyclobacillus acidocaldarius belongs to the subfamily E1 of family 9 glycoside hydrolases, many members of which have an N-terminal Ig-like domain followed by the catalytic domain. The Ig-like domain is not directly involved in either carbohydrate binding or biocatalysis; however, deletion of the Ig-domain promotes loss of enzymatic activity. We have investigated the functional role of the Ig-like domain using molecular dynamics simulations. Our simulations indicate that residues within the Ig-like domain are dynamically correlated with residues in the carbohydrate-binding pocket and with key catalytic residues of Cel9A. Free energy perturbation simulations indicate that the Ig-like domain stabilizes the catalytic domain and may be responsible for the enhanced thermostability of Cel9A.

Structure of endoglucanase Cel9A from the thermoacidophilic Alicyclobacillus acidocaldarius.

The production of biofuels using biomass is an alternative route to support the growing global demand for energy and to also reduce the environmental problems caused by the burning of fossil fuels. Cellulases are likely to play an important role in the degradation of biomass and the production of sugars for subsequent fermentation to fuel. Here, the crystal structure of an endoglucanase, Cel9A, from Alicyclobacillus acidocaldarius (Aa_Cel9A) is reported which displays a modular architecture composed of an N-terminal Ig-like domain connected to the catalytic domain. This paper describes the overall structure and the detailed contacts between the two modules. Analysis suggests that the interaction involving the residues Gln13 (from the Ig-like module) and Phe439 (from the catalytic module) is important in maintaining the correct conformation of the catalytic module required for protein activity. Moreover, the Aa_Cel9A structure shows three metal-binding sites that are associated with the thermostability and/or substrate affinity of the enzyme.

Purification of human transcription factors Nanog and Sox2, each in complex with Skp, an Escherichia coli periplasmic chaperone.

Nanog and Sox2 are key transcriptional factors involved in self-renewal and pluripotency of stem cells in human and other mammals. Nanog and Sox2 contain homeodomain (HD) and high-mobility group (HMG) DNA-binding domain, respectively, for targeting them to their regulatory regions and the other regions with transactivation function by providing sites for recruiting other transcriptional regulators. To gain insights in the biochemical and biophysical characteristics of the other regions of Nanog and Sox2, we have tried to overproduce and purify full length wild-type human Nanog and Sox2 expressed in Escherichia coli. Interestingly, we found that Nanog and Sox2 were individually stabilized by tight interaction with Skp, an E. coli periplasmic chaperone, thereby enabling stable over-expression and purification of Nanog and Sox2, each in complex with Skp. Purified Skp complexes of Nanog and Sox maintained DNA-binding activity toward its cognate DNA sequence. A similar approach may be applicable for some other mammalian proteins that are unstable or difficult to over-express in E. coli.

Structure of human Brn-5 transcription factor in complex with CRH gene promoter.

The Brn-5 protein, highly expressed in human brain, belongs to the POU family; a class of transcription factors involved in a wide variety of biological processes ranging from programming of embryonic stem cells to cellular housekeeping. This functional diversity is conferred by two DNA-binding subdomains that can assume several configurations due to a bipartite arrangement of POU-specific (POU(S)) and POU-homeo (POU(H)) subdomains separated by a linker region. The crystal structure of human Brn-5 transcription factor in complex with corticotrophin-releasing hormone (CRH) gene promoter reveals an unexpected recognition mode of the protein to its cognate DNA. Moreover, the structure also shows the role of the linker in allowing diverse configurations that can be assumed by the two subdomains.

Crystallization and preliminary X-ray analysis of human Brn-5 transcription factor in complex with DNA.

The Brn-5 protein plays an important role in the control of cellular development and belongs to a class of transcription factors that usually contain two domains: the POU homeodomain (POU(HD)) and the POU-specific domain (POU(S)). Since high-quality crystals suitable for crystallographic studies of the proteins of this class are difficult to obtain, all the known structural information available is for POU(HD) and/or POU(S). This paper describes several critical steps that allowed the production of high-quality crystals of the full-length Brn-5 protein complexed with its cognate DNA.

The inhibition of 5-enolpyruvylshikimate-3-phosphate synthase as a model for development of novel antimicrobials.

EPSP synthase (EPSPS) is an essential enzyme in the shikimate pathway, transferring the enolpyruvyl group of phosphoenolpyruvate to shikimate-3-phosphate to form 5-enolpyruvyl-3-shikimate phosphate and inorganic phosphate. This enzyme is composed of two domains, which are formed by three copies of betaalphabetaalphabetabeta-folding units; in between there are two crossover chain segments hinging the nearly topologically symmetrical domains together and allowing conformational changes necessary for substrate conversion. The reaction is ordered with shikimate-3-phosphate binding first, followed by phosphoenolpyruvate, and then by the subsequent release of phosphate and EPSP. N-[phosphomethyl]glycine (glyphosate) is the commercial inhibitor of this enzyme. Apparently, the binding of shikimate-3-phosphate is necessary for glyphosate binding, since it induces the closure of the two domains to form the active site in the interdomain cleft. However, it is somehow controversial whether binding of shikimate-3-phosphate alone is enough to induce the complete conversion to the closed state. The phosphoenolpyruvate binding site seems to be located mainly on the C-terminal domain, while the binding site of shikimate-3-phosphate is located primarily in the N-terminal domain residues. However, recent results demonstrate that the active site of the enzyme undergoes structural changes upon inhibitor binding on a scale that cannot be predicted by conventional computational methods. Studies of molecular docking based on the interaction of known EPSPS structures with (R)- phosphonate TI analogue reveal that more experimental data on the structure and dynamics of various EPSPS-ligand complexes are needed to more effectively apply structure-based drug design of this enzyme in the future.