The Small β-Barrel Domain: A Survey-Based Structural Analysis

Unveiling the structural and functional implications of uncharacterized NSPs and variations in the molecular toolkit across arteriviruses

Despite considerable scrutiny of mammalian arterivirus genomes, their genomic architecture remains incomplete, with several unannotated non-structural proteins (NSPs) and the enigmatic absence of methyltransferase (MTase) domains. Additionally, the host range of arteriviruses has expanded to include seven newly sequenced genomes from non-mammalian hosts, which remain largely unannotated and await detailed comparisons alongside mammalian isolates. Utilizing comparative genomics approaches and comprehensive sequence-structure analysis, we provide enhanced genomic architecture and annotations for arterivirus genomes. We identified the previously unannotated C-terminal domain of NSP3 as a winged helix-turn-helix domain and classified NSP7 as a new small β-barrel domain, both likely involved in interactions with viral RNA. NSP12 is identified as a derived variant of the N7-MTase-like Rossmann fold domain that retains core structural alignment with N7-MTases in Nidovirales but likely lacks enzymatic functionality due to the erosion of catalytic residues, indicating a unique role specific to mammalian arteriviruses. In contrast, non-mammalian arteriviruses sporadically retain a 2'-O-MTase and an exonuclease (ExoN) domain, which are typically absent in mammalian arteriviruses, highlighting contrasting evolutionary trends and variations in their molecular toolkit. Similar lineage-specific patterns are observed in the diversification of papain-like proteases and structural proteins. Overall, the study extends our knowledge of arterivirus genomic diversity and evolution.

A building blocks perspective on protein emergence and evolution

Recent findings increasingly suggest the emergence of proteins by mix and match of short peptides, or 'building blocks'. What are these building blocks, and how did they evolve into contemporary proteins? We review two complementary approaches to tackling these questions. First, a bottom-up approach that involves identifying putative components of primordial peptides, and the synthetic routes through which these peptides may have emerged. Second, searches in protein space to reveal building blocks that make up the contemporary protein repertoire; proteins that are not closely related to one another may nevertheless have certain parts in common, suggesting common ancestry. Identifying such shared building blocks, and characterizing their functions, can shed light on the ancient molecules from which proteins emerged, and hint at the mechanisms that govern their evolution. A key challenge lies in merging these two approaches to create a cohesive narrative of how proteins emerged and continue to evolve.

Unusual traits shape the architecture of the Ig ancestor molecule

Understanding the ancestral Ig domain's molecular structure and tracing the evolution of Ig-like proteins are fundamental components missing from our comprehension of their evolutionary trajectory and function. We have determined high-resolution structures of two Ig-like proteins from the evolutionary most ancestral phylum, Porifera. The structures reveal N-terminal Ig-like domains with an unconventional configuration of features that set them apart from canonical Ig domains. These findings prompted us to call this novel domain as Ig "Early Variable" (EV)-set. Remarkably, the EV-sets are linked to C1-set domains. To the best of our knowledge, the C1-set has not been previously reported in non-vertebrates. The IgV and IgC1 tandems and their combination into functional Ig-like receptors are part of the adaptive immune system in higher vertebrates, which allows for highly specific immune responses. By unveiling important clues into the molecular configuration of ancestral Ig domains, these findings challenge and expand our understanding of how immunity has evolved within its current landscape.

Massively parallel assessment of designed protein solution properties using mass spectrometry and peptide barcoding

Library screening and selection methods can determine the binding activities of individual members of large protein libraries given a physical link between protein and nucleotide sequence, which enables identification of functional molecules by DNA sequencing. However, the solution properties of individual protein molecules cannot be probed using such approaches because they are completely altered by DNA attachment. Mass spectrometry enables parallel evaluation of protein properties amenable to physical fractionation such as solubility and oligomeric state, but current approaches are limited to libraries of 1,000 or fewer proteins. Here, we improved mass spectrometry barcoding by co-synthesizing proteins with barcodes optimized to be highly multiplexable and minimally perturbative, scaling to libraries of >5,000 proteins. We use these barcodes together with mass spectrometry to assay the solution behavior of libraries of de novo-designed monomeric scaffolds, oligomers, binding proteins and nanocages, rapidly identifying design failure modes and successes.

Deep generative models of protein structure uncover distant relationships across a continuous fold space

Our views of fold space implicitly rest upon many assumptions that impact how we analyze, interpret and understand protein structure, function and evolution. For instance, is there an optimal granularity in viewing protein structural similarities (e.g., architecture, topology or some other level)? Similarly, the discrete/continuous dichotomy of fold space is central, but remains unresolved. Discrete views of fold space bin similar folds into distinct, non-overlapping groups; unfortunately, such binning can miss remote relationships. While hierarchical systems like CATH are indispensable resources, less heuristic and more conceptually flexible approaches could enable more nuanced explorations of fold space. Building upon an Urfold model of protein structure, here we present a deep generative modeling framework, termed DeepUrfold, for analyzing protein relationships at scale. DeepUrfold's learned embeddings occupy high-dimensional latent spaces that can be distilled for a given protein in terms of an amalgamated representation uniting sequence, structure and biophysical properties. This approach is structure-guided, versus being purely structure-based, and DeepUrfold learns representations that, in a sense, define superfamilies. Deploying DeepUrfold with CATH reveals evolutionarily-remote relationships that evade existing methodologies, and suggests a mostly-continuous view of fold space-a view that extends beyond simple geometric similarity, towards the realm of integrated sequence ↔ structure ↔ function properties.

Molecular insights into the initiation step of the Rcs signaling pathway

