Comparative functional genomics identifies an iron-limited bottleneck in a Saccharomyces cerevisiae strain with a cytosolic-localized isobutanol pathway.

Metabolic engineering strategies have been successfully implemented to improve the production of isobutanol, a next-generation biofuel, in Saccharomyces cerevisiae. Here, we explore how two of these strategies, pathway re-localization and redox cofactor-balancing, affect the performance and physiology of isobutanol producing strains. We equipped yeast with isobutanol cassettes which had either a mitochondrial or cytosolic localized isobutanol pathway and used either a redox-imbalanced (NADPH-dependent) or redox-balanced (NADH-dependent) ketol-acid reductoisomerase enzyme. We then conducted transcriptomic, proteomic and metabolomic analyses to elucidate molecular differences between the engineered strains. Pathway localization had a large effect on isobutanol production with the strain expressing the mitochondrial-localized enzymes producing 3.8-fold more isobutanol than strains expressing the cytosolic enzymes. Cofactor-balancing did not improve isobutanol titers and instead the strain with the redox-imbalanced pathway produced 1.5-fold more isobutanol than the balanced version, albeit at low overall pathway flux. Functional genomic analyses suggested that the poor performances of the cytosolic pathway strains were in part due to a shortage in cytosolic Fe-S clusters, which are required cofactors for the dihydroxyacid dehydratase enzyme. We then demonstrated that this cofactor limitation may be partially recovered by disrupting iron homeostasis with a fra2 mutation, thereby increasing cellular iron levels. The resulting isobutanol titer of the fra2 null strain harboring a cytosolic-localized isobutanol pathway outperformed the strain with the mitochondrial-localized pathway by 1.3-fold, demonstrating that both localizations can support flux to isobutanol.

CRISpy-Pop: A Web Tool for Designing CRISPR/Cas9-Driven Genetic Modifications in Diverse Populations.

