Unification and extensive diversification of M/Orf3-related ion channel proteins in coronaviruses and other nidoviruses

Mechanism, structural and functional insights into nidovirus-induced double-membrane vesicles

During infection, positive-stranded RNA causes a rearrangement of the host cell membrane, resulting in specialized membrane structure formation aiding viral genome replication. Double-membrane vesicles (DMVs), typical structures produced by virus-induced membrane rearrangements, are platforms for viral replication. Nidoviruses, one of the most complex positive-strand RNA viruses, have the ability to infect not only mammals and a few birds but also invertebrates. Nidoviruses possess a distinctive replication mechanism, wherein their nonstructural proteins (nsps) play a crucial role in DMV biogenesis. With the participation of host factors related to autophagy and lipid synthesis pathways, several viral nsps hijack the membrane rearrangement process of host endoplasmic reticulum (ER), Golgi apparatus, and other organelles to induce DMV formation. An understanding of the mechanisms of DMV formation and its structure and function in the infectious cycle of nidovirus may be essential for the development of new and effective antiviral strategies in the future.

Ribonuclease T2 represents a distinct circularly permutated version of the BECR RNases

Detection of homologous relationships among proteins and understanding their mechanisms of diversification are major topics in the fields of protein science, bioinformatics, and phylogenetics. Recent developments in sequence/profile-based and structural similarity-based methods have greatly facilitated the unification and classification of many protein families into superfamilies or folds, yet many proteins remain unclassified in current protein databases. As one of the three earliest identified RNases in biology, ribonuclease T2, also known as RNase I in Escherichia coli, RNase Rh in fungi, or S-RNase in plant, is thought to be an ancient RNase family due to its widespread distribution and distinct structure. In this study, we present evidence that RNase T2 represents a circularly permutated version of the BECR (Barnase-EndoU-Colicin E5/D-RelE) fold RNases. This subtle relationship cannot be detected by traditional methods such as sequence/profile-based comparisons, structure-similarity searches, and circular permutation detections. However, we were able to identify the structural similarity using rational reconstruction of a theoretical RNase T2 ancestor via a reverse circular permutation process, followed by structural modeling using AlphaFold2, and structural comparisons. This relationship is further supported by the fact that RNase T2 and other typical BECR RNases, namely Colicin D, RNase A, and BrnT, share similar catalytic site configurations, all involving an analogous set of conserved residues on the α0 helix and the β4 strand of the BECR fold. This study revealed a hidden root of RNase T2 in bacterial toxin systems and demonstrated that reconstruction and modeling of ancestral topology is an effective strategy to identify remote relationship between proteins.

Using Alphafold2 to Predict the Structure of the Gp5/M Dimer of Porcine Respiratory and Reproductive Syndrome Virus

Porcine reproductive and respiratory syndrome virus is a positive-stranded RNA virus of the family Arteriviridae. The Gp5/M dimer, the major component of the viral envelope, is required for virus budding and is an antibody target. We used alphafold2, an artificial-intelligence-based system, to predict a credible structure of Gp5/M. The short disulfide-linked ectodomains lie flat on the membrane, with the exception of the erected N-terminal helix of Gp5, which contains the antibody epitopes and a hypervariable region with a changing number of carbohydrates. The core of the dimer consists of six curved and tilted transmembrane helices, and three are from each protein. The third transmembrane regions extend into the cytoplasm as amphiphilic helices containing the acylation sites. The endodomains of Gp5 and M are composed of seven β-strands from each protein, which interact via β-strand seven. The area under the membrane forms an open cavity with a positive surface charge. The M and Orf3a proteins of coronaviruses have a similar structure, suggesting that all four proteins are derived from the same ancestral gene. Orf3a, like Gp5/M, is acylated at membrane-proximal cysteines. The role of Gp5/M during virus replication is discussed, in particular the mechanisms of virus budding and models of antibody-dependent virus neutralization.

