Dehydrated Biomimetic Membranes with Skinlike Structure and Function.

Novel vapor-permeable materials are sought after for applications in protective wear, energy generation, and water treatment. Current impermeable protective materials effectively block harmful agents but trap heat due to poor water vapor transfer. Here we present a new class of materials, vapor permeable dehydrated nanoporous biomimetic membranes (DBMs), based on channel proteins. This application for biomimetic membranes is unexpected as channel proteins and biomimetic membranes were assumed to be unstable under dry conditions. DBMs mimic human skin's structure to offer both high vapor transport and small molecule exclusion under dry conditions. DBMs feature highly organized pores resembling sweat pores in human skin, but at super high densities (>1012 pores/cm2). These DBMs achieved exceptional water vapor transport rates, surpassing commercial breathable fabrics by up to 6.2 times, despite containing >2 orders of magnitude smaller pores (1 nm vs >700 nm). These DBMs effectively excluded model biological agents and harmful chemicals both in liquid and vapor phases, again in contrast with the commercial breathable fabrics. Remarkably, while hydrated biomimetic membranes were highly permeable to liquid water, they exhibited higher water resistances after dehydration at values >38 times that of commercial breathable fabrics. Molecular dynamics simulations support our hypothesis that dehydration induced protein hydrophobicity increases which enhanced DBM performance. DBMs hold promise for various applications, including membrane distillation, dehumidification, and protective barriers for atmospheric water harvesting materials.

Mitochondrial ATP generation is more proteome efficient than glycolysis.

Metabolic efficiency profoundly influences organismal fitness. Nonphotosynthetic organisms, from yeast to mammals, derive usable energy primarily through glycolysis and respiration. Although respiration is more energy efficient, some cells favor glycolysis even when oxygen is available (aerobic glycolysis, Warburg effect). A leading explanation is that glycolysis is more efficient in terms of ATP production per unit mass of protein (that is, faster). Through quantitative flux analysis and proteomics, we find, however, that mitochondrial respiration is actually more proteome efficient than aerobic glycolysis. This is shown across yeast strains, T cells, cancer cells, and tissues and tumors in vivo. Instead of aerobic glycolysis being valuable for fast ATP production, it correlates with high glycolytic protein expression, which promotes hypoxic growth. Aerobic glycolytic yeasts do not excel at aerobic growth but outgrow respiratory cells during oxygen limitation. We accordingly propose that aerobic glycolysis emerges from cells maintaining a proteome conducive to both aerobic and hypoxic growth.

Genome streamlining to improve performance of a fast-growing cyanobacterium Synechococcus elongatus UTEX 2973.

Cyanobacteria are photosynthetic organisms that have garnered significant recognition as potential hosts for sustainable bioproduction. However, their complex regulatory networks pose significant challenges to major metabolic engineering efforts, thereby limiting their feasibility as production hosts. Genome streamlining has been demonstrated to be a successful approach for improving productivity and fitness in heterotrophs but is yet to be explored to its full potential in phototrophs. Here, we present the systematic reduction of the genome of the cyanobacterium exhibiting the fastest exponential growth, Synechococcus elongatus UTEX 2973. This work, the first of its kind in a photoautotroph, involved an iterative process using state-of-the-art genome-editing technology guided by experimental analysis and computational tools. CRISPR-Cas3 enabled large, progressive deletions of predicted dispensable regions and aided in the identification of essential genes. The large deletions were combined to obtain a strain with 55-kb genome reduction. The strains with streamlined genome showed improvement in growth (up to 23%) and productivity (by 22.7%) as compared to the wild type (WT). This streamlining strategy not only has the potential to develop cyanobacterial strains with improved growth and productivity traits but can also facilitate a better understanding of their genome-to-phenome relationships.IMPORTANCEGenome streamlining is an evolutionary strategy used by natural living systems to dispense unnecessary genes from their genome as a mechanism to adapt and evolve. While this strategy has been successfully borrowed to develop synthetic heterotrophic microbial systems with desired phenotype, it has not been extensively explored in photoautotrophs. Genome streamlining strategy incorporates both computational predictions to identify the dispensable regions and experimental validation using genome-editing tool, and in this study, we have employed a modified strategy with the goal to minimize the genome size to an extent that allows optimal cellular fitness under specified conditions. Our strategy has explored a novel genome-editing tool in photoautotrophs, which, unlike other existing tools, enables large, spontaneous optimal deletions from the genome. Our findings demonstrate the effectiveness of this modified strategy in obtaining strains with streamlined genome, exhibiting improved fitness and productivity.