CRISPR/Cas9 is a powerful tool for editing genomes, but design decisions are generally made with respect to a single reference genome. With population genomic data becoming available for an increasing number of model organisms, researchers are interested in manipulating multiple strains and lines. CRISpy-pop is a web application that generates and filters guide RNA sequences for CRISPR/Cas9 genome editing for diverse yeast and bacterial strains. The current implementation designs and predicts the activity of guide RNAs against more than 1000 Saccharomyces cerevisiae genomes, including 167 strains frequently used in bioenergy research. Zymomonas mobilis, an increasingly popular bacterial bioenergy research model, is also supported. CRISpy-pop is available as a web application (https://CRISpy-pop.glbrc.org/) with an intuitive graphical user interface. CRISpy-pop also cross-references the human genome to allow users to avoid the selection of guide RNAs with potential biosafety concerns. Additionally, CRISpy-pop predicts the strain coverage of each guide RNA within the supported strain sets, which aids in functional population genetic studies. Finally, we validate how CRISpy-pop can accurately predict the activity of guide RNAs across strains using population genomic data.

Synthetic hybrids of six yeast species.

Allopolyploidy generates diversity by increasing the number of copies and sources of chromosomes. Many of the best-known evolutionary radiations, crops, and industrial organisms are ancient or recent allopolyploids. Allopolyploidy promotes differentiation and facilitates adaptation to new environments, but the tools to test its limits are lacking. Here we develop an iterative method of Hybrid Production (iHyPr) to combine the genomes of multiple budding yeast species, generating Saccharomyces allopolyploids of at least six species. When making synthetic hybrids, chromosomal instability and cell size increase dramatically as additional copies of the genome are added. The six-species hybrids initially grow slowly, but they rapidly regain fitness and adapt, even as they retain traits from multiple species. These new synthetic yeast hybrids and the iHyPr method have potential applications for the study of polyploidy, genome stability, chromosome segregation, and bioenergy.

Natural Variation in the Multidrug Efflux Pump SGE1 Underlies Ionic Liquid Tolerance in Yeast.

Imidazolium ionic liquids (IILs) have a range of biotechnological applications, including as pretreatment solvents that extract cellulose from plant biomass for microbial fermentation into sustainable bioenergy. However, residual levels of IILs, such as 1-ethyl-3-methylimidazolium chloride ([C2C1im]Cl), are toxic to biofuel-producing microbes, including the yeast Saccharomyces cerevisiae. S. cerevisiae strains isolated from diverse ecological niches differ in genomic sequence and in phenotypes potentially beneficial for industrial applications, including tolerance to inhibitory compounds present in hydrolyzed plant feedstocks. We evaluated >100 genome-sequenced S. cerevisiae strains for tolerance to [C2C1im]Cl and identified one strain with exceptional tolerance. By screening a library of genomic DNA fragments from the [C2C1im]Cl-tolerant strain for improved IIL tolerance, we identified SGE1, which encodes a plasma membrane multidrug efflux pump, and a previously uncharacterized gene that we named ionic liquid tolerance 1 (ILT1), which encodes a predicted membrane protein. Analyses of SGE1 sequences from our panel of S. cerevisiae strains together with growth phenotypes implicated two single nucleotide polymorphisms (SNPs) that associated with IIL tolerance and sensitivity. We confirmed these phenotypic effects by transferring the SGE1 SNPs into a [C2C1im]Cl-sensitive yeast strain using CRISPR/Cas9 genome editing. Further studies indicated that these SNPs affect Sge1 protein stability and cell surface localization, influencing the amount of toxic IILs that cells can pump out of the cytoplasm. Our results highlight the general potential for discovering useful biotechnological functions from untapped natural sequence variation and provide functional insight into emergent SGE1 alleles with reduced capacities to protect against IIL toxicity.

Repeated Cis-Regulatory Tuning of a Metabolic Bottleneck Gene during Evolution.

Repeated evolutionary events imply underlying genetic constraints that can make evolutionary mechanisms predictable. Morphological traits are thought to evolve frequently through cis-regulatory changes because these mechanisms bypass constraints in pleiotropic genes that are reused during development. In contrast, the constraints acting on metabolic traits during evolution are less well studied. Here we show how a metabolic bottleneck gene has repeatedly adopted similar cis-regulatory solutions during evolution, likely due to its pleiotropic role integrating flux from multiple metabolic pathways. Specifically, the genes encoding phosphoglucomutase activity (PGM1/PGM2), which connect GALactose catabolism to glycolysis, have gained and lost direct regulation by the transcription factor Gal4 several times during yeast evolution. Through targeted mutations of predicted Gal4-binding sites in yeast genomes, we show this galactose-mediated regulation of PGM1/2 supports vigorous growth on galactose in multiple yeast species, including Saccharomyces uvarum and Lachancea kluyveri. Furthermore, the addition of galactose-inducible PGM1 alone is sufficient to improve the growth on galactose of multiple species that lack this regulation, including Saccharomyces cerevisiae. The strong association between regulation of PGM1/2 by Gal4 even enables remarkably accurate predictions of galactose growth phenotypes between closely related species. This repeated mode of evolution suggests that this specific cis-regulatory connection is a common way that diverse yeasts can govern flux through the pathway, likely due to the constraints imposed by this pleiotropic bottleneck gene. Since metabolic pathways are highly interconnected, we argue that cis-regulatory evolution might be widespread at pleiotropic genes that control metabolic bottlenecks and intersections.

A cell-free method for expressing and reconstituting membrane proteins enables functional characterization of the plant receptor-like protein kinase FERONIA.

There are more than 600 receptor-like kinases (RLKs) in Arabidopsis, but due to challenges associated with the characterization of membrane proteins, only a few have known biological functions. The plant RLK FERONIA is a peptide receptor and has been implicated in plant growth regulation, but little is known about its molecular mechanism of action. To investigate the properties of this enzyme, we used a cell-free wheat germ-based expression system in which mRNA encoding FERONIA was co-expressed with mRNA encoding the membrane scaffold protein variant MSP1D1. With the addition of the lipid cardiolipin, assembly of these proteins into nanodiscs was initiated. FERONIA protein kinase activity in nanodiscs was higher than that of soluble protein and comparable with other heterologously expressed protein kinases. Truncation experiments revealed that the cytoplasmic juxtamembrane domain is necessary for maximal FERONIA activity, whereas the transmembrane domain is inhibitory. An ATP analogue that reacts with lysine residues inhibited catalytic activity and labeled four lysines; mutagenesis demonstrated that two of these, Lys-565 and Lys-663, coordinate ATP in the active site. Mass spectrometric phosphoproteomic measurements further identified phosphorylation sites that were examined using phosphomimetic mutagenesis. The results of these experiments are consistent with a model in which kinase-mediated phosphorylation within the C-terminal region is inhibitory and regulates catalytic activity. These data represent a step further toward understanding the molecular basis for the protein kinase catalytic activity of FERONIA and show promise for future characterization of eukaryotic membrane proteins.

Mitochondrial Protein Interaction Mapping Identifies Regulators of Respiratory Chain Function.

Mitochondria are essential for numerous cellular processes, yet hundreds of their proteins lack robust functional annotation. To reveal functions for these proteins (termed MXPs), we assessed condition-specific protein-protein interactions for 50 select MXPs using affinity enrichment mass spectrometry. Our data connect MXPs to diverse mitochondrial processes, including multiple aspects of respiratory chain function. Building upon these observations, we validated C17orf89 as a complex I (CI) assembly factor. Disruption of C17orf89 markedly reduced CI activity, and its depletion is found in an unresolved case of CI deficiency. We likewise discovered that LYRM5 interacts with and deflavinates the electron-transferring flavoprotein that shuttles electrons to coenzyme Q (CoQ). Finally, we identified a dynamic human CoQ biosynthetic complex involving multiple MXPs whose topology we map using purified components. Collectively, our data lend mechanistic insight into respiratory chain-related activities and prioritize hundreds of additional interactions for further exploration of mitochondrial protein function.

Expression platforms for producing eukaryotic proteins: a comparison of E. coli cell-based and wheat germ cell-free synthesis, affinity and solubility tags, and cloning strategies.

Vectors designed for protein production in Escherichia coli and by wheat germ cell-free translation were tested using 21 well-characterized eukaryotic proteins chosen to serve as controls within the context of a structural genomics pipeline. The controls were carried through cloning, small-scale expression trials, large-scale growth or synthesis, and purification. Successfully purified proteins were also subjected to either crystallization trials or (1)H-(15)N HSQC NMR analyses. Experiments evaluated: (1) the relative efficacy of restriction/ligation and recombinational cloning systems; (2) the value of maltose-binding protein (MBP) as a solubility enhancement tag; (3) the consequences of in vivo proteolysis of the MBP fusion as an alternative to post-purification proteolysis; (4) the effect of the level of LacI repressor on the yields of protein obtained from E. coli using autoinduction; (5) the consequences of removing the His tag from proteins produced by the cell-free system; and (6) the comparative performance of E. coli cells or wheat germ cell-free translation. Optimal promoter/repressor and fusion tag configurations for each expression system are discussed.

Mitochondrial ADCK3 employs an atypical protein kinase-like fold to enable coenzyme Q biosynthesis.

The ancient UbiB protein kinase-like family is involved in isoprenoid lipid biosynthesis and is implicated in human diseases, but demonstration of UbiB kinase activity has remained elusive for unknown reasons. Here, we quantitatively define UbiB-specific sequence motifs and reveal their positions within the crystal structure of a UbiB protein, ADCK3. We find that multiple UbiB-specific features are poised to inhibit protein kinase activity, including an N-terminal domain that occupies the typical substrate binding pocket and a unique A-rich loop that limits ATP binding by establishing an unusual selectivity for ADP. A single alanine-to-glycine mutation of this loop flips this coenzyme selectivity and enables autophosphorylation but inhibits coenzyme Q biosynthesis in vivo, demonstrating functional relevance for this unique feature. Our work provides mechanistic insight into UbiB enzyme activity and establishes a molecular foundation for further investigation of how UbiB family proteins affect diseases and diverse biological pathways.

Mitochondrial COQ9 is a lipid-binding protein that associates with COQ7 to enable coenzyme Q biosynthesis.

Coenzyme Q (CoQ) is an isoprenylated quinone that is essential for cellular respiration and is synthesized in mitochondria by the combined action of at least nine proteins (COQ1-9). Although most COQ proteins are known to catalyze modifications to CoQ precursors, the biochemical role of COQ9 remains unclear. Here, we report that a disease-related COQ9 mutation leads to extensive disruption of the CoQ protein biosynthetic complex in a mouse model, and that COQ9 specifically interacts with COQ7 through a series of conserved residues. Toward understanding how COQ9 can perform these functions, we solved the crystal structure of Homo sapiens COQ9 at 2.4 Å. Unexpectedly, our structure reveals that COQ9 has structural homology to the TFR family of bacterial transcriptional regulators, but that it adopts an atypical TFR dimer orientation and is not predicted to bind DNA. Our structure also reveals a lipid-binding site, and mass spectrometry-based analyses of purified COQ9 demonstrate that it associates with multiple lipid species, including CoQ itself. The conserved COQ9 residues necessary for its interaction with COQ7 comprise a surface patch around the lipid-binding site, suggesting that COQ9 might serve to present its bound lipid to COQ7. Collectively, our data define COQ9 as the first, to our knowledge, mammalian TFR structural homolog and suggest that its lipid-binding capacity and association with COQ7 are key features for enabling CoQ biosynthesis.

The oligomeric states of the purified sigma-1 receptor are stabilized by ligands.

Sigma-1 receptor (S1R) is a mammalian member of the ERG2 and sigma-1 receptor-like protein family (pfam04622). It has been implicated in drug addiction and many human neurological disorders, including Alzheimer and Parkinson diseases and amyotrophic lateral sclerosis. A broad range of synthetic small molecules, including cocaine, (+)-pentazocine, haloperidol, and small endogenous molecules such as N,N-dimethyltryptamine, sphingosine, and steroids, have been identified as regulators of S1R. However, the mechanism of activation of S1R remains obscure. Here, we provide evidence in vitro that S1R has ligand binding activity only in an oligomeric state. The oligomeric state is prone to decay into an apparent monomeric form when exposed to elevated temperature, with loss of ligand binding activity. This decay is suppressed in the presence of the known S1R ligands such as haloperidol, BD-1047, and sphingosine. S1R has a GXXXG motif in its second transmembrane region, and these motifs are often involved in oligomerization of membrane proteins. Disrupting mutations within the GXXXG motif shifted the fraction of the higher oligomeric states toward smaller states and resulted in a significant decrease in specific (+)-[(3)H]pentazocine binding. Results presented here support the proposal that S1R function may be regulated by its oligomeric state. Possible mechanisms of molecular regulation of interacting protein partners by S1R in the presence of small molecule ligands are discussed.

Structural and functional characterization of CalS11, a TDP-rhamnose 3'-O-methyltransferase involved in calicheamicin biosynthesis.

Sugar methyltransferases (MTs) are an important class of tailoring enzymes that catalyze the transfer of a methyl group from S-adenosyl-l-methionine to sugar-based N-, C- and O-nucleophiles. While sugar N- and C-MTs involved in natural product biosynthesis have been found to act on sugar nucleotide substrates prior to a subsequent glycosyltransferase reaction, corresponding sugar O-methylation reactions studied thus far occur after the glycosyltransfer reaction. Herein we report the first in vitro characterization using (1)H-(13)C-gHSQC with isotopically labeled substrates and the X-ray structure determination at 1.55 Å resolution of the TDP-3'-O-rhamnose-methyltransferase CalS11 from Micromonospora echinospora. This study highlights a unique NMR-based methyltransferase assay, implicates CalS11 to be a metal- and general acid/base-dependent O-methyltransferase, and as a first crystal structure for a TDP-hexose-O-methyltransferase, presents a new template for mechanistic studies and/or engineering.

Improved expression and purification of sigma 1 receptor fused to maltose binding protein by alteration of linker sequence.

Sigma 1 receptor (S1R) is a eukaryotic membrane protein that functions as an inter-organelle signaling modulator and chaperone. Here we report an improved expression of S1R in Escherichia coli as a fusion to maltose binding protein (MBP) and a high-yield purification. Variants with linking amino acid sequences consisting of 0-5 alanine residues between MBP and S1R were created and tested in several E. coli expression strains in order to determine the best combination of construct and host for production of active MBP-S1R. Among the linker variations, the protein containing a 4-Ala linker exhibited superior expression characteristics (MBP-4A-S1R); this construct was most productively paired with E. coli B834-pRARE2 and a chemically defined growth and expression medium. A 3-step purification was developed, including extraction from the E. coli membrane fraction using a mixture of Triton X-100 and n-dodecyl-beta-D-maltopyranoside identified by screening constrainted by retention of binding function, and purification by amylose affinity and gel filtration chromatographies. This procedure yields ∼3.5mg of purified fusion protein per L of bacterial culture medium. Purified MBP-4A-S1R showed a 175-fold purification from the starting cellular lysate with respect to specific ligand binding activity, and is stable during concentration and freeze-thaw cycling.

Cell-free translation of integral membrane proteins into unilamelar liposomes.

Wheat germ cell-free translation is shown to be an effective method to produce integral membrane proteins in the presence of unilamelar liposomes. In this chapter, we describe the expression vectors, preparation of mRNA, two types of cell-free translation reactions performed in the presence of liposomes, a simple and highly efficient purification of intact proteoliposomes using density gradient ultracentrifugation, and some of the types of characterization studies that are facilitated by this facile preparative approach. The in vitro transfer of newly translated, membrane proteins into liposomes compatible with direct measurements of their catalytic function is contrasted with existing approaches to extract membrane proteins from biological membranes using detergents and subsequently transfer them back to liposomes for functional studies.

The Center for Eukaryotic Structural Genomics.

The Center for Eukaryotic Structural Genomics (CESG) is a "specialized" or "technology development" center supported by the Protein Structure Initiative (PSI). CESG's mission is to develop improved methods for the high-throughput solution of structures from eukaryotic proteins, with a very strong weighting toward human proteins of biomedical relevance. During the first three years of PSI-2, CESG selected targets representing 601 proteins from Homo sapiens, 33 from mouse, 10 from rat, 139 from Galdieria sulphuraria, 35 from Arabidopsis thaliana, 96 from Cyanidioschyzon merolae, 80 from Plasmodium falciparum, 24 from yeast, and about 25 from other eukaryotes. Notably, 30% of all structures of human proteins solved by the PSI Centers were determined at CESG. Whereas eukaryotic proteins generally are considered to be much more challenging targets than prokaryotic proteins, the technology now in place at CESG yields success rates that are comparable to those of the large production centers that work primarily on prokaryotic proteins. We describe here the technological innovations that underlie CESG's platforms for bioinformatics and laboratory information management, target selection, protein production, and structure determination by X-ray crystallography or NMR spectroscopy.

Flexi vector cloning.

A protocol for ligation-dependent cloning using the Flexi Vector method in a 96-well format is described. The complete protocol includes PCR amplification of the desired gene to append Flexi Vector cloning sequences, restriction digestion of the PCR products, ligation of the digested PCR products into a similarly digested acceptor vector, transformation and growth of host cells, analysis of the transformed clones, and storage of a sequence-verified clone. The protocol also includes transfer of the sequence-verified clones into another Flexi Vector plasmid backbone. Smaller numbers of cloning reactions can be undertaken by appropriate scaling of the indicated reaction volumes.

Structural and functional characterization of a novel phosphatase from the Arabidopsis thaliana gene locus At1g05000.

The crystal structure of the protein product of the gene locus At1g05000, a hypothetical protein from A. thaliana, was determined by the multiple-wavelength anomalous diffraction method and was refined to an R factor of 20.4% (R(free) = 24.9%) at 3.3 A. The protein adopts the alpha/beta fold found in cysteine phosphatases, a superfamily of phosphatases that possess a catalytic cysteine and form a covalent thiol-phosphate intermediate during the catalytic cycle. In At1g05000, the analogous cysteine (Cys(150)) is located at the bottom of a positively-charged pocket formed by residues that include the conserved arginine (Arg(156)) of the signature active site motif, HCxxGxxRT. Of 74 model phosphatase substrates tested, purified recombinant At1g05000 showed highest activity toward polyphosphate (poly-P(12-13)) and deoxyribo- and ribonucleoside triphosphates, and less activity toward phosphoenolpyruvate, phosphotyrosine, phosphotyrosine-containing peptides, and phosphatidyl inositols. Divalent metal cations were not required for activity and had little effect on the reaction.

High efficiency single step production of expression plasmids from cDNA clones using the Flexi Vector cloning system.

The success of structural genomics and proteomics initiatives is dependent on the availability of target genes in vectors suitable for protein production. Here, we compare two high-throughput methods for producing expression vectors from plasmid-derived cDNA fragments. Expression vectors were constructed for compatibility with the Gateway recombination cloning system and the Flexi Vector restriction-based cloning system. Cloning protocols for each system were conducted in parallel for 96 different target genes from PCR through the production of sequence-verified expression clones. The short nucleotide sequences required to prepare the target open reading frames for Flexi Vector cloning allowed a single-step PCR protocol, resulting in fewer mutations relative to the Gateway protocol. Furthermore, through initial cloning of the target open reading frames directly into an expression vector, the Flexi Vector system gave time and cost savings compared to the protocol required for the Gateway system. Within the Flexi Vector system, genes were transferred between four different expression vectors. The efficiency of gene transfer between Flexi Vectors depended on including a region of sequence identity adjacent to one of the restriction sites. With the proper construction in the flanking sequence of the vector, gene transfer efficiencies of 95-98% were demonstrated.

Comparison of cell-based and cell-free protocols for producing target proteins from the Arabidopsis thaliana genome for structural studies.

We describe a comparative study of protein production from 96 Arabidopsis thaliana open reading frames (ORFs) by cell-based and cell-free protocols. Each target was carried through four pipeline protocols used by the Center for Eukaryotic Structural Genomics (CESG), one for the production of unlabeled protein to be used in crystallization trials and three for the production of 15N-labeled proteins to be analyzed by 1H-15N NMR correlation spectroscopy. Two of the protocols involved Escherichia coli cell-based and two involved wheat germ cell-free technology. The progress of each target through each of the protocols was followed with all failures and successes noted. Failures were of the following types: ORF not cloned, protein not expressed, low protein yield, no cleavage of fusion protein, insoluble protein, protein not purified, NMR sample too dilute. Those targets that reached the goal of analysis by 1H-15N NMR correlation spectroscopy were scored as HSQC+ (protein folded and suitable for NMR structural analysis), HSQC+/- (protein partially disordered or not in a single stable conformational state), HSQC- (protein unfolded, misfolded, or aggregated and thus unsuitable for NMR structural analysis). Targets were also scored as X- for failing to crystallize and X+ for successful crystallization. The results constitute a rich database for understanding differences between targets and protocols. In general, the wheat germ cell-free platform offers the advantage of greater genome coverage for NMR-based structural proteomics whereas the E. coli platform when successful yields more protein, as currently needed for crystallization trials for X-ray structure determination.