1. Biochemistry and Chemical Biology
  2. Evolutionary Biology

Nitrogenase resurrection and the evolution of a singular enzymatic mechanism

  1. Amanda K Garcia
  2. Derek F Harris
  3. Alex J Rivier
  4. Brooke M Carruthers
  5. Azul Pinochet-Barros
  6. Lance C Seefeldt
  7. Betül Kaçar  Is a corresponding author
  1. Department of Bacteriology, University of Wisconsin–Madison, United States
  2. Department of Chemistry and Biochemistry, Utah State University, United States
https://doi.org/10.7554/eLife.85003
Share this article
Cite this article
  1. Amanda K Garcia
  2. Derek F Harris
  3. Alex J Rivier
  4. Brooke M Carruthers
  5. Azul Pinochet-Barros
  6. Lance C Seefeldt
  7. Betül Kaçar
(2023)
Nitrogenase resurrection and the evolution of a singular enzymatic mechanism
eLife 12:e85003.
https://doi.org/10.7554/eLife.85003
  2. Download BibTeX
  3. Download .RIS

Abstract

The planetary biosphere is powered by a suite of key metabolic innovations that emerged early in the history of life. However, it is unknown whether life has always followed the same set of strategies for performing these critical tasks. Today, microbes access atmospheric sources of bioessential nitrogen through the activities of just one family of enzymes, nitrogenases. Here, we show that the only dinitrogen reduction mechanism known to date is an ancient feature conserved from nitrogenase ancestors. We designed a paleomolecular engineering approach wherein ancestral nitrogenase genes were phylogenetically reconstructed and inserted into the genome of the diazotrophic bacterial model, Azotobacter vinelandii, enabling an integrated assessment of both in vivo functionality and purified nitrogenase biochemistry. Nitrogenase ancestors are active and robust to variable incorporation of one or more ancestral protein subunits. Further, we find that all ancestors exhibit the reversible enzymatic mechanism for dinitrogen reduction, specifically evidenced by hydrogen inhibition, which is also exhibited by extant A. vinelandii nitrogenase isozymes. Our results suggest that life may have been constrained in its sampling of protein sequence space to catalyze one of the most energetically challenging biochemical reactions in nature. The experimental framework established here is essential for probing how nitrogenase functionality has been shaped within a dynamic, cellular context to sustain a globally consequential metabolism.

Editor's evaluation

This manuscript reports valuable findings regarding the evolution of nitrogenases through ancestral sequence reconstruction and resurrection. The results are convincing and support the conclusions of the study, and highlight the historical constraints that have been acting on this enzyme. The findings will be of interest to people interested in enzyme evolution in general and particularly to those interested in the evolution of nitrogenases.

https://doi.org/10.7554/eLife.85003.sa0

Introduction

The evolutionary history of life on Earth has generated tremendous ecosystem diversity, the sum of which is orders of magnitude larger than that which exists at present (Jablonski, 2004). Life’s historical diversity provides a measure of its ability to solve adaptive problems within an integrated planetary system. Accessing these solutions requires a deeper understanding of the selective forces that have shaped the evolution of molecular-scale, metabolic innovations. However, the early histories of many of life’s key metabolic pathways and the enzymes that catalyze them remain coarsely resolved.

Important efforts to advance understanding of early metabolic innovations have included phylogenetic inference (Gold et al., 2017; Sánchez-Baracaldo and Cardona, 2020), systems-level network reconstructions (Goldford et al., 2017), and the leveraging of extant or mutant biological models as proxies for their ancient counterparts (Zerkle et al., 2006; Soboh et al., 2010; Hurley et al., 2021). However, methods to directly study these metabolic evolutionary histories across past environmental and cellular transitions remain underexplored. A unified, experimental strategy that integrates historical changes to enzymes, which serve as the primary interface between metabolism and environment, and clarifies their impact within specific cellular and physiochemical contexts is necessary. To address this, phylogenetic reconstructions of enzymes can be directly integrated within laboratory microbial model systems (Garcia and Kaçar, 2019; Kacar et al., 2017b; Kędzior et al., 2022). In this paleomolecular framework, predicted ancestral enzymes can be ‘resurrected’ within a compatible host organism for functional characterization. These experimental systems can ultimately integrate multiple levels of historical analysis by interrogating critical features of ancient enzymes as well as dynamic interactions between enzymes, their broader metabolic networks, and the external environment.

The study of biological nitrogen fixation offers a promising testbed to thread these investigations of early metabolic evolution. Both phylogenetic and geological evidence (Garcia et al., 2020; Raymond et al., 2004; Boyd et al., 2011a; Stüeken et al., 2015; Parsons et al., 2021) indicate that the origin of biological nitrogen fixation was a singular and ancient evolutionary event on which the modern biosphere has since been built. The only known nitrogen fixation pathway (compared to, for example, at least seven carbon-fixation pathways [Garcia et al., 2021a]) is catalyzed by an early-evolved family of metalloenzymes called nitrogenases that reduce highly inert, atmospheric dinitrogen (N2) to bioavailable ammonia (NH3). The nitrogenase family comprises three isozymes that vary in their metal dependence (i.e. molybdenum, vanadium, and iron) and, in certain cases, all coexist within the same host organism (Mus et al., 2018). Many diazotrophs depend on genetic strategies for coordinating the biosynthesis and expression of multiple nitrogenase isozymes and their respective metalloclusters (Pérez-González et al., 2021; Burén et al., 2020), and, in oxic environments, protecting the oxygen-sensitive metalloclusters from degradation (Gallon, 1992). Thus, nitrogenase enzymes are a central component of a broader, co-evolving nitrogen fixation machinery. These features that create significant experimental challenges for nitrogen fixation engineering (Smanski et al., 2014; Bennett et al., 2023; Burén and Rubio, 2018) nevertheless also make this metabolism an ideal candidate for systems-level, paleomolecular study.

How biological nitrogen fixation emerged and evolved under past environmental conditions is still poorly constrained relative to its importance in Earth’s planetary and biological history. Because nitrogen has been a limiting nutrient over geological timescales (Falkowski, 1997; Allen et al., 2019), nitrogenase has long been a key constituent of the expanding Earth’s biosphere. The impact of nitrogen limitation is underscored by human reliance on the industrial Haber-Bosch process, an energetically and environmentally costly workaround for nitrogen fertilizer production (Vicente and Dean, 2017) designed to supplement a remarkable molecular innovation that biology has tinkered with for more than three billion years. How the structural domains and regulatory network of nitrogenase were recruited (Boyd et al., 2015; Mus et al., 2019) and under what selective pressures the metal dependence of nitrogenases was shaped (Garcia et al., 2020; Boyd et al., 2011b) remain open questions. Importantly, it is not known how the enzymatic mechanism for dinitrogen reduction has been tuned by both peptide and metallocluster to achieve one of the most difficult reactions in nature (Seefeldt et al., 2020; Stripp et al., 2022; Harris et al., 2019). At the enzyme level, previous insights into nitrogenase sequence-function relationships have primarily derived from single or dual substitution studies. These have often yielded diminished or abolished nitrogenase activity (Seefeldt et al., 2020; Stripp et al., 2022), though in certain cases improved reactivity toward alternate, industrially relevant substrates (Seefeldt et al., 2020). Despite illuminating key features of extant nitrogenase mechanisms in select model organisms, the combination of detailed functional studies within an explicit evolutionary scheme has not previously been accomplished for the nitrogen fixation system.

Here, we seek guidance from the Earth’s evolutionary past to reconstruct the history of the key metabolic enzyme, nitrogenase. We establish an evolutionary systems biology approach for the cellular- and molecular-level characterization of ancestral nitrogenases resurrected within the model diazotrophic bacterium, A. vinelandii. We find that variably replacing different protein subunits of the nitrogenase complex with inferred ancestral counterparts enables nitrogen fixation in A. vinelandii. Purified ancestral enzymes exhibit the specific N2 reduction mechanism retained by their studied, extant counterparts, and maintain the same catalytic selectivity between N2 and protons. Thus, the core strategy for biological nitrogen fixation is conserved across the investigated timeline. Our paleomolecular approach opens a new route to study the ancient functionality and evolution of nitrogenases both deeper in its ancestry and within the broader context of its supporting, cellular machinery.

Results

A resurrection strategy for ancestral nitrogenases

We designed an engineering pipeline for the resurrection and experimental characterization of ancestral nitrogenases (Figure 1A). In this scheme, phylogenetically inferred, ancestral nitrogenase genes are synthesized and engineered into the genome of a modern diazotrophic bacterium, enabling the assessment of in vivo nitrogenase activity and expression in parallel with biochemical analysis of the purified enzyme. The engineering and functional assessment of ancestral nitrogenases required a suitable diazotrophic microbial host, owing to challenges associated with nitrogenase heterologous expression (Vicente and Dean, 2017). We selected the obligately aerobic gammaproteobacterium, A. vinelandii (strain ‘DJ’), an ideal experimental model due to its genetic tractability and the availability of detailed studies on the genetics and biochemistry of its nitrogen fixation machinery (Noar and Bruno-Bárcena, 2018). We specifically targeted the extant A. vinelandii molybdenum-dependent (‘Mo-;) nitrogenase (hereafter referred to as wild-type, ‘WT’), which is the best-studied isozyme (Seefeldt et al., 2020) relative to the A. vinelandii vanadium (‘V-’) and iron (‘Fe-’) dependent nitrogenases. The WT Mo-nitrogenase complex comprises multiple subunits, NifH, NifD, and NifK (encoded by nifHDK genes), which are arranged into two catalytic components: a NifH homodimer and a NifDK heterotetramer (Figure 1B). During catalysis, both components transiently associate to transfer one electron from NifH to NifDK and subsequently dissociate. Transferred electrons accumulate at the active-site Mo-containing metallocluster (‘FeMoco’) housed within the NifD subunits for reduction of the N2 substrate to NH3.

Figure 1
Download asset Open asset
Engineering strategy for ancestral nitrogenase resurrection.