KETCHUP: Parameterizing of large-scale kinetic models using multiple datasets with different reference states.

Large-scale kinetic models provide the computational means to dynamically link metabolic reaction fluxes to metabolite concentrations and enzyme levels while also conforming to substrate level regulation. However, the development of broadly applicable frameworks for efficiently and robustly parameterizing models remains a challenge. Challenges arise due to both the heterogeneity, paucity, and difficulty in obtaining flux and/or concentration data but also due to the computational difficulties of the underlying parameter identification problem. Both the computational demands for parameterization, degeneracy of obtained parameter solutions and interpretability of results has so far limited widespread adoption of large-scale kinetic models despite their potential. Herein, we introduce the Kinetic Estimation Tool Capturing Heterogeneous Datasets Using Pyomo (KETCHUP), a flexible parameter estimation tool that leverages a primal-dual interior-point algorithm to solve a nonlinear programming (NLP) problem that identifies a set of parameters capable of recapitulating the (non)steady-state fluxes and concentrations in wild-type and perturbed metabolic networks. KETCHUP is benchmarked against previously parameterized large-scale kinetic models demonstrating an at least an order of magnitude faster convergence than the tool K-FIT while at the same time attaining better data fits. This versatile toolbox accepts different kinetic descriptions, metabolic fluxes, enzyme levels and metabolite concentrations, under either steady-state or instationary conditions to enable robust kinetic model construction and parameterization. KETCHUP supports the SBML format and can be accessed at https://github.com/maranasgroup/KETCHUP.

A multi-organ maize metabolic model connects temperature stress with energy production and reducing power generation.

Climate change has adversely affected maize productivity. Thereby, a holistic understanding of metabolic crosstalk among its organs is important to address this issue. Thus, we reconstructed the first multi-organ maize metabolic model, iZMA6517, and contextualized it with heat and cold stress transcriptomics data using expression distributed reaction flux measurement (EXTREAM) algorithm. Furthermore, implementing metabolic bottleneck analysis on contextualized models revealed differences between these stresses. While both stresses had reducing power bottlenecks, heat stress had additional energy generation bottlenecks. We also performed thermodynamic driving force analysis, revealing thermodynamics-reducing power-energy generation axis dictating the nature of temperature stress responses. Thus, a temperature-tolerant maize ideotype can be engineered by leveraging the proposed thermodynamics-reducing power-energy generation axis. We experimentally inoculated maize root with a beneficial mycorrhizal fungus, Rhizophagus irregularis, and as a proof-of-concept demonstrated its efficacy in alleviating temperature stress. Overall, this study will guide the engineering effort of temperature stress-tolerant maize ideotypes.

Deuterated water as a substrate-agnostic isotope tracer for investigating reversibility and thermodynamics of reactions in central carbon metabolism.

