Functional and evolutionary analysis of viral proteins containing a Rossmann-like fold

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.

The cluster model of energy transduction in biological systems

The central problem in transduction is to explain how the energy caught from sunlight by chloroplasts becomes biological work. Or to express it in different terms: how does the energy remain trapped in the biological network and not get lost through thermalization into the environment? The pathway consists of an immensely large number of steps crossing hierarchical levels - some upwards, to larger assemblies, others downwards into energy rich molecules - before fuelling an action potential or a contracting cell. Accepting the assumption that steps are executed by protein domains, we expect that transduction mechanisms are the result of conformational changes, which in turn involve rearrangements of the bonds responsible for the protein fold. But why are these essential changes so difficult to detect? In this presentation, the metabolic pathway is viewed as equivalent to an energy conduit composed of equally sized units - the protein domains - rather than a row of catalysts. The flow of energy through them occurs by the same mechanism as through the cytoplasmic medium (water). This mechanism is based on the cluster-wave model of water structure, which successfully explains the transfer of energy through the liquid medium responsible for the build up of osmotic pressure. The analogy to the line of balls called "Newton's cradle" provides a useful comparison, since there the transfer is also invisible to us because the intermediate balls are motionless. It is further proposed that the spatial arrangements of the H-bonds of the α and β secondary structures support wave motion, with the linear and lateral forms of the groups of bonds belonging to the helices and sheets executing the longitudinal and transverse modes, respectively.

Ongoing shuffling of protein fragments diversifies core viral functions linked to interactions with bacterial hosts

Biological modularity enhances evolutionary adaptability. This principle is vividly exemplified by bacterial viruses (phages), which display extensive genomic modularity. Phage genomes are composed of independent functional modules that evolve separately and recombine in various configurations. While genomic modularity in phages has been extensively studied, less attention has been paid to protein modularity-proteins consisting of distinct building blocks that can evolve and recombine, enhancing functional and genetic diversity. Here, we use a set of 133,574 representative phage proteins and highly sensitive homology detection to capture instances of domain mosaicism, defined as fragment sharing between two otherwise unrelated proteins, and to understand its relationship with functional diversity in phage genomes. We discover that unrelated proteins from diverse functional classes frequently share homologous domains. This phenomenon is particularly pronounced within receptor-binding proteins, endolysins, and DNA polymerases. We also identify multiple instances of recent diversification via domain shuffling in receptor-binding proteins, neck passage structures, endolysins and some members of the core replication machinery, often transcending distant taxonomic and ecological boundaries. Our findings suggest that ongoing diversification via domain shuffling is reflective of a co-evolutionary arms race, driven by the need to overcome various bacterial resistance mechanisms against phages.

Insights into peculiar fungal LPMO family members holding a short C-terminal sequence reminiscent of phosphate binding motifs

Lytic polysaccharide monooxygenases (LPMOs) are taxonomically widespread copper-enzymes boosting biopolymers conversion (e.g. cellulose, chitin) in Nature. White-rot Polyporales, which are major fungal wood decayers, may possess up to 60 LPMO-encoding genes belonging to the auxiliary activities family 9 (AA9). Yet, the functional relevance of such multiplicity remains to be uncovered. Previous comparative transcriptomic studies of six Polyporales fungi grown on cellulosic substrates had shown the overexpression of numerous AA9-encoding genes, including some holding a C-terminal domain of unknown function ("X282"). Here, after carrying out structural predictions and phylogenetic analyses, we selected and characterized six AA9-X282s with different C-term modularities and atypical features hitherto unreported. Unexpectedly, after screening a large array of conditions, these AA9-X282s showed only weak binding properties to cellulose, and low to no cellulolytic oxidative activity. Strikingly, proteomic analysis revealed the presence of multiple phosphorylated residues at the surface of these AA9-X282s, including a conserved residue next to the copper site. Further analyses focusing on a 9 residues glycine-rich C-term extension suggested that it could hold phosphate-binding properties. Our results question the involvement of these AA9 proteins in the degradation of plant cell wall and open new avenues as to the divergence of function of some AA9 members.

