Understanding how transmembrane domains regulate interactions between human BST-2 and the SARS-CoV-2 accessory protein ORF7a

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the causative agent of COVID, replicates at intracellular membranes. Bone marrow stromal antigen 2 (BST-2; tetherin) is an antiviral response protein that inhibits transport of viral particles after budding within infected cells. RNA viruses such as SARS-CoV-2 use various strategies to disable BST-2, including use of transmembrane ‘accessory’ proteins that interfere with BST-2 oligomerization. ORF7a is a small, transmembrane protein present in SARS-CoV-2 shown previously to alter BST-2 glycosylation and function. In this study, we investigated the structural basis for BST-2 ORF7a interactions, with a particular focus on transmembrane and juxtamembrane interactions. Our results indicate that transmembrane domains play an important role in BST-2 ORF7a interactions and mutations to the transmembrane domain of BST-2 can alter these interactions, particularly single-nucleotide polymorphisms in BST-2 that result in mutations such as I28S. Using molecular dynamics simulations, we identified specific interfaces and interactions between BST-2 and ORF7a to develop a structural basis for the transmembrane interactions. Differences in glycosylation are observed for BST-2 transmembrane mutants interacting with ORF7a, consistent with the idea that transmembrane domains play a key role in their heterooligomerization. Overall, our results indicate that ORF7a transmembrane domain interactions play a key role along with extracellular and juxtamembrane domains in modulating BST-2 function.     

Sequence to structure analysis of the ORF4 protein from Hepatitis E virus

Hepatitis E virus (HEV) is the main cause of acute hepatitis worldwide. HEV accounts for up to 30% mortality rate in pregnant women, with highest incidences reported for genotype 1 (G1) HEV. The contributing factors in adverse cases during pregnancy in women due to HEV infection is still debated. The mechanism underlying the pathogenesis of viral infection is attributed to different genomic component of HEV, i.e., open reading frames (ORFs): ORF1, ORF2, ORF3 and ORF4. Recently, ORF4 has been discovered in enhancing the replication of GI isolates of HEV through regulation of an IRES-like RNA element. However, its characterization through computational methodologies remains unexplored. In this novel study, we provide comprehensive overview of ORF4 protein''s genetic and molecular characteristics through analyzing its sequence and different structural levels. A total of three different datasets (Human, Rat and Ferret) of ORF4 genomes were built and comparatively analyzed. Several non-synonymous mutations in conjunction with higher entropy values were observed in rat and ferret datasets, however, limited variation was observed in human ORF4 genomes. Higher transition to tranversion ratio was observed in the ORF4 genomes. Studies have reported the association of intrinsic disordered proteins (IDP) with drug discovery due to its role in several signaling and regulatory processes through protein-protein interactions (PPIs). As PPIs are potent drug target sources, thus the ORF4 protein was explored by analyzing its polypeptide structure in order to shed light on its intrinsic disorder. Pressures that lead towards preponderance of disordered-promoting amino acid residues shaped the evolution of ORF4. The intrinsic disorder propensity analysis revealed ORF4 protein (Human) as a highly disordered protein (IDP). Predominance of coils and lack of secondary structure further substantiated our findings suggesting its involvement in binding to ligand molecules. Thus, ORF4 contributes to cellular signaling processes through protein-protein interactions, as IDPs are targets for regulation to accelerate the process of drug designing strategies against HEV infections.     

Variations in Orf3a protein of SARS-CoV-2 alter its structure and function

Severe acquired respiratory syndrome coronavirus 2 (SARS-CoV-2) rapidly spread worldwide and acquired multiple mutations in its genome. Orf3a, an accessory protein encoded by the genome of SARS-CoV-2, plays a significant role in viral infection and pathogenesis. In the present in-silico study, 15,928 sequences of Orf3a reported worldwide were compared to identify variations in this protein. Our analysis revealed the occurrence of mutations at 173 residues of Orf3a protein. Subsequently, protein modelling was performed that revealed twelve mutations which can considerably affect the stability of Orf3a. Among the 12 mutations, three mutations (Y160H, D210Y and S171L) also lead to alterations in secondary structure and protein disorder parameters of the Orf3a protein. Further, we used predictive tools to identify five promising epitopes of B-cells, which resides in the mutated regions of Orf3a. Altogether, our study sheds light on the variations occurring in Orf3a that might contribute to alteration in protein structure and function.     

Host E3 ligase HUWE1 attenuates the proapoptotic activity of the MERS-CoV accessory protein ORF3 by promoting its ubiquitin-dependent degradation