The Rcs pathway is repressed by the inner membrane protein IgaA under non-stressed conditions. This repression is hypothesized to be relieved by the binding of the outer membrane-anchored RcsF to IgaA. However, the precise mechanism by which RcsF binding triggers the signaling remains unclear. Here, we present the 1.8 Å resolution crystal structure capturing the interaction between IgaA and RcsF. Our comparative structural analysis, examining both the bound and unbound states of the periplasmic domain of IgaA (IgaAp), highlights rotational flexibility within IgaAp. Conversely, the conformation of RcsF remains unchanged upon binding. Our in vivo and in vitro studies do not support the model of a stable complex involving RcsF, IgaAp, and RcsDp. Instead, we demonstrate that the elements beyond IgaAp play a role in the interaction between IgaA and RcsD. These findings collectively allow us to propose a potential mechanism for the signaling across the inner membrane through IgaA.

Protein Fold Usages in Ribosomes: Another Glance to the Past

The analysis of protein fold usage, similar to codon usage, offers profound insights into the evolution of biological systems and the origins of modern proteomes. While previous studies have examined fold distribution in modern genomes, our study focuses on the comparative distribution and usage of protein folds in ribosomes across bacteria, archaea, and eukaryotes. We identify the prevalence of certain 'super-ribosome folds,' such as the OB fold in bacteria and the SH3 domain in archaea and eukaryotes. The observed protein fold distribution in the ribosomes announces the future power-law distribution where only a few folds are highly prevalent, and most are rare. Additionally, we highlight the presence of three copies of proto-Rossmann folds in ribosomes across all kingdoms, showing its ancient and fundamental role in ribosomal structure and function. Our study also explores early mechanisms of molecular convergence, where different protein folds bind equivalent ribosomal RNA structures in ribosomes across different kingdoms. This comparative analysis enhances our understanding of ribosomal evolution, particularly the distinct evolutionary paths of the large and small subunits, and underscores the complex interplay between RNA and protein components in the transition from the RNA world to modern cellular life. Transcending the concept of folds also makes it possible to group a large number of ribosomal proteins into five categories of urfolds or metafolds, which could attest to their ancestral character and common origins. This work also demonstrates that the gradual acquisition of extensions by simple but ordered folds constitutes an inexorable evolutionary mechanism. This observation supports the idea that simple but structured ribosomal proteins preceded the development of their disordered extensions.

HSV-1 ICP0 dimer domain adopts a novel β-barrel fold

Infected cell protein 0 (ICP0) is an immediate-early regulatory protein of herpes simplex virus 1 (HSV-1) that possesses E3 ubiquitin ligase activity. ICP0 transactivates viral genes, in part, through its C-terminal dimer domain (residues 555-767). Deletion of this dimer domain results in reduced viral gene expression, lytic infection, and reactivation from latency. Since ICP0's dimer domain is associated with its transactivation activity and efficient viral replication, we wanted to determine the structure of this specific domain. The C-terminus of ICP0 was purified from bacteria and analyzed by X-ray crystallography to solve its structure. Each subunit or monomer in the ICP0 dimer is composed of nine β-strands and two α-helices. Interestingly, two adjacent β-strands from one monomer "reach" into the adjacent subunit during dimer formation, generating two β-barrel-like structures. Additionally, crystallographic analyses indicate a tetramer structure is formed from two β-strands of each dimer, creating a "stacking" of the β-barrels. The structural protein database searches indicate the fold or structure adopted by the ICP0 dimer is novel. The dimer is held together by an extensive network of hydrogen bonds. Computational analyses reveal that ICP0 can either form a dimer or bind to SUMO1 via its C-terminal SUMO-interacting motifs but not both. Understanding the structure of the dimer domain will provide insights into the activities of ICP0 and, ultimately, the HSV-1 life cycle.

Fatty acid binding to the human transport proteins FABP3, FABP4, and FABP5 from a Ligand's perspective

Fatty acid binding proteins (FABPs) are a family of amphiphilic transport proteins with high diversity in terms of their amino acid sequences and binding preferences. Beyond their main biological role as cytosolic fatty acid transporters, many aspects regarding their binding mechanism and functional specializations in human cells remain unclear. In this work, the binding properties and thermodynamics of FABP3, FABP4, and FABP5 were analyzed under various physical conditions. For this purpose, the FABPs were loaded with fatty acids bearing fluorescence or spin probes as model ligands, comparing their binding affinities via microscale thermophoresis (MST) and continuous-wave electron paramagnetic resonance (CW EPR) spectroscopy. The CW EPR spectra of non-covalently bound 5- and 16-DOXYL stearic acid (5/16-DSA) deliver in-depth information about the dynamics and chemical environments of ligands inside the binding pockets of the FABPs. EPR spectral simulations allow the construction of binding curves, revealing two different binding states ('intermediately' and 'strongly' bound). The proportion of bound 5/16-DSA depends strongly on the FABP concentration and the temperature but with remarkable differences between the three isoforms. Additionally, the more dynamic state ('intermediately bound') seems to dominate at body temperature with thermodynamic preference. The ligand binding studies were supplemented by aggregation studies via dynamic light scattering and bioinformatic analyses. Beyond the remarkably fine-tuned binding properties exhibited by each FABP, which were discernible with our EPR-centered approach, the results of this work attest to the power of simple spectroscopic experiments to provide new insights into the ligand binding mechanisms of proteins in general on a molecular level.

The pathogen-encoded signalling receptor Tir exploits host-like intrinsic disorder for infection