Insertions and deletions mediated functional divergence of Rossmann fold enzymes

Nucleobase-containing coenzymes are hypothesized to be relics of an early RNA-based world that preceded the emergence of proteins. Despite the importance of coenzyme-protein synergisms, their emergence and evolution remain understudied. An excellent target to address this issue is the Rossmann fold, the most catalytically diverse and abundant protein architecture in nature. We investigated two main Rossmann lineages: the nicotinamide adenine dinucleotide phosphate (NAD(P)) and the S-adenosyl methionine (SAM)- binding superfamilies. To identify the evolutionary changes that lead to a coenzyme specificity switch on these superfamilies, we performed structural and sequence-based Hidden Markov model analysis to systematically search for key motifs in their coenzyme-binding pockets. Our analyses revealed that through insertions and deletions (InDels) and a residue substitution, the ancient β1-loop-α1 coenzyme-binding structure of NAD(P) could be reshaped into the SAM-binding β1-loop-α1 structure. To experimentally prove this obsevation, we removed three amino acids from the NAD(P)-binding pocket and solved the structure of the resulting mutant, revealing the characteristic loop features of the SAM-binding pocket. To confirm the binding to SAM, we performed isothermal titration calorimetry measurements. Molecular dynamics simulations also corroborated the role of InDels in abolishing NAD binding and acquiring SAM binding. Our results uncovered how nature may have utilized insertions and deletions to optimize the different coenzyme-binding pockets and the distinct functionalities observed for Rossmann superfamilies. This work also proposes a general mechanism by which protein templates could have been recycled through the course of evolution to adopt different coenzymes and confer distinct chemistries.

Methyltransferases of Riboviria

Many viruses from the realm Riboviria infecting eukaryotic hosts encode protein domains with sequence similarity to S-adenosylmethionine-dependent methyltransferases. These protein domains are thought to be involved in methylation of the 5'-terminal cap structures in virus mRNAs. Some methyltransferase-like domains of Riboviria are homologous to the widespread cellular FtsJ/RrmJ-like methyltransferases involved in modification of cellular RNAs; other methyltransferases, found in a subset of positive-strand RNA viruses, have been assigned to a separate "Sindbis-like" family; and coronavirus-specific Nsp13/14-like methyltransferases appeared to be different from both those classes. The representative structures of proteins from all three groups belong to a specific variety of the Rossmann fold with a seven-stranded β-sheet, but it was unclear whether this structural similarity extends to the level of conserved sequence signatures. Here I survey methyltransferases in Riboviria and derive a joint sequence alignment model that covers all groups of virus methyltransferases and subsumes the previously defined conserved sequence motifs. Analysis of the spatial structures indicates that two highly conserved residues, a lysine and an aspartate, frequently contact a water molecule, which is located in the enzyme active center next to the methyl group of S-adenosylmethionine cofactor and could play a key role in the catalytic mechanism of the enzyme. Phylogenetic evidence indicates a likely origin of all methyltransferases of Riboviria from cellular RrmJ-like enzymes and their rapid divergence with infrequent horizontal transfer between distantly related viruses.

Structure of the Core Postfusion Porcine Endogenous Retrovirus Fusion Protein