(A) Experimental pipeline for nitrogenase resurrection in A. vinelandii and subsequent characterization, as described in the main text. (B) Structural overview of ancestral nitrogenases reconstructed in this study. Homology models (template PDB 1M34) of Anc1B NifH and NifDK proteins are shown with ancestral substitutions (relative to WT) highlighted in red. Select substitutions at relatively conserved sites in proximity to FeMoco (NifD, I355V) and the NifD:NifK interface (NifD, F429Y; NifK, R108K) are displayed in the insets. (C) Parallel genome engineering strategies were executed in this study, involving both ancestral replacement of only nifD (Anc1A and Anc2) and replacement of nifHDK (Anc1B). ‘PnifH’: nifH promoter, ‘KanR’: kanamycin resistance cassette. *Anc1A and Anc1B were each reconstructed from equivalent nodes of alternate phylogenies (see Materials and methods).

To infer Mo-nitrogenase ancestors, we built a maximum-likelihood nitrogenase phylogeny from a concatenated alignment of NifHDK amino acid sequences (Figure 2A; Figure 2—figure supplement 1). The phylogeny contains 385 sets of homologs representative of known nitrogenase molecular sequence diversity (including Mo-, V-, and Fe-nitrogenases), and is rooted by dark-operative protochlorophyllide oxidoreductase proteins classified within the nitrogenase superfamily (Ghebreamlak and Mansoorabadi, 2020). For this study, we selected ancestors that fall within the direct evolutionary lineage of A. vinelandii WT (Figure 2B), ‘Anc1’ and ‘Anc2’ (listed in order of increasing age), having ~90% and~85% amino acid sequence identity to WT across the full length of their concatenated NifHDK proteins, respectively (Figure 2C; Supplementary file 1a). A relatively conservative percentage identity threshold was chosen based on prior studies benchmarking the functional expression of ancestral elongation factor proteins in Escherichia coli (Kacar et al., 2017a) and Synechococcus elongatus (Kędzior et al., 2022). The high-dimensional, nitrogenase protein sequence space occupied by both extant and ancestral homologs is visualized in two dimensions in Figure 2D by machine-learning embeddings (see Materials and methods). This analysis highlights the swath of sequence space targeted here, as well as that made accessible by the resurrection of nitrogenase ancestors more broadly. WT, Anc1, and Anc2 lie within a Mo-nitrogenase clade (previously termed ‘Group I’ [Raymond et al., 2004]) that contains homologs from diverse aerobic and facultatively anaerobic taxa, including proteobacteria and cyanobacteria (Figure 2A). A maximum age constraint of ~2.5 Ga for Group I nitrogenases (and thus for both Anc1 and Anc2) can be reasoned based on the timing of the Great Oxidation Event (Lyons et al., 2014) and downstream emergence of aerobic taxa represented nearly exclusively within this clade.

Figure 2 with 3 supplements see all
Download asset Open asset
Phylogenetic and genomic context of resurrected ancestral nitrogenases.

(A) Maximum-likelihood phylogenetic tree from which ancestral nitrogenases were inferred. Extant nodes are colored by microbial host taxonomic diversity. Red box highlights the clade targeted in this study and depicted in (B). Tree shown was used to infer Anc1A and Anc2 sequences (an alternate tree was used for Anc1B inference; see Materials and methods). (B) nif gene cluster complexity within the targeted nitrogenase clade. Presence and absence of each nif gene are indicated by gray and white colors, respectively. Because some homologs for phylogenetic analysis were obtained from organisms lacking fully assembled genomes, the absence of accessory nif genes may result from missing genomic information. (C) Amino acid sequence identity matrix of nitrogenase protein subunits harbored by WT and engineered A. vinelandii strains. (D) Extant and ancestral nitrogenase protein sequence space visualized by machine-learning embeddings, with the resulting dimensionality reduced to two-dimensional space. UMAP dimension axes are in arbitrary units. The field demarcated by dashed lines in the left plot is expanded on the right plot.

Residue-level differences between ancestral and WT nitrogenases (‘ancestral substitutions’) are broadly distributed along the length of each ancestral sequence (Figure 1B; Figure 2—figure supplement 2). An ancestral substitution proximal to the active-site FeMoco metallocluster lies within a loop considered important for FeMoco insertion (Dos Santos et al., 2004) (NifD I355V; residue numbering from WT) and is observed across all targeted NifD ancestors (Figure 1B). Other ancestral substitutions are notable for their location at relatively conserved residue sites (assessed by ConSurf Ashkenazy et al., 2016; Figure 2—figure supplement 3; see Materials and methods) and/or within subunit interfaces, including two at the NifD:NifK interface that are proximal to one another, F429Y (NifD) and R108K (NifK). The Anc2 NifD protein contains five more ancestral substitutions at conserved sites than the younger Anc1 NifD protein. In all studied ancestors, the C275 and H442 FeMoco ligands, as well as other strictly conserved nitrogenase residues, are retained.

Phylogenetic analysis informs the compatibility of selected ancestors in extant microbial hosts. Extant nitrogenases within the Group I nitrogenase clade (which include Anc1 and Anc2 descendants) are associated with numerous accessory genes likely recruited to optimize the synthesis and regulation of the oxygen-sensitive nitrogenase for aerobic or facultative metabolisms (Boyd et al., 2015). For example, in addition to the structural nifHDK genes, A. vinelandii WT is assembled and regulated with the help of >15 additional nif genes. Likewise, the extant descendants of Anc1 and Anc2 are primarily aerobic or facultative proteobacteria and are thus associated with higher complexity nif gene clusters (Figure 2B). We hypothesized that the likely oxygen-tolerant, ancient proteobacterium harboring these ancestral nitrogenases were similar in nif cluster complexity to extant A. vinelandii. Thus, we predicted that the nif accessory genes present in A. vinelandii would support the functional expression of resurrected nitrogenase ancestors.

To gauge the degree of compatibility between ancestral and extant nitrogenase proteins and maximize the chance of recovering functional nitrogenase variants, we executed two parallel genome engineering strategies. First, we constructed A. vinelandii strains harboring only the ancestral nifD gene from both targeted ancestral nodes (‘Anc1A’ and ‘Anc2’), thereby expressing ‘hybrid ancestral-WT’ nitrogenase complexes (Figure 1C). This strategy is similar to in vitro ‘cross-reaction’ studies that have evaluated the compatibility of nitrogenase protein components from differing host taxa (Smith et al., 1976). Second, we constructed a strain harboring all Anc1 nifHDK genes, expressing a fully ancestral nitrogenase complex (‘Anc1B’; sequence reconstructed from a node equivalent to Anc1A from an alternate phylogeny, see Materials and methods). A. vinelandii strains were constructed by markerless genomic integration of ancestral nitrogenase genes, as described in Materials and methods.

Ancestral nitrogenases enable diazotrophic microbial growth

All A. vinelandii constructs harboring ancestral genes enabled diazotrophic growth in molybdenum-containing, nitrogen-free media. All strains had comparable doubling times to WT during the exponential phase (p>0.05; Figure 3A and B). The only significant difference among strains was a~14 hr increase in the lag phase of strain Anc2 relative to WT, harboring the oldest nitrogenase ancestor (p ≈ 2e-7). We did not detect growth under the same conditions for a control ΔnifD strain (DJ2278, see Supplementary file 1c). This result confirmed that the growth observed for ancestral strains did not stem from leaky expression of the alternative, V- or Fe-dependent nitrogen fixation genes in A. vinelandii, which were left intact.

Figure 3
Download asset Open asset
Cellular-level characterization of ancestral nitrogenase activity and expression.

(A) Diazotrophic growth curves of A. vinelandii strains measured by the optical density at 600 nm (‘OD600’). A smoothed curve is shown alongside individual data points obtained from five biological replicates per strain. The non-diazotrophic DJ2278 (ΔnifD) strain was used as a negative control. (B) Mean doubling and midpoint times of A. vinelandii strains, calculated from data in (A). (C) In vivo acetylene (C2H2) reduction rates quantified by the production of ethylene (C2H4). Bars represent the mean of biological replicates (n=3) per strain. (D) Immunodetection and protein quantification of Strep-II-tagged WT (‘WT-Strep,’ strain DJ2102) and ancestral NifD. Top gel image shows Strep-II-tagged NifD proteins detected by anti-Strep antibody and bottom gel image shows total protein stain. Plot displays relative immunodetected NifD signal intensity normalized to total protein intensity and expressed relative to WT. Bars in the plot represent the mean of biological replicates (n=3) per strain. (B–D) Error bars indicate ±1 SD and asterisks indicate p<.01 (one-way ANOVA, post-hoc Tukey HSD) compared to WT or WT-Strep.

Figure 3—source data 1

Source Excel file for diazotrophic growth curve data and statistical analyses.

https://cdn.elifesciences.org/articles/85003/elife-85003-fig3-data1-v2.xlsx
Figure 3—source data 2

Source Excel file for in vivo acetylene reduction assay data and statistical analyses.

https://cdn.elifesciences.org/articles/85003/elife-85003-fig3-data2-v2.xlsx
Figure 3—source data 3

Source Excel file for NifD protein densitometry data and statistical analyses.

https://cdn.elifesciences.org/articles/85003/elife-85003-fig3-data3-v2.xlsx
Figure 3—source data 4

Zip archive of Western blot image data (total protein stain, all strains, replicate 1), containing labeled and unlabeled image files.

https://cdn.elifesciences.org/articles/85003/elife-85003-fig3-data4-v2.zip
Figure 3—source data 5

Zip archive of Western blot image data (all strains, replicate 1), containing labeled and unlabeled image files.

https://cdn.elifesciences.org/articles/85003/elife-85003-fig3-data5-v2.zip
Figure 3—source data 6

Zip archive of Western blot image data (total protein stain, all strains, replicate 2), containing labeled and unlabeled image files.

https://cdn.elifesciences.org/articles/85003/elife-85003-fig3-data6-v2.zip
Figure 3—source data 7

Zip archive of Western blot image data (all strains, replicate 2), containing labeled and unlabeled image files.

https://cdn.elifesciences.org/articles/85003/elife-85003-fig3-data7-v2.zip
Figure 3—source data 8

Zip archive of Western blot image data (total protein stain, all strains, replicate 3), containing labeled and unlabeled image files.

https://cdn.elifesciences.org/articles/85003/elife-85003-fig3-data8-v2.zip
Figure 3—source data 9

Zip archive of Western blot image data (all strains, replicate 3), containing labeled and unlabeled image files.

https://cdn.elifesciences.org/articles/85003/elife-85003-fig3-data9-v2.zip