The translocated intimin receptor (Tir) is an essential type III secretion system (T3SS) effector of attaching and effacing pathogens contributing to the global foodborne disease burden. Tir acts as a cell-surface receptor in host cells, rewiring intracellular processes by targeting multiple host proteins. We investigated the molecular basis for Tir's binding diversity in signalling, finding that Tir is a disordered protein with host-like binding motifs. Unexpectedly, also are several other T3SS effectors. By an integrative approach, we reveal that Tir dimerises via an antiparallel OB-fold within a highly disordered N-terminal cytosolic domain. Also, it has a long disordered C-terminal cytosolic domain partially structured at host-like motifs that bind lipids. Membrane affinity depends on lipid composition and phosphorylation, highlighting a previously unrecognised host interaction impacting Tir-induced actin polymerisation and cell death. Furthermore, multi-site tyrosine phosphorylation enables Tir to engage host SH2 domains in a multivalent fuzzy complex, consistent with Tir's scaffolding role and binding promiscuity. Our findings provide insights into the intracellular Tir domains, highlighting the ability of T3SS effectors to exploit host-like protein disorder as a strategy for host evasion.

HSV-1 ICP0 Dimer Domain Adopts a Novel β-barrel Fold

Infected cell protein 0 (ICP0) is an immediate-early regulatory protein of herpes simplex virus 1 (HSV-1) that possesses E3 ubiquitin ligase activity. ICP0 transactivates viral genes, in part, through its C-terminal dimer domain (residues 555-767). Deletion of this dimer domain results in reduced viral gene expression, lytic infection, and reactivation from latency. Since ICP0's dimer domain is associated with its transactivation activity and efficient viral replication, we wanted to determine the structure of this specific domain. The C-terminus of ICP0 was purified from bacteria and analyzed by X-ray crystallography to solve its structure. Each subunit or monomer in the ICP0 dimer is composed of nine β-strands and two α-helices. Interestingly, two adjacent β-strands from one monomer "reach" into the adjacent subunit during dimer formation, generating two β-barrel-like structures. Additionally, crystallographic analyses indicate a tetramer structure is formed from two β-strands of each dimer, creating a "stacking" of the β-barrels. The structural protein database searches indicate the fold or structure adopted by the ICP0 dimer is novel. The dimer is held together by an extensive network of hydrogen bonds. Computational analyses reveal that ICP0 can either form a dimer or bind to SUMO1 via its C-terminal SUMO-interacting motifs but not both. Understanding the structure of the dimer domain will provide insights into the activities of ICP0 and, ultimately, the HSV-1 life cycle.

Revisiting evolutionary trajectories and the organization of the Pleolipoviridae family

Archaeal pleomorphic viruses belonging to the Pleolipoviridae family represent an enigmatic group as they exhibit unique genomic features and are thought to have evolved through recombination with different archaeal plasmids. However, most of our understanding of the diversity and evolutionary trajectories of this clade comes from a handful of isolated representatives. Here we present 164 new genomes of pleolipoviruses obtained from metagenomic data of Australian hypersaline lakes and publicly available metagenomic data. We perform a comprehensive analysis on the diversity and evolutionary relationships of the newly discovered viruses and previously described pleolipoviruses. We propose to classify the viruses into five genera within the Pleolipoviridae family, with one new genus represented only by virus genomes retrieved in this study. Our data support the current hypothesis that pleolipoviruses reshaped their genomes through recombining with multiple different groups of plasmids, which is reflected in the diversity of their predicted replication strategies. We show that the proposed genus Epsilonpleolipovirus has evolutionary ties to pRN1-like plasmids from Sulfolobus, suggesting that this group could be infecting other archaeal phyla. Interestingly, we observed that the genome size of pleolipoviruses is correlated to the presence or absence of an integrase. Analyses of the host range revealed that all but one virus exhibit an extremely narrow range, and we show that the predicted tertiary structure of the spike protein is strongly associated with the host family, suggesting a specific adaptation to the host S-layer glycoprotein organization.

Cryo-EM structure of the chain-elongating E3 ubiquitin ligase UBR5

UBR5 is a nuclear E3 ligase that ubiquitinates a vast range of substrates for proteasomal degradation. This HECT domain-containing ubiquitin ligase has recently been identified as an important regulator of oncogenes, e.g., MYC, but little is known about its structure or mechanisms of substrate engagement and ubiquitination. Here, we present the cryo-EM structure of human UBR5, revealing an α-solenoid scaffold with numerous protein-protein interacting motifs, assembled into an antiparallel dimer that adopts further oligomeric states. Using cryo-EM processing tools, we observe the dynamic nature of the UBR5 catalytic domain, which we postulate is important for its enzymatic activity. We characterise the proteasomal nuclear import factor AKIRIN2 as an interacting protein and propose UBR5 as an efficient ubiquitin chain elongator. This preference for ubiquitinated substrates and several distinct domains for protein-protein interactions may explain how UBR5 is linked to several different signalling pathways and cancers. Together, our data expand on the limited knowledge of the structure and function of HECT E3 ligases.

Protective Effects of Caffeine on Chikungunya and Zika Virus Infections: An in Vitro and in Silico Study

Infection by viruses Chikungunya (CHIKV) and Zika (ZIKV) continue to be serious problems in tropical and subtropical areas of the world. Here, we evaluated the antiviral and virucidal activity of caffeine against CHIKV and ZIKV in Vero, A549, and Huh-7 cell lines. Results showed that caffeine displays antiviral properties against both viruses. By pre-and post-infection treatment, caffeine significantly inhibited CHIKV and ZIKV replication in a dose-dependent manner. Furthermore, caffeine showed a virucidal effect against ZIKV. Molecular docking suggests the possible binding of caffeine with envelope protein and RNA-dependent RNA polymerase of CHIKV and ZIKV. This is the first study that showed an antiviral effect of caffeine against CHIKV and ZIKV. Although further studies are needed to better understand the mechanism of caffeine-mediated repression of viral replication, caffeine appears to be a promising compound that could be used for in vivo studies, perhaps in synergy with other compounds present in daily beverages.