Retroviral elements from endogenous retroviruses have functions in mammalian physiology. The best-known examples are the envelope proteins that function in placenta development and immune suppression. Porcine endogenous retroviruses (PERVs) are an understudied class of endogenous retroviruses that infect cultured human cells, raising concern regarding porcine xenografts. The PERV envelope glycoprotein has also been proposed as a possible swine syncytin with a role in placental development. Despite the growing interest in PERVs, their envelope glycoproteins remain poorly characterized. Here, we successfully determined the postfusion crystal structure of the PERV core fusion ectodomain. The PERV fusion protein structure reveals a conserved class I viral fusion protein six-helix bundle. Biophysical experiments demonstrated that the thermodynamic stability of the PERV fusion protein secondary structure was the same at physiological and acidic pHs. A conserved surface analysis highlights the high degree of sequence conservation among retroviral fusogens in the chain reversal region that facilitates the large-scale conformational change required for membrane fusion. Further structural alignment of class I viral fusogens revealed a phylogenetic clustering that shows evolution into various lineages that correlate with virus mechanisms of cell entry. Our work indicates that structural dendrograms can be used to qualitatively infer insights into the fusion mechanisms of newly discovered class I viral fusogen structures. IMPORTANCE Class I viral fusion proteins represent a diverse group of fusogens that catalyze membrane fusion. Although structural studies have focused on those from exogenous viruses, ancient retroviral infections of germ line cells have immortalized ancient fusogens in eukaryotic genomes. These "fossilized" glycoproteins are poorly defined compared to modern fusogens. In this study, we characterized and determined the structure of the porcine endogenous retrovirus fusogen, an ancient retroviral element captured by swine. This fusion protein revealed remarkable alignment to exogenous retroviral fusion proteins, suggesting that fossil fusogens utilize similar structural determinants to perform membrane fusion. Moreover, structural phylogenetic analysis demonstrates that class I viral fusogens cluster into distinct lineages defined by mechanism of membrane fusion. Our results suggest that structural dendrograms can be used to infer mechanistic insights for uncharacterized fusion proteins.

Structure-function analysis of the nsp14 N7-guanine methyltransferase reveals an essential role in Betacoronavirus replication

As coronaviruses (CoVs) replicate in the host cell cytoplasm, they rely on their own capping machinery to ensure the efficient translation of their messenger RNAs (mRNAs), protect them from degradation by cellular 5' exoribonucleases (ExoNs), and escape innate immune sensing. The CoV nonstructural protein 14 (nsp14) is a bifunctional replicase subunit harboring an N-terminal 3'-to-5' ExoN domain and a C-terminal (N7-guanine)-methyltransferase (N7-MTase) domain that is presumably involved in viral mRNA capping. Here, we aimed to integrate structural, biochemical, and virological data to assess the importance of conserved N7-MTase residues for nsp14's enzymatic activities and virus viability. We revisited the crystal structure of severe acute respiratory syndrome (SARS)-CoV nsp14 to perform an in silico comparative analysis between betacoronaviruses. We identified several residues likely involved in the formation of the N7-MTase catalytic pocket, which presents a fold distinct from the Rossmann fold observed in most known MTases. Next, for SARS-CoV and Middle East respiratory syndrome CoV, site-directed mutagenesis of selected residues was used to assess their importance for in vitro enzymatic activity. Most of the engineered mutations abolished N7-MTase activity, while not affecting nsp14-ExoN activity. Upon reverse engineering of these mutations into different betacoronavirus genomes, we identified two substitutions (R310A and F426A in SARS-CoV nsp14) abrogating virus viability and one mutation (H424A) yielding a crippled phenotype across all viruses tested. Our results identify the N7-MTase as a critical enzyme for betacoronavirus replication and define key residues of its catalytic pocket that can be targeted to design inhibitors with a potential pan-coronaviral activity spectrum.

Structural bioinformatic analysis of DsbA proteins and their pathogenicity associated substrates

The disulfide bond (DSB) forming system and in particular DsbA, is a key bacterial oxidative folding catalyst. Due to its role in promoting the correct assembly of a wide range of virulence factors required at different stages of the infection process, DsbA is a master virulence rheostat, making it an attractive target for the development of new virulence blockers. Although DSB systems have been extensively studied across different bacterial species, to date, little is known about how DsbA oxidoreductases are able to recognize and interact with such a wide range of substrates. This review summarizes the current knowledge on the DsbA enzymes, with special attention on their interaction with the partner oxidase DsbB and substrates associated with bacterial virulence. The structurally and functionally diverse set of bacterial proteins that rely on DsbA-mediated disulfide bond formation are summarized. Local sequence and secondary structure elements of these substrates are analyzed to identify common elements recognized by DsbA enzymes. This not only provides information on protein folding systems in bacteria but also offers tools for identifying new DsbA substrates and informs current efforts aimed at developing DsbA targeted anti-microbials.