An acetylene reduction assay was performed to measure cellular nitrogenase activity in engineered strains. This assay quantifies the reduction rate of the non-physiological substrate acetylene (C2H2) to ethylene (C2H4) (Hardy et al., 1968), here normalized to total protein content. A. vinelandii strains harboring only ancestral nifD (Anc1A, Anc2) exhibited mean C2H2 reduction rates of ~5–6 μmol C2H4/mg total protein/hr,~40–45% that of WT (p ≈ 6e-4 and p ≈ 3e-4, respectively) (Figure 3C). Strain Anc1B, harboring ancestral nifHDK, exhibited a mean acetylene reduction rate of ~9 μmol C2H4/mg total protein/hr,~70% that of WT (p ≈ 3e-2).

The phenotypic variability we observed among engineered and WT A. vinelandii strains might result both from differences in nitrogenase expression and nitrogenase activity. To provide insights into these disparate effects, we quantified nitrogenase protein expression in engineered strains by immunodetection of ancestral and WT Strep-tagged NifD proteins (the latter from strain DJ2102, see Supplementary file 1c) and did not conclusively detect significant differences in protein quantity relative to WT (p>0.05; Figure 3D).

Purified ancestral nitrogenases conserve extant N2 reduction mechanisms and efficiency

Ancestral nitrogenase NifDK protein components were expressed and purified for biochemical characterization. All ancestral NifDK proteins were assayed with WT NifH protein (ancestral NifH proteins were not purified) for reduction of H+, N2, and C2H2. Ancestors were found to reduce all three substrates in vitro, supporting the cellular-level evidence of ancestral nitrogenase activity (Figure 4A).

Figure 4 with 1 supplement see all
Download asset Open asset
In vitro analyses of ancestral nitrogenase activity profiles and mechanism.

All measurements were obtained from assays using purified NifDK assayed with WT NifH. (A) Specific activities were measured for H+, N2, and C2H2 substrates. (B) Partial schematic of the reductive-elimination N2-reduction mechanism of nitrogenase is shown above, centering on the N2-binding E4(4 H) state of FeMoco (see main text for discussion) (Harris et al., 2019). (C) Inhibition of N2 reduction by H2, evidencing the mechanism illustrated in (B). (D) Catalytic efficiencies of ancestral nitrogenases, described by the ratio of formed H2 to reduced N2 (H2/N2), mapped across the targeted phylogenetic clade. NifD homology models (PDB 1M34 template) are displayed with ancestral substitutions highlighted in red. (A,C) Bars represent the mean of independent experiments (n=2) with individual data points shown as black circles.

Figure 4—source data 1

Source Excel file for nitrogenase in vitro activity data.

https://cdn.elifesciences.org/articles/85003/elife-85003-fig4-data1-v2.xlsx
Figure 4—source data 2

Source Excel file for nitrogenase in vitro H2 inhibition data.

https://cdn.elifesciences.org/articles/85003/elife-85003-fig4-data2-v2.xlsx

We investigated whether the ancestral nitrogenases studied here would exhibit the general mechanism for N2 binding and reduction that has been observed for the studied, extant nitrogenase isozymes of A. vinelandii (Mo, V, and Fe) (Harris et al., 2019; Harris et al., 2022). This mechanism involves the accumulation of four electrons/protons on the active-site cofactor as metal-bound hydrides, generating the E4(4 H) state (Figure 4B). Once generated, N2 can bind to the E4(4 H) state through a reversible reductive elimination/oxidative addition (re/oa) mechanism, which results in the release (re) of a single molecule of hydrogen gas (H2). N2 binding is reversible in the presence of sufficient H2, which displaces bound N2 and results in the reformation of E4(4 H) with two hydrides (oa). Thus, a classic test of the (re/oa) mechanism is the ability of H2 to inhibit N2 reduction. We observed that the reduction of N2 to NH3 for all nitrogenase ancestors was inhibited in the presence of H2, indicating that the ancestors follow the same mechanism of N2 binding determined for extant enzymes (Figure 4C).

In the event the E4(4 H) state fails to capture N2, nitrogenases will simply produce H2 from the E4(4 H) state to generate the E2(2 H) state. The ratio of H2 formed to N2 reduced (H2/N2) can be used as a measure of the efficiency of nitrogenases in using ATP and reducing equivalents for N2 reduction. The stoichiometric minimum of the mechanism is H2/N2=1. Experimentally (under 1 atm N2), a ratio of ~2 is seen for Mo-nitrogenase and ~5 and~7 for V- and Fe-nitrogenase, respectively (Harris et al., 2019). H2/N2 values for all ancestors under 1 atm N2 was ~2, similar to extant Mo-nitrogenase (Figure 4D).

Discussion

In this study, we leverage a new approach to investigate ancient nitrogen fixation by the resurrection and functional assessment of ancestral nitrogenase enzymes. We demonstrate that engineered A. vinelandii cells can reduce N2 and C2H2 and exhibit diazotrophic growth rates comparable to WT, though we observe that the oldest ancestor, Anc2, has a significantly longer lag phase. Purified nitrogenase ancestors are active for the reduction of H+, N2, and C2H2, while maintaining the catalytic efficiency (described by the H2/N2 ratio) of WT enzymes. Our results also show that ancestral N2 reduction is inhibited by H2, indicating an early emergence of the reductive-elimination N2-reduction mechanism preserved by characterized, extant nitrogenases (Harris et al., 2019; Harris et al., 2022). These properties are maintained despite substantial residue-level changes to the peripheral nitrogenase structure (including relatively conserved sites), as well as a handful within the active-site or protein-interface regions within the enzyme complex.

It is important to consider that the nitrogenase ancestors resurrected here represent hypotheses regarding the true ancestral state. Uncertainty underlying ancestral reconstructions might derive from incomplete extant molecular sequence data, as well as incorrect assumptions associated with the implemented evolutionary models (Garcia and Kaçar, 2019). For example, a complicating feature of nitrogenase evolution is that it has been shaped significantly by horizontal gene transfer (Raymond et al., 2004; Parsons et al., 2021), which in certain cases has led to different evolutionary trajectories of individual nitrogenase structural genes. Specifically, certain H-subunit genes of the V-nitrogenase appear to have different evolutionary histories relative to their DK components (Raymond et al., 2004). However, since our study targets a Mo-nitrogenase lineage, we do not expect horizontal transfer to be a significant source of uncertainty in our reconstructions.

The N2-reduction activity of nitrogenase ancestors suggests that the required protein-protein interactions—both between subunits that comprise the nitrogenase complex as well as those required for nitrogenase assembly in A. vinelandii—and metallocluster interactions are sufficiently maintained for primary function. Still, our results reveal the degree to which the organism-level phenotype of host strains can be perturbed by varying both the number and age of ancestral subunits. Importantly, these changes appear to impact phenotypic properties in complex ways, representative of the type of cellular constraints on nitrogenase evolution that would be unobservable through an in vitro study alone. For example, we observed comparable growth characteristics of strains harboring single (Anc1A, ancestral NifD) versus multiple (Anc1B, ancestral NifHDK) ancestral subunits of equivalent age, whereas the lag phase of Anc2 hosting a single, older subunit (ancestral NifD) was increased. Here, growth is sensitive to older, ancestral substitutions in a single subunit while permissive of more recent ancestral substitutions across one or more subunits within the nitrogenase complex. However, a different pattern is observed across in vivo acetylene reduction rates. These are most negatively impacted relative to WT in strains with a single NifD ancestor (Anc1A, Anc2), whereas rates are more modestly decreased in a strain with a complete, ancestral NifHDK complex (Anc1A). These results suggest that ancestral subunits of equivalent age have greater compatibility and yield greater in vivo activity compared to subunits of disparate ages, perhaps owing to modified protein interactions within the nitrogenase complexes The discrepancy between in vivo activity and growth characteristics may also be attributable to impacted cellular processes external to the biochemical properties of the nitrogenase complex itself, and yet nevertheless vital in determining the overall fitness of the host organism. Finally, though we do not detect significant differences in ancestral protein expression here, it is possible that phenotypic outcomes of future reconstructions might be impacted by perturbed expression levels (e.g. Kędzior et al., 2022; Garcia et al., 2021b). To what degree these expression levels are representative of the ancestral state and impact the phenotypic property of interest should be considered in future work.

That nitrogenase ancestors perform the reductive-elimination N2-reduction mechanism—as the distantly related (Garcia et al., 2020), extant Mo-, V-, and Fe-nitrogenases of A. vinelandii do today (Harris et al., 2019)—likely indicates that this enzymatic characteristic was set early in nitrogenase evolutionary history and sustained through significant past environmental change (Lyons et al., 2014; Som et al., 2016; Catling and Zahnle, 2020) and ecological diversification (Boyd et al., 2015; Zehr et al., 2003). It is possible that life’s available strategies for achieving N2 reduction may be fundamentally limited, and that a defining constraint of nitrogenase evolution has been the preservation of the same N2 reduction mechanism across shifting selective pressures. For example, in the acquisition of V- and Fe-dependence from Mo-dependent ancestors (Garcia et al., 2020), nitrogenases may have required substantial sequence and structural changes (Sippel and Einsle, 2017; Eady, 1996) in order to facilitate reductive elimination given a different active-site metallocluster. It is also possible that alternate strategies for biological nitrogen fixation evolved early in the history of life and were subsequently outcompeted, leaving no trace of their existence in extant microbial genomes. Why these alternate possibilities were evidently not explored by nature to the same degree remains an open question, particularly given the several abiotic mechanisms for nitrogen fixation (Cherkasov et al., 2015; Dörr et al., 2003; Yung and McElroy, 1979) and the multiple biological pathways for another, globally significant metabolism, carbon fixation (Garcia et al., 2021a). Because our paleomolecular approach is ultimately informed by extant sequence data, it cannot directly evaluate extinct sequences that, for instance, due to contingency or entrenchment, did not persist and become preserved in extant microbial genomes. Nevertheless, evolutionarily informed studies of nitrogenase functionality that define the sequence-function space of this enzyme family will provide a foundation for laboratory efforts aimed toward exploring alternate scenarios. Future work that explores deeper into nitrogenase evolutionary history (and across extant and ancestral nitrogenase sequence space, as charted here (Figure 2D)) will clarify the degree of functional constraint exhibited by the nitrogenase family, both past and present.

Conclusion