Evolutionary analysis and structure modelling of the Rcs-repressor IgaA unveil a functional role of two cytoplasmic small β-barrel (SBB) domains

The Rcs sensor system, comprising the RcsB/RcsC/RcsD and RcsF proteins, is used by bacteria of the order Enterobacterales to withstand envelope damage. In non-stress conditions, Rcs is repressed by IgaA, a membrane protein with three cytoplasmic regions (cyt-1, cyt-2 and cyt-3). How the Rcs-IgaA axis evolved within Enterobacterales has not been yet explored. Here, we report phylogenetic data supporting co-evolution of IgaA with RcsC/RcsD. Functional exchange assays showed that IgaA from Shigella and Dickeya, but not from Yersinia or the endosymbionts Photorhabdus and Sodalis, repress the Rcs system of Salmonella. IgaA from Dickeya, however, repress only partially the Rcs system despite being produced at high levels in the complementation assay. The modelled structures of these IgaA variants uncovered one periplasmic and two cytoplasmic conserved β-rich architectures forming partially closed small β-barrel (SBB) domains. Conserved residues map in a connector linking cytoplasmic SSB-1 and SBB-2 domains (E180-R265); a region of cyt-1 facing cyt-2 (R188-E194-D309 and T191-H326); and between cyt-2 and cyt-3 (H293-E328-R686). These structures validated early in vivo studies in Salmonella that assigned a role in function to R188, T191 and G262, and in addition revealed a previously unnoticed "hybrid" SBB-2 domain to which cyt-1 and cyt-2 contribute. IgaA variants not functional or partially functional in Salmonella lack H192-P249 and R255-D313 interactions. Among these variants, only IgaA from Dickeya conserves the helix α6 in SSB-1 that is present in IgaA from Salmonella and Shigella. RcsF and RcsD, which interact directly with IgaA, failed to show structural features linked to specific IgaA variants. Altogether, our data provide new insights into IgaA by mapping residues selected differently during evolution and involved in function. Our data also infer contrasting lifestyles of Enterobacterales bacteria as source of variability in the IgaA-RcsD/IgaA-RcsF interactions.

c-di-GMP regulates activity of the PlzA RNA chaperone from the Lyme disease spirochete

PlzA is a c-di-GMP-binding protein crucial for adaptation of the Lyme disease spirochete Borrelia (Borreliella) burgdorferi during its enzootic life cycle. Unliganded apo-PlzA is important for vertebrate infection, while liganded holo-PlzA is important for survival in the tick; however, the biological function of PlzA has remained enigmatic. Here, we report that PlzA has RNA chaperone activity that is inhibited by c-di-GMP binding. Holo- and apo-PlzA bind RNA and accelerate RNA annealing, while only apo-PlzA can strand displace and unwind double-stranded RNA. Guided by the crystal structure of PlzA, we identified several key aromatic amino acids protruding from the N- and C-terminal domains that are required for RNA-binding and unwinding activity. Our findings illuminate c-di-GMP as a switch controlling the RNA chaperone activity of PlzA, and we propose that complex RNA-mediated modulatory mechanisms allow PlzA to regulate gene expression during both the vector and host phases of the B. burgdorferi life cycle.

Characterization of a fold in TANGO1 evolved from SH3 domains for the export of bulky cargos

Bulky cargos like procollagens, apolipoproteins, and mucins exceed the size of conventional COPII vesicles. During evolution a process emerged in metazoans, predominantly governed by the TANGO1 protein family, that organizes cargo at the exit sites of the endoplasmic reticulum and facilitates export by the formation of tunnel-like connections between the ER and Golgi. Hitherto, cargo-recognition appeared to be mediated by an SH3-like domain. Based on structural and dynamic data as well as interaction studies from NMR spectroscopy and microscale thermophoresis presented here, we show that the luminal cargo-recognition domain of TANGO1 adopts a new functional fold for which we suggest the term MOTH (MIA, Otoraplin, TALI/TANGO1 homology) domain. These MOTH domains, as well as an evolutionary intermediate found in invertebrates, constitute a distinct domain family that emerged from SH3 domains and acquired the ability to bind collagen.

Structure of the human UBR5 E3 ubiquitin ligase

The human UBR5 is a single polypeptide chain homology to E6AP C terminus (HECT)-type E3 ubiquitin ligase essential for embryonic development in mammals. Dysregulated UBR5 functions like an oncoprotein to promote cancer growth and metastasis. Here, we report that UBR5 assembles into a dimer and a tetramer. Our cryoelectron microscopy (cryo-EM) structures reveal that two crescent-shaped UBR5 monomers assemble head to tail to form the dimer, and two dimers bind face to face to form the cage-like tetramer with all four catalytic HECT domains facing the central cavity. Importantly, the N-terminal region of one subunit and the HECT of the other form an "intermolecular jaw" in the dimer. We show the jaw-lining residues are important for function, suggesting that the intermolecular jaw functions to recruit ubiquitin-loaded E2 to UBR5. Further work is needed to understand how oligomerization regulates UBR5 ligase activity. This work provides a framework for structure-based anticancer drug development and contributes to a growing appreciation of E3 ligase diversity.

Cryo-EM structure of ssDNA bacteriophage ΦCjT23 provides insight into early virus evolution