Genetic and structural analyses of ssRNA viruses pave the way for the discovery of novel antiviral pharmacological targets

In the era of big data and artificial intelligence, a lot of new discoveries have influenced the fields of antiviral drug design and pharmacophore identification. Viruses have always been a threat to society in terms of public health and economic stability. Viruses not only affect humans but also livestock and agriculture with a direct impact on food safety, economy and environmental imprint. Most recently, with the pandemic of COVID-19, it was made clear that a single virus can have a devastating impact on global well-being and economy. In this direction, there is an emerging need for the identification of promising pharmacological targets in viruses. Herein, an effort has been made to discuss the current knowledge, state-of-the-art applications and future implications for the main pharmacological targets of single-stranded RNA viruses.

A Fifth of the Protein World: Rossmann-like Proteins as an Evolutionarily Successful Structural unit

The Rossmann-like fold is the most prevalent and diversified doubly-wound superfold of ancient evolutionary origin. Rossmann-like domains are present in a variety of metabolic enzymes and are capable of binding diverse ligands. Discerning evolutionary relationships among these domains is challenging because of their diverse functions and ancient origin. We defined a minimal Rossmann-like structural motif (RLM), identified RLM-containing domains among known 3D structures (20%) and classified them according to their homologous relationships. New classifications were incorporated into our Evolutionary Classification of protein Domains (ECOD) database. We defined 156 homology groups (H-groups), which were further clustered into 123 possible homology groups (X-groups). Our analysis revealed that RLM-containing proteins constitute approximately 15% of the human proteome. We found that disease-causing mutations are more frequent within RLM domains than within non-RLM domains of these proteins, highlighting the importance of RLM-containing proteins for human health.

β-Strand-mediated interactions of protein domains

Protein domains exist by themselves or in combination with other domains to form complex multidomain proteins. Defining domain boundaries in proteins is essential for understanding their evolution and function but is not trivial. More specifically, partitioning domains that interact by forming a single β-sheet is known to be particularly troublesome for automatic structure-based domain decomposition pipelines. Here, we study edge-to-edge β-strand interactions between domains in a protein chain, to help define the boundaries for some more difficult cases where a single β-sheet spanning over two domains gives an appearance of one. We give a number of examples where β-strands belonging to a single β-sheet do not belong to a single domain and highlight the difficulties of automatic domain parsers on these examples. This work can be used as a baseline for defining domain boundaries in homologous proteins or proteins with similar domain interactions in the future.

Ca2+-based allosteric switches and shape shifting in RGLG1 VWA domain

RGLG1 is an E3 ubiquitin ligase in Arabidopsis thaliana that participates in ABA signaling and regulates apical dominance. Here, we present crystal structures of RGLG1 VWA domain, revealing two novel calcium ions binding sites (NCBS1 and NCBS2). Furthermore, the structures with guided mutagenesis in NCBS1 prove that Ca²⁺ ions play important roles in controlling conformational change of VWA, which is stabilized in open state with Ca²⁺ bound and converted to closed state after Ca²⁺ removal. This allosteric regulation mechanism is distinct from the ever reported one involving the C-terminal helix in integrin α and β I domains. The mutation of a key residue in NCBS2 do not abolish its Ca²⁺-binding potential, with no conformational change. MD simulations reveals that open state of RGLG1 VWA has higher ligand affinity than its closed state, consisting with integrin. Structural comparison of ion-free-MIDAS with Mg²⁺-MIDAS reveals that Mg²⁺ binding to MIDAS does not induce conformational change. With acquisition of first structure of plant VWA domain in both open state and closed state, we carefully analyze the conformational change and propose a totally new paradigm for its transition of open-closed states, which will be of great value for guiding future researches on VWA proteins and their important biological significance.

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.

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.