Broadening the historical level of analysis beyond a single enzyme to the organism level is necessary to generate comprehensive insights into the evolutionary history and engineering potential of nitrogen fixation. Paleomolecular work that has expanded toward the systems-level investigation of early-evolved, crucial metabolic pathways remains in its infancy, despite the potential for provocative connections between molecular-scale innovations and planetary history (Garcia and Kaçar, 2019; Kędzior et al., 2022). Our results highlight the evolutionary conservation of a critical metabolic pathway that has shaped the biosphere over billions of years, as well as establish the tractability of leveraging phylogenetic models to carry out extensive, empirical manipulations of challenging enzymatic systems and their microbial hosts. Building on the empirical framework presented here will illuminate the evolutionary design principles behind ancient metabolic systems more broadly as well as leverage these histories to understand how key enzymes that allowed organisms to access nitrogen from the atmosphere evolved.

Materials and methods

Key resources table
Reagent type (species) or resourceDesignationSource or referenceIdentifiersAdditional information
strain, strain background (A. vinelandii)DJDOI:10.1128/JB.00504–09n/aDennis Dean, Virginia Tech; Wild-type (WT); Nif+
genetic reagent (A. vinelandii)DJ2102DOI:10.1016/bs.mie.2018.10.007n/aDennis Dean, Virginia Tech; Strep-tagged WT NifD; Nif+
genetic reagent (A. vinelandii)DJ2278Othern/aDennis Dean, Virginia Tech; ΔnifD::KanR; Nif-
genetic reagent (A. vinelandii)DJ884Othern/aDennis Dean, Virginia Tech; nifDR187I mutant; Nif+(slow); overexpresses NifH
genetic reagent (A. vinelandii)AK022This papern/aΔnifHDK::KanR; Nif-
genetic reagent (A. vinelandii)AK013This papern/a‘Anc1A’; ΔnifD::nifDAnc1A; Nif+
genetic reagent (A. vinelandii)AK023This papern/a‘Anc1B’;
ΔnifHDK::nifHDKAnc1B; Nif+
genetic reagent (A. vinelandii)AK014This papern/a‘Anc2’;
ΔnifD::nifDAnc2; Nif+
antibodyStrepMAB-Classic (Mouse monoclonal)IBA LifesciencesCat# 2-1507-001, RRID: AB_513133WB (1:5000)
recombinant DNA reagentpAG25This papern/aKanR cassette (APH(3’)-I gene)+400 bp nifHDK flanking homology regions, synthesized into XbaI/KpnI sites in pUC19; used to construct strain AK022 from DJ
recombinant DNA reagentpAG13This papern/anifDAnc1A + 400-bp nifD flanking homology regions, synthesized into XbaI/KpnI sites in pUC19; used to construct strain Anc1A from AK022
recombinant DNA reagentpAG19This papern/anifHDKAnc1B + 400-bp nifHDK flanking homology regions, synthesized into XbaI/KpnI sites in pUC19; used to construct strain Anc1B from AK022
recombinant DNA reagentpAG14This papern/anifDAnc2 +400 bp nifD flanking homology regions, synthesized into XbaI/KpnI sites in pUC19; used to construct strain Anc2 from AK022
sequence-based reagent306_nifH_FThis paperPCR primersGCCGAACGTTCAAGTGGAAA
sequence-based reagent307_nifH_RThis paperPCR primersAGAGCCAATCTGCCCTGTC
sequence-based reagent308_nifD_FThis paperPCR primersCACCCGTTACCCGCATATGA
sequence-based reagent309_nifD_RThis paperPCR primersACTCATCTGTGAACGGCGTT
sequence-based reagent310_nifK_FThis paperPCR primersGCTAACGCCGTTCACAGATG
sequence-based reagent311_nifK_RThis paperPCR primersTCAGTTGGCCTTCGTCGTTG
software, algorithmMAFFTMAFFTRRID:SCR_011811
software, algorithmtrimAltrimAlRRID:SCR_017334
software, algorithmIQ-TREEIQ-TREERRID:SCR_017254
software, algorithmRAxMLRAxMLRRID:SCR_006086
software, algorithmPAMLPAMLRRID:SCR_014932
software, algorithmMODELLERMODELLERRRID:SCR_008395
software, algorithmChimeraXChimeraXRRID:SCR_015872
software, algorithmGrowthcurverGrowthcurvern/aR package

Nitrogenase ancestral sequence reconstruction and selection