Stable isotope tracers are a powerful tool for the quantitative analysis of microbial metabolism, enabling pathway elucidation, metabolic flux quantification, and assessment of reaction and pathway thermodynamics. 13C and 2H metabolic flux analysis commonly relies on isotopically labeled carbon substrates, such as glucose. However, the use of 2H-labeled nutrient substrates faces limitations due to their high cost and limited availability in comparison to 13C-tracers. Furthermore, isotope tracer studies in industrially relevant bacteria that metabolize complex substrates such as cellulose, hemicellulose, or lignocellulosic biomass, are challenging given the difficulty in obtaining these as isotopically labeled substrates. In this study, we examine the potential of deuterated water (2H2O) as an affordable, substrate-neutral isotope tracer for studying central carbon metabolism. We apply 2H2O labeling to investigate the reversibility of glycolytic reactions across three industrially relevant bacterial species -C. thermocellum, Z. mobilis, and E. coli-harboring distinct glycolytic pathways with unique thermodynamics. We demonstrate that 2H2O labeling recapitulates previous reversibility and thermodynamic findings obtained with established 13C and 2H labeled nutrient substrates. Furthermore, we exemplify the utility of this 2H2O labeling approach by applying it to high-substrate C. thermocellum fermentations -a setting in which the use of conventional tracers is impractical-thereby identifying the glycolytic enzyme phosphofructokinase as a major bottleneck during high-substrate fermentations and unveiling critical insights that will steer future engineering efforts to enhance ethanol production in this cellulolytic organism. This study demonstrates the utility of deuterated water as a substrate-agnostic isotope tracer for examining flux and reversibility of central carbon metabolic reactions, which yields biological insights comparable to those obtained using costly 2H-labeled nutrient substrates.

A single amino acid change led to structural and functional differentiation of PvHd1 to control flowering in switchgrass.

Switchgrass, a forage and bioenergy crop, occurs as two main ecotypes with different but overlapping ranges of adaptation. The two ecotypes differ in a range of characteristics, including flowering time. Flowering time determines the duration of vegetative development and therefore biomass accumulation, a key trait in bioenergy crops. No causal variants for flowering time differences between switchgrass ecotypes have, as yet, been identified. In this study, we mapped a robust flowering time quantitative trait locus (QTL) on chromosome 4K in a biparental F2 population and characterized the flowering-associated transcription factor gene PvHd1, an ortholog of CONSTANS in Arabidopsis and Heading date 1 in rice, as the underlying causal gene. Protein modeling predicted that a serine to glycine substitution at position 35 (p.S35G) in B-Box domain 1 greatly altered the global structure of the PvHd1 protein. The predicted variation in protein compactness was supported in vitro by a 4 °C shift in denaturation temperature. Overexpressing the PvHd1-p.35S allele in a late-flowering CONSTANS-null Arabidopsis mutant rescued earlier flowering, whereas PvHd1-p.35G had a reduced ability to promote flowering, demonstrating that the structural variation led to functional divergence. Our findings provide us with a tool to manipulate the timing of floral transition in switchgrass cultivars and, potentially, expand their cultivation range.

Rank-ordering of known enzymes as starting points for re-engineering novel substrate activity using a convolutional neural network.

Retro-biosynthetic approaches have made significant advances in predicting synthesis routes of target biofuel, bio-renewable or bio-active molecules. The use of only cataloged enzymatic activities limits the discovery of new production routes. Recent retro-biosynthetic algorithms increasingly use novel conversions that require altering the substrate or cofactor specificities of existing enzymes while connecting pathways leading to a target metabolite. However, identifying and re-engineering enzymes for desired novel conversions are currently the bottlenecks in implementing such designed pathways. Herein, we present EnzRank, a convolutional neural network (CNN) based approach, to rank-order existing enzymes in terms of their suitability to undergo successful protein engineering through directed evolution or de novo design towards a desired specific substrate activity. We train the CNN model on 11,800 known active enzyme-substrate pairs from the BRENDA database as positive samples and data generated by scrambling these pairs as negative samples using substrate dissimilarity between an enzyme's native substrate and all other molecules present in the dataset using Tanimoto similarity score. EnzRank achieves an average recovery rate of 80.72% and 73.08% for positive and negative pairs on test data after using a 10-fold holdout method for training and cross-validation. We further developed a web-based user interface (available at https://huggingface.co/spaces/vuu10/EnzRank) to predict enzyme-substrate activity using SMILES strings of substrates and enzyme sequence as input to allow convenient and easy-to-use access to EnzRank. In summary, this effort can aid de novo pathway design tools to prioritize starting enzyme re-engineering candidates for novel reactions as well as in predicting the potential secondary activity of enzymes in cell metabolism.

Little Brown Bats (Myotis lucifugus) Support the Binding of SARS-CoV-2 Spike and Are Likely Susceptible to SARS-CoV-2 Infection.

Severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2), believed to have originated from a bat species, can infect a wide range of non-human hosts. Bats are known to harbor hundreds of coronaviruses capable of spillover into human populations. Recent studies have shown a significant variation in the susceptibility among bat species to SARS-CoV-2 infection. We show that little brown bats (LBB) express angiotensin-converting enzyme 2 receptor and the transmembrane serine protease 2, which are accessible to and support SARS-CoV-2 binding. All-atom molecular dynamics (MD) simulations revealed that LBB ACE2 formed strong electrostatic interactions with the RBD similar to human and cat ACE2 proteins. In summary, LBBs, a widely distributed North American bat species, could be at risk of SARS-CoV-2 infection and potentially serve as a natural reservoir. Finally, our framework, combining in vitro and in silico methods, is a useful tool to assess the SARS-CoV-2 susceptibility of bats and other animal species.

Evaluating proteome allocation of Saccharomyces cerevisiae phenotypes with resource balance analysis.

Saccharomyces cerevisiae is an important model organism and a workhorse in bioproduction. Here, we reconstructed a compact and tractable genome-scale resource balance analysis (RBA) model (i.e., named scRBA) to analyze metabolic fluxes and proteome allocation in a computationally efficient manner. Resource capacity models such as scRBA provide the quantitative means to identify bottlenecks in biosynthetic pathways due to enzyme, compartment size, and/or ribosome availability limitations. ATP maintenance rate and in vivo apparent turnover numbers (kapp) were regressed from metabolic flux and protein concentration data to capture observed physiological growth yield and proteome efficiency and allocation, respectively. Estimated parameter values were found to vary with oxygen and nutrient availability. Overall, this work (i) provides condition-specific model parameters to recapitulate phenotypes corresponding to different extracellular environments, (ii) alludes to the enhancing effect of substrate channeling and post-translational activation on in vivo enzyme efficiency in glycolysis and electron transport chain, and (iii) reveals that the Crabtree effect is underpinned by specific limitations in mitochondrial proteome capacity and secondarily ribosome availability rather than overall proteome capacity.

A detailed genome-scale metabolic model of Clostridium thermocellum investigates sources of pyrophosphate for driving glycolysis.

Lignocellulosic biomass is an abundant and renewable source of carbon for chemical manufacturing, yet it is cumbersome in conventional processes. A promising, and increasingly studied, candidate for lignocellulose bioprocessing is the thermophilic anaerobe Clostridium thermocellum given its potential to produce ethanol, organic acids, and hydrogen gas from lignocellulosic biomass under high substrate loading. Possessing an atypical glycolytic pathway which substitutes GTP or pyrophosphate (PPi) for ATP in some steps, including in the energy-investment phase, identification, and manipulation of PPi sources are key to engineering its metabolism. Previous efforts to identify the primary pyrophosphate have been unsuccessful. Here, we explore pyrophosphate metabolism through reconstructing, updating, and analyzing a new genome-scale stoichiometric model for C. thermocellum, iCTH669. Hundreds of changes to the former GEM, iCBI655, including correcting cofactor usages, addressing charge and elemental balance, standardizing biomass composition, and incorporating the latest experimental evidence led to a MEMOTE score improvement to 94%. We found agreement of iCTH669 model predictions across all available fermentation and biomass yield datasets. The feasibility of hundreds of PPi synthesis routes, newly identified and previously proposed, were assessed through the lens of the iCTH669 model including biomass synthesis, tRNA synthesis, newly identified sources, and previously proposed PPi-generating cycles. In all cases, the metabolic cost of PPi synthesis is at best equivalent to investment of one ATP suggesting no direct energetic advantage for the cofactor substitution in C. thermocellum. Even though no unique source of PPi could be gleaned by the model, by combining with gene expression data two most likely scenarios emerge. First, previously investigated PPi sources likely account for most PPi production in wild-type strains. Second, alternate metabolic routes as encoded by iCTH669 can collectively maintain PPi levels even when previously investigated synthesis cycles are disrupted. Model iCTH669 is available at github.com/maranasgroup/iCTH669.