The origin of viruses remains an open question. While lack of detectable sequence similarity hampers the analysis of distantly related viruses, structural biology investigations of conserved capsid protein structures facilitate the study of distant evolutionary relationships. Here we characterize the lipid-containing ssDNA temperate bacteriophage ΦCjT23, which infects Flavobacterium sp. (Bacteroidetes). We report ΦCjT23-like sequences in the genome of strains belonging to several Flavobacterium species. The virion structure determined by cryogenic electron microscopy reveals similarities to members of the viral kingdom Bamfordvirae that currently consists solely of dsDNA viruses with a major capsid protein composed of two upright β-sandwiches. The minimalistic structure of ΦCjT23 suggests that this phage serves as a model for the last common ancestor between ssDNA and dsDNA viruses in the Bamfordvirae. Both ΦCjT23 and the related phage FLiP infect Flavobacterium species found in several environments, suggesting that these types of viruses have a global distribution and a shared evolutionary origin. Detailed comparisons to related, more complex viruses not only expand our knowledge about this group of viruses but also provide a rare glimpse into early virus evolution.

3β-Corner Stability by Comparative Molecular Dynamics Simulations

This study explored the mechanisms by which the stability of super-secondary structures of the 3β-corner type autonomously outside the protein globule are maintained in an aqueous environment. A molecular dynamic (MD) study determined the behavioral diversity of a large set of non-homologous 3β-corner structures of various origins. We focused on geometric parameters such as change in gyration radius, solvent-accessible area, major conformer lifetime and torsion angles, and the number of hydrogen bonds. Ultimately, a set of 3β-corners from 330 structures was characterized by a root mean square deviation (RMSD) of less than 5 Å, a change in the gyration radius of no more than 5%, and the preservation of amino acid residues positioned within the allowed regions on the Ramachandran map. The studied structures retained their topologies throughout the MD experiments. Thus, the 3β-corner structure was found to be rather stable per se in a water environment, i.e., without the rest of a protein molecule, and can act as the nucleus or “ready-made” building block in protein folding. The 3β-corner can also be considered as an independent object for study in field of structural biology.

A widespread family of WYL-domain transcriptional regulators co-localizes with diverse phage defence systems and islands

Bacteria are under constant assault by bacteriophages and other mobile genetic elements. As a result, bacteria have evolved a multitude of systems that protect from attack. Genes encoding bacterial defence mechanisms can be clustered into 'defence islands', providing a potentially synergistic level of protection against a wider range of assailants. However, there is a comparative paucity of information on how expression of these defence systems is controlled. Here, we functionally characterize a transcriptional regulator, BrxR, encoded within a recently described phage defence island from a multidrug resistant plasmid of the emerging pathogen Escherichia fergusonii. Using a combination of reporters and electrophoretic mobility shift assays, we discovered that BrxR acts as a repressor. We present the structure of BrxR to 2.15 Å, the first structure of this family of transcription factors, and pinpoint a likely binding site for ligands within the WYL-domain. Bioinformatic analyses demonstrated that BrxR-family homologues are widespread amongst bacteria. About half (48%) of identified BrxR homologues were co-localized with a diverse array of known phage defence systems, either alone or clustered into defence islands. BrxR is a novel regulator that reveals a common mechanism for controlling the expression of the bacterial phage defence arsenal.

BioS2Net: Holistic Structural and Sequential Analysis of Biomolecules Using a Deep Neural Network

BACKGROUND

For decades, the rate of solving new biomolecular structures has been exceeding that at which their manual classification and feature characterisation can be carried out efficiently. Therefore, a new comprehensive and holistic tool for their examination is needed.

METHODS

Here we propose the Biological Sequence and Structure Network (BioS2Net), which is a novel deep neural network architecture that extracts both sequential and structural information of biomolecules. Our architecture consists of four main parts: (i) a sequence convolutional extractor, (ii) a 3D structure extractor, (iii) a 3D structure-aware sequence temporal network, as well as (iv) a fusion and classification network.

RESULTS

We have evaluated our approach using two protein fold classification datasets. BioS2Net achieved a 95.4% mean class accuracy on the eDD dataset and a 76% mean class accuracy on the F184 dataset. The accuracy of BioS2Net obtained on the eDD dataset was comparable to results achieved by previously published methods, confirming that the algorithm described in this article is a top-class solution for protein fold recognition.

CONCLUSIONS

BioS2Net is a novel tool for the holistic examination of biomolecules of known structure and sequence. It is a reliable tool for protein analysis and their unified representation as feature vectors.

OB-fold Families of Genome Guardians: A Universal Theme Constructed From the Small β-barrel Building Block

The maintenance of genome stability requires the coordinated actions of multiple proteins and protein complexes, that are collectively known as genome guardians. Within this broadly defined family is a subset of proteins that contain oligonucleotide/oligosaccharide-binding folds (OB-fold). While OB-folds are widely associated with binding to single-stranded DNA this view is no longer an accurate depiction of how these domains are utilized. Instead, the core of the OB-fold is modified and adapted to facilitate binding to a variety of DNA substrates (both single- and double-stranded), phospholipids, and proteins, as well as enabling catalytic function to a multi-subunit complex. The flexibility accompanied by distinctive oligomerization states and quaternary structures enables OB-fold genome guardians to maintain the integrity of the genome via a myriad of complex and dynamic, protein-protein; protein-DNA, and protein-lipid interactions in both prokaryotes and eukaryotes.

Structural basis of Ty3 retrotransposon integration at RNA Polymerase III-transcribed genes

