Structure and interactions of HIV-1 gp41 CHR-NHR reverse hairpin constructs reveal molecular determinants of antiviral activity.

Engineered reverse hairpin constructs containing a partial C-heptad repeat (CHR) sequence followed by a short loop and full-length N-heptad repeat (NHR) were previously shown to form trimers in solution and to be nanomolar inhibitors of HIV-1 Env mediated fusion. Their target is the in situ gp41 fusion intermediate, and they have similar potency to other previously reported NHR trimers. However, their design implies that the NHR is partially covered by CHR, which would be expected to limit potency. An exposed hydrophobic pocket in the folded structure may be sufficient to confer the observed potency, or they may exist in a partially unfolded state exposing full length NHR. Here we examined their structure by crystallography, CD and fluorescence, establishing that the proteins are folded hairpins both in crystal form and in solution. We examined unfolding in the milieu of the fusion reaction by conducting experiments in the presence of a membrane mimetic solvent and by engineering a disulfide bond into the structure to prevent partial unfolding. We further examined the role of the hydrophobic pocket, using a hairpin-small molecule adduct that occluded the pocket, as confirmed by X-ray footprinting. The results demonstrated that the NHR region nominally covered by CHR in the engineered constructs and the hydrophobic pocket region that is exposed by design were both essential for nanomolar potency and that interaction with membrane is likely to play a role in promoting the required inhibitor structure. The design concepts can be applied to other Class 1 viral fusion proteins.

Systematic engineering for production of anti-aging sunscreen compound in Pseudomonas putida.

Sunscreen has been used for thousands of years to protect skin from ultraviolet radiation. However, the use of modern commercial sunscreen containing oxybenzone, ZnO, and TiO2 has raised concerns due to their negative effects on human health and the environment. In this study, we aim to establish an efficient microbial platform for production of shinorine, a UV light absorbing compound with anti-aging properties. First, we methodically selected an appropriate host for shinorine production by analyzing central carbon flux distribution data from prior studies alongside predictions from genome-scale metabolic models (GEMs). We enhanced shinorine productivity through CRISPRi-mediated downregulation and utilized shotgun proteomics to pinpoint potential competing pathways. Simultaneously, we improved the shinorine biosynthetic pathway by refining its design, optimizing promoter usage, and altering the strength of ribosome binding sites. Finally, we conducted amino acid feeding experiments under various conditions to identify the key limiting factors in shinorine production. The study combines meta-analysis of 13C-metabolic flux analysis, GEMs, synthetic biology, CRISPRi-mediated gene downregulation, and omics analysis to improve shinorine production, demonstrating the potential of Pseudomonas putida KT2440 as platform for shinorine production.

Expression in poplar of dehydroshikimate dehydratase induces transcriptional and metabolic changes in the phenylpropanoid pathway.

Modification of lignin in feedstocks via genetic engineering aims to reduce biomass recalcitrance to facilitate efficient conversion processes. These improvements can be achieved by expressing exogenous enzymes that interfere with native biosynthetic pathways responsible for the production of the lignin precursors. In-planta expression of a 3-dehydroshikimate dehydratase (QsuB) in poplar trees reduced lignin content and altered their monomer composition, which enabled higher yields of sugars after cell wall polysaccharide hydrolysis. Understanding how plants respond to such genetic modifications at the transcriptional and metabolic levels is needed to facilitate further improvement and field deployment. In this work, we amassed fundamental knowledge on lignin-modified QsuB poplar using RNA-seq and metabolomics. The data clearly demonstrate that changes in gene expression and metabolite abundance can occur in a strict spatiotemporal fashion, revealing tissue-specific responses in the xylem, phloem, or periderm. In the poplar line that exhibits the strongest reduction in lignin, we found that 3% of the transcripts had altered expression levels and ~19% of the detected metabolites had differential abundance in the xylem from older stems. Changes affect predominantly the shikimate and phenylpropanoid pathways as wells as secondary cell wall metabolism, and result in significant accumulation of hydroxybenzoates derived from protocatechuate and salicylate.

Complete biosynthesis of QS-21 in engineered yeast.