Structure of SARS-CoV-2 M protein in lipid nanodiscs

SARS-CoV-2 encodes four structural proteins incorporated into virions, spike (S), envelope (E), nucleocapsid (N), and membrane (M). M plays an essential role in viral assembly by organizing other structural proteins through physical interactions and directing them to sites of viral budding. As the most abundant protein in the viral envelope and a target of patient antibodies, M is a compelling target for vaccines and therapeutics. Still, the structure of M and molecular basis for its role in virion formation are unknown. Here, we present the cryo-EM structure of SARS-CoV-2 M in lipid nanodiscs to 3.5 Å resolution. M forms a 50 kDa homodimer that is structurally related to the SARS-CoV-2 ORF3a viroporin, suggesting a shared ancestral origin. Structural comparisons reveal how intersubunit gaps create a small, enclosed pocket in M and large open cavity in ORF3a, consistent with a structural role and ion channel activity, respectively. M displays a strikingly electropositive cytosolic surface that may be important for interactions with N, S, and viral RNA. Molecular dynamics simulations show a high degree of structural rigidity in a simple lipid bilayer and support a role for M homodimers in scaffolding viral assembly. Together, these results provide insight into roles for M in coronavirus assembly and structure.

Structure of SARS-CoV-2 membrane protein essential for virus assembly

The coronavirus membrane protein (M) is the most abundant viral structural protein and plays a central role in virus assembly and morphogenesis. However, the process of M protein-driven virus assembly are largely unknown. Here, we report the cryo-electron microscopy structure of the SARS-CoV-2 M protein in two different conformations. M protein forms a mushroom-shaped dimer, composed of two transmembrane domain-swapped three-helix bundles and two intravirion domains. M protein further assembles into higher-order oligomers. A highly conserved hinge region is key for conformational changes. The M protein dimer is unexpectedly similar to SARS-CoV-2 ORF3a, a viral ion channel. Moreover, the interaction analyses of M protein with nucleocapsid protein (N) and RNA suggest that the M protein mediates the concerted recruitment of these components through the positively charged intravirion domain. Our data shed light on the M protein-driven virus assembly mechanism and provide a structural basis for therapeutic intervention targeting M protein.

M protein plays essential roles in virus assembly and morphogenesis. Here, authors reveal two cryo-EM structures of M protein from SARS-CoV-2 that suggest conformational dynamics of M protein and its role in virus assembly.

Discovering Biological Conflict Systems Through Genome Analysis: Evolutionary Principles and Biochemical Novelty

Biological replicators, from genes within a genome to whole organisms, are locked in conflicts. Comparative genomics has revealed a staggering diversity of molecular armaments and mechanisms regulating their deployment, collectively termed biological conflict systems. These encompass toxins used in inter- and intraspecific interactions, self/nonself discrimination, antiviral immune mechanisms, and counter-host effectors deployed by viruses and intragenomic selfish elements. These systems possess shared syntactical features in their organizational logic and a set of effectors targeting genetic information flow through the Central Dogma, certain membranes, and key molecules like NAD⁺. These principles can be exploited to discover new conflict systems through sensitive computational analyses. This has led to significant advances in our understanding of the biology of these systems and furnished new biotechnological reagents for genome editing, sequencing, and beyond. We discuss these advances using specific examples of toxins, restriction-modification, apoptosis, CRISPR/second messenger-regulated systems, and other enigmatic nucleic acid-targeting systems. Expected final online publication date for the Annual Review of Biomedical Data Science, Volume 5 is August 2022. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

Homology-based classification of accessory proteins in coronavirus genomes uncovers extremely dynamic evolution of gene content