Retrotransposons are endogenous elements that have the ability to mobilise their DNA between different locations in the host genome. The Ty3 retrotransposon integrates with an exquisite specificity in a narrow window upstream of RNA Polymerase (Pol) III-transcribed genes, representing a paradigm for harmless targeted integration. Here we present the cryo-EM reconstruction at 4.0 Å of an active Ty3 strand transfer complex bound to TFIIIB transcription factor and a tRNA gene. The structure unravels the molecular mechanisms underlying Ty3 targeting specificity at Pol III-transcribed genes and sheds light into the architecture of retrotransposon machinery during integration. Ty3 intasome contacts a region of TBP, a subunit of TFIIIB, which is blocked by NC2 transcription regulator in RNA Pol II-transcribed genes. A newly-identified chromodomain on Ty3 integrase interacts with TFIIIB and the tRNA gene, defining with extreme precision the integration site position.

Fold Evolution before LUCA: Common Ancestry of SH3 Domains and OB Domains

SH3 and OB are the simplest, oldest, and most common protein domains within the translation system. SH3 and OB domains are β-barrels that are structurally similar but are topologically distinct. To transform an OB domain to a SH3 domain, β-strands must be permuted in a multistep and evolutionarily implausible mechanism. Here, we explored relationships between SH3 and OB domains of ribosomal proteins, initiation, and elongation factors using a combined sequence- and structure-based approach. We detect a common core of SH3 and OB domains, as a region of significant structure and sequence similarity. The common core contains four β-strands and a loop, but omits the fifth β-strand, which is variable and is absent from some OB and SH3 domain proteins. The structure of the common core immediately suggests a simple permutation mechanism for interconversion between SH3 and OB domains, which appear to share an ancestor. The OB domain was formed by duplication and adaptation of the SH3 domain core, or vice versa, in a simple and probable transformation. By employing the folding algorithm AlphaFold2, we demonstrated that an ancestral reconstruction of a permuted SH3 sequence folds into an OB structure, and an ancestral reconstruction of a permuted OB sequence folds into a SH3 structure. The tandem SH3 and OB domains in the universal ribosomal protein uL2 share a common ancestor, suggesting that the divergence of these two domains occurred before the last universal common ancestor.

Topological and Structural Plasticity of the Single Ig Fold and the Double Ig Fold Present in CD19

The Ig fold has had a remarkable success in vertebrate evolution, with a presence in over 2% of human genes. The Ig fold is not just the elementary structural domain of antibodies and TCRs, it is also at the heart of a staggering 30% of immunologic cell surface receptors, making it a major orchestrator of cell–cell interactions. While BCRs, TCRs, and numerous Ig-based cell surface receptors form homo- or heterodimers on the same cell surface (in cis), many of them interface as ligand-receptors (checkpoints) on interacting cells (in trans) through their Ig domains. New Ig-Ig interfaces are still being discovered between Ig-based cell surface receptors, even in well-known families such as B7. What is largely ignored, however, is that the Ig fold itself is pseudosymmetric, a property that makes the Ig domain a versatile self-associative 3D structure and may, in part, explain its success in evolution, especially through its ability to bind in cis or in trans in the context of cell surface receptor–ligand interactions. In this paper, we review the Ig domains’ tertiary and quaternary pseudosymmetries, with particular attention to the newly identified double Ig fold in the solved CD19 molecular structure to highlight the underlying fundamental folding elements of Ig domains, i.e., Ig protodomains. This pseudosymmetric property of Ig domains gives us a decoding frame of reference to understand the fold, relate all Ig domain forms, single or double, and suggest new protein engineering avenues.

Comparative genomics provides structural and functional insights into Bacteroides RNA biology

Bacteria employ noncoding RNA molecules for a wide range of biological processes, including scaffolding large molecular complexes, catalyzing chemical reactions, defending against phages, and controlling gene expression. Secondary structures, binding partners, and molecular mechanisms have been determined for numerous small noncoding RNAs (sRNAs) in model aerobic bacteria. However, technical hurdles have largely prevented analogous analyses in the anaerobic gut microbiota. While experimental techniques are being developed to investigate the sRNAs of gut commensals, computational tools and comparative genomics can provide immediate functional insight. Here, using Bacteroides thetaiotaomicron as a representative microbiota member, we illustrate how comparative genomics improves our understanding of the RNA biology in an understudied gut bacterium. We investigate putative RNA-binding proteins and predict a Bacteroides cold-shock protein homologue to have an RNA-related function. We apply an in-silico protocol incorporating both sequence and structural analysis to determine the consensus structures and conservation of nine Bacteroides noncoding RNA families. Using structure probing, we validate and refine these predictions, and deposit them in the Rfam database. Through synteny analyses, we illustrate how genomic co-conservation can serve as a predictor of sRNA function. Altogether, this work showcases the power of RNA informatics for investigating the RNA biology of anaerobic microbiota members.

The β-link motif in protein architecture

The β-link is a composite protein motif consisting of a G1β β-bulge and a type II β-turn, and is generally found at the end of two adjacent strands of antiparallel β-sheet. The 1,2-positions of the β-bulge are also the 3,4-positions of the β-turn, with the result that the N-terminal portion of the polypeptide chain is orientated at right angles to the β-sheet. Here, it is reported that the β-link is frequently found in certain protein folds of the SCOPe structural classification at specific locations where it connects a β-sheet to another area of a protein. It is found at locations where it connects one β-sheet to another in the β-sandwich and related structures, and in small (four-, five- or six-stranded) β-barrels, where it connects two β-strands through the polypeptide chain that crosses an open end of the barrel. It is not found in larger (eight-stranded or more) β-barrels that are straightforward β-meanders. In some cases it initiates a connection between a single β-sheet and an α-helix. The β-link also provides a framework for catalysis in serine proteases, where the catalytic serine is part of a conserved β-link, and in cysteine proteases, including M^pro of human SARS-CoV-2, in which two residues of the active site are located in a conserved β-link.