Metabolic engineering of low-pH-tolerant non-model yeast, Issatchenkia orientalis, for production of citramalate.

Methyl methacrylate (MMA) is an important petrochemical with many applications. However, its manufacture has a large environmental footprint. Combined biological and chemical synthesis (semisynthesis) may be a promising alternative to reduce both cost and environmental impact, but strains that can produce the MMA precursor (citramalate) at low pH are required. A non-conventional yeast, Issatchenkia orientalis, may prove ideal, as it can survive extremely low pH. Here, we demonstrate the engineering of I. orientalis for citramalate production. Using sequence similarity network analysis and subsequent DNA synthesis, we selected a more active citramalate synthase gene (cimA) variant for expression in I. orientalis. We then adapted a piggyBac transposon system for I. orientalis that allowed us to simultaneously explore the effects of different cimA gene copy numbers and integration locations. A batch fermentation showed the genome-integrated-cimA strains produced 2.0 g/L citramalate in 48 h and a yield of up to 7% mol citramalate/mol consumed glucose. These results demonstrate the potential of I. orientalis as a chassis for citramalate production.

Comparative study of two Saccharomyces cerevisiae strains with kinetic models at genome-scale.

The parameterization of kinetic models requires measurement of fluxes and/or metabolite levels for a base strain and a few genetic perturbations thereof. Unlike stoichiometric models that are mostly invariant to the specific strain, it remains unclear whether kinetic models constructed for different strains of the same species have similar or significantly different kinetic parameters. This important question underpins the applicability range and prediction limits of kinetic reconstructions. To this end, herein we parameterize two separate large-scale kinetic models using K-FIT with genome-wide coverage corresponding to two distinct strains of Saccharomyces cerevisiae: CEN.PK 113-7D strain (model k-sacce306-CENPK), and growth-deficient BY4741 (isogenic to S288c; model k-sacce306-BY4741). The metabolic network for each model contains 306 reactions, 230 metabolites, and 119 substrate-level regulatory interactions. The two models (for CEN.PK and BY4741) recapitulate, within one standard deviation, 77% and 75% of the fitted dataset fluxes, respectively, determined by 13C metabolic flux analysis for wild-type and eight single-gene knockout mutants of each strain. Strain-specific kinetic parameterization results indicate that key enzymes in the TCA cycle, glycolysis, and arginine and proline metabolism drive the metabolic differences between these two strains of S. cerevisiae. Our results suggest that although kinetic models cannot be readily used across strains as stoichiometric models, they can capture species-specific information through the kinetic parameterization process.

De novo design and Rosetta-based assessment of high-affinity antibody variable regions (Fv) against the SARS-CoV-2 spike receptor binding domain (RBD).

The continued emergence of new SARS-CoV-2 variants has accentuated the growing need for fast and reliable methods for the design of potentially neutralizing antibodies (Abs) to counter immune evasion by the virus. Here, we report on the de novo computational design of high-affinity Ab variable regions (Fv) through the recombination of VDJ genes targeting the most solvent-exposed hACE2-binding residues of the SARS-CoV-2 spike receptor binding domain (RBD) protein using the software tool OptMAVEn-2.0. Subsequently, we carried out computational affinity maturation of the designed variable regions through amino acid substitutions for improved binding with the target epitope. Immunogenicity of designs was restricted by preferring designs that match sequences from a 9-mer library of "human Abs" based on a human string content score. We generated 106 different antibody designs and reported in detail on the top five that trade-off the greatest computational binding affinity for the RBD with human string content scores. We further describe computational evaluation of the top five designs produced by OptMAVEn-2.0 using a Rosetta-based approach. We used Rosetta SnugDock for local docking of the designs to evaluate their potential to bind the spike RBD and performed "forward folding" with DeepAb to assess their potential to fold into the designed structures. Ultimately, our results identified one designed Ab variable region, P1.D1, as a particularly promising candidate for experimental testing. This effort puts forth a computational workflow for the de novo design and evaluation of Abs that can quickly be adapted to target spike epitopes of emerging SARS-CoV-2 variants or other antigenic targets.

