Emergence of fractal geometries in the evolution of a metabolic enzyme.

Fractals are patterns that are self-similar across multiple length-scales1. Macroscopic fractals are common in nature2-4; however, so far, molecular assembly into fractals is restricted to synthetic systems5-12. Here we report the discovery of a natural protein, citrate synthase from the cyanobacterium Synechococcus elongatus, which self-assembles into Sierpinski triangles. Using cryo-electron microscopy, we reveal how the fractal assembles from a hexameric building block. Although different stimuli modulate the formation of fractal complexes and these complexes can regulate the enzymatic activity of citrate synthase in vitro, the fractal may not serve a physiological function in vivo. We use ancestral sequence reconstruction to retrace how the citrate synthase fractal evolved from non-fractal precursors, and the results suggest it may have emerged as a harmless evolutionary accident. Our findings expand the space of possible protein complexes and demonstrate that intricate and regulatable assemblies can evolve in a single substitution.

Evolution-inspired engineering of nonribosomal peptide synthetases.

Many clinically used drugs are derived from or inspired by bacterial natural products that often are produced through nonribosomal peptide synthetases (NRPSs), megasynthetases that activate and join individual amino acids in an assembly line fashion. In this work, we describe a detailed phylogenetic analysis of several bacterial NRPSs that led to the identification of yet undescribed recombination sites within the thiolation (T) domain that can be used for NRPS engineering. We then developed an evolution-inspired "eXchange Unit between T domains" (XUT) approach, which allows the assembly of NRPS fragments over a broad range of GC contents, protein similarities, and extender unit specificities, as demonstrated for the specific production of a proteasome inhibitor designed and assembled from five different NRPS fragments.

Two distinct ferredoxins are essential for nitrogen fixation by the iron nitrogenase in Rhodobacter capsulatus.

Nitrogenases are the only enzymes able to fix gaseous nitrogen into bioavailable ammonia and hence are essential for sustaining life. Catalysis by nitrogenases requires both a large amount of ATP and electrons donated by strongly reducing ferredoxins or flavodoxins. Our knowledge about the mechanisms of electron transfer to nitrogenase enzymes is limited: The electron transport to the iron (Fe)-nitrogenase has hardly been investigated. Here, we characterized the electron transfer pathway to the Fe-nitrogenase in Rhodobacter capsulatus via proteome analyses, genetic deletions, complementation studies, and phylogenetics. Proteome analyses revealed an upregulation of four ferredoxins under nitrogen-fixing conditions reliant on the Fe-nitrogenase in a molybdenum nitrogenase knockout strain, compared to non-nitrogen-fixing conditions. Based on these findings, R. capsulatus strains with deletions of ferredoxin (fdx) and flavodoxin (fld, nifF) genes were constructed to investigate their roles in nitrogen fixation by the Fe-nitrogenase. R. capsulatus deletion strains were characterized by monitoring diazotrophic growth and Fe-nitrogenase activity in vivo. Only deletions of fdxC or fdxN resulted in slower growth and reduced Fe-nitrogenase activity, whereas the double deletion of both fdxC and fdxN abolished diazotrophic growth. Differences in the proteomes of ∆fdxC and ∆fdxN strains, in conjunction with differing plasmid complementation behaviors of fdxC and fdxN, indicate that the two Fds likely possess different roles and functions. These findings will guide future engineering of the electron transport systems to nitrogenase enzymes, with the aim of increased electron flux and product formation.IMPORTANCENitrogenases are essential for biological nitrogen fixation, converting atmospheric nitrogen gas to bioavailable ammonia. The production of ammonia by diazotrophic organisms, harboring nitrogenases, is essential for sustaining plant growth. Hence, there is a large scientific interest in understanding the cellular mechanisms for nitrogen fixation via nitrogenases. Nitrogenases rely on highly reduced electrons to power catalysis, although we lack knowledge as to which proteins shuttle the electrons to nitrogenases within cells. Here, we characterized the electron transport to the iron (Fe)-nitrogenase in the model diazotroph Rhodobacter capsulatus, showing that two distinct ferredoxins are very important for nitrogen fixation despite having different redox centers. In addition, our research expands upon the debate on whether ferredoxins have functional redundancy or perform distinct roles within cells. Here, we observe that both essential ferredoxins likely have distinct roles based on differential proteome shifts of deletion strains and different complementation behaviors.