QS-21 is a potent vaccine adjuvant and remains the only saponin-based adjuvant that has been clinically approved for use in humans1,2. However, owing to the complex structure of QS-21, its availability is limited. Today, the supply depends on laborious extraction from the Chilean soapbark tree or on low-yielding total chemical synthesis3,4. Here we demonstrate the complete biosynthesis of QS-21 and its precursors, as well as structural derivatives, in engineered yeast strains. The successful biosynthesis in yeast requires fine-tuning of the host's native pathway fluxes, as well as the functional and balanced expression of 38 heterologous enzymes. The required biosynthetic pathway spans seven enzyme families-a terpene synthase, P450s, nucleotide sugar synthases, glycosyltransferases, a coenzyme A ligase, acyl transferases and polyketide synthases-from six organisms, and mimics in yeast the subcellular compartmentalization of plants from the endoplasmic reticulum membrane to the cytosol. Finally, by taking advantage of the promiscuity of certain pathway enzymes, we produced structural analogues of QS-21 using this biosynthetic platform. This microbial production scheme will allow for the future establishment of a structure-activity relationship, and will thus enable the rational design of potent vaccine adjuvants.

Comparative Pore Structure and Dynamics for Bacterial Microcompartment Shell Protein Assemblies in Sheets or Shells.

Bacterial microcompartments (BMCs) are protein-bound organelles found in some bacteria which encapsulate enzymes for enhanced catalytic activity. These compartments spatially sequester enzymes within semi-permeable shell proteins, analogous to many membrane-bound organelles. The shell proteins assemble into multimeric tiles; hexamers, trimers, and pentamers, and these tiles self-assemble into larger assemblies with icosahedral symmetry. While icosahedral shells are the predominant form in vivo, the tiles can also form nanoscale cylinders or sheets. The individual multimeric tiles feature central pores that are key to regulating transport across the protein shell. Our primary interest is to quantify pore shape changes in response to alternative component morphologies at the nanoscale. We use molecular modeling tools to develop atomically detailed models for both planar sheets of tiles and curved structures representative of the complete shells found in vivo. Subsequently, these models were animated using classical molecular dynamics simulations. From the resulting trajectories, we analyzed overall structural stability, water accessibility to individual residues, water residence time, and pore geometry for the hexameric and trimeric protein tiles from the Haliangium ochraceum model BMC shell. These exhaustive analyses suggest no substantial variation in pore structure or solvent accessibility between the flat and curved shell geometries. We additionally compare our analysis to hydroxyl radical footprinting data to serve as a check against our simulation results, highlighting specific residues where water molecules are bound for a long time. Although with little variation in morphology or water interaction, we propose that the planar and capsular morphology can be used interchangeably when studying permeability through BMC pores.

Edible mycelium bioengineered for enhanced nutritional value and sensory appeal using a modular synthetic biology toolkit.

Filamentous fungi are critical in the transition to a more sustainable food system. While genetic modification of these organisms has promise for enhancing the nutritional value, sensory appeal, and scalability of fungal foods, genetic tools and demonstrated use cases for bioengineered food production by edible strains are lacking. Here, we develop a modular synthetic biology toolkit for Aspergillus oryzae, an edible fungus used in fermented foods, protein production, and meat alternatives. Our toolkit includes a CRISPR-Cas9 method for gene integration, neutral loci, and tunable promoters. We use these tools to elevate intracellular levels of the nutraceutical ergothioneine and the flavor-and color molecule heme in the edible biomass. The strain overproducing heme is red in color and is readily formulated into imitation meat patties with minimal processing. These findings highlight the promise of synthetic biology to enhance fungal foods and provide useful genetic tools for applications in food production and beyond.

Genome-scale and pathway engineering for the sustainable aviation fuel precursor isoprenol production in Pseudomonas putida.

Sustainable aviation fuel (SAF) will significantly impact global warming in the aviation sector, and important SAF targets are emerging. Isoprenol is a precursor for a promising SAF compound DMCO (1,4-dimethylcyclooctane) and has been produced in several engineered microorganisms. Recently, Pseudomonas putida has gained interest as a future host for isoprenol bioproduction as it can utilize carbon sources from inexpensive plant biomass. Here, we engineer metabolically versatile host P. putida for isoprenol production. We employ two computational modeling approaches (Bilevel optimization and Constrained Minimal Cut Sets) to predict gene knockout targets and optimize the "IPP-bypass" pathway in P. putida to maximize isoprenol production. Altogether, the highest isoprenol production titer from P. putida was achieved at 3.5 g/L under fed-batch conditions. This combination of computational modeling and strain engineering on P. putida for an advanced biofuels production has vital significance in enabling a bioproduction process that can use renewable carbon streams.

