The growth of metagenomics-derived amino acid sequence data has transformed our understanding of protein function, microbial diversity, and evolutionary relationships. However, the vast majority of these proteins remain functionally uncharacterized. Grouping the millions of such uncharacterized sequences with the few experimentally characterized ones allows the transfer of annotations, while the inspection of conserved residues with multiple sequence alignments can provide clues to function, even in the absence of existing functional information. To address the challenges associated with this data surge and the need to group sequences, we present a scalable, open-source, parametrizable Nextflow pipeline (nf-core/proteinfamilies) that generates nascent protein families or assigns new proteins to existing families. The computational benchmarks demonstrated that resource usage scales approximately linearly with input size, and the biological benchmarks showed that the generated protein families closely resemble manually curated families in widely used databases.

Cytochrome bd is a terminal oxidase expressed under low oxygen conditions and central for the survival of many pathogens. Here, we characterize the cyt bd-II from Mycobacterium smegmatis, a member of a hitherto uncharacterized evolutionary group (qOR-2) of bd oxidases, by combining biochemical studies with cryo-electron microscopy (cryo-EM), and multiscale simulations. Overexpressing the appCB operon in its native host led to production of a highly active bd-II (k_obs = 30 e⁻ s⁻¹) that together with a high-resolution (2.8 Å) cryo-EM structure and multiscale simulations reveal unique proton pathways and oxygen channels responsible for its function. We propose that a pH-dependent molecular switch, involving coordination changes of heme d and surrounding bulky residues regulate substrate access into the active site. Taken together, our findings provide detailed mechanistic insight of qOR-2 type bd oxidases, and a basis for understanding the evolution of the superfamily.

NrdR is a bacterial transcriptional repressor consisting of a zinc (Zn)-ribbon domain followed by an ATP-cone domain. Understanding its mechanism of action could aid the design of novel antibacterials. NrdR binds specifically to two “NrdR boxes” upstream of ribonucleotide reductase operons, of which Escherichia coli has three: nrdHIEF, nrdDG and nrdAB, in the last of which we identified a new box. We show that E. coli NrdR (EcoNrdR) has similar binding strength to all three sites when loaded with ATP plus deoxyadenosine triphosphate (dATP) or equivalent diphosphate combinations. No other combination of adenine nucleotides promotes binding to DNA. We present crystal structures of EcoNrdR–ATP–dATP and EcoNrdR–ADP–dATP, which are the first high-resolution crystal structures of an NrdR. We have also determined cryo-electron microscopy structures of DNA-bound EcoNrdR–ATP–dATP and novel filaments of EcoNrdR–ATP. Tetrameric forms of EcoNrdR involve alternating interactions between pairs of Zn-ribbon domains and ATP-cones. The structures reveal considerable flexibility in relative orientation of ATP-cones vs Zn-ribbon domains. The structure of DNA-bound EcoNrdR–ATP–dATP shows that significant conformational rearrangements between ATP-cones and Zn-ribbons accompany DNA binding while the ATP-cones retain the same relative orientation. In contrast, ATP-loaded EcoNrdR filaments show rearrangements of the ATP-cone pairs and sequester the DNA-binding residues of NrdR such that they are unable to bind to DNA. Our results, in combination with a previous structural and biochemical study, point to highly flexible EcoNrdR structures that, when loaded with the correct nucleotides, adapt to an optimal promoter-binding conformation.

Class I ribonucleotide reductases (RNRs) convert ribonucleotides into deoxyribonucleotides under oxic conditions. The R2 subunit provides a radical required for catalysis conducted by the R1 subunit. In most R2s the radical is generated on a tyrosine via oxidation by an adjacent metal site. The discovery of a metal-free R2 defined the new RNR subclass Ie. In R2e, three of the otherwise strictly conserved metal-binding glutamates in the active site are substituted. Two variants have been found, VPK and QSK. To date, the VPK version has been the focus of biochemical characterization. Here we characterize a QSK variant of R2e. We analyse the organismal distribution of the two R2e versions and find dozens of organisms relying solely on the QSK RNR for deoxyribonucleotide production. We demonstrate that the R2e_QSK of the human pathogen Gardnerella vaginalis (Bifidobacterium vaginale) modifies the active site-adjacent tyrosine to DOPA. The amount of modified protein is shown to be dependent on coexpression with the other proteins encoded in the RNR operon. The DOPA containing R2e_QSK can support ribonucleotide reduction in vitro while the unmodified protein cannot. Finally, we determined the first structures of R2e_QSK in the unmodified and DOPA states.