With the outbreak of severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2), coronaviruses have begun to attract great attention across the world. Of the known human coronaviruses, however, Middle East respiratory syndrome coronavirus (MERS-CoV) is the most lethal. Coronavirus proteins can be divided into three groups: nonstructural proteins, structural proteins, and accessory proteins. While the number of each of these proteins varies greatly among different coronaviruses, accessory proteins are most closely related to the pathogenicity of the virus. We found for the first time that the ORF3 accessory protein of MERS-CoV, which closely resembles the ORF3a proteins of severe acute respiratory syndrome coronavirus and SARS-CoV-2, has the ability to induce apoptosis in cells in a dose-dependent manner. Through bioinformatics analysis and validation, we revealed that ORF3 is an unstable protein and has a shorter half-life in cells compared to that of severe acute respiratory syndrome coronavirus and SARS-CoV-2 ORF3a proteins. After screening, we identified a host E3 ligase, HUWE1, that specifically induces MERS-CoV ORF3 protein ubiquitination and degradation through the ubiquitin–proteasome system. This results in the diminished ability of ORF3 to induce apoptosis, which might partially explain the lower spread of MERS-CoV compared to other coronaviruses. In summary, this study reveals a pathological function of MERS-CoV ORF3 protein and identifies a potential host antiviral protein, HUWE1, with an ability to antagonize MERS-CoV pathogenesis by inducing ORF3 degradation, thus enriching our knowledge of the pathogenesis of MERS-CoV and suggesting new targets and strategies for clinical development of drugs for MERS-CoV treatment.     

A conservation and rigidity based method for detecting critical protein residues

Background

Certain amino acids in proteins play a critical role in determining their structural stability and function. Examples include flexible regions such as hinges which allow domain motion, and highly conserved residues on functional interfaces which allow interactions with other proteins. Detecting these regions can aid in the analysis and simulation of protein rigidity and conformational changes, and helps characterizing protein binding and docking. We present an analysis of critical residues in proteins using a combination of two complementary techniques. One method performs in-silico mutations and analyzes the protein's rigidity to infer the role of a point substitution to Glycine or Alanine. The other method uses evolutionary conservation to find functional interfaces in proteins.

Results

We applied the two methods to a dataset of proteins, including biomolecules with experimentally known critical residues as determined by the free energy of unfolding. Our results show that the combination of the two methods can detect the vast majority of critical residues in tested proteins.

Conclusions

Our results show that the combination of the two methods has the potential to detect more information than each method separately. Future work will provide a confidence level for the criticalness of a residue to improve the accuracy of our method and eliminate false positives. Once the combined methods are integrated into one scoring function, it can be applied to other domains such as estimating functional interfaces.

     

Differential binding of SARS-CoV-2 Spike protein variants to its cognate receptor hACE2 using molecular modeling based binding analysis

The current emergence of novel coronavirus, SARS-CoV-2 and its ceaseless expansion worldwide has posed a global health emergency that has adversely affected the humans. With the entire world striving to understand the newly emerged virus, differences in morbidity and infection rate of SARS-CoV-2 have been observed across varied geographic areas, which have been ascribed to viral mutation and evolution over time. The homotrimeric Spike (S) glycoprotein on the viral envelope surface is responsible for binding, priming, and initiating infection in the host. Our phylogeny analysis of 1947 sequences of S proteins indicated there is a change in amino acid (aa) from aspartate (Group-A) to glycine (Group-B) at position 614, near the receptor- binding domain (RBD; aa positions 331-524). The two variants are reported to be in circulation, disproportionately across the world, with Group-A dominant in Asia and Group-B in North America. The trimeric, monomeric, and RBD of S protein of both the variant groups (A & B) were modeled using the Swiss-Model server and were docked with the human receptor angiotensin-converting enzyme 2 (hACE2) employing the PatchDock server and visualized in PyMol. Group-A S protein''s RBD bound imperceptibly to the two binding clefts of the hACE2 protein, on the other hand, Group-B S protein''s RBD perfectly interacted inside the binding clefts of hACE2, with higher number of hydrogen and hydrophobic interactions. This implies that the S protein''s amino acid at position 614 near the core RBD influences its interaction with the cognate hACE2 receptor, which may induce its infectivity that should be explored further with molecular and biochemical studies.     