The World of Stable Ribonucleoproteins and Its Mapping With Grad-Seq and Related Approaches

Macromolecular complexes of proteins and RNAs are essential building blocks of cells. These stable supramolecular particles can be viewed as minimal biochemical units whose structural organization, i.e., the way the RNA and the protein interact with each other, is directly linked to their biological function. Whether those are dynamic regulatory ribonucleoproteins (RNPs) or integrated molecular machines involved in gene expression, the comprehensive knowledge of these units is critical to our understanding of key molecular mechanisms and cell physiology phenomena. Such is the goal of diverse complexomic approaches and in particular of the recently developed gradient profiling by sequencing (Grad-seq). By separating cellular protein and RNA complexes on a density gradient and quantifying their distributions genome-wide by mass spectrometry and deep sequencing, Grad-seq charts global landscapes of native macromolecular assemblies. In this review, we propose a function-based ontology of stable RNPs and discuss how Grad-seq and related approaches transformed our perspective of bacterial and eukaryotic ribonucleoproteins by guiding the discovery of new RNA-binding proteins and unusual classes of noncoding RNAs. We highlight some methodological aspects and developments that permit to further boost the power of this technique and to look for exciting new biology in understudied and challenging biological models.

Structure of a bacterial OapB protein with its OLE RNA target gives insights into the architecture of the OLE ribonucleoprotein complex

Bacterial noncoding RNAs (ncRNAs) play key roles in many biological processes including gene regulation, RNA processing and modification, and protein synthesis and translocation. OLE RNAs, found in many Gram-positive species, are one of the largest highly structured ncRNA classes whose biochemical functions remain unknown. In Bacillus halodurans, OLE RNAs interact with at least two proteins, OapA and OapB, which are required to assemble a functional OLE ribonucleoprotein (RNP) complex contributing to cellular responses to certain environmental stresses. We established X-ray structural models that reveal the sequence elements and tertiary structural features of OLE RNA that are critical for its specific recognition by OapB, which will aid future exploration of the biological and biochemical functions of the unusual OLE RNP complex.

The OLE (ornate, large, and extremophilic) RNA class is one of the most complex and well-conserved bacterial noncoding RNAs known to exist. This RNA is known to be important for bacterial responses to stress caused by short-chain alcohols, cold, and elevated Mg²⁺ concentrations. These biological functions have been shown to require the formation of a ribonucleoprotein (RNP) complex including at least two protein partners: OLE-associated protein A (OapA) and OLE-associated protein B (OapB). OapB directly binds OLE RNA with high-affinity and specificity and is believed to assist in assembling the functional OLE RNP complex. To provide the atomic details of OapB–OLE RNA interaction and to potentially reveal previously uncharacterized protein–RNA interfaces, we determined the structure of OapB from Bacillus halodurans alone and in complex with an OLE RNA fragment at resolutions of 1.0 Å and 2.0 Å, respectively. The structure of OapB exhibits a K-shaped overall architecture wherein its conserved KOW motif and additional unique structural elements of OapB form a bipartite RNA-binding surface that docks to the P13 hairpin and P12.2 helix of OLE RNA. These high-resolution structures elucidate the molecular contacts used by OapB to form a stable RNP complex and explain the high conservation of sequences and structural features at the OapB–OLE RNA-binding interface. These findings provide insight into the role of OapB in the assembly and biological function of OLE RNP complex and can guide the exploration of additional possible OLE RNA-binding interactions present in OapB.

Outer membrane protein evolution

Outer membrane proteins have remarkably homogeneous structure. They are all up down β-barrels. Up down barrels themselves are composed of repeated sets of β-hairpins. The consistency of the usage of the β-hairpin throughout the outer membrane milieu allows for interrogation of the evolution of these repetitive structures. Here we describe recent investigations of outer membrane protein evolution and how evolutionary precepts have been used for novel outer membrane protein design.

Functional analysis of Rossmann-like domains reveals convergent evolution of topology and reaction pathways

Rossmann folds are ancient, frequently diverged domains found in many biological reaction pathways where they have adapted for different functions. Consequently, discernment and classification of their homologous relations and function can be complicated. We define a minimal Rossmann-like structure motif (RLM) that corresponds for the common core of known Rossmann domains and use this motif to identify all RLM domains in the Protein Data Bank (PDB), thus finding they constitute about 20% of all known 3D structures. The Evolutionary Classification of protein structure Domains (ECOD) classifies RLM domains in a number of groups that lack evidence for homology (X-groups), which suggests that they could have evolved independently multiple times. Closely related, homologous RLM enzyme families can diverge to bind different ligands using similar binding sites and to catalyze different reactions. Conversely, non-homologous RLM domains can converge to catalyze the same reactions or to bind the same ligand with alternate binding modes. We discuss a special case of such convergent evolution that is relevant to the polypharmacology paradigm, wherein the same drug (methotrexate) binds to multiple non-homologous RLM drug targets with different topologies. Finally, assigning proteins with RLM domain to the Enzyme Commission classification suggest that RLM enzymes function mainly in metabolism (and comprise 38% of reference metabolic pathways) and are overrepresented in extant pathways that represent ancient biosynthetic routes such as nucleotide metabolism, energy metabolism, and metabolism of amino acids. In fact, RLM enzymes take part in five out of eight enzymatic reactions of the Wood-Ljungdahl metabolic pathway thought to be used by the last universal common ancestor (LUCA). The prevalence of RLM domains in this ancient metabolism might explain their wide distribution among enzymes.