The virulence regulator VirB from Shigella flexneri uses a CTP-dependent switch mechanism to activate gene expression.

The transcriptional antisilencer VirB acts as a master regulator of virulence gene expression in the human pathogen Shigella flexneri. It binds DNA sequences (virS) upstream of VirB-dependent promoters and counteracts their silencing by the nucleoid-organizing protein H-NS. However, its precise mode of action remains unclear. Notably, VirB is not a classical transcription factor but related to ParB-type DNA-partitioning proteins, which have recently been recognized as DNA-sliding clamps using CTP binding and hydrolysis to control their DNA entry gate. Here, we show that VirB binds CTP, embraces DNA in a clamp-like fashion upon its CTP-dependent loading at virS sites and slides laterally on DNA after clamp closure. Mutations that prevent CTP-binding block VirB loading in vitro and abolish the formation of VirB nucleoprotein complexes as well as virulence gene expression in vivo. Thus, VirB represents a CTP-dependent molecular switch that uses a loading-and-sliding mechanism to control transcription during bacterial pathogenesis.

Ion binding with charge inversion combined with screening modulates DEAD box helicase phase transitions.

Membraneless organelles, or biomolecular condensates, enable cells to compartmentalize material and processes into unique biochemical environments. While specific, attractive molecular interactions are known to stabilize biomolecular condensates, repulsive interactions, and the balance between these opposing forces, are largely unexplored. Here, we demonstrate that repulsive and attractive electrostatic interactions regulate condensate stability, internal mobility, interfaces, and selective partitioning of molecules both in vitro and in cells. We find that signaling ions, such as calcium, alter repulsions between model Ddx3 and Ddx4 condensate proteins by directly binding to negatively charged amino acid sidechains and effectively inverting their charge, in a manner fundamentally dissimilar to electrostatic screening. Using a polymerization model combined with generalized stickers and spacers, we accurately quantify and predict condensate stability over a wide range of pH, salt concentrations, and amino acid sequences. Our model provides a general quantitative treatment for understanding how charge and ions reversibly control condensate stability.

Molecular basis and design principles of switchable front-rear polarity and directional migration in Myxococcus xanthus.

During cell migration, front-rear polarity is spatiotemporally regulated; however, the underlying design of regulatory interactions varies. In rod-shaped Myxococcus xanthus cells, a spatial toggle switch dynamically regulates front-rear polarity. The polarity module establishes front-rear polarity by guaranteeing front pole-localization of the small GTPase MglA. Conversely, the Frz chemosensory system, by acting on the polarity module, causes polarity inversions. MglA localization depends on the RomR/RomX GEF and MglB/RomY GAP complexes that localize asymmetrically to the poles by unknown mechanisms. Here, we show that RomR and the MglB and MglC roadblock domain proteins generate a positive feedback by forming a RomR/MglC/MglB complex, thereby establishing the rear pole with high GAP activity that is non-permissive to MglA. MglA at the front engages in negative feedback that breaks the RomR/MglC/MglB positive feedback allosterically, thus ensuring low GAP activity at this pole. These findings unravel the design principles of a system for switchable front-rear polarity.

Functional Degeneracy in Paracoccus denitrificans Pd1222 Is Coordinated via RamB, Which Links Expression of the Glyoxylate Cycle to Activity of the Ethylmalonyl-CoA Pathway.