Coronaviruses (CoVs) have complex genomes that encode a fixed array of structural and nonstructural components, as well as a variety of accessory proteins that differ even among closely related viruses. Accessory proteins often play a role in the suppression of immune responses and may represent virulence factors. Despite their relevance for CoV phenotypic variability, information on accessory proteins is fragmentary. We applied a systematic approach based on homology detection to create a comprehensive catalogue of accessory proteins encoded by CoVs. Our analyses grouped accessory proteins into 379 orthogroups and 12 super‐groups. No orthogroup was shared by the four CoV genera and very few were present in all or most viruses in the same genus, reflecting the dynamic evolution of CoV genomes. We observed differences in the distribution of accessory proteins in CoV genera. Alphacoronaviruses harboured the largest diversity of accessory open reading frames (ORFs), deltacoronaviruses the smallest. However, the average number of accessory proteins per genome was highest in betacoronaviruses. Analysis of the evolutionary history of some orthogroups indicated that the different CoV genera adopted similar evolutionary strategies. Thus, alphacoronaviruses and betacoronaviruses acquired phosphodiesterases and spike‐like accessory proteins independently, whereas horizontal gene transfer from reoviruses endowed betacoronaviruses and deltacoronaviruses with fusion‐associated small transmembrane (FAST) proteins. Finally, analysis of accessory ORFs in annotated CoV genomes indicated ambiguity in their naming. This complicates cross‐communication among researchers and hinders automated searches of large data sets (e.g., PubMed, GenBank). We suggest that orthogroup membership is used together with a naming system to provide information on protein function.

An issue of concern: unique truncated ORF8 protein variants of SARS-CoV-2

Open reading frame 8 (ORF8) shows one of the highest levels of variability among accessory proteins in Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2), the causative agent of Coronavirus Disease 2019 (COVID-19). It was previously reported that the ORF8 protein inhibits the presentation of viral antigens by the major histocompatibility complex class I (MHC-I), which interacts with host factors involved in pulmonary inflammation. The ORF8 protein assists SARS-CoV-2 in evading immunity and plays a role in SARS-CoV-2 replication. Among many contributing mutations, Q27STOP, a mutation in the ORF8 protein, defines the B.1.1.7 lineage of SARS-CoV-2, engendering the second wave of COVID-19. In the present study, 47 unique truncated ORF8 proteins (T-ORF8) with the Q27STOP mutations were identified among 49,055 available B.1.1.7 SARS-CoV-2 sequences. The results show that only one of the 47 T-ORF8 variants spread to over 57 geo-locations in North America, and other continents, which include Africa, Asia, Europe and South America. Based on various quantitative features, such as amino acid homology, polar/non-polar sequence homology, Shannon entropy conservation, and other physicochemical properties of all specific 47 T-ORF8 protein variants, nine possible T-ORF8 unique variants were defined. The question as to whether T-ORF8 variants function similarly to the wild type ORF8 is yet to be investigated. A positive response to the question could exacerbate future COVID-19 waves, necessitating severe containment measures.

Understanding the Role of SARS-CoV-2 ORF3a in Viral Pathogenesis and COVID-19

The ongoing SARS-CoV-2 pandemic has shocked the world due to its persistence, COVID-19-related morbidity and mortality, and the high mutability of the virus. One of the major concerns is the emergence of new viral variants that may increase viral transmission and disease severity. In addition to mutations of spike protein, mutations of viral proteins that affect virulence, such as ORF3a, also must be considered. The purpose of this article is to review the current literature on ORF3a, to summarize the molecular actions of SARS-CoV-2 ORF3a, and its role in viral pathogenesis and COVID-19. ORF3a is a polymorphic, multifunctional viral protein that is specific to SARS-CoV/SARS-CoV-2. It was acquired from β-CoV lineage and likely originated from bats through viral evolution. SARS-CoV-2 ORF3a is a viroporin that interferes with ion channel activities in host plasma and endomembranes. It is likely a virion-associated protein that exerts its effect on the viral life cycle during viral entry through endocytosis, endomembrane-associated viral transcription and replication, and viral release through exocytosis. ORF3a induces cellular innate and pro-inflammatory immune responses that can trigger a cytokine storm, especially under hypoxic conditions, by activating NLRP3 inflammasomes, HMGB1, and HIF-1α to promote the production of pro-inflammatory cytokines and chemokines. ORF3a induces cell death through apoptosis, necrosis, and pyroptosis, which leads to tissue damage that affects the severity of COVID-19. ORF3a continues to evolve along with spike and other viral proteins to adapt in the human cellular environment. How the emerging ORF3a mutations alter the function of SARS-CoV-2 ORF3a and its role in viral pathogenesis and COVID-19 is largely unknown. This review provides an in-depth analysis of ORF3a protein's structure, origin, evolution, and mutant variants, and how these characteristics affect its functional role in viral pathogenesis and COVID-19.