NrdR is a universal transcriptional repressor of bacterial genes coding for ribonucleotide reductases (RNRs), essential enzymes that provide DNA building blocks in all living cells. Despite its bacterial prevalence, the NrdR mechanism has been scarcely studied. We report the biochemical, biophysical, and bioinformatical characterization of NrdR and its binding sites from two major bacterial pathogens of the phylum Bacillota Listeria monocytogenes and Streptococcus pneumoniae. NrdR consists of a Zn-ribbon domain followed by an ATP-cone domain. We show that it forms tetramers that bind to DNA when loaded with ATP and dATP, but if loaded with only ATP, NrdR forms various oligomeric complexes unable to bind DNA. The DNA-binding site in L. monocytogenes is a pair of NrdR boxes separated by 15–16 bp, whereas in S. pneumoniae, the NrdR boxes are separated by unusually long spacers of 25–26 bp. This observation triggered a comprehensive binding study of four NrdRs from L. monocytogenes, S. pneumoniae, Escherichia coli, and Streptomyces coelicolor to a series of dsDNA fragments where the NrdR boxes were separated by 12–27 bp. The in vitro results were confirmed in vivo in E. coli and revealed that NrdR binds most efficiently when there is an integer number of DNA turns between the center of the two NrdR boxes. The study facilitates the prediction of NrdR binding sites in bacterial genomes and suggests that the NrdR mechanism is conserved throughout the bacterial domain. It sheds light on RNR regulation in Listeria and Streptococcus, and since NrdR does not occur in eukaryotes, opens a way to the development of novel antibiotics.

All nitrogen-fixing root nodule symbioses of angiosperms–legume and actinorhizal symbioses–possess a common ancestor. Molecular processes for the induction of root nodules are modulated by phytohormones, as is the case of the first nodulation-related transcription factor NODULE INCEPTION (NIN), whose expression can be induced by exogenous cytokinin in legumes. The process of actinorhizal nodule organogenesis is less well understood. To study the changes exerted by phytohormones on the expression of the orthologs of CYCLOPS, NIN, and NF-YA1 in the actinorhizal host Datisca glomerata, an axenic hydroponic system was established and used to examine the transcriptional responses (RTqPCR) in roots treated with the synthetic cytokinin 6-Benzylaminopurine (BAP), the natural auxin Phenylacetic acid (PAA), and the synthetic auxin 1-Naphthaleneacetic acid (NAA). The model legume Lotus japonicus was used as positive control. Molecular readouts for auxins and cytokinin were established: DgSAUR1 for PAA, DgGH3.1. for NAA, and DgARR9 for BAP. L. japonicus NIN was induced by BAP, PAA, and NAA in a dosage- and time-dependent manner. While expression of D. glomerata NIN2 could not be induced in roots, D. glomerata NIN1 was induced by PAA; this induction was abolished in the presence of exogenous BAP. Furthermore, the induction of DgNIN1 expression by PAA required ethylene and gibberellic acid. This study suggests that while cytokinin signaling is central for cortex-induced nodules of L. japonicus, it acts antagonistically to the induction of nodule primordia of D. glomerata by PAA in the root pericycle.