Metabolic degeneracy describes the phenomenon that cells can use one substrate through different metabolic routes, while metabolic plasticity, refers to the ability of an organism to dynamically rewire its metabolism in response to changing physiological needs. A prime example for both phenomena is the dynamic switch between two alternative and seemingly degenerate acetyl-CoA assimilation routes in the alphaproteobacterium Paracoccus denitrificans Pd1222: the ethylmalonyl-CoA pathway (EMCP) and the glyoxylate cycle (GC). The EMCP and the GC each tightly control the balance between catabolism and anabolism by shifting flux away from the oxidation of acetyl-CoA in the tricarboxylic acid (TCA) cycle toward biomass formation. However, the simultaneous presence of both the EMCP and GC in P. denitrificans Pd1222 raises the question of how this apparent functional degeneracy is globally coordinated during growth. Here, we show that RamB, a transcription factor of the ScfR family, controls expression of the GC in P. denitrificans Pd1222. Combining genetic, molecular biological and biochemical approaches, we identify the binding motif of RamB and demonstrate that CoA-thioester intermediates of the EMCP directly bind to the protein. Overall, our study shows that the EMCP and the GC are metabolically and genetically linked with each other, demonstrating a thus far undescribed bacterial strategy to achieve metabolic plasticity, in which one seemingly degenerate metabolic pathway directly drives expression of the other. IMPORTANCE Carbon metabolism provides organisms with energy and building blocks for cellular functions and growth. The tight regulation between degradation and assimilation of carbon substrates is central for optimal growth. Understanding the underlying mechanisms of metabolic control in bacteria is of importance for applications in health (e.g., targeting of metabolic pathways with new antibiotics, development of resistances) and biotechnology (e.g., metabolic engineering, introduction of new-to-nature pathways). In this study, we use the alphaproteobacterium P. denitrificans as model organism to study functional degeneracy, a well-known phenomenon of bacteria to use the same carbon source through two different (competing) metabolic routes. We demonstrate that two seemingly degenerate central carbon metabolic pathways are metabolically and genetically linked with each other, which allows the organism to control the switch between them in a coordinated manner during growth. Our study elucidates the molecular basis of metabolic plasticity in central carbon metabolism, which improves our understanding of how bacterial metabolism is able to partition fluxes between anabolism and catabolism.

A bacterial tungsten-containing aldehyde oxidoreductase forms an enzymatic decorated protein nanowire.

Aldehyde oxidoreductases (AORs) are tungsten enzymes catalyzing the oxidation of many different aldehydes to the corresponding carboxylic acids. In contrast to other known AORs, the enzyme from the denitrifying betaproteobacterium Aromatoleum aromaticum (AORAa) consists of three different subunits (AorABC) and uses nicotinamide adenine dinucleotide (NAD) as an electron acceptor. Here, we reveal that the enzyme forms filaments of repeating AorAB protomers that are capped by a single NAD-binding AorC subunit, based on solving its structure via cryo-electron microscopy. The polyferredoxin-like subunit AorA oligomerizes to an electron-conducting nanowire that is decorated with enzymatically active and W-cofactor (W-co) containing AorB subunits. Our structure further reveals the binding mode of the native substrate benzoate in the AorB active site. This, together with quantum mechanics:molecular mechanics (QM:MM)-based modeling for the coordination of the W-co, enables formulation of a hypothetical catalytic mechanism that paves the way to further engineering for applications in synthetic biology and biotechnology.

Ancestral reconstruction of the MotA stator subunit reveals that conserved residues far from the pore are required to drive flagellar motility.

The bacterial flagellar motor (BFM) is a rotary nanomachine powered by the translocation of ions across the inner membrane through the stator complex. The stator complex consists of two membrane proteins: MotA and MotB (in H+-powered motors), or PomA and PomB (in Na+-powered motors). In this study, we used ancestral sequence reconstruction (ASR) to probe which residues of MotA correlate with function and may have been conserved to preserve motor function. We reconstructed 10 ancestral sequences of MotA and found four of them were motile in combination with contemporary Escherichia coli MotB and in combination with our previously published functional ancestral MotBs. Sequence comparison between wild-type (WT) E. coli MotA and MotA-ASRs revealed 30 critical residues across multiple domains of MotA that were conserved among all motile stator units. These conserved residues included pore-facing, cytoplasm-facing, and MotA-MotA intermolecular facing sites. Overall, this work demonstrates the role of ASR in assessing conserved variable residues in a subunit of a molecular complex.

Fortuitously compatible protein surfaces primed allosteric control in cyanobacterial photoprotection.

Highly specific interactions between proteins are a fundamental prerequisite for life, but how they evolve remains an unsolved problem. In particular, interactions between initially unrelated proteins require that they evolve matching surfaces. It is unclear whether such surface compatibilities can only be built by selection in small incremental steps, or whether they can also emerge fortuitously. Here, we used molecular phylogenetics, ancestral sequence reconstruction and biophysical characterization of resurrected proteins to retrace the evolution of an allosteric interaction between two proteins that act in the cyanobacterial photoprotection system. We show that this interaction between the orange carotenoid protein (OCP) and its unrelated regulator, the fluorescence recovery protein (FRP), evolved when a precursor of FRP was horizontally acquired by cyanobacteria. FRP's precursors could already interact with and regulate OCP even before these proteins first encountered each other in an ancestral cyanobacterium. The OCP-FRP interaction exploits an ancient dimer interface in OCP, which also predates the recruitment of FRP into the photoprotection system. Together, our work shows how evolution can fashion complex regulatory systems easily out of pre-existing components.