Development of Corynebacterium glutamicum as a monoterpene production platform.

Monoterpenes are commonly known for their role in the flavors and fragrances industry and are also gaining attention for other uses like insect repellant and as potential renewable fuels for aviation. Corynebacterium glutamicum, a Generally Recognized as Safe microbe, has been a choice organism in industry for the annual million ton-scale bioproduction of amino acids for more than 50 years; however, efforts to produce monoterpenes in C. glutamicum have remained relatively limited. In this study, we report a further expansion of the C. glutamicum biosynthetic repertoire through the development and optimization of a mevalonate-based monoterpene platform. In the course of our plasmid design iterations, we increased flux through the mevalonate-based bypass pathway, measuring isoprenol production as a proxy for monoterpene precursor abundance and demonstrating the highest reported titers in C. glutamicum to date at 1504.6 mg/L. Our designs also evaluated the effects of backbone, promoter, and GPP synthase homolog origin on monoterpene product titers. Monoterpene production was further improved by disrupting competing pathways for isoprenoid precursor supply and by implementing a biphasic production system to prevent volatilization. With this platform, we achieved 321.1 mg/L of geranoids, 723.6 mg/L of 1,8-cineole, and 227.8 mg/L of linalool. Furthermore, we determined that C. glutamicum first oxidizes geraniol through an aldehyde intermediate before it is asymmetrically reduced to citronellol. Additionally, we demonstrate that the aldehyde reductase, AdhC, possesses additional substrate promiscuity for acyclic monoterpene aldehydes.

3D reconstructions of parasite development and the intracellular niche of the microsporidian pathogen Encephalitozoon intestinalis.

Microsporidia are an early-diverging group of fungal pathogens with a wide host range. Several microsporidian species cause opportunistic infections in humans that can be fatal. As obligate intracellular parasites with highly reduced genomes, microsporidia are dependent on host metabolites for successful replication and development. Our knowledge of microsporidian intracellular development remains rudimentary, and our understanding of the intracellular niche occupied by microsporidia has relied on 2D TEM images and light microscopy. Here, we use serial block-face scanning electron microscopy (SBF-SEM) to capture 3D snapshots of the human-infecting species, Encephalitozoon intestinalis, within host cells. We track E. intestinalis development through its life cycle, which allows us to propose a model for how its infection organelle, the polar tube, is assembled de novo in developing spores. 3D reconstructions of parasite-infected cells provide insights into the physical interactions between host cell organelles and parasitophorous vacuoles, which contain the developing parasites. The host cell mitochondrial network is substantially remodeled during E. intestinalis infection, leading to mitochondrial fragmentation. SBF-SEM analysis shows changes in mitochondrial morphology in infected cells, and live-cell imaging provides insights into mitochondrial dynamics during infection. Our data provide insights into parasite development, polar tube assembly, and microsporidia-induced host mitochondria remodeling.

Maximizing Heterologous Expression of Engineered Type I Polyketide Synthases: Investigating Codon Optimization Strategies.