The nitrogenase protein sequence dataset was assembled by BLASTp (Camacho et al., 2009) search of the NCBI non-redundant protein database (accessed August 2020) with A. vinelandii NifH (WP_012698831.1), NifD (WP_012698832.1), and NifK (WP_012698833.1) queries and a 1e-5 Expect value threshold (Supplementary file 1b). BLASTp hits were manually curated to remove partially sequenced, misannotated, and taxonomically overrepresented homologs. BLASTp hits included protein sequences from homologous Mo-, V-, and Fe-nitrogenase isozymes (Garcia et al., 2020). H-, D, and K-subunit sequences from these isozymes were individually aligned by MAFFT v7.450 (Katoh and Standley, 2013) and concatenated along with outgroup dark-operative protochlorophyllide oxidoreductase sequences (Bch/ChlLNB). The final dataset included 385 nitrogenase sequences and 385 outgroup sequences. For sequences used to construct Anc1A and Anc2 (internal nodes #960 and #929, respectively), tree reconstruction (using a trimmed alignment generated by trimAl v1.2 [Capella-Gutiérrez et al., 2009]), and ancestral sequence inference (using the initial untrimmed alignment) were both performed by RAxML v8.2.10 (Stamatakis, 2014) with the LG +G + F evolutionary model (model testing performed by the ModelFinder Kalyaanamoorthy et al., 2017 in the IQ-TREE v.1.6.12 package [Nguyen et al., 2015]).

Due to concerns that RAxML v.8.2 does not implement full, marginal ancestral sequence reconstruction as described by Yang et al., 1995, we performed a second phylogenetic analysis as follows. The extant sequence dataset described above was realigned by MAFFT (untrimmed) and tree reconstruction was again performed by RAxML. Ancestral sequence reconstruction was instead performed by PAML v4.9j (Yang, 2007) using the LG +G + F model. From this second reconstruction, Anc1B (internal node #1312), equivalent to Anc1A, was selected for experimental analysis. Anc1B and Anc1A have identical sets of descendent homologs, and their NifD proteins are 95% identical (Figure 2B).

Only the ancestral sequences inferred with the most probable residue at each protein site were considered for this study (mean posterior probabilities of targeted nitrogenase subunits range from 0.95 to 0.99; see Supplementary file 1a). All ancestral sequences were reconstructed from well-supported clades (SH-like aLRT = 99–100 [Anisimova and Gascuel, 2006]).

Ancestral nitrogenase structural modeling and sequence analysis

Structural homology models of ancestral sequences were generated by MODELLER v10.2 (Webb and Sali, 2016) using PDB 1M34 as a template for all nitrogenase protein subunits and visualized by ChimeraX v1.3 (Pettersen et al., 2021).

Extant and ancestral protein sequence space was visualized by machine-learning embeddings, where each protein embedding represents protein features in a fixed-size, multidimensional vector space. The analysis was conducted on concatenated (HDK) nitrogenase protein sequences in our phylogenetic dataset. The embeddings were obtained using the pre-trained language model ESM2 (Lin et al., 2022; Rives et al., 2021), a transformer architecture trained to reproduce correlations at the sequence level in a dataset containing hundreds of millions of protein sequences. Layer 33 of this transformer was used, as recommended by the authors. The resulting 1024 dimensions were reduced by UMAP (McInnes et al., 2020) for visualization in a two-dimensional space.

Protein site-wise conservation analysis was performed using the Consurf server (Ashkenazy et al., 2016). An input alignment containing only extant, Group I Mo-nitrogenases was submitted for analysis under default parameters. Conserved sites were defined by a Consurf conservation score >7.

A. vinelandii strain engineering

Nucleotide sequences of targeted ancestral nitrogenase proteins were codon-optimized for A. vinelandii by a semi-randomized strategy that maximized ancestral nucleotide sequence identity to WT genes. Ancestral and WT protein sequences were compared using the alignment output of ancestral sequence reconstruction (RAxML or PAML). For sites where the ancestral and WT residues were identical, the WT codon was assigned. At sites where the residues were different, the codon was assigned randomly, weighted by A. vinelandii codon frequencies (Codon Usage Database, https://www.kazusa.or.jp/codon/). Nucleotide sequences were synthesized into XbaI/KpnI sites of pUC19 vectors (unable to replicate in A. vinelandii) (Twist Bioscience; GenScript). Inserts were designed with 400-base-pair flanking regions for homology-directed recombination at the relevant A. vinelandii nif locus. An ‘ASWSHPQFEK’ Strep-II-tag was included at the N-terminus of each synthetic nifD gene for downstream NifD immunodetection and NifDK affinity purification. See Supplementary file 1c for a list of strains and plasmids used in this study.

Engineering of A. vinelandii strains used established methods, following Dos Santos, 2019. A. vinelandii WT (‘DJ’), DJ2278 (ΔnifD::KanR), DJ2102 (Strep-II-tagged WT NifD), and DJ884 (NifH-overexpression mutant) strains were generously provided by Dennis Dean (Virginia Tech) (Supplementary file 1c). Strains Anc1A and Anc2 were constructed from the DJ2278 parent strain via transformation with plasmids pAG13 and pAG19, respectively (Supplementary file 1c). For the construction of strain Anc2, we first generated a ΔnifHDK strain, AK022, by transforming the DJ strain with pAG25. Genetic competency was induced by subculturing relevant parent strains in Mo- and Fe-free Burk’s medium (see below). Competent cells were transformed with at least 1 μg of donor plasmid. Transformants were screened on solid Burk’s medium for the rescue of the diazotrophic phenotype (‘Nif+’) and loss of kanamycin resistance, followed by Sanger sequencing of the PCR-amplified nifHDK cluster (see Supplementary file 1d for a list of primers). Transformants were passaged at least three times to ensure phenotypic stability prior to storage at –80 °C in phosphate buffer containing 7% DMSO.

A. vinelandii culturing and growth analysis

A. vinelandii strains were grown diazotrophically in nitrogen-free Burk’s medium (containing 1 μM Na2MoO4) at 30 °C and agitated at 300 rpm. To induce genetic competency for transformation experiments, Mo and Fe salts were excluded. For transformant screening, kanamycin antibiotic was added to solid Burk’s medium at a final concentration of 0.6 μg/mL. 50 mL seed cultures for growth rate and acetylene reduction rate quantification were grown non-diazotrophically in flasks with Burk’s medium containing 13 mM ammonium acetate.

For growth rate quantification, seed cultures were inoculated into 100 mL nitrogen-free Burk’s medium to an optical density of ~0.01 at 600 nm (OD600), after Carruthers et al., 2021, and monitored for 72 hr. Growth parameters were modeled using the R package Growthcurver (Sprouffske and Wagner, 2016).

Microbial acetylene reduction assays

A. vinelandii seed cultures representing independent biological replicates were prepared as described above and used to inoculate 100 mL of nitrogen-free Burk’s medium to an OD600 ≈ 0.01. Cells were grown diazotrophically to an OD600 ≈ 0.5, at which point a rubber septum cap was affixed to the mouth of each flask. 25 mL of headspace was removed and replaced by injecting an equivalent volume of acetylene gas. The cultures were subsequently shaken at 30 °C and agitated at 300 rpm. Headspace samples were taken after 15, 30, 45, and 60 min of incubation for ethylene quantification by a Nexis GC-2030 gas chromatograph (Shimadzu). After the 60 min incubation period, cells were pelleted at 4700 rpm for 10 min, washed once with 4 mL of phosphate buffer, and pelleted once more under the same conditions prior to storage at –80 °C. Total protein was quantified using the Quick Start Bradford Protein Assay kit (Bio-Rad) according to manufacturer instructions and a CLARIOstar Plus plate reader (BMG Labtech). Acetylene reduction rates for each replicate were normalized to total protein.

Nitrogenase expression analysis

Strep-II-tagged NifD protein quantification was performed on all ancestral strains (Anc1A, Anc1B, Anc2) and DJ2102 (harboring Strep-II-tagged WT NifD). Diazotrophic A. vinelandii cultures (100 mL) representing three independent biological replicates were prepared as described above and harvested at an OD600 ≈ 1. Cell pellets were resuspended in TE lysis buffer (10 mM Tris, 1 mM EDTA, 1 mg/mL lysozyme) and heated at 95 °C for 10 min. Cell lysates were centrifuged at 5000 rpm for 15 min. Total protein in the resulting supernatant was quantified using the Pierce BCA Protein Assay kit (ThermoFisher) following manufacturer instructions. Normalized protein samples were diluted in 2×Laemmli buffer at a 1:1 (v/v) ratio prior to SDS-PAGE analysis. Proteins were transferred to nitrocellulose membranes (ThermoFisher), stained with Revert 700 Total Protein Stain (LI-COR), and imaged on an Odyssey Fc Imager (LI-COR). Membranes were then destained with Revert Destaining Solution (LI-COR) and blocked with 5% non-fat milk in PBS solution (137 mM NaCl, 2.7 mM KCl, 10 mM Na2HPO4, 1.8 mM KH2PO4.) for 1 hr at room temperature. Membranes were rinsed once with PBS-T (PBS with 0.01% Tween-20) and incubated with primary Strep-tag II antibody (Strep-MAB-Classic, IBA Lifesciences, Cat# 2-1507-001, RRID: AB_513133; 1:5000 in 0.2% BSA) for 2 hr at room temperature. Membranes were then incubated in LI-COR blocking buffer containing 1:15,000 IRDye 680RD Goat anti-Mouse (LI-COR) for 2 hr at room temperature and subsequently imaged with an Odyssey Fc Imager (LI-COR). Densitometry analysis was performed with ImageJ (Schneider et al., 2012), with Strep-II-tagged NifD signal intensity normalized to that of the total protein stain.

Nitrogenase expression, purification, and biochemical characterization

Ancestral nitrogenase NifDK proteins were expressed from relevant A. vinelandii strains (Anc1A, Anc1B, Anc2) and purified according to previously published methods (Jiménez-Vicente et al., 2018) with the following modifications: cells were grown diazotrophically in nitrogen-free Burk’s medium and no derepression step to a sufficient OD600 (~1.8) before harvesting. WT NifH was expressed in A. vinelandii strain DJ884 and purified by previously published methods (Christiansen et al., 1998). Protein purity was assessed at ≥95% by SDS-PAGE gel with Coomassie blue staining (Figure 4—figure supplement 1).

Assays were performed in 9.4 mL vials with a MgATP regeneration buffer (6.7 mM MgCl2, 30 mM phosphocreatine, 5 mM ATP, 0.2 mg/mL creatine phosphokinase, 1.2 mg/mL BSA) and 10 mM sodium dithionite in 100 mM MOPS buffer at pH 7.0. Reaction vials were made anaerobic and relevant gases (N2, C2H2, H2) were added to appropriate concentrations with the headspace balanced by argon. NifDK proteins (~240 kDa) were added to 0.42 µM, the vial vented to atmospheric pressure, and the reaction initiated by the addition of NifH (~60 kDa) protein to 8.4 µM. Reactions were run, shaking, at 30 °C for 8 min and stopped by the addition of 500 µL of 400 mM EDTA pH 8.0. NH3 was quantified using a fluorescence protocol (Corbin, 1984) with the following modifications: an aliquot of the sample was added to a solution containing 200 mM potassium phosphate pH 7.3, 20 mM o-phthalaldehyde, and 3.5 mM 2-mercaptoethanol, and incubated for 30 min in the dark. Fluorescence was measured at λexcitation of 410 nm and λemission of 472 nm and NH3 was quantified using a standard generated with NH4Cl. H2 and C2H4 were quantified by gas chromatography with a thermal conductivity detector (GC-TCD) and gas chromatography with a flame ionization detector (GC-FID) respectively, according to published methods (Khadka et al., 2016; Yang et al., 2011).

Statistical analyses

Experimental data were statistically analyzed by one-way ANOVA with the post-hoc Tukey HSD test.

Data availability

Materials including bacterial strains and plasmids are available to the scientific community upon request. Phylogenetic data, including sequence alignments and phylogenetic trees, and the script for ancestral gene codon-optimization are publicly available at https://github.com/kacarlab/garcia_nif2023, (copy archived at swh:1:rev:c9b3cf5021e50b4a0995b3972ad81d5cedea4ed5). All other data are included as source data and supplementary files.

References

    1. Gallon JR
    (1992)
    Reconciling the incompatible - N2 fixation and O2
    The New Phytologist 122:571–609.
    1. Kacar B
    2. Guy L
    3. Smith E
    4. Baross J
    (2017b) Resurrecting ancestral genes in bacteria to interpret ancient biosignatures
    Philosophical Transactions. Series A, Mathematical, Physical, and Engineering Sciences 375:20160352.
    https://doi.org/10.1098/rsta.2016.0352

Decision letter

  1. Christian R Landry
    Senior and Reviewing Editor; Université Laval, Canada
  2. Matilda Newton
    Reviewer; University of Colorado Boulder, United States
  3. Christian B Macdonald
    Reviewer; University of California, San Francisco, San Francisco, United States

Our editorial process produces two outputs: (i) public reviews designed to be posted alongside the preprint for the benefit of readers; (ii) feedback on the manuscript for the authors, including requests for revisions, shown below. We also include an acceptance summary that explains what the editors found interesting or important about the work.

Decision letter after peer review:

Thank you for submitting your article "Nitrogenase resurrection and the evolution of a singular enzymatic mechanism" for consideration by eLife. Your article has been reviewed by 3 peer reviewers, and the evaluation has been overseen by a Reviewing Editor and Christian Landry as the Senior Editor. The following individuals involved in the review of your submission have agreed to reveal their identity: Matilda Newton (Reviewer #2); Christian B Macdonald (Reviewer #3).

The reviewers have discussed their reviews with one another, and the Reviewing Editor has drafted this to help you prepare a revised submission.

Essential revisions:

– The reviewers raised questions about the impact of horizontal gene transfers on phylogenetic reconstructions. Further discussions led to the conclusion that this is probably not a major issue but it would be important to address this point in the paper.

– Reviewer 1 raised the issue that two distinct phylogenies were obtained with the same dataset as well as issues with the reconstruction methods implemented in the computational tools used. These would be important points to verify and clarify as needed.

Reviewer #1 (Recommendations for the authors):

I would strongly suggest the RaxML issue is explained in the methods. As a community, we should make sure that rigorous methods are used, especially since there is virtually no significant computational cost to the full algorithm anymore. As the Kacar lab is one of the leading groups in this field, it would be great to be clear about this.

If I misunderstood, and the authors' version of RaxML did in fact use the correct algorithm, this part of my public review should be deleted. It would still raise the question as to why the PAML and RaxML reconstructions differ, which would still need to be explained in the methods.

I did not understand exactly what the language model analysis shows. It is only mentioned in one sentence in the paper that does not explain the meaning of this analysis.

Reviewer #2 (Recommendations for the authors):

• Include a caveat that the reconstructed enzymes can only ever be hypotheses.

• I appreciated the insertion of the gene into the genome as opposed to plasmid-based expression

• ASR can only ever reconstruct the ancestor to extant enzymes, we can't rule out that there were other N2 fixation strategies competing pre-LUCA but they have left no trace.

• Why did you choose to use maximum likelihood inference instead of Bayesian?

• Figure 3A – please show the y-axis in the log.

• Could you please represent catalysis in units for kcat and KM? This would make the kinetics easier to compare to other enzymes, both extant and reconstructed.

• 4C please include WT.

• In light of the conclusion about a highly conserved mechanism, I would like to see more nuanced mechanistic studies. Pre-steady state kinetics; MD; pH studies.

• Please infer the age of the hypothetical ancestors expressing anc1 and anc2.

• Is there a notable difference between how the all-anc complexes (Anc1B) interact as opposed to the hybrids (Anc1A, Anc2)? Is there a notable difference in melting temperature or oligomeric state? This is discussed to a degree in the paragraph beginning line 323, but many of the statements are general and do not posit molecular explanations. What do the authors mean by "historical" amino acid substitutions?

• The discussion notes the surprising conservation of inhibition by H+ despite "substantial residue-level changes to the peripheral nitrogenase structure, as well as a handful within relatively conserved, active-site or protein-interface regions within the enzyme complex" Please elaborate on this, with specific attention to the active site. Are the residues involved/mechanisms known? Specify the mutations, how chemically conserved are they?

• I would be interested to see a discussion of potential ancestral promoters and expression levels. Expression levels are mentioned briefly in the results. I know this experiment must be compatible with the biological system used for the experiment (i.e. it would be impractical to also reconstruct ancestral promoters), but do the authors speculate that an ancestral nitrogenase would have been overexpressed to compensate for lower efficiency or that N2 fixation would merely have been rate-limited?

• For the uninitiated, please clearly introduce the evolution of metal ion dependence in nitrogenases.

• How do your results compare to enzyme reconstructions of a similar "age"?

• I find these results unsurprising. There is sufficient ASR literature to predict that a reconstructed enzyme will have comparable activity to the one or two extant enzymes it is compared to. When using a clade of conserved enzymes using a conserved mechanism, it is not surprising the ancestor conforms. What have we learned? It would be more interesting to probe these reconstructions for promiscuous activities or additional inhibitors. Are they easier to evolve than extant enzymes? If you reconstruct the whole pathway, does it behave differently? Does it act inefficiently and leak metabolites? Are other methods of fixation conceivable?

Reviewer #3 (Recommendations for the authors):

I have several questions and suggestions for the phylogenetic analysis that I do not believe will alter any of the results but may help with presentation or for a better understanding of the uncertainty with them.

• Why were these particular nodes picked for reconstruction? Are they the highest confidence reconstructions?

• Why was LG+G+F used? Was any model testing done?

• What tool was used to align the sequences?

• Why was ASR performed with RAxML and PAML both?

• What is the difference between the forak023 and forak013-14 files in the github repository? The topology of tree_forak013-14_branchsupport.tre seems to be the one in the manuscript, but there are no branch support values in that file.

As mentioned in the public comments, I worry whether HGT may cause issues during tree inference. I believe the simplest way to find out would be to reconstruct a gene tree for each individual nif gene and see how the trees differ. It could also be worth examining whether these or the tree in the manuscript agree with bacterial phylogenies.

The randomly-selected codon (de)optimization process is a nice inclusion. I was a little unsure about how it was performed, though – were codons swapped randomly until some metric was reached, or a fixed number, or some other procedure? Given the other controls, I do not expect this to change any results, but it would be a nice method for other groups to potentially use.

Given that the WT and ancestral sequences are 83% identical or higher (roughly akin to humans and mice), is the result that function is conserved surprising? Is it possible that these results say more about sequence (and functional) conservation, rather than a constraint? The UMAP embedding is an interesting approach, but makes this point in a different way, as the ancestral and WT sequences are extremely close in UMAP space. I believe some discussion of sequence.

https://doi.org/10.7554/eLife.85003.sa1

Author response

Essential revisions:

– The reviewers raised questions about the impact of horizontal gene transfers on phylogenetic reconstructions. Further discussions led to the conclusion that this is probably not a major issue but it would be important to address this point in the paper.

– Reviewer 1 raised the issue that two distinct phylogenies were obtained with the same dataset as well as issues with the reconstruction methods implemented in the computational tools used. These would be important points to verify and clarify as needed.

The essential revision recommendations relate to (1) the impact of horizontal gene transfer on our reconstructions, and (2) our use of two distinct phylogenies.

Regarding horizontal gene transfer, we have addressed the bulk of these concerns in our response to Reviewer #3. Briefly, we do acknowledge that horizontal gene transfer has significantly shaped the evolution of nitrogenases, and we would not expect agreement between protein and species trees. However, these do not impact our central conclusions since we do not make inferences of species-level divergence. A potential impact is if the individual proteins in our concatenated alignments follow different evolutionary trajectories. For the reasons described in more detail in the following responses, we do not expect this to contribute significantly to the uncertainty associated with our specific reconstructions.

We thank Reviewer #1 for outlining the issue with RAxML, which we became aware of in the process of carrying out this study. We therefore repeated our study with PAML, which, as Reviewer #1 points out, appears to yield comparable experimental results. We have followed the reviewer’s recommendation to highlight the discrepancies in ancestral sequence reconstruction algorithms, so the broader community does not inadvertently use the incorrect algorithm.

Reviewer #1 (Recommendations for the authors):

I would strongly suggest the RaxML issue is explained in the methods. As a community, we should make sure that rigorous methods are used, especially since there is virtually no significant computational cost to the full algorithm anymore. As the Kacar lab is one of the leading groups in this field, it would be great to be clear about this.

We fully welcome the reviewer’s recommendation to explain the RAxML v.8 issue in our Materials and methods. We initially used RAxML v.8 for ancestral sequence reconstruction (ASR) (e.g., following Aadland and Kolaczkowski, Genome Biol Evol, 2020). Despite RAxML v.8 documentation describing the algorithm as marginal reconstruction (as the reviewer notes), we were subsequently made aware that it is not the correct implementation. We therefore repeated our analysis using RAxML v.8 for tree reconstruction and PAML for ASR. As the reviewer also points out, ancestral sequences reconstructed from equivalent nodes in the RAxML and PAML analyses were very similar (~95% identical) and both exhibit the core N2 reduction mechanism described in the main text. Therefore, uncertainty associated with our use of the incorrect ASR algorithm does not impact the central findings of our study.

In our revised text, we clarify that RAxML v.8 does not implement full marginal ancestral sequence reconstruction and justify our repeated ASR analysis in more detail.

If I misunderstood, and the authors' version of RaxML did in fact use the correct algorithm, this part of my public review should be deleted. It would still raise the question as to why the PAML and RaxML reconstructions differ, which would still need to be explained in the methods.

The reviewer is correct that the RAxML version used did not implement the full algorithm, as outlined above.

I did not understand exactly what the language model analysis shows. It is only mentioned in one sentence in the paper that does not explain the meaning of this analysis.

Our aim with the language model analysis was not to test a specific hypothesis, but to visualize the protein sequence space occupied both by extant and ancestral nitrogenases. On one hand, it places the studied ancestors in the context of this available diversity, and also charts a “roadmap” for future studies that will navigate a broader swath of this sequence space. We have now clarified the text accordingly.

Reviewer #2 (Recommendations for the authors):

• Include a caveat that the reconstructed enzymes can only ever be hypotheses.

We have included a new paragraph in our Discussion that addresses uncertainty in ancestral sequence reconstruction, including the fact that reconstructed ancestral enzymes represent hypotheses regarding the true ancestral state.

• I appreciated the insertion of the gene into the genome as opposed to plasmid-based expression

We’re glad the reviewer appreciated this crucial feature of our study.

• ASR can only ever reconstruct the ancestor to extant enzymes, we can't rule out that there were other N2 fixation strategies competing pre-LUCA but they have left no trace.

We agree that ASR is fundamentally limited to reconstructing the ancestors of whatever descendent enzymes have survived to the present. It is true that we cannot exclude the possibility that other, early-evolved nitrogen-fixing enzymes followed different strategies that were subsequently outcompeted. We address these possibilities in additional text in our Discussion, and have clarified text that discusses evolutionary constraints on nitrogen fixation strategies so as to capture these possibilities.

• Why did you choose to use maximum likelihood inference instead of Bayesian?

Both maximum-likelihood (ML) and Bayesian phylogenetic methods have previously been applied to nitrogenase evolutionary studies, for example in our (authors Garcia, Kaçar) earlier work (Garcia et al., Geobiology, 2020; Garcia et al., Genome Biol Evol, 2021) and others (Boyd et al., Front Microbiol, 2011). With either method, we observe that the general topology of nitrogenase trees is maintained. Additionally, previous work has demonstrated that Bayesian methods don’t necessarily generate more accurate ancestral sequences than ML methods (Hanson-Smith, Mol Biol Evol, 2010). Finally, given the greater computational expense of Bayesian methods – potentially weeks of computation using our resources with an alignment containing several hundred concatenated sequences – we elected to prioritize broader sequence sampling and used ML.

• Figure 3A – please show the y-axis in the log.

The plot for Figure 3A has been edited to show the y-axis in log.

• Could you please represent catalysis in units for kcat and KM? This would make the kinetics easier to compare to other enzymes, both extant and reconstructed.

We have used “specific activity,” (units of nmol product/nmol protein/s), which is kcat. Nitrogenase is a complex system for which Km is not an appropriate metric. A brief description is in the text and further details can be found in Harris et al., Biochemistry, 2019; 2022.

• 4C please include WT.

We have included data for WT in a revised Figure 4C.

• In light of the conclusion about a highly conserved mechanism, I would like to see more nuanced mechanistic studies. Pre-steady state kinetics; MD; pH studies.

Our central finding concerns a major aspect of nitrogenase mechanism (i.e., the reductive elimination/oxidative addition, “re/oa”, model for N2 binding and reduction). These conclusions are not dependent on a more exhaustive investigation of other properties relating to nitrogenase mechanism. We direct the reviewer to Harris et al., Biochemistry, 2019; 2022 which provide a more in-depth description of the re/oa model and kinetics, which is also cited in the main text.

• Please infer the age of the hypothetical ancestors expressing anc1 and anc2.

We did not perform time calibrations for our nitrogenase phylogeny, and age estimates based on species divergence are challenged by horizontal gene transfer (e.g., Parsons et al., Geobiology, 2021). Previous studies have estimated the timing of nitrogenase emergence (Parsons et al., Geobiology, 2021, Boyd et al., Geobiology, 2011), but more detailed constraints for the timeline targeted here are not presently available. We look forward to performing our own analysis in a future study.

• Is there a notable difference between how the all-anc complexes (Anc1B) interact as opposed to the hybrids (Anc1A, Anc2)? Is there a notable difference in melting temperature or oligomeric state? This is discussed to a degree in the paragraph beginning line 323, but many of the statements are general and do not posit molecular explanations. What do the authors mean by "historical" amino acid substitutions?

Our ability to copurify NifDK complexes from all ancestors (Anc1A, Anc1B, Anc2) and their exhibited activity in vitro suggests that there is no substantial distinction in their oligomeric states (i.e., they all form NifDK heterotetramers that interact with the NifH homodimer during the catalytic cycle). We did not investigate their melting temperatures.

Nevertheless, we do highlight the organism-level phenotypic differences observed between Anc1A and Anc1B strains, which likely stem from the additional substitutions (relative to WT) in the NifH and NifK proteins of Anc1B, which are not present in Anc1A. We interpret that these phenotypic outcomes might result from perturbed interactions within the complex (though evidently not enough to change the core oligomeric state), which can be more deeply explored in future work. We have now clarified these inferences in the Discussion.

By historical substitutions, we refer to differences between WT and ancestral amino acids at a given site. Upon reviewer’s feedback, we think that “ancestral substitution” might better capture this intended meaning. We have replaced this terminology and defined it at its first occurrence in the article.

• The discussion notes the surprising conservation of inhibition by H+ despite "substantial residue-level changes to the peripheral nitrogenase structure, as well as a handful within relatively conserved, active-site or protein-interface regions within the enzyme complex" Please elaborate on this, with specific attention to the active site. Are the residues involved/mechanisms known? Specify the mutations, how chemically conserved are they?

Our original Results text included a description of the sequence level differences observed between WT and Anc1/Anc2 ancestors. There, we highlight specific residues in functionally significant regions (e.g., I355V in the active site, within a loop considered important for cofactor insertion; F429Y and R108K within the NifD:NifK interface). A global visualization of sequence level differences relative to their conservation is shown in SI Appendix, Figure S3, and a listing of substitutions at relatively conserved positions is found in SI Appendix, Table S1. To our knowledge, the functional significance of many of these amino acid sites are not well characterized.

• I would be interested to see a discussion of potential ancestral promoters and expression levels. Expression levels are mentioned briefly in the results. I know this experiment must be compatible with the biological system used for the experiment (i.e. it would be impractical to also reconstruct ancestral promoters), but do the authors speculate that an ancestral nitrogenase would have been overexpressed to compensate for lower efficiency or that N2 fixation would merely have been rate-limited?

We agree that ancestral protein expression is worthy of study, though, to date, not well explored. As the reviewer suggests, we (authors Garcia, Kaçar) recently reported that expression levels of an ancestral and less active RuBisCO enzyme are increased relative to the WT enzyme (Kedzior et al., Cell Reports, 2022). However, we don’t see strong evidence of nitrogenase overexpression in the present study (Figure 3D).

Since expression levels in our experiments are dictated by regulatory mechanisms possessed by our model, Azotobacter, it’s challenging to infer whether our specific results would represent ancestral protein expression levels in an ancient host organism. We can speculate that expression of an ancestral, Mo-dependent nitrogenase might itself be ultimately limited by Mo availability, particularly early in Earth history when bulk marine Mo concentrations were extremely low (e.g., Anbar, Science, 2008). We envision that these possibilities can be tested in future genome engineering studies building off our presented experimental system.

We have acknowledged the role of protein expression in shaping ancestral phenotypes in additional text within the Discussion.

• For the uninitiated, please clearly introduce the evolution of metal ion dependence in nitrogenases.

We assume the reviewer is referring to a statement mentioning evolution of metal dependence in nitrogenases in our Introduction and agree that diversity of nitrogenase metal dependence is not well introduced. We now include additional text earlier in the introduction that describes the variability of nitrogenase metal dependence (i.e., relying on Mo, V, and Fe).

• How do your results compare to enzyme reconstructions of a similar "age"?

To our knowledge, our study represents the first laboratory reconstruction of nitrogenase enzymes. Therefore, we cannot yet compare our nitrogenase reconstructions to others of similar age.

• I find these results unsurprising. There is sufficient ASR literature to predict that a reconstructed enzyme will have comparable activity to the one or two extant enzymes it is compared to. When using a clade of conserved enzymes using a conserved mechanism, it is not surprising the ancestor conforms. What have we learned? It would be more interesting to probe these reconstructions for promiscuous activities or additional inhibitors. Are they easier to evolve than extant enzymes? If you reconstruct the whole pathway, does it behave differently? Does it act inefficiently and leak metabolites? Are other methods of fixation conceivable?

(Some of this text was also provided in response to similar comments from other reviewers).

We appreciate the reviewers’ comment and have revised relevant sections in our manuscript to clarify the novelty our study and include additional nuance into our discussions of conservation and constraints.

Nitrogenases are deep time enzymes and are a challenging target for engineering and functional study due to the number of protein components involved, their interactions with a broader cellular network, and their oxygen sensitivity (we now elaborate on these points in our Introduction). To date, only three, modern nitrogenase enzymes have recently been characterized with respect to their specific mechanism for N2 binding and reduction (the “reversible reductive elimination/oxidative addition” mechanism described in our article) (Harris et al., Biochemistry, 2019). Our study is therefore the first demonstration of this mechanism in nitrogenase ancestors and effectively doubles the number of nitrogenases that have been characterized to this degree.

Our broader point centers on the implications of this conservation in the evolution of biological nitrogen fixation strategies. Only one family of nitrogenase enzymes has evolved and survived to the present day. The comparable mechanistic features of extant nitrogenases and, now, ancestral nitrogenases, suggests that this one family has not only catalyzed N2 reduction for billions of years, but has achieved this incredibly challenging reaction in the same, specific manner.

The ecological importance of biological nitrogen fixation is on par with carbon fixation, though there are at least seven known pathways for achieving the latter. How life had become constrained to this particular N2 reduction mechanism remains an open question (particularly given several strategies for abiotic nitrogen fixation e.g., Cherkasov et al., Chem Engineer Process 2015; Dorr, Angew Chem Int Ed, 2003; Yung and McElroy, Science, 1979), but is one that can be further explored with the experimental approach presented here.

We agree that other aspects of nitrogenase functionality, including promiscuous activities, evolvability, and its specific interactions with other proteins involved in the nitrogen fixation pathway, are all excellent research targets that can also leverage our paleomolecular approach.

We have expanded our Discussion to include these key points.

Reviewer #3 (Recommendations for the authors):

I have several questions and suggestions for the phylogenetic analysis that I do not believe will alter any of the results but may help with presentation or for a better understanding of the uncertainty with them.

• Why were these particular nodes picked for reconstruction? Are they the highest confidence reconstructions?

The selected nodes are indeed well-supported (SH-like aLRT values ≈ 98-100). We have included this information in the Materials and methods. In addition, we chose a relatively conservative percentage identify threshold based on our (Garcia, Kaçar) prior resurrection studies to ensure we could recover active nitrogenase ancestors in our experimental system. We look forward to publishing the results of our ongoing work on older ancestral nodes across the nitrogenase tree.

We have updated the Results and Materials and methods with additional node selection rationale.

• Why was LG+G+F used? Was any model testing done?

Model testing performed by ModelFinder in IQ-TREE, now specified in Materials and methods.

• What tool was used to align the sequences?

MAFFT v7.450, now specified in Materials and methods.

• Why was ASR performed with RAxML and PAML both?

As we described in a response to Reviewer 1, we initially performed ASR with RAxML v.8, and constructed strains Anc1A and Anc2. However, due to concerns that this version of RAxML does not implement full marginal ancestral reconstruction, we repeated the analysis with PAML, constructing strains Anc1B (equivalent node to Anc1A). As we describe in our main text, Anc1A and Anc1B have the same set of descendant proteins and their NifD amino acid sequences are 95% identical.

• What is the difference between the forak023 and forak013-14 files in the github repository? The topology of tree_forak013-14_branchsupport.tre seems to be the one in the manuscript, but there are no branch support values in that file.

We thank the reviewer for noting this issue. We have updated the Github repository files and provided filenames that are more easily identifiable.

As mentioned in the public comments, I worry whether HGT may cause issues during tree inference. I believe the simplest way to find out would be to reconstruct a gene tree for each individual nif gene and see how the trees differ. It could also be worth examining whether these or the tree in the manuscript agree with bacterial phylogenies.

Indeed, as the reviewer points out, nitrogenase evolution is known to have been affected by significant HGT (e.g., Raymond et al., MBE, 2004; Parsons et al., Geobiology, 2021). We would therefore not expect protein trees and bacterial/archaeal species trees to agree, given our dataset. Discordance between the protein and species trees shouldn’t impact our main conclusions since we are not drawing inferences concerning species evolution from our protein tree. Previous work has also demonstrated that the individual nitrogenase genes have followed similar evolutionary trajectories (Raymond et al., MBE, 2004; Garcia et al., Genome Biol Evol, 2022), with the exception of H-subunit genes from V-nitrogenases (VnfH), some of which have diverged recently from molybdenum-dependent genes (NifH). However, given that (1) we are targeting a specific Nif lineage, (2) only one of our reconstructions includes an ancestral NifH gene, and (3) we nevertheless observe phenotypic consistency across multiple reconstructions, we do not expect issues stemming from HGT to have significantly impacted the present results.

Nevertheless, we concede that HGT is an important aspect of nitrogenase evolution and should be thoughtfully considered in ours and others’ future work, particularly for reconstructions that extend deeper in the nitrogenase phylogeny.

The randomly-selected codon (de)optimization process is a nice inclusion. I was a little unsure about how it was performed, though – were codons swapped randomly until some metric was reached, or a fixed number, or some other procedure? Given the other controls, I do not expect this to change any results, but it would be a nice method for other groups to potentially use.

We thank the reviewer for their comment and agree that providing more details on our codon optimization process would be helpful for the community. We used a semi-randomized optimization process to maximize identify between WT and ancestral nucleotide sequences. Ancestral and WT A. vinelandii protein sequences were compared using the alignment output of ancestral sequence reconstruction (RAxML or PAML). For sites where the ancestral and WT residues were identical, the WT codon was assigned. At sites where the residues were different, the codon was assigned randomly, weighted by A. vinelandii codon frequencies (Codon Usage Database, https://www.kazusa.or.jp/codon/). We have elaborated on our strategy in the Materials and methods and have added the relevant script to the public Kaçar Lab GitHub repository.

Given that the WT and ancestral sequences are 83% identical or higher (roughly akin to humans and mice), is the result that function is conserved surprising? Is it possible that these results say more about sequence (and functional) conservation, rather than a constraint? The UMAP embedding is an interesting approach, but makes this point in a different way, as the ancestral and WT sequences are extremely close in UMAP space. I believe some discussion of sequence.

(Some of this text was also provided in response to similar comments from other reviewers. It also seems that the reviewer comment here was truncated at the end. We’ve done our best to respond to what we expect were the reviewer’s intended recommendations).

We appreciate the reviewers’ comment and have revised relevant sections in our manuscript to clarify the novelty our study and include additional nuance into our discussions of conservation and constraints.

As the reviewer notes in their public comments, nitrogenases are indeed a challenging target for engineering and functional study due to the number of protein components involved, their interactions with a broader cellular network, and their oxygen sensitivity (we now elaborate on these points in our Introduction). To date, only three, modern nitrogenase enzymes have recently been characterized with respect to their specific mechanism for N2 binding and reduction (the “reversible reductive elimination/oxidative addition” mechanism described in our article) (Harris et al., Biochemistry, 2019). Our study is therefore the first demonstration of this mechanism in nitrogenase ancestors and effectively doubles the number of nitrogenases that have been characterized to this degree.

We agree that mechanistic conservation is likely an outcome of underlying sequence conservation. Our broader point centers on the implications of this conservation in the evolution of biological nitrogen fixation strategies. Only one family of nitrogenase enzymes has evolved and survived to the present day. The comparable mechanistic features of extant nitrogenases and, now, ancestral nitrogenases, suggests that this one family has not only catalyzed N2 reduction for billions of years, but has achieved this incredibly challenging reaction in the same, specific manner.

The ecological importance of biological nitrogen fixation is on par with carbon fixation, though there are at least seven known pathways for achieving the latter. How life had become constrained to this particular N2 reduction mechanism remains an open question (particularly given several strategies for abiotic nitrogen fixation e.g., Cherkasov et al., Chem Engineer Process 2015; Dorr, Angew Chem Int Ed, 2003; Yung and McElroy, Science, 1979), but is one that can be further explored with the experimental approach presented here.

We have expanded our Discussion to include these key points.

https://doi.org/10.7554/eLife.85003.sa2

Article and author information

Author details

  1. Amanda K Garcia

    Department of Bacteriology, University of Wisconsin–Madison, Madison, United States
    Contribution
    Conceptualization, Data curation, Formal analysis, Validation, Investigation, Visualization, Writing – original draft, Writing – review and editing
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-1936-2568

  2. Derek F Harris

    Department of Chemistry and Biochemistry, Utah State University, Logan, United States
    Contribution
    Data curation, Formal analysis, Investigation, Visualization, Writing – review and editing
    Competing interests
    No competing interests declared

  3. Alex J Rivier

    Department of Bacteriology, University of Wisconsin–Madison, Madison, United States
    Contribution
    Data curation, Formal analysis, Investigation
    Competing interests
    No competing interests declared

  4. Brooke M Carruthers

    Department of Bacteriology, University of Wisconsin–Madison, Madison, United States
    Contribution
    Data curation, Formal analysis, Investigation
    Competing interests
    No competing interests declared

  5. Azul Pinochet-Barros

    Department of Bacteriology, University of Wisconsin–Madison, Madison, United States
    Contribution
    Data curation
    Competing interests
    No competing interests declared

  6. Lance C Seefeldt

    Department of Chemistry and Biochemistry, Utah State University, Logan, United States
    Contribution
    Resources, Supervision, Writing – review and editing
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-6457-9504

  7. Betül Kaçar

    Department of Bacteriology, University of Wisconsin–Madison, Madison, United States
    Contribution
    Conceptualization, Resources, Data curation, Formal analysis, Supervision, Funding acquisition, Investigation, Methodology, Writing – original draft, Project administration, Writing – review and editing
    For correspondence
    bkacar@wisc.edu
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-0482-2357

Funding

National Aeronautics and Space Administration (19- ICAR19_2-0007)

  • Amanda K Garcia
  • Derek F Harris
  • Alex J Rivier
  • Brooke M Carruthers
  • Lance Seefeldt
  • Betül Kaçar
  • Azul Pinochet-Barros

John Templeton Foundation (61926)

  • Betül Kaçar

National Science Foundation (2228495)

  • Betül Kaçar
  • Alex J Rivier
  • Brooke M Carruthers
  • Amanda K Garcia

University of Wisconsin-Madison

  • Betül Kaçar

Arizona Space Grant Consortium

  • Brooke M Carruthers

National Aeronautics and Space Administration (80NSSC19K1617)

  • Betül Kaçar

The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.

Acknowledgements

We thank Dennis Dean and Valerie Cash for providing A. vinelandii strains DJ, DJ2278, DJ2102, and DJ884 and for guidance in genomic manipulations; Jean-Michel Ané, April MacIntyre, and Junko Maeda for guidance and instrumentation support in performing the in vivo acetylene reduction assays; Bruno Cuevas for assistance with visualizing nitrogenase sequence space; and the members of the Metal Selection and Utilization Across Eons (MUSE) Consortium for helpful suggestions and discussions. This research was supported by the National Aeronautics and Space Administration (NASA) Interdisciplinary Consortium for Astrobiology Research: Metal Utilization and Selection Across Eons, MUSE (19- ICAR19_2–0007), the University of Wisconsin-Madison College of Agricultural and Life Sciences, the NASA Postdoctoral Program (AKG), NASA Arizona Space Grant (BMC), the John Templeton Foundation (BK; 61926), the National Science Foundation (BK; 2228495), the NASA Early Career Faculty Award (BK), and the Hypothesis Fund Award (BK).

Senior and Reviewing Editor

  1. Christian R Landry, Université Laval, Canada

Reviewers

  1. Matilda Newton, University of Colorado Boulder, United States
  2. Christian B Macdonald, University of California, San Francisco, San Francisco, United States

Publication history

  1. Preprint posted: May 19, 2022 (view preprint)
  2. Received: November 18, 2022
  3. Accepted: February 16, 2023
  4. Accepted Manuscript published: February 17, 2023 (version 1)
  5. Version of Record published: March 1, 2023 (version 2)

Copyright

© 2023, Garcia et al.

This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.

Metrics

  • 1,060
    Page views
  • 163
    Downloads
  • 0
    Citations

Article citation count generated by polling the highest count across the following sources: Crossref, PubMed Central, Scopus.

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

  1. Amanda K Garcia
  2. Derek F Harris
  3. Alex J Rivier
  4. Brooke M Carruthers
  5. Azul Pinochet-Barros
  6. Lance C Seefeldt
  7. Betül Kaçar
(2023)
Nitrogenase resurrection and the evolution of a singular enzymatic mechanism
eLife 12:e85003.
https://doi.org/10.7554/eLife.85003

Categories and tags

Research organism

Further reading

    1. Biochemistry and Chemical Biology
    2. Cancer Biology

    Multi-targeted therapy resistance via drug-induced secretome fucosylation

    Mark Borris D Aldonza, Junghwa Cha ... Yoosik Kim
    Research Article

    Cancer secretome is a reservoir for aberrant glycosylation. How therapies alter this post- translational cancer hallmark and the consequences thereof remain elusive. Here we show that an elevated secretome fucosylation is a pan-cancer signature of both response and resistance to multiple targeted therapies. Large-scale pharmacogenomics revealed that fucosylation genes display widespread association with resistance to these therapies. In cancer cell cultures, xenograft mouse models, and patients, targeted kinase inhibitors distinctively induced core fucosylation of secreted proteins less than 60 kDa. Label-free proteomics of N-glycoproteomes identified fucosylation of the antioxidant PON1 as a critical component of the therapy-induced secretome (TIS). N-glycosylation of TIS and target core fucosylation of PON1 are mediated by the fucose salvage-FUT8-SLC35C1 axis with PON3 directly modulating GDP-Fuc transfer on PON1 scaffolds. Core fucosylation in the Golgi impacts PON1 stability and folding prior to secretion, promoting a more degradation-resistant PON1. Global and PON1-specific secretome de-N-glycosylation both limited the expansion of resistant clones in a tumor regression model. We defined the resistance-associated transcription factors (TFs) and genes modulated by the N-glycosylated TIS via a focused and transcriptome-wide analyses. These genes characterize the oxidative stress, inflammatory niche, and unfolded protein response as important factors for this modulation. Our findings demonstrate that core fucosylation is a common modification indirectly induced by targeted therapies that paradoxically promotes resistance.

    1. Biochemistry and Chemical Biology
    2. Developmental Biology

    O-GlcNAc glycosylation orchestrates fate decision and niche function of bone marrow stromal progenitors

    Zengdi Zhang, Zan Huang ... Hai-Bin Ruan
    Research Article Updated

    In mammals, interactions between the bone marrow (BM) stroma and hematopoietic progenitors contribute to bone-BM homeostasis. Perinatal bone growth and ossification provide a microenvironment for the transition to definitive hematopoiesis; however, mechanisms and interactions orchestrating the development of skeletal and hematopoietic systems remain largely unknown. Here, we establish intracellular O-linked β-N-acetylglucosamine (O-GlcNAc) modification as a posttranslational switch that dictates the differentiation fate and niche function of early BM stromal cells (BMSCs). By modifying and activating RUNX2, O-GlcNAcylation promotes osteogenic differentiation of BMSCs and stromal IL-7 expression to support lymphopoiesis. In contrast, C/EBPβ-dependent marrow adipogenesis and expression of myelopoietic stem cell factor (SCF) is inhibited by O-GlcNAcylation. Ablating O-GlcNAc transferase (OGT) in BMSCs leads to impaired bone formation, increased marrow adiposity, as well as defective B-cell lymphopoiesis and myeloid overproduction in mice. Thus, the balance of osteogenic and adipogenic differentiation of BMSCs is determined by reciprocal O-GlcNAc regulation of transcription factors, which simultaneously shapes the hematopoietic niche.

    1. Biochemistry and Chemical Biology
    2. Neuroscience

    Glia-neuron coupling via a bipartite sialylation pathway promotes neural transmission and stress tolerance in Drosophila

    Hilary Scott, Boris Novikov ... Vladislav Panin
    Research Article

    Modification by sialylated glycans can affect protein functions, underlying mechanisms that control animal development and physiology. Sialylation relies on a dedicated pathway involving evolutionarily conserved enzymes, including CMP-sialic acid synthetase (CSAS) and sialyltransferase (SiaT) that mediate the activation of sialic acid and its transfer onto glycan termini, respectively. In Drosophila, CSAS and DSiaT genes function in the nervous system, affecting neural transmission and excitability. We found that these genes function in different cells: the function of CSAS is restricted to glia, while DSiaT functions in neurons. This partition of the sialylation pathway allows for regulation of neural functions via a glia-mediated control of neural sialylation. The sialylation genes were shown to be required for tolerance to heat and oxidative stress and for maintenance of the normal level of voltage-gated sodium channels. Our results uncovered a unique bipartite sialylation pathway that mediates glia-neuron coupling and regulates neural excitability and stress tolerance.