A small, nucleotide-binding domain, the ATP-cone, is found at the N-terminus of most ribonucleotide reductase (RNR) catalytic subunits. By binding adenosine triphosphate (ATP) or deoxyadenosine triphosphate (dATP) it regulates the enzyme activity of all classes of RNR. Functional and structural work on aerobic RNRs has revealed a plethora of ways in which dATP inhibits activity by inducing oligomerisation and preventing a productive radical transfer from one subunit to the active site in the other. Anaerobic RNRs, on the other hand, store a stable glycyl radical next to the active site and the basis for their dATP-dependent inhibition is completely unknown. We present biochemical, biophysical, and structural information on the effects of ATP and dATP binding to the anaerobic RNR from Prevotella copri. The enzyme exists in a dimer-tetramer equilibrium biased towards dimers when two ATP molecules are bound to the ATP-cone and tetramers when two dATP molecules are bound. In the presence of ATP, P. copri NrdD is active and has a fully ordered glycyl radical domain (GRD) in one monomer of the dimer. Binding of dATP to the ATP-cone results in loss of activity and increased dynamics of the GRD, such that it cannot be detected in the cryo-EM structures. The glycyl radical is formed even in the dATP-bound form, but the substrate does not bind. The structures implicate a complex network of interactions in activity regulation that involve the GRD more than 30 Å away from the dATP molecules, the allosteric substrate specificity site and a conserved but previously unseen flap over the active site. Taken together, the results suggest that dATP inhibition in anaerobic RNRs acts by increasing the flexibility of the flap and GRD, thereby preventing both substrate binding and radical mobilisation.

Aerobic ribonucleotide reductases (RNRs) initiate synthesis of DNA building blocks by generating a free radical within the R2 subunit; the radical is subsequently shuttled to the catalytic R1 subunit through proton-coupled electron transfer (PCET). We present a high-resolution room temperature structure of the class Ie R2 protein radical captured by x-ray free electron laser serial femtosecond crystallography. The structure reveals conformational reorganization to shield the radical and connect it to the translocation path, with structural changes propagating to the surface where the protein interacts with the catalytic R1 subunit. Restructuring of the hydrogen bond network, including a notably short O–O interaction of 2.41 angstroms, likely tunes and gates the radical during PCET. These structural results help explain radical handling and mobilization in RNR and have general implications for radical transfer in proteins. 

Ribonucleotide reductase (RNR) is an essential enzyme that catalyzes the synthesis of DNA building blocks in virtually all living cells. NrdR, an RNR-specific repressor, controls the transcription of RNR genes and, often, its own, in most bacteria and some archaea. NrdR senses the concentration of nucleotides through its ATP-cone, an evolutionarily mobile domain that also regulates the enzymatic activity of many RNRs, while a Zn-ribbon domain mediates binding to NrdR boxes upstream of and overlapping the transcription start site of RNR genes. Here, we combine biochemical and cryo-EM studies of NrdR from Streptomyces coelicolor to show, at atomic resolution, how NrdR binds to DNA. The suggested mechanism involves an initial dodecamer loaded with two ATP molecules that cannot bind to DNA. When dATP concentrations increase, an octamer forms that is loaded with one molecule each of dATP and ATP per monomer. A tetramer derived from this octamer then binds to DNA and represses transcription of RNR. In many bacteria — including well-known pathogens such as Mycobacterium tuberculosis — NrdR simultaneously controls multiple RNRs and hence DNA synthesis, making it an excellent target for novel antibiotics development.

Ribonucleotide reductases (RNRs) catalyze the reduction of ribonucleotides to the corresponding deoxyribonucleotides. The catalytic activity of most RNRs depends on the formation of a dimer of the catalytic subunits. The active site is located at the interface, and part of the substrate binding site and regulatory mechanisms work across the subunit in the dimer. In this study, we describe and characterize a novel domain responsible for forming the catalytic dimer in several class II RNRs. The 3D structure of the class II RNR from Rhodobacter sphaeroides reveals a so far undescribed α-helical domain in the dimer interface, which is embracing the other subunit. Genetic removal of this HUG domain leads to a severe reduction of activity paired with reduced dimerization capability. In comparison with other described RNRs, the enzyme with this domain is less dependent on the presence of nucleotides to act as allosteric effectors in the formation of dimers. The HUG domain appears to serve as an interlock to keep the dimer intact and functional even at low enzyme and/or effector concentrations.