Type I polyketide synthases (T1PKSs) hold enormous potential as a rational production platform for the biosynthesis of specialty chemicals. However, despite great progress in this field, the heterologous expression of PKSs remains a major challenge. One of the first measures to improve heterologous gene expression can be codon optimization. Although controversial, choosing the wrong codon optimization strategy can have detrimental effects on the protein and product levels. In this study, we analyzed 11 different codon variants of an engineered T1PKS and investigated in a systematic approach their influence on heterologous expression in Corynebacterium glutamicum, Escherichia coli, and Pseudomonas putida. Our best performing codon variants exhibited a minimum 50-fold increase in PKS protein levels, which also enabled the production of an unnatural polyketide in each of these hosts. Furthermore, we developed a free online tool (https://basebuddy.lbl.gov) that offers transparent and highly customizable codon optimization with up-to-date codon usage tables. In this work, we not only highlight the significance of codon optimization but also establish the groundwork for the high-throughput assembly and characterization of PKS pathways in alternative hosts.

Simultaneous carbon catabolite repression governs sugar and aromatic co-utilization in Pseudomonas putida M2.

Pseudomonas putida have emerged as promising biocatalysts for the conversion of sugars and aromatic compounds obtained from lignocellulosic biomass. Understanding the role of carbon catabolite repression (CCR) in these strains is critical to optimize biomass conversion to fuels and chemicals. The CCR functioning in P. putida M2, a strain capable of consuming both hexose and pentose sugars as well as aromatic compounds, was investigated by cultivation experiments, proteomics, and CRISPRi-based gene repression. Strain M2 co-utilized sugars and aromatic compounds simultaneously; however, during cultivation with glucose and aromatic compounds (p-coumarate and ferulate) mixture, intermediates (4-hydroxybenzoate and vanillate) accumulated, and substrate consumption was incomplete. In contrast, xylose-aromatic consumption resulted in transient intermediate accumulation and complete aromatic consumption, while xylose was incompletely consumed. Proteomics analysis revealed that glucose exerted stronger repression than xylose on the aromatic catabolic proteins. Key glucose (Eda) and xylose (XylX) catabolic proteins were also identified at lower abundance during cultivation with aromatic compounds implying simultaneous catabolite repression by sugars and aromatic compounds. Reduction of crc expression via CRISPRi led to faster growth and glucose and p-coumarate uptake in the CRISPRi strains compared to the control, while no difference was observed on xylose+p-coumarate. The increased abundances of Eda and amino acid biosynthesis proteins in the CRISPRi strain further supported these observations. Lastly, small RNAs (sRNAs) sequencing results showed that CrcY and CrcZ homologues levels in M2, previously identified in P. putida strains, were lower under strong CCR (glucose+p-coumarate) condition compared to when repression was absent (p-coumarate or glucose only).IMPORTANCEA newly isolated Pseudomonas putida strain, P. putida M2, can utilize both hexose and pentose sugars as well as aromatic compounds making it a promising host for the valorization of lignocellulosic biomass. Pseudomonads have developed a regulatory strategy, carbon catabolite repression, to control the assimilation of carbon sources in the environment. Carbon catabolite repression may impede the simultaneous and complete metabolism of sugars and aromatic compounds present in lignocellulosic biomass and hinder the development of an efficient industrial biocatalyst. This study provides insight into the cellular physiology and proteome during mixed-substrate utilization in P. putida M2. The phenotypic and proteomics results demonstrated simultaneous catabolite repression in the sugar-aromatic mixtures, while the CRISPRi and sRNA sequencing demonstrated the potential role of the crc gene and small RNAs in carbon catabolite repression.

Maximizing microbial bioproduction from sustainable carbon sources using iterative systems engineering.

Maximizing the production of heterologous biomolecules is a complex problem that can be addressed with a systems-level understanding of cellular metabolism and regulation. Specifically, growth-coupling approaches can increase product titers and yields and also enhance production rates. However, implementing these methods for non-canonical carbon streams is challenging due to gaps in metabolic models. Over four design-build-test-learn cycles, we rewire Pseudomonas putida KT2440 for growth-coupled production of indigoidine from para-coumarate. We explore 4,114 potential growth-coupling solutions and refine one design through laboratory evolution and ensemble data-driven methods. The final growth-coupled strain produces 7.3 g/L indigoidine at 77% maximum theoretical yield in para-coumarate minimal medium. The iterative use of growth-coupling designs and functional genomics with experimental validation was highly effective and agnostic to specific hosts, carbon streams, and final products and thus generalizable across many systems.

3D reconstructions of parasite development and the intracellular niche of the microsporidian pathogen E. intestinalis.

Microsporidia are an early-diverging group of fungal pathogens that infect a wide range of hosts. Several microsporidian species infect humans, and infections can lead to fatal disease in immunocompromised individuals. As obligate intracellular parasites with highly reduced genomes, microsporidia are dependent on metabolites from their hosts for successful replication and development. Our knowledge of how microsporidian parasites develop inside the host remains rudimentary, and our understanding of the intracellular niche occupied by microsporidia has thus far relied largely on 2D TEM images and light microscopy. Here, we use serial block face scanning electron microscopy (SBF-SEM) to capture 3D snapshots of the human-infecting microsporidian, Encephalitozoon intestinalis , within host cells. We track the development of E. intestinalis through its life cycle, which allows us to propose a model for how its infection organelle, the polar tube, is assembled de novo in each developing spore. 3D reconstructions of parasite-infected cells provide insights into the physical interactions between host cell organelles and parasitophorous vacuoles, which contain the developing parasites. The host cell mitochondrial network is substantially remodeled during E. intestinalis infection, leading to mitochondrial fragmentation. SBF-SEM analysis shows changes in mitochondrial morphology in infected cells, and live-cell imaging provides insights into mitochondrial dynamics during infection. Together, our data provide insights into parasite development, polar tube assembly, and microsporidia-induced mitochondrial remodeling in the host cell.

Alkaline-SDS cell lysis of microbes with acetone protein precipitation for proteomic sample preparation in 96-well plate format.

Plate-based proteomic sample preparation offers a solution to the large sample throughput demands in the biotechnology field where hundreds or thousands of engineered microbes are constructed for testing is routine. Meanwhile, sample preparation methods that work efficiently on broader microbial groups are desirable for new applications of proteomics in other fields, such as microbial communities. Here, we detail a step-by-step protocol that consists of cell lysis in an alkaline chemical buffer (NaOH/SDS) followed by protein precipitation with high-ionic strength acetone in 96-well format. The protocol works for a broad range of microbes (e.g., Gram-negative bacteria, Gram-positive bacteria, non-filamentous fungi) and the resulting proteins are ready for tryptic digestion for bottom-up quantitative proteomic analysis without the need for desalting column cleanup. The yield of protein using this protocol increases linearly with respect to the amount of starting biomass from 0.5-2.0 OD*mL of cells. By using a bench-top automated liquid dispenser, a cost-effective and environmentally-friendly option to eliminating pipette tips and reducing reagent waste, the protocol takes approximately 30 minutes to extract protein from 96 samples. Tests on mock mixtures showed expected results that the biomass composition structure is in close agreement with the experimental design. Lastly, we applied the protocol for the composition analysis of a synthetic community of environmental isolates grown on two different media. This protocol has been developed to facilitate rapid, low-variance sample preparation of hundreds of samples and allow flexibility for future protocol development.

Expression of blue pigment synthetase a from Streptomyces lavenduale reveals insights on the effects of refactoring biosynthetic megasynthases for heterologous expression in Escherichia coli.

High GC bacteria from the genus Streptomyces harbor expansive secondary metabolism. The expression of biosynthetic proteins and the characterization and identification of biological "parts" for synthetic biology purposes from such pathways are of interest. However, the high GC content of proteins from actinomycetes in addition to the large size and multi-domain architecture of many biosynthetic proteins (such as non-ribosomal peptide synthetases; NRPSs, and polyketide synthases; PKSs often called "megasynthases") often presents issues with full-length translation and folding. Here we evaluate a non-ribosomal peptide synthetase (NRPS) from Streptomyces lavenduale, a multidomain "megasynthase" gene that comes from a high GC (72.5%) genome. While a preliminary step in revealing differences, to our knowledge this presents the first head-to-head comparison of codon-optimized sequences versus a native sequence of proteins of streptomycete origin heterologously expressed in E. coli. We found that any disruption in co-translational folding from codon mismatch that reduces the titer of indigoidine is explainable via the formation of more inclusion bodies as opposed to compromising folding or posttranslational modification in the soluble fraction. This result supports that one could apply any refactoring strategies that improve soluble expression in E. coli without concern that the protein that reaches the soluble fraction is differentially folded.

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.

PeakDecoder enables machine learning-based metabolite annotation and accurate profiling in multidimensional mass spectrometry measurements.

Multidimensional measurements using state-of-the-art separations and mass spectrometry provide advantages in untargeted metabolomics analyses for studying biological and environmental bio-chemical processes. However, the lack of rapid analytical methods and robust algorithms for these heterogeneous data has limited its application. Here, we develop and evaluate a sensitive and high-throughput analytical and computational workflow to enable accurate metabolite profiling. Our workflow combines liquid chromatography, ion mobility spectrometry and data-independent acquisition mass spectrometry with PeakDecoder, a machine learning-based algorithm that learns to distinguish true co-elution and co-mobility from raw data and calculates metabolite identification error rates. We apply PeakDecoder for metabolite profiling of various engineered strains of Aspergillus pseudoterreus, Aspergillus niger, Pseudomonas putida and Rhodosporidium toruloides. Results, validated manually and against selected reaction monitoring and gas-chromatography platforms, show that 2683 features could be confidently annotated and quantified across 116 microbial sample runs using a library built from 64 standards.

Expanding Extender Substrate Selection for Unnatural Polyketide Biosynthesis by Acyltransferase Domain Exchange within a Modular Polyketide Synthase.

Modular polyketide synthases (PKSs) are polymerases that employ α-carboxyacyl-CoAs as extender substrates. This enzyme family contains several catalytic modules, where each module is responsible for a single round of polyketide chain extension. Although PKS modules typically use malonyl-CoA or methylmalonyl-CoA for chain elongation, many other malonyl-CoA analogues are used to diversify polyketide structures in nature. Previously, we developed a method to alter an extension substrate of a given module by exchanging an acyltransferase (AT) domain while maintaining protein folding. Here, we report in vitro polyketide biosynthesis by 13 PKSs (the wild-type PKS and 12 AT-exchanged PKSs with unusual ATs) and 14 extender substrates. Our ∼200 in vitro reactions resulted in 13 structurally different polyketides, including several polyketides that have not been reported. In some cases, AT-exchanged PKSs produced target polyketides by >100-fold compared to the wild-type PKS. These data also indicate that most unusual AT domains do not incorporate malonyl-CoA and methylmalonyl-CoA but incorporate various rare extender substrates that are equal to in size or slightly larger than natural substrates. We developed a computational workflow to predict the approximate AT substrate range based on active site volumes to support the selection of ATs. These results greatly enhance our understanding of rare AT domains and demonstrate the benefit of using the proposed PKS engineering strategy to produce novel chemicals in vitro.

Revealing oxidative pentose metabolism in new Pseudomonas putida isolates.

The Pseudomonas putida group in the Gammaproteobacteria has been intensively studied for bioremediation and plant growth promotion. Members of this group have recently emerged as promising hosts to convert intermediates derived from plant biomass to biofuels and biochemicals. However, most strains of P. putida cannot metabolize pentose sugars derived from hemicellulose. Here, we describe three isolates that provide a broader view of the pentose sugar catabolism in the P. putida group. One of these isolates clusters with the well-characterized P. alloputida KT2440 (Strain BP6); the second isolate clustered with plant growth-promoting strain P. putida W619 (Strain M2), while the third isolate represents a new species in the group (Strain BP8). Each of these isolates possessed homologous genes for oxidative xylose catabolism (xylDXA) and a potential xylonate transporter. Strain M2 grew on arabinose and had genes for oxidative arabinose catabolism (araDXA). A CRISPR interference (CRISPRi) system was developed for strain M2 and identified conditionally essential genes for xylose growth. A glucose dehydrogenase was found to be responsible for initial oxidation of xylose and arabinose in strain M2. These isolates have illuminated inherent diversity in pentose catabolism in the P. putida group and may provide alternative hosts for biomass conversion.

Mixed heavy metal stress induces global iron starvation response.

Multiple heavy metal contamination is an increasingly common global problem. Heavy metals have the potential to disrupt microbially mediated biogeochemical cycling. However, systems-level studies on the effects of combinations of heavy metals on bacteria are lacking. For this study, we focused on the Oak Ridge Reservation (ORR; Oak Ridge, TN, USA) subsurface which is contaminated with several heavy metals and high concentrations of nitrate. Using a native Bacillus cereus isolate that represents a dominant species at this site, we assessed the combined impact of eight metal contaminants, all at site-relevant concentrations, on cell processes through an integrated multi-omics approach that included discovery proteomics, targeted metabolomics, and targeted gene-expression profiling. The combination of eight metals impacted cell physiology in a manner that could not have been predicted from summing phenotypic responses to the individual metals. Exposure to the metal mixture elicited a global iron starvation response not observed during individual metal exposures. This disruption of iron homeostasis resulted in decreased activity of the iron-cofactor-containing nitrate and nitrite reductases, both of which are important in biological nitrate removal at the site. We propose that the combinatorial effects of simultaneous exposure to multiple heavy metals is an underappreciated yet significant form of cell stress in the environment with the potential to disrupt global nutrient cycles and to impede bioremediation efforts at mixed waste sites. Our work underscores the need to shift from single- to multi-metal studies for assessing and predicting the impacts of complex contaminants on microbial systems.