Towards the application of Tc toxins as a universal protein translocation system

Tc toxins are bacterial protein complexes that inject cytotoxic enzymes into target cells using a syringe-like mechanism. Tc toxins are composed of a membrane translocator and a cocoon that encapsulates a toxic enzyme. The toxic enzyme varies between Tc toxins from different species and is not conserved. Here, we investigate whether the toxic enzyme can be replaced by other small proteins of different origin and properties, namely Cdc42, herpes simplex virus ICP47, Arabidopsis thaliana iLOV, Escherichia coli DHFR, Ras-binding domain of CRAF kinase, and TEV protease. Using a combination of electron microscopy, X-ray crystallography and in vitro translocation assays, we demonstrate that it is possible to turn Tc toxins into customizable molecular syringes for delivering proteins of interest across membranes. We also infer the guidelines that protein cargos must obey in terms of size, charge, and fold in order to apply Tc toxins as a universal protein translocation system.

Tc toxins are a major class of bacterial toxin translocation systems that inject toxic enzymes into target cells. Here the authors present functional and structural data showing that the toxic enzyme can be replaced by other small proteins and identify prerequisites required for successful translocation, which could facilitate the development of functional Tc-based protein injection devices.

The Urfold: Structural similarity just above the superfold level?

We suspect that there is a level of granularity of protein structure intermediate between the classical levels of "architecture" and "topology," as reflected in such phenomena as extensive three-dimensional structural similarity above the level of (super)folds. Here, we examine this notion of architectural identity despite topological variability, starting with a concept that we call the "Urfold." We believe that this model could offer a new conceptual approach for protein structural analysis and classification: indeed, the Urfold concept may help reconcile various phenomena that have been frequently recognized or debated for years, such as the precise meaning of "significant" structural overlap and the degree of continuity of fold space. More broadly, the role of structural similarity in sequence↔structure↔function evolution has been studied via many models over the years; by addressing a conceptual gap that we believe exists between the architecture and topology levels of structural classification schemes, the Urfold eventually may help synthesize these models into a generalized, consistent framework. Here, we begin by qualitatively introducing the concept.

Composition of the Survival Motor Neuron (SMN) Complex in Drosophila melanogaster

Spinal Muscular Atrophy (SMA) is caused by homozygous mutations in the human survival motor neuron 1 (SMN1) gene. SMN protein has a well-characterized role in the biogenesis of small nuclear ribonucleoproteins (snRNPs), core components of the spliceosome. SMN is part of an oligomeric complex with core binding partners, collectively called Gemins. Biochemical and cell biological studies demonstrate that certain Gemins are required for proper snRNP assembly and transport. However, the precise functions of most Gemins are unknown. To gain a deeper understanding of the SMN complex in the context of metazoan evolution, we investigated its composition in Drosophila melanogaster Using transgenic flies that exclusively express Flag-tagged SMN from its native promoter, we previously found that Gemin2, Gemin3, Gemin5, and all nine classical Sm proteins, including Lsm10 and Lsm11, co-purify with SMN. Here, we show that CG2941 is also highly enriched in the pulldown. Reciprocal co-immunoprecipitation reveals that epitope-tagged CG2941 interacts with endogenous SMN in Schneider2 cells. Bioinformatic comparisons show that CG2941 shares sequence and structural similarity with metazoan Gemin4. Additional analysis shows that three other genes (CG14164, CG31950 and CG2371) are not orthologous to Gemins 6-7-8, respectively, as previously suggested. In D.melanogaster, CG2941 is located within an evolutionarily recent genomic triplication with two other nearly identical paralogous genes (CG32783 and CG32786). RNAi-mediated knockdown of CG2941 and its two close paralogs reveals that Gemin4 is essential for organismal viability.

Protodomains: Symmetry-Related Supersecondary Structures in Proteins and Self-Complementarity

We will consider in this chapter supersecondary structures (SSS) as a set of secondary structure elements (SSEs) found in protein domains. Some SSS arrangements/topologies have been consistently observed within known tertiary structural domains. We use them in the context of repeating supersecondary structures that self-assemble in a symmetric arrangement to form a domain. We call them protodomains (or protofolds). Protodomains are some of the most interesting and insightful SSSs. Within a given 3D protein domain/fold, recognizing such sets may give insights into a possible evolutionary process of duplication, fusion, and coevolution of these protodomains, pointing to possible original protogenes. On protein folding itself, pseudosymmetric domains may point to a "directed" assembly of pseudosymmetric protodomains, directed by the only fact that they are tethered together in a protein chain. On function, tertiary functional sites often occur at protodomain interfaces, as they often occur at domain-domain interfaces in quaternary arrangements.First, we will briefly review some lessons learned from a previously published census of pseudosymmetry in protein domains (Myers-Turnbull, D. et al., J Mol Biol. 426:2255-2268, 2014) to introduce protodomains/protofolds. We will observe that the most abundant and diversified folds, or superfolds, in the currently known protein structure universe are indeed pseudosymmetric. Then, we will learn by example and select a few domain representatives of important pseudosymmetric folds and chief among them the immunoglobulin (Ig) fold and go over a pseudosymmetry supersecondary structure (protodomain) analysis in tertiary and quaternary structures. We will point to currently available software tools to help in identifying pseudosymmetry, delineating protodomains, and see how the study of pseudosymmetry and the underlying supersecondary structures can enrich a structural analysis. This should potentially help in protein engineering, especially in the development of biologics and immunoengineering.