Parameterization of a single H-bond in Orange Carotenoid Protein by atomic mutation reveals principles of evolutionary design of complex chemical photosystems.

Introduction: Dissecting the intricate networks of covalent and non-covalent interactions that stabilize complex protein structures is notoriously difficult and requires subtle atomic-level exchanges to precisely affect local chemical functionality. The function of the Orange Carotenoid Protein (OCP), a light-driven photoswitch involved in cyanobacterial photoprotection, depends strongly on two H-bonds between the 4-ketolated xanthophyll cofactor and two highly conserved residues in the C-terminal domain (Trp288 and Tyr201). Method: By orthogonal translation, we replaced Trp288 in Synechocystis OCP with 3-benzothienyl-L-alanine (BTA), thereby exchanging the imino nitrogen for a sulphur atom. Results: Although the high-resolution (1.8 Å) crystal structure of the fully photoactive OCP-W288_BTA protein showed perfect isomorphism to the native structure, the spectroscopic and kinetic properties changed distinctly. We accurately parameterized the effects of the absence of a single H-bond on the spectroscopic and thermodynamic properties of OCP photoconversion and reveal general principles underlying the design of photoreceptors by natural evolution. Discussion: Such "molecular surgery" is superior over trial-and-error methods in hypothesis-driven research of complex chemical systems.

Evidence for an early green/red photocycle that precedes the diversification of GAF domain photoreceptor cyanobacteriochromes.

Phytochromes are linear tetrapyrrole-binding photoreceptors in eukaryotes and bacteria, primarily responding to red and far-red light signals reversibly. Among the GAF domain-based phytochrome superfamily, cyanobacteria-specific cyanobacteriochromes show various optical properties covering the entire visible region. It is unknown what physiological demands drove the evolution of cyanobacteriochromes in cyanobacteria. Here, we utilize ancestral sequence reconstruction and biochemical verification to show that the resurrected ancestral cyanobacteriochrome proteins reversibly respond to green and red light signals. pH titration analyses indicate that the deprotonation of the bound phycocyanobilin chromophore is crucial to perceive green light. The ancestral cyanobacteriochromes show only modest thermal reversion to the green light-absorbing form, suggesting that they evolved to sense the incident green/red light ratio. Many cyanobacteria can utilize green light for photosynthesis using phycobilisome light-harvesting complexes. The green/red sensing cyanobacteriochromes may have allowed better acclimation to changing light environments by rearranging the absorption capacity of the phycobilisome through chromatic acclimation.

Structural basis for lipid and copper regulation of the ABC transporter MsbA.

A critical step in lipopolysaccharide (LPS) biogenesis involves flipping lipooligosaccharide, an LPS precursor, from the cytoplasmic to the periplasmic leaflet of the inner membrane, an operation carried out by the ATP-binding cassette transporter MsbA. Although LPS binding to the inner cavity of MsbA is well established, the selectivity of MsbA-lipid interactions at other site(s) remains poorly understood. Here we use native mass spectrometry (MS) to characterize MsbA-lipid interactions and guide structural studies. We show the transporter co-purifies with copper(II) and metal binding modulates protein-lipid interactions. A 2.15 Å resolution structure of an N-terminal region of MsbA in complex with copper(II) is presented, revealing a structure reminiscent of the GHK peptide, a high-affinity copper(II) chelator. Our results demonstrate conformation-dependent lipid binding affinities, particularly for the LPS-precursor, 3-deoxy-D-manno-oct-2-ulosonic acid (Kdo)2-lipid A (KDL). We report a 3.6 Å-resolution structure of MsbA trapped in an open, outward-facing conformation with adenosine 5'-diphosphate and vanadate, revealing a distinct KDL binding site, wherein the lipid forms extensive interactions with the transporter. Additional studies provide evidence that the exterior KDL binding site is conserved and a positive allosteric modulator of ATPase activity, serving as a feedforward activation mechanism to couple transporter activity with LPS biosynthesis.