SARS-CoV-2-Specific Immune Response and the Pathogenesis of COVID-19

The review aims to consolidate research findings on the molecular mechanisms and virulence and pathogenicity characteristics of coronavirus disease (COVID-19) causative agent, severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), and their relevance to four typical stages in the development of acute viral infection. These four stages are invasion; primary blockade of antiviral innate immunity; engagement of the virus's protection mechanisms against the factors of adaptive immunity; and acute, long-term complications of COVID-19. The invasion stage entails the recognition of the spike protein (S) of SARS-CoV-2 target cell receptors, namely, the main receptor (angiotensin-converting enzyme 2, ACE2), its coreceptors, and potential alternative receptors. The presence of a diverse repertoire of receptors allows SARS-CoV-2 to infect various types of cells, including those not expressing ACE2. During the second stage, the majority of the polyfunctional structural, non-structural, and extra proteins SARS-CoV-2 synthesizes in infected cells are involved in the primary blockage of antiviral innate immunity. A high degree of redundancy and systemic action characterizing these pathogenic factors allows SARS-CoV-2 to overcome antiviral mechanisms at the initial stages of invasion. The third stage includes passive and active protection of the virus from factors of adaptive immunity, overcoming of the barrier function at the focus of inflammation, and generalization of SARS-CoV-2 in the body. The fourth stage is associated with the deployment of variants of acute and long-term complications of COVID-19. SARS-CoV-2's ability to induce autoimmune and autoinflammatory pathways of tissue invasion and development of both immunosuppressive and hyperergic mechanisms of systemic inflammation is critical at this stage of infection.

Photoglobin, a distinct family of non-heme binding globins, defines a potential photosensor in prokaryotic signal transduction systems

Globins constitute an ancient superfamily of proteins, exhibiting enormous structural and functional diversity, as demonstrated by many heme-binding families and two non-heme binding families that were discovered in bacterial stressosome component RsbR and in light-harvesting phycobiliproteins (phycocyanin) in cyanobacteria and red algae. By comprehensively exploring the globin repertoire using sensitive computational analyses of sequences, structures, and genomes, we present the identification of the third family of non-heme binding globins—the photoglobin. By conducting profile-based comparisons, clustering analyses, and structural modeling, we demonstrate that photoglobin is related to, but distinct from, the phycocyanin family. Photoglobin preserves a potential ligand-binding pocket, whose residue configuration closely resembles that of phycocyanin, indicating that photoglobin potentially binds to a comparable linear tetrapyrrole. By exploring the contextual information provided by the photoglobin’s domain architectures and gene-neighborhoods, we found that photoglobin is frequently associated with the B12-binding light sensor domain and many domains typical of prokaryotic signal transduction systems. Structural modeling using AlphaFold2 demonstrated that photoglobin and B12-binding domains form a structurally conserved hub among different domain architecture contexts. Based on these strong associations, we predict that the coupled photoglobin and B12-binding domains act as a light-sensing regulatory bundle, with each domain sensing different wavelengths of light resulting in switch-like regulation of downstream signaling effectors. Thus, based on the above lines of evidence, we present a distinct non-heme binding globin family and propose that it may define a new type of light sensor, by means of a linear tetrapyrrole, in complex prokaryotic signal transduction systems.