SARS-CoV-2 Omicron (B.1.1.529) Infection of Wild White-Tailed Deer in New York City.

There is mounting evidence of SARS-CoV-2 spillover from humans into many domestic, companion, and wild animal species. Research indicates that humans have infected white-tailed deer, and that deer-to-deer transmission has occurred, indicating that deer could be a wildlife reservoir and a source of novel SARS-CoV-2 variants. We examined the hypothesis that the Omicron variant is actively and asymptomatically infecting the free-ranging deer of New York City. Between December 2021 and February 2022, 155 deer on Staten Island, New York, were anesthetized and examined for gross abnormalities and illnesses. Paired nasopharyngeal swabs and blood samples were collected and analyzed for the presence of SARS-CoV-2 RNA and antibodies. Of 135 serum samples, 19 (14.1%) indicated SARS-CoV-2 exposure, and 11 reacted most strongly to the wild-type B.1 lineage. Of the 71 swabs, 8 were positive for SARS-CoV-2 RNA (4 Omicron and 4 Delta). Two of the animals had active infections and robust neutralizing antibodies, revealing evidence of reinfection or early seroconversion in deer. Variants of concern continue to circulate among and may reinfect US deer populations, and establish enzootic transmission cycles in the wild: this warrants a coordinated One Health response, to proactively surveil, identify, and curtail variants of concern before they can spill back into humans.

Development and Validation of Indirect Enzyme-Linked Immunosorbent Assays for Detecting Antibodies to SARS-CoV-2 in Cattle, Swine, and Chicken.

Multiple domestic and wild animal species are susceptible to SARS-CoV-2 infection. Cattle and swine are susceptible to experimental SARS-CoV-2 infection. The unchecked transmission of SARS-CoV-2 in animal hosts could lead to virus adaptation and the emergence of novel variants. In addition, the spillover and subsequent adaptation of SARS-CoV-2 in livestock could significantly impact food security as well as animal and public health. Therefore, it is essential to monitor livestock species for SARS-CoV-2 spillover. We developed and optimized species-specific indirect ELISAs (iELISAs) to detect anti-SARS-CoV-2 antibodies in cattle, swine, and chickens using the spike protein receptor-binding domain (RBD) antigen. Serum samples collected prior to the COVID-19 pandemic were used to determine the cut-off threshold. RBD hyperimmunized sera from cattle (n = 3), swine (n = 6), and chicken (n = 3) were used as the positive controls. The iELISAs were evaluated compared to a live virus neutralization test using cattle (n = 150), swine (n = 150), and chicken (n = 150) serum samples collected during the COVID-19 pandemic. The iELISAs for cattle, swine, and chicken were found to have 100% sensitivity and specificity. These tools facilitate the surveillance that is necessary to quickly identify spillovers into the three most important agricultural species worldwide.

Examining organic acid production potential and growth-coupled strategies in Issatchenkia orientalis using constraint-based modeling.

Growth-coupling product formation can facilitate strain stability by aligning industrial objectives with biological fitness. Organic acids make up many building block chemicals that can be produced from sugars obtainable from renewable biomass. Issatchenkia orientalis is a yeast strain tolerant to acidic conditions and is thus a promising host for industrial production of organic acids. Here, we use constraint-based methods to assess the potential of computationally designing growth-coupled production strains for I. orientalis that produce 22 different organic acids under aerobic or microaerobic conditions. We explore native and engineered pathways using glucose or xylose as the carbon substrates as proxy constituents of hydrolyzed biomass. We identified growth-coupled production strategies for 37 of the substrate-product pairs, with 15 pairs achieving production for any growth rate. We systematically assess the strain design solutions and categorize the underlying principles involved.

A CNN model for predicting binding affinity changes between SARS-CoV-2 spike RBD variants and ACE2 homologues.