ATPase Activity of Bacillus subtilis RecA Affects the Dynamic Formation of RecA Filaments at DNA Double Strand Breaks.

RecA plays a central role in DNA repair and is a main actor involved in homologous recombination (HR). In vivo, RecA forms filamentous structures termed "threads," which are essential for HR, but whose nature is still ill defined. We show that RecA from Bacillus subtilis having lower ATP binding activity can still form nucleoprotein filaments in vitro, features lower dsDNA binding activity, but still retains most of wild type RecA activity in vivo. Contrarily, loss of ATPase activity strongly reduced formation of nucleoprotein filaments in vitro, and effectivity to repair double strand breaks (DSBs) in vivo. In the presence of wild type RecA protein, additionally expressed RecA with lowered ATPbinding activity only moderately affected RecA dynamics, while loss of ATPase activity leads to a large reduction of the formation of threads, as well as of their dynamic changes observed in a seconds-scale. Single molecule tracking of RecA revealed incorporation of freely diffusing and nonspecifically DNA-bound molecules into threads upon induction of a single DSB. This change of dynamics was highly perturbed in the absence of ATPase activity, revealing that filamentous forms of RecA as well as their dynamics depend on ATPase activity. Based on the idea that ATPase activity of RecA is most important for DNA strand exchange activity, our data suggest that extension and retraction of threads due is to many local strand invasion events during the search for sequences homologous to the induced DNA break site. IMPORTANCE Single-strand (ss) DNA binding ATPase RecA is the central recombinase in homologous recombination, and therefore essential for DNA repair pathways involving DNA strand exchange reactions. In several bacterial, RecA forms filamentous structures along the long axis of cells after induction of double strand breaks (DSBs) in the chromosome. These striking assemblies likely reflect RecA/ssDNA nucleoprotein filaments, which can extend and remodel within a time frame of few minutes. We show that ATPase activity of RecA is pivotal for these dynamic rearrangements, which include recruitment of freely diffusing molecules into low-mobile molecules within filaments. Our data suggest that ssDNA binding- and unbinding reactions are at the heart of RecA dynamics that power the dynamics of subcellular filamentous assemblies, leading to strand exchange reactions over a distance of several micrometers.

Evolution of increased complexity and specificity at the dawn of form I Rubiscos.

The evolution of ribulose-1,5-bisphosphate carboxylase/oxygenases (Rubiscos) that discriminate strongly between their substrate carbon dioxide and the undesired side substrate dioxygen was an important event for photosynthetic organisms adapting to an oxygenated environment. We use ancestral sequence reconstruction to recapitulate this event. We show that Rubisco increased its specificity and carboxylation efficiency through the gain of an accessory subunit before atmospheric oxygen was present. Using structural and biochemical approaches, we retrace how this subunit was gained and became essential. Our work illuminates the emergence of an adaptation to rising ambient oxygen levels, provides a template for investigating the function of interactions that have remained elusive because of their essentiality, and sheds light on the determinants of specificity in Rubisco.

Simple mechanisms for the evolution of protein complexity.

Proteins are tiny models of biological complexity: specific interactions among their many amino acids cause proteins to fold into elaborate structures, assemble with other proteins into higher-order complexes, and change their functions and structures upon binding other molecules. These complex features are classically thought to evolve via long and gradual trajectories driven by persistent natural selection. But a growing body of evidence from biochemistry, protein engineering, and molecular evolution shows that naturally occurring proteins often exist at or near the genetic edge of multimerization, allostery, and even new folds, so just one or a few mutations can trigger acquisition of these properties. These sudden transitions can occur because many of the physical properties that underlie these features are present in simpler proteins as fortuitous by-products of their architecture. Moreover, complex features of proteins can be encoded by huge arrays of sequences, so they are accessible from many different starting points via many possible paths. Because the bridges to these features are both short and numerous, random chance can join selection as a key factor in explaining the evolution of molecular complexity.

Diadenosine tetraphosphate regulates biosynthesis of GTP in Bacillus subtilis.