An issue of concern: unique truncated ORF8 protein variants of SARS-CoV-2

Open reading frame 8 (ORF8) shows one of the highest levels of variability among accessory proteins in Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2), the causative agent of Coronavirus Disease 2019 (COVID-19). It was previously reported that the ORF8 protein inhibits the presentation of viral antigens by the major histocompatibility complex class I (MHC-I), which interacts with host factors involved in pulmonary inflammation. The ORF8 protein assists SARS-CoV-2 in evading immunity and plays a role in SARS-CoV-2 replication. Among many contributing mutations, Q27STOP, a mutation in the ORF8 protein, defines the B.1.1.7 lineage of SARS-CoV-2, engendering the second wave of COVID-19. In the present study, 47 unique truncated ORF8 proteins (T-ORF8) with the Q27STOP mutations were identified among 49,055 available B.1.1.7 SARS-CoV-2 sequences. The results show that only one of the 47 T-ORF8 variants spread to over 57 geo-locations in North America, and other continents, which include Africa, Asia, Europe and South America. Based on various quantitative features, such as amino acid homology, polar/non-polar sequence homology, Shannon entropy conservation, and other physicochemical properties of all specific 47 T-ORF8 protein variants, nine possible T-ORF8 unique variants were defined. The question as to whether T-ORF8 variants function similarly to the wild type ORF8 is yet to be investigated. A positive response to the question could exacerbate future COVID-19 waves, necessitating severe containment measures.

Comparative Phylogenomic Analysis Reveals Evolutionary Genomic Changes and Novel Toxin Families in Endophytic Liberibacter Pathogens

Liberibacter pathogens are the causative agents of several severe crop diseases worldwide, including citrus Huanglongbing and potato zebra chip. These bacteria are endophytic and nonculturable, which makes experimental approaches challenging and highlights the need for bioinformatic analysis in advancing our understanding about Liberibacter pathogenesis. Here, we performed an in-depth comparative phylogenomic analysis of the Liberibacter pathogens and their free-living, nonpathogenic, ancestral species, aiming to identify major genomic changes and determinants associated with their evolutionary transitions in living habitats and pathogenicity. Using gene neighborhood analysis and phylogenetic classification, we systematically uncovered, annotated, and classified all prophage loci into four types, including one previously unrecognized group. We showed that these prophages originated through independent gene transfers at different evolutionary stages of Liberibacter and only the SC-type prophage was associated with the emergence of the pathogens. Using ortholog clustering, we vigorously identified two additional sets of genomic genes, which were either lost or gained in the ancestor of the pathogens. Consistent with the habitat change, the lost genes were enriched for biosynthesis of cellular building blocks. Importantly, among the gained genes, we uncovered several previously unrecognized toxins, including new toxins homologous to the EspG/VirA effectors, a YdjM phospholipase toxin, and a secreted endonuclease/exonuclease/phosphatase (EEP) protein. Our results substantially extend the knowledge of the evolutionary events and potential determinants leading to the emergence of endophytic, pathogenic Liberibacter species, which will facilitate the design of functional experiments and the development of new methods for detection and blockage of these pathogens. IMPORTANCELiberibacter pathogens are associated with several severe crop diseases, including citrus Huanglongbing, the most destructive disease to the citrus industry. Currently, no effective cure or treatments are available, and no resistant citrus variety has been found. The fact that these obligate endophytic pathogens are not culturable has made it extremely challenging to experimentally uncover the genes/proteins important to Liberibacter pathogenesis. Further, earlier bioinformatics studies failed to identify key genomic determinants, such as toxins and effector proteins, that underlie the pathogenicity of the bacteria. In this study, an in-depth comparative genomic analysis of Liberibacter pathogens along with their ancestral nonpathogenic species identified the prophage loci and several novel toxins that are evolutionarily associated with the emergence of the pathogens. These results shed new light on the disease mechanism of Liberibacter pathogens and will facilitate the development of new detection and blockage methods targeting the toxins.