The cellular entry of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) involves the association of its receptor binding domain (RBD) with human angiotensin converting enzyme 2 (hACE2) as the first crucial step. Efficient and reliable prediction of RBD-hACE2 binding affinity changes upon amino acid substitutions can be valuable for public health surveillance and monitoring potential spillover and adaptation into non-human species. Here, we introduce a convolutional neural network (CNN) model trained on protein sequence and structural features to predict experimental RBD-hACE2 binding affinities of 8,440 variants upon single and multiple amino acid substitutions in the RBD or ACE2. The model achieves a classification accuracy of 83.28% and a Pearson correlation coefficient of 0.85 between predicted and experimentally calculated binding affinities in five-fold cross-validation tests and predicts improved binding affinity for most circulating variants. We pro-actively used the CNN model to exhaustively screen for novel RBD variants with combinations of up to four single amino acid substitutions and suggested candidates with the highest improvements in RBD-ACE2 binding affinity for human and animal ACE2 receptors. We found that the binding affinity of RBD variants against animal ACE2s follows similar trends as those against human ACE2. White-tailed deer ACE2 binds to RBD almost as tightly as human ACE2 while cattle, pig, and chicken ACE2s bind weakly. The model allows testing whether adaptation of the virus for increased binding with other animals would cause concomitant increases in binding with hACE2 or decreased fitness due to adaptation to other hosts.

Developmental changes in lignin composition are driven by both monolignol supply and laccase specificity.

The factors controlling lignin composition remain unclear. Catechyl (C)-lignin is a homopolymer of caffeyl alcohol with unique properties as a biomaterial and precursor of industrial chemicals. The lignin synthesized in the seed coat of Cleome hassleriana switches from guaiacyl (G)- to C-lignin at around 12 to 14 days after pollination (DAP), associated with a rerouting of the monolignol pathway. Lack of synthesis of caffeyl alcohol limits C-lignin formation before around 12 DAP, but coniferyl alcohol is still synthesized and highly accumulated after 14 DAP. We propose a model in which, during C-lignin biosynthesis, caffeyl alcohol noncompetitively inhibits oxidation of coniferyl alcohol by cell wall laccases, a process that might limit movement of coniferyl alcohol to the apoplast. Developmental changes in both substrate availability and laccase specificity together account for the metabolic fates of G- and C-monolignols in the Cleome seed coat.

Metabolic engineering of Rhodotorula toruloides IFO0880 improves C16 and C18 fatty alcohol production from synthetic media.

BACKGROUND: The oleaginous, carotenogenic yeast Rhodotorula toruloides has been increasingly explored as a platform organism for the production of terpenoids and fatty acid derivatives. Fatty alcohols, a fatty acid derivative widely used in the production of detergents and surfactants, can be produced microbially with the expression of a heterologous fatty acyl-CoA reductase. Due to its high lipid production, R. toruloides has high potential for fatty alcohol production, and in this study several metabolic engineering approaches were investigated to improve the titer of this product. RESULTS: Fatty acyl-CoA reductase from Marinobacter aqueolei was co-expressed with SpCas9 in R. toruloides IFO0880 and a panel of gene overexpressions and Cas9-mediated gene deletions were explored to increase the fatty alcohol production. Two overexpression targets (ACL1 and ACC1, improving cytosolic acetyl-CoA and malonyl-CoA production, respectively) and two deletion targets (the acyltransferases DGA1 and LRO1) resulted in significant (1.8 to 4.4-fold) increases to the fatty alcohol titer in culture tubes. Combinatorial exploration of these modifications in bioreactor fermentation culminated in a 3.7 g/L fatty alcohol titer in the LRO1Δ mutant. As LRO1 deletion was not found to be beneficial for fatty alcohol production in other yeasts, a lipidomic comparison of the DGA1 and LRO1 knockout mutants was performed, finding that DGA1 is the primary acyltransferase responsible for triacylglyceride production in R. toruloides, while LRO1 disruption simultaneously improved fatty alcohol production, increased diacylglyceride and triacylglyceride production, and increased glucose consumption. CONCLUSIONS: The fatty alcohol titer of fatty acyl-CoA reductase-expressing R. toruloides was significantly improved through the deletion of LRO1, or the deletion of DGA1 combined with overexpression of ACC1 and ACL1. Disruption of LRO1 surprisingly increased both lipid and fatty alcohol production, creating a possible avenue for future study of the lipid metabolism of this yeast.