Diadenosine tetraphosphate (Ap4A) is a putative second messenger molecule that is conserved from bacteria to humans. Nevertheless, its physiological role and the underlying molecular mechanisms are poorly characterized. We investigated the molecular mechanism by which Ap4A regulates inosine-5'-monophosphate dehydrogenase (IMPDH, a key branching point enzyme for the biosynthesis of adenosine or guanosine nucleotides) in Bacillus subtilis. We solved the crystal structure of BsIMPDH bound to Ap4A at a resolution of 2.45 Å to show that Ap4A binds to the interface between two IMPDH subunits, acting as the glue that switches active IMPDH tetramers into less active octamers. Guided by these insights, we engineered mutant strains of B. subtilis that bypass Ap4A-dependent IMPDH regulation without perturbing intracellular Ap4A pools themselves. We used metabolomics, which suggests that these mutants have a dysregulated purine, and in particular GTP, metabolome and phenotypic analysis, which shows increased sensitivity of B. subtilis IMPDH mutant strains to heat compared with wild-type strains. Our study identifies a central role for IMPDH in remodelling metabolism and heat resistance, and provides evidence that Ap4A can function as an alarmone.

GGDEF domain as spatial on-switch for a phosphodiesterase by interaction with landmark protein HubP.

In bacteria, the monopolar localization of enzymes and protein complexes can result in a bimodal distribution of enzyme activity between the dividing cells and heterogeneity of cellular behaviors. In Shewanella putrefaciens, the multidomain hybrid diguanylate cyclase/phosphodiesterase PdeB, which degrades the secondary messenger c-di-GMP, is located at the flagellated cell pole. Here, we show that direct interaction between the inactive diguanylate cyclase (GGDEF) domain of PdeB and the FimV domain of the polar landmark protein HubP is crucial for full function of PdeB as a phosphodiesterase. Thus, the GGDEF domain serves as a spatially controlled on-switch that effectively restricts PdeBs activity to the flagellated cell pole. PdeB regulates abundance and activity of at least two crucial surface-interaction factors, the BpfA surface-adhesion protein and the MSHA type IV pilus. The heterogeneity in c-di-GMP concentrations, generated by differences in abundance and timing of polar appearance of PdeB, orchestrates the population behavior with respect to cell-surface interaction and environmental spreading.

Non-adaptive complexity and biochemical function.

Intricate biochemical structures are usually thought to be useful, because natural selection preserves them from degradation by a constant hail of destructive mutations. Biochemists therefore often deliberately disrupt them to understand how complexity improves protein function or fitness. However, evolutionary theory suggests that even useless complexity that never improved fitness can become completely essential if a simple set of evolutionary conditions is fulfilled. We review evidence that stable protein complexes, protein-chaperone interactions, and complexes consisting of several paralogs all fulfill these conditions. This makes reverse genetics or destructive mutagenesis unsuitable for assigning functions to these kinds of complexity. Instead, we advocate that incorporating evolutionary approaches into biochemistry overcomes this difficulty and allows us to distinguish useless from useful biochemical complexity.

Structure and mechanistic features of the prokaryotic minimal RNase P.

Endonucleolytic removal of 5'-leader sequences from tRNA precursor transcripts (pre-tRNAs) by ribonuclease P (RNase P) is essential for protein synthesis. Beyond RNA-based RNase P enzymes, protein-only versions of the enzyme exert this function in various eukarya (there termed PRORPs) and in some bacteria (Aquifex aeolicus and close relatives); both enzyme types belong to distinct subgroups of the PIN domain metallonuclease superfamily. Homologs of Aquifex RNase P (HARPs) are also expressed in some other bacteria and many archaea, where they coexist with RNA-based RNase P and do not represent the main RNase P activity. Here, we solved the structure of the bacterial HARP from Halorhodospira halophila by cryo-electron microscopy, revealing a novel screw-like dodecameric assembly. Biochemical experiments demonstrate that oligomerization is required for RNase P activity of HARPs. We propose that the tRNA substrate binds to an extended spike-helix (SH) domain that protrudes from the screw-like assembly to position the 5'-end in close proximity to the active site of the neighboring dimer. The structure suggests that eukaryotic PRORPs and prokaryotic HARPs recognize the same structural elements of pre-tRNAs (tRNA elbow region and cleavage site). Our analysis thus delivers the structural and mechanistic basis for pre-tRNA processing by the prokaryotic HARP system.