Insights on cross-species transmission of SARS-CoV-2 from structural modeling

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is responsible for the ongoing global pandemic that has infected more than 31 million people in more than 180 countries worldwide. Like other coronaviruses, SARS-CoV-2 is thought to have been transmitted to humans from wild animals. Given the scale and widespread geographical distribution of the current pandemic and confirmed cases of cross-species transmission, the question of the extent to which this transmission is possible emerges, as well as what molecular features distinguish susceptible from non-susceptible animal species. Here, we investigated the structural properties of several ACE2 orthologs bound to the SARS-CoV-2 spike protein. We found that species known not to be susceptible to SARS-CoV-2 infection have non-conservative mutations in several ACE2 amino acid residues that disrupt key polar and charged contacts with the viral spike protein. Our models also allow us to predict affinity-enhancing mutations that could be used to design ACE2 variants for therapeutic purposes. Finally, our study provides a blueprint for modeling viral-host protein interactions and highlights several important considerations when designing these computational studies and analyzing their results.     

The landscape of antibody binding in SARS-CoV-2 infection

The search for potential antibody-based diagnostics, vaccines, and therapeutics for pandemic severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has focused almost exclusively on the spike (S) and nucleocapsid (N) proteins. Coronavirus membrane (M), ORF3a, and ORF8 proteins are humoral immunogens in other coronaviruses (CoVs) but remain largely uninvestigated for SARS-CoV-2. Here, we use ultradense peptide microarray mapping to show that SARS-CoV-2 infection induces robust antibody responses to epitopes throughout the SARS-CoV-2 proteome, particularly in M, in which 1 epitope achieved excellent diagnostic accuracy. We map 79 B cell epitopes throughout the SARS-CoV-2 proteome and demonstrate that antibodies that develop in response to SARS-CoV-2 infection bind homologous peptide sequences in the 6 other known human CoVs. We also confirm reactivity against 4 of our top-ranking epitopes by enzyme-linked immunosorbent assay (ELISA). Illness severity correlated with increased reactivity to 9 SARS-CoV-2 epitopes in S, M, N, and ORF3a in our population. Our results demonstrate previously unknown, highly reactive B cell epitopes throughout the full proteome of SARS-CoV-2 and other CoV proteins.

Profiling of antibody binding from naïve and COVID-19 convalescent human sera to the entire proteome of SARS-CoV-2 and other human, bat and pangolin coronaviruses identifies 79 B cell epitopes throughout the SARS-CoV-2 proteome, finding that the most sensitive and specific binding occurred in the membrane (M) protein, and revealing cross-reactivity patterns.     

Using affinity chromatography to engineer and characterize pH‐dependent protein switches

Conformational changes play important roles in the regulation of many enzymatic reactions. Specific motions of side chains, secondary structures, or entire protein domains facilitate the precise control of substrate selection, binding, and catalysis. Likewise, the engineering of allostery into proteins is envisioned to enable unprecedented control of chemical reactions and molecular assembly processes. We here study the structural effects of engineered ionizable residues in the core of the glutathione‐S‐transferase to convert this protein into a pH‐dependent allosteric protein. The underlying rational of these substitutions is that in the neutral state, an uncharged residue is compatible with the hydrophobic environment. In the charged state, however, the residue will invoke unfavorable interactions, which are likely to induce conformational changes that will affect the function of the enzyme. To test this hypothesis, we have engineered a single aspartate, cysteine, or histidine residue at a distance from the active site into the protein. All of the mutations exhibit a dramatic effect on the protein's affinity to bind glutathione. Whereas the aspartate or histidine mutations result in permanently nonbinding or binding versions of the protein, respectively, mutant GST50C exhibits distinct pH‐dependent GSH‐binding affinity. The crystal structures of the mutant protein GST50C under ionizing and nonionizing conditions reveal the recruitment of water molecules into the hydrophobic core to produce conformational changes that influence the protein's active site. The methodology described here to create and characterize engineered allosteric proteins through affinity chromatography may lead to a general approach to engineer effector‐specific allostery into a protein structure.     

Identification and expression localization of a Ran homologue in mollusc abalone, Haliotis diversicolor supertexta

Ran protein is a central molecular in several key nuclear functions, including nucleocytoplasmic transport, cell-cycle progression and nuclear envelope assembly. In this study, we have isolated and characterized a Ran homologue from a gastropod abalone which we named ab-Ran. The full-length cDNA consists of 1239 bp with an ORF encoding a 220 amino acid protein. The deduced amino acid sequence of ab-Ran shows highly similar to that of other Ran members (84-88%). Moreover, the ab-Ran contains five conserved regions and four carboxy-terminal residues CAAX-box. RT-PCR analysis showed that the ab-Ran was ubiquitously expressed in abalone tissues. The intracellular localization examined by immunofluorescence and immunohistochemistry staining displayed that ab-Ran was largely concentrated in the nuclei and partially in the cytoplasm. To the best of our knowledge, this is the first identification and characterization of a Ran homologue in mollusk.     

The role of S477N mutation in the molecular behavior of SARS-CoV-2 spike protein: An in-silico perspective

The attachment of SARA-CoV-2 happens between ACE2 and the receptor binding domain (RBD) on the spike protein. Mutations in this domain can affect the binding affinity of the spike protein for ACE2. S477N, one of the most common mutations reported in the recent variants, is located in the RBD. Today's computational approaches in biology, especially during the SARS-CoV-2 pandemic, assist researchers in predicting a protein's behavior in contact with other proteins in more detail. In this study, we investigated the interactions of the S477N-hACE2 in silico to find the impact of this mutation on its binding affinity for ACE2 and immunity responses using dynamics simulation, protein–protein docking, and immunoinformatics methods. Our computational analysis revealed an increased binding affinity of N477 for ACE2. Four new hydrogen and hydrophobic bonds in the mutant RBD-ACE2 were formed (with S19 and Q24 of ACE2), which do not exist in the wild type. Also, the protein spike structure in this mutation was associated with an increase in stabilization and a decrease in its fluctuations at the atomic level. N477 mutation can be considered as the cause of increased escape from the immune system through MHC-II.     

Great interactions: How binding incorrect partners can teach us about protein recognition and function

Protein–protein interactions play a key part in most biological processes and understanding their mechanism is a fundamental problem leading to numerous practical applications. The prediction of protein binding sites in particular is of paramount importance since proteins now represent a major class of therapeutic targets. Amongst others methods, docking simulations between two proteins known to interact can be a useful tool for the prediction of likely binding patches on a protein surface. From the analysis of the protein interfaces generated by a massive cross‐docking experiment using the 168 proteins of the Docking Benchmark 2.0, where all possible protein pairs, and not only experimental ones, have been docked together, we show that it is also possible to predict a protein's binding residues without having any prior knowledge regarding its potential interaction partners. Evaluating the performance of cross‐docking predictions using the area under the specificity‐sensitivity ROC curve (AUC) leads to an AUC value of 0.77 for the complete benchmark (compared to the 0.5 AUC value obtained for random predictions). Furthermore, a new clustering analysis performed on the binding patches that are scattered on the protein surface show that their distribution and growth will depend on the protein's functional group. Finally, in several cases, the binding‐site predictions resulting from the cross‐docking simulations will lead to the identification of an alternate interface, which corresponds to the interaction with a biomolecular partner that is not included in the original benchmark. Proteins 2016; 84:1408–1421. © 2016 The Authors Proteins: Structure, Function, and Bioinformatics Published by Wiley Periodicals, Inc.     

The DbpA catalytic core unwinds double-helix substrates by directly loading on them

DbpA is a DEAD-box RNA helicase implicated in the assembly of the large ribosomal subunit. Similar to all the members of the DEAD-box family, the DbpA protein has two N-terminal RecA-like domains, which perform the RNA unwinding. However, unlike other members of this family, the DbpA protein also possesses a structured C-terminal RNA-binding domain that mediates specific tethering of DbpA to hairpin 92 of the Escherichia coli 23S ribosomal RNA. Previous studies using model RNA molecules containing hairpin 92 show that the RNA molecules support the DbpA protein''s double-helix unwinding activity, provided that the double helix has a 3′ single-stranded region. The 3′ single-stranded region was suggested to be the start site of the DbpA protein''s catalytic unwinding activity. The data presented here demonstrate that the single-stranded region 3′ of the double-helix substrate is not required for the DbpA protein''s unwinding activity and the DbpA protein unwinds the double-helix substrates by directly loading on them.     

Modification of protein crystal packing by systematic mutations of surface residues: Implications on biotemplating and crystal porosity

Bioinspired nano‐scale biotemplating for the development of novel composite materials has recently culminated in several demonstrations of nano‐structured hybrid materials. Protein crystals, routinely prepared for the elucidation of protein 3D structures by X‐ray crystallography, present an ordered and highly accurate 3D array of protein molecules. Inherent to the 3D arrangement of the protein “building blocks” in the crystal, a complementary 3D array of interconnected cavities—voids array, exhibiting highly ordered porosity is formed. The porous arrays of protein crystal may serve as a nano‐structured, accurate biotemplate by a “filling” process. These cavities arrays are shaped by the mode of protein packing throughout the crystallization process. Here we propose and demonstrate feasibility of targeting site specific mutations to modify protein's surface to affect protein crystal packing, enabling the generation of a series of protein crystal “biotemplates” all originating from same parent protein. The selection of these modification sites was based on in silico analysis of protein–protein interface contact areas in the parent crystal. The model protein selected for this study was the N‐terminal type II cohesin from the cellulosomal scaffold in ScaB subunit of Acetivibrio cellulolyticus and mutations were focused on lysine residues involved in protein packing as prime target. The impact of systematically mutating these lysine residues on protein packing and its resulting interconnected cavities array were found to be most significant when surface lysine residues were substituted to tryptophan residues. Our results demonstrate the feasibility of using pre‐designed site directed mutations for the generation of a series of protein crystal biotemplates from a “parent” protein. Biotechnol. Bioeng. 2009; 104: 444–457 © 2009 Wiley Periodicals, Inc.     

Self‐association of TPR domains: Lessons learned from a designed,consensus‐based TPR oligomer

The tetratricopeptide repeat (TPR) motif is a protein–protein interaction module that acts as an organizing centre for complexes regulating a multitude of biological processes. Despite accumulating evidence for the formation of TPR oligomers as an additional level of regulation there is a lack of structural and solution data explaining TPR self‐association. In the present work we characterize the trimeric TPR‐containing protein YbgF, which is linked to the Tol system in Gram‐negative bacteria. By subtracting previously identified TPR consensus residues required for stability of the fold from residues conserved across YbgF homologs, we identified residues involved in oligomerization of the C‐terminal YbgF TPR domain. Crafting these residues, which are located in loop regions between TPR motifs, onto the monomeric consensus TPR protein CTPR3 induced the formation of oligomers. The crystal structure of this engineered oligomer shows an asymmetric trimer where stacking interactions between the introduced tyrosines and displacement of the C‐terminal hydrophilic capping helix, present in most TPR domains, are key to oligomerization. Asymmetric trimerization of the YbgF TPR domain and CTPR3Y3 leads to the formation of higher order oligomers both in the crystal and in solution. However, such open‐ended self‐association does not occur in full‐length YbgF suggesting that the protein's N‐terminal coiled‐coil domain restricts further oligomerization. This interpretation is borne out in experiments where the coiled‐coil domain of YbgF was engineered onto the N‐terminus of CTPR3Y3 and shown to block self‐association beyond trimerization. Our study lays the foundations for understanding the structural basis for TPR domain self‐association and how such self‐association can be regulated in TPR domain‐containing proteins. Proteins 2010. © 2010 Wiley‐Liss, Inc.     

Sequence,biophysical, and structural analyses of the PstS lipoprotein (BB0215) from Borrelia burgdorferi reveal a likely binding component of an ABC‐type phosphate transporter

The Lyme disease agent Borrelia burgdorferi, which is transmitted via a tick vector, is dependent on its tick and mammalian hosts for a number of essential nutrients. Like other bacterial diderms, it must transport these biochemicals from the extracellular milieu across two membranes, ultimately to the B. burgdorferi cytoplasm. In the current study, we established that a gene cluster comprising genes bb0215 through bb0218 is cotranscribed and is therefore an operon. Sequence analysis of these proteins suggested that they are the components of an ABC‐type transporter responsible for translocating phosphate anions from the B. burgdorferi periplasm to the cytoplasm. Biophysical experiments established that the putative ligand‐binding protein of this system, BbPstS (BB0215), binds to phosphate in solution. We determined the high‐resolution (1.3 Å) crystal structure of the protein in the absence of phosphate, revealing that the protein's fold is similar to other phosphate‐binding proteins, and residues that are implicated in phosphate binding in other such proteins are conserved in BbPstS. Taken together, the gene products of bb0215‐0218 function as a phosphate transporter for B. burgdorferi.     

Dynamics in an unusual acyl carrier protein from a ladderane lipid-synthesizing organism

Anaerobic ammonium-oxidizing (anammox) bacteria express a distinct acyl carrier protein implicated in the biosynthesis of the highly unusual “ladderane” lipids these organisms produce. This “anammox-specific” ACP, or amxACP, shows several unique features such as a conserved FF motif and an unusual sequence in the functionally important helix III. Investigation of the protein's structure and dynamics, both in the crystal by ensemble refinement and by MD simulations, reveals that helix III adopts a rare six-residue-long 3₁₀-helical conformation that confers a large degree of conformational and positional variability on this part of the protein. This way of introducing structural flexibility by using the inherent properties of 3₁₀-helices appears unique among ACPs. Moreover, the structure suggests a role for the FF motif in shielding the thioester linkage between the protein's prosthetic group and its acyl cargo from hydrolysis.