A protein coevolution method uncovers critical features of the Hepatitis C Virus fusion mechanism

A sequence-based evolutionary distance method for Phylogenetic analysis of highly divergent proteins

Because of the limited effectiveness of prevailing phylogenetic methods when applied to highly divergent protein sequences, the phylogenetic analysis problem remains challenging. Here, we propose a sequence-based evolutionary distance algorithm termed sequence distance (SD), which innovatively incorporates site-to-site correlation within protein sequences into the distance estimation. In protein superfamilies, SD can effectively distinguish evolutionary relationships both within and between protein families, producing phylogenetic trees that closely align with those based on structural information, even with sequence identity less than 20%. SD is highly correlated with the similarity of the protein structure, and can calculate evolutionary distances for thousands of protein pairs within seconds using a single CPU, which is significantly faster than most protein structure prediction methods that demand high computational resources and long run times. The development of SD will significantly advance phylogenetics, providing researchers with a more accurate and reliable tool for exploring evolutionary relationships.

Insights into Hepatitis C Virus E2core Interactions with Human Cellular Receptor CD81 at Different pHs from Molecular Simulations

Hepatitis C virus (HCV) is the second viral agent that causes the majority of chronic hepatic infections worldwide, following Hepatitis B virus (HBV) infection. HCV infection comprises several steps, from the attachment to the receptors to the delivery of the viral genetic material and replication inside the cells. Tetraspanin CD81 is a key entry factor for HCV as it accompanies the virus during attachment and internalization through clathrin-mediated endocytosis. HCV-CD81 binding takes place through the viral glycoprotein E2. We performed full-atom molecular dynamics simulations reproducing the pH conditions that occur during the viral attachment to the hepatocytes (pH 7.4) and internalization (pH 6.2-4.6). We observed that changing the pH from 7.4 to 6.2 triggers a large conformational change in the binding orientation between E2_core (E2_core corresponds to residues 412-645 of the viral glycoprotein E2) and CD81_LEL (CD81_LEL corresponds to residues 112-204 of CD81) that occurs even more rapidly at low pH 4.6. This pH-induced switching mechanism has never been observed before and could allow the virus particles to sense the right moment during the maturation of the endosome to start fusion.

Can next-generation humanized mice that reconstituted with both functional human immune system and hepatocytes model the progression of viral hepatitis to hepatocarcinogenesis?

Hepatitis B virus (HBV) and Hepatitis C virus (HCV) chronic infections cause liver immunopathological diseases such as hepatitis, fibrosis, cirrhosis, and hepatocellular carcinomas, which are difficult to treat and continue to be major health problems globally. Due to the species-specific hepato-tropism of HBV and HCV, conventional rodent models are limited in their utility for studying the infection and associated liver immunopathogenesis. Humanized mice reconstituted with both functional human immune system and hepatocytes (HIS-HuHEP mice) have been extremely instrumental for in vivo studies of HBV or HCV infection and human-specific aspects of the progression of liver immunopathogenesis. However, none of the current HIS-HuHEP mice can model the progression of viral hepatitis to hepatocarcinogenesis which may be a notorious result of HBV or HCV chronic infection in patients, suggesting that they were functionally compromised and that there is still significant space to improve and establish next-generation of HIS-HuHEP mice with more sophisticated functions. In this review, we first summarize the principal requirements to establish HIS-HuHEP mice. We then discuss the respective protocols for current HIS-HuHEP mice and their applications, as well as their advantages and disadvantages. We also raise perspectives for further improving and establishing next-generation HIS-HuHEP mice.

Web-Based Protein Interactions Calculator Identifies Likely Proteome Coevolution with Alzheimer’s Disease-Associated Proteins

Protein–protein functional interactions arise from either transitory or permanent biomolecular associations and often lead to the coevolution of the interacting residues. Although mutual information has traditionally been used to identify coevolving residues within the same protein, its application between coevolving proteins remains largely uncharacterized. Therefore, we developed the Protein Interactions Calculator (PIC) to efficiently identify coevolving residues between two protein sequences using mutual information. We verified the algorithm using 2102 known human protein interactions and 233 known bacterial protein interactions, with a respective 1975 and 252 non-interacting protein controls. The average PIC score for known human protein interactions was 4.5 times higher than non-interacting proteins (p = 1.03 × 10⁻¹⁰⁸) and 1.94 times higher in bacteria (p = 1.22 × 10⁻³⁵). We then used the PIC scores to determine the probability that two proteins interact. Using those probabilities, we paired 37 Alzheimer’s disease-associated proteins with 8608 other proteins and determined the likelihood that each pair interacts, which we report through a web interface. The PIC had significantly higher sensitivity and residue-specific resolution not available in other algorithms. Therefore, we propose that the PIC can be used to prioritize potential protein interactions, which can lead to a better understanding of biological processes and additional therapeutic targets belonging to protein interaction groups.

iBIS2Analyzer: a web server for a phylogeny-driven coevolution analysis of protein families

Residue coevolution within and between proteins is used as a marker of physical interaction and/or residue functional cooperation. Pairs or groups of coevolving residues are extracted from multiple sequence alignments based on a variety of computational approaches. However, coevolution signals emerging in subsets of sequences might be lost if the full alignment is considered. iBIS2Analyzer is a web server dedicated to a phylogeny-driven coevolution analysis of protein families with different evolutionary pressure. It is based on the iterative version, iBIS², of the coevolution analysis method BIS, Blocks in Sequences. iBIS² is designed to iteratively select and analyse subtrees in phylogenetic trees, possibly large and comprising thousands of sequences. With iBIS2Analyzer, openly accessible at http://ibis2analyzer.lcqb.upmc.fr/, the user visualizes, compares and inspects clusters of coevolving residues by mapping them onto sequences, alignments or structures of choice, greatly simplifying downstream analysis steps. A rich and interactive graphic interface facilitates the biological interpretation of the results.

Secretory glycoprotein NS1 plays a crucial role in the particle formation of flaviviruses

Flaviviruses, which are globally distributed and cause a spectrum of potentially severe illnesses, pose a major threat to public health. Although Flaviviridae viruses, including flaviviruses, possess similar genome structures, only the flaviviruses encode the non-structural protein NS1, which resides in the endoplasmic reticulum (ER) and is secreted from cells after oligomerization. The ER-resident NS1 is known to be involved in viral genome replication, but the essential roles of secretory NS1 in the virus life cycle are not fully understood. Here we characterized the roles of secretory NS1 in the particle formation of flaviviruses. We first identified an amino acid residue essential for the NS1 secretion but not for viral genome replication by using protein-protein interaction network analyses and mutagenesis scanning. By using the recombinant flaviviruses carrying the identified NS1 mutation, we clarified that the mutant flaviviruses employed viral genome replication. We then constructed a recombinant NS1 with the identified mutation and demonstrated by physicochemical assays that the mutant NS1 was unable to form a proper oligomer or associate with liposomes. Finally, we showed that the functions of NS1 that were lost by the identified mutation could be compensated for by the in trans-expression of E^rns of pestiviruses and host exchangeable apolipoproteins, which participate in the infectious particle formation of pestiviruses and hepaciviruses in the family Flaviviridae, respectively. Collectively, our study suggests that secretory NS1 plays a role in the particle formation of flaviviruses through its interaction with the lipid membrane.

It is difficult to characterize the function of NS1 in the post-genome replication stages in the virus life cycle of flaviviruses. Here, by means of protein-protein interaction network analyses and mutagenesis scanning, we identified a unique mutation in NS1 by which the protein loses its secretory capacity while retaining its genome replication activity. Physicochemical assays using the mutant NS1 revealed that oligomerization of NS1 is responsible for the lipid association and secretion of NS1. In addition, we established a complementation assay that can evaluate the particle formation of Flaviviridae viruses. By using recombinant flaviviruses possessing the identified mutation in NS1, we clarified that NS1 is involved in particle formation. Our findings reveal that the flavivirus NS1 has at least two roles in the virus life cycles—namely, a role in infectious particle formation and a role in viral genome replication.

Sites of vulnerability in HCV E1E2 identified by comprehensive functional screening

The E1 and E2 envelope proteins of hepatitis C virus (HCV) form a heterodimer that drives virus-host membrane fusion. Here, we analyze the role of each amino acid in E1E2 function, expressing 545 individual alanine mutants of E1E2 in human cells, incorporating them into infectious viral pseudoparticles, and testing them against 37 different monoclonal antibodies (MAbs) to ascertain full-length translation, folding, heterodimer assembly, CD81 binding, viral pseudoparticle incorporation, and infectivity. We propose a model describing the role of each critical residue in E1E2 functionality and use it to examine how MAbs neutralize infection by exploiting functionally critical sites of vulnerability on E1E2. Our results suggest that E1E2 is a surprisingly fragile protein complex where even a single alanine mutation at 92% of positions disrupts its function. The amino-acid-level targets identified are highly conserved and functionally critical and can be exploited for improved therapies and vaccines.

Viral Causality of Human Cancer and Potential Roles of Human Endogenous Retroviruses in the Multi-Omics Era: An Evolutionary Epidemiology Review

Being responsible for almost 12% of cancers worldwide, viruses are among the oldest known and most prevalent oncogenic agents. The quality of the evidence for the in vivo tumorigenic potential of microorganisms varies, thus accordingly, viruses were classified in 4 evidence-based categories by the International Agency for Research on Cancer in 2009. Since then, our understanding of the role of viruses in cancer has significantly improved, firstly due to the emergence of high throughput sequencing technologies that allowed the “brute-force” recovery of unknown viral genomes. At the same time, multi-omics approaches unravelled novel virus-host interactions in stem-cell biology. We now know that viral elements, either exogenous or endogenous, have multiple sometimes conflicting roles in human pathophysiology and the development of cancer. Here we integrate emerging evidence on viral causality in human cancer from basic mechanisms to clinical studies. We analyze viral tumorigenesis under the scope of deep-in-time human-virus evolutionary relationships and critically comment on the evidence through the eyes of clinical epidemiology, firstly by reviewing recognized oncoviruses and their mechanisms of inducing tumorigenesis, and then by examining the potential role of integrated viruses in our genome in the process of carcinogenesis.

Identification and Characteristics of Fusion Peptides Derived From Enveloped Viruses

Membrane fusion events allow enveloped viruses to enter and infect cells. The study of these processes has led to the identification of a number of proteins that mediate this process. These proteins are classified according to their structure, which vary according to the viral genealogy. To date, three classes of fusion proteins have been defined, but current evidence points to the existence of additional classes. Despite their structural differences, viral fusion processes follow a common mechanism through which they exert their actions. Additional studies of the viral fusion proteins have demonstrated the key role of specific proteinogenic subsequences within these proteins, termed fusion peptides. Such peptides are able to interact and insert into membranes for which they hold interest from a pharmacological or therapeutic viewpoint. Here, the different characteristics of fusion peptides derived from viral fusion proteins are described. These criteria are useful to identify new fusion peptides. Moreover, this review describes the requirements of synthetic fusion peptides derived from fusion proteins to induce fusion by themselves. Several sequences of the viral glycoproteins E1 and E2 of HCV were, for example, identified to be able to induce fusion, which are reviewed here.

COVTree: Coevolution in OVerlapped sequences by Tree analysis server

Overlapping genes are commonplace in viruses and play an important role in their function and evolution. For these genes, molecular coevolution may be seen as a mechanism to decrease the evolutionary constraints of amino acid positions in the overlapping regions and to tolerate or compensate unfavorable mutations. Tracing these mutational sites, could help to gain insight on the direct or indirect effect of the mutations in the corresponding overlapping proteins. In the past, coevolution analysis has been used to identify residue pairs and coevolutionary signatures within or between proteins that served as markers of physical interactions and/or functional relationships. Coevolution in OVerlapped sequences by Tree analysis (COVTree) is a web server providing the online analysis of coevolving amino-acid pairs in overlapping genes, where residues might be located inside or outside the overlapping region. COVTree is designed to handle protein families with various characteristics, among which those that typically display a small number of highly conserved sequences. It is based on BIS2, a fast version of the coevolution analysis tool Blocks in Sequences (BIS). COVTree provides a rich and interactive graphical interface to ease biological interpretation of the results and it is openly accessible at http://www.lcqb.upmc.fr/COVTree/.

Chlorcyclizine Inhibits Viral Fusion of Hepatitis C Virus Entry by Directly Targeting HCV Envelope Glycoprotein 1

Chlorcyclizine (CCZ) is a potent hepatitis C virus (HCV) entry inhibitor, but its molecular mechanism is unknown. Here, we show that CCZ directly targets the fusion peptide of HCV E1 and interferes with the fusion process. Generation of CCZ resistance-associated substitutions of HCV in vitro revealed six missense mutations in the HCV E1 protein, five being in the putative fusion peptide. A viral fusion assay demonstrated that CCZ blocked HCV entry at the membrane fusion step and that the mutant viruses acquired resistance to CCZ's action in blocking membrane fusion. UV cross-linking of photoactivatable CCZ-diazirine-biotin in both HCV-infected cells and recombinant HCV E1/E2 protein demonstrated direct binding to HCV E1 glycoprotein. Mass spectrometry analysis revealed that CCZ cross-linked to an E1 sequence adjacent to the putative fusion peptide. Docking simulations demonstrate a putative binding model, wherein CCZ binds to a hydrophobic pocket of HCV E1 and forms extensive interactions with the fusion peptide.

CoevDB: a database of intramolecular coevolution among protein-coding genes of the bony vertebrates

The study of molecular coevolution, due to its potential to identify gene regions under functional or structural constraints, has recently been subject to numerous scientific inquiries. Particular efforts have been conducted to develop methods predicting the presence of coevolution in molecular sequences. Among these methods, a few aim to model the underlying evolutionary process of coevolution, which enable to differentiate the shared history of genes to coevolution and thus improve their accuracy. However, the usage of such methods remains sparse due to their expensive computational cost and the lack of resources alleviating this issue. Here we present CoevDB (http://phylodb.unil.ch/CoevDB), a database containing the result of a large-scale analysis of intramolecular coevolution of 8201 protein-coding genes of bony vertebrates. The web interface of CoevDB gives access to the results to 800 millions of statistical tests corresponding to all the pairs of sites analyzed. Several type of queries enable users to explore the database by either targeting specific genes or by discovering genes having promising estimations of coevolution.

Hepatitis C Virus Entry: An Intriguingly Complex and Highly Regulated Process

Hepatitis C virus (HCV) is a major cause of chronic hepatitis and liver disease worldwide. Its tissue and species tropism are largely defined by the viral entry process that is required for subsequent productive viral infection and establishment of chronic infection. This review provides an overview of the viral and host factors involved in HCV entry into hepatocytes, summarizes our understanding of the molecular mechanisms governing this process and highlights the therapeutic potential of host-targeting entry inhibitors.

Flexibility and intrinsic disorder are conserved features of hepatitis C virus E2 glycoprotein

The glycoproteins of hepatitis C virus, E1E2, are unlike any other viral fusion machinery yet described, and are the current focus of immunogen design in HCV vaccine development; thus, making E1E2 both scientifically and medically important. We used pre-existing, but fragmentary, structures to model a complete ectodomain of the major glycoprotein E2 from three strains of HCV. We then performed molecular dynamic simulations to explore the conformational landscape of E2, revealing a number of important features. Despite high sequence divergence, and subtle differences in the models, E2 from different strains behave similarly, possessing a stable core flanked by highly flexible regions, some of which perform essential functions such as receptor binding. Comparison with sequence data suggest that this consistent behaviour is conferred by a network of conserved residues that act as hinge and anchor points throughout E2. The variable regions (HVR-1, HVR-2 and VR-3) exhibit particularly high flexibility, and bioinformatic analysis suggests that HVR-1 is a putative intrinsically disordered protein region. Dynamic cross-correlation analyses demonstrate intramolecular communication and suggest that specific regions, such as HVR-1, can exert influence throughout E2. To support our computational approach we performed small-angle X-ray scattering with purified E2 ectodomain; this data was consistent with our MD experiments, suggesting a compact globular core with peripheral flexible regions. This work captures the dynamic behaviour of E2 and has direct relevance to the interaction of HCV with cell-surface receptors and neutralising antibodies.

Disease-relevant mutations alter amino acid co-evolution networks in the second nucleotide binding domain of CFTR

Cystic Fibrosis (CF) is an inherited disease caused by mutations in the cystic fibrosis transmembrane conductance regulator (CFTR) ion channel. Mutations in CFTR cause impaired chloride ion transport in the epithelial tissues of patients leading to cardiopulmonary decline and pancreatic insufficiency in the most severely affected patients. CFTR is composed of twelve membrane-spanning domains, two nucleotide-binding domains (NBDs), and a regulatory domain. The most common mutation in CFTR is a deletion of phenylalanine at position 508 (ΔF508) in NBD1. Previous research has primarily concentrated on the structure and dynamics of the NBD1 domain; However numerous pathological mutations have also been found in the lesser-studied NBD2 domain. We have investigated the amino acid co-evolved network of interactions in NBD2, and the changes that occur in that network upon the introduction of CF and CF-related mutations (S1251N(T), S1235R, D1270N, N1303K(T)). Extensive coupling between the α- and β-subdomains were identified with residues in, or near Walker A, Walker B, H-loop and C-loop motifs. Alterations in the predicted residue network varied from moderate for the S1251T perturbation to more severe for N1303T. The S1235R and D1270N networks varied greatly compared to the wildtype, but these CF mutations only affect ion transport preference and do not severely disrupt CFTR function, suggesting dynamic flexibility in the network of interactions in NBD2. Our results also suggest that inappropriate interactions between the β-subdomain and Q-loop could be detrimental. We also identified mutations predicted to stabilize the NBD2 residue network upon introduction of the CF and CF-related mutations, and these predicted mutations are scored as benign by the MUTPRED2 algorithm. Our results suggest the level of disruption of the co-evolution predictions of the amino acid networks in NBD2 does not have a straightforward correlation with the severity of the CF phenotypes observed.

Coevolution analysis of amino-acids reveals diversified drug-resistance solutions in viral sequences: a case study of hepatitis B virus

The study of mutational landscapes of viral proteins is fundamental for the understanding of the mechanisms of cross-resistance to drugs and the design of effective therapeutic strategies based on several drugs. Antiviral therapy with nucleos(t)ide analogues targeting the hepatitis B virus (HBV) polymerase protein (Pol) can inhibit disease progression by suppression of HBV replication and makes it an important case study. In HBV, treatment may fail due to the emergence of drug-resistant mutants. Primary and compensatory mutations have been associated with lamivudine resistance, whereas more complex mutational patterns are responsible for resistance to other HBV antiviral drugs. So far, all known drug-resistance mutations are located in one of the four Pol domains, called reverse transcriptase. We demonstrate that sequence covariation identifies drug-resistance mutations in viral sequences. A new algorithmic strategy, BIS2TreeAnalyzer, is designed to apply the coevolution analysis method BIS2, successfully used in the past on small sets of conserved sequences, to large sets of evolutionary related sequences. When applied to HBV, BIS2TreeAnalyzer highlights diversified viral solutions by discovering thirty-seven positions coevolving with residues known to be associated with drug resistance and located on the four Pol domains. These results suggest a sequential mechanism of emergence for some mutational patterns. They reveal complex combinations of positions involved in HBV drug resistance and contribute with new information to the landscape of HBV evolutionary solutions. The computational approach is general and can be applied to other viral sequences when compensatory mutations are presumed.

Functional expression and characterization of the envelope glycoprotein E1E2 heterodimer of hepatitis C virus

Hepatitis C virus (HCV) is a member of Hepacivirus and belongs to the family of Flaviviridae. HCV infects millions of people worldwide and may lead to cirrhosis and hepatocellular carcinoma. HCV envelope proteins, E1 and E2, play critical roles in viral cell entry and act as major epitopes for neutralizing antibodies. However, unlike other known flaviviruses, it has been challenging to study HCV envelope proteins E1E2 in the past decades as the in vitro expressed E1E2 heterodimers are usually of poor quality, making the structural and functional characterization difficult. Here we express the ectodomains of HCV E1E2 heterodimer with either an Fc-tag or a de novo designed heterodimeric tag and are able to isolate soluble E1E2 heterodimer suitable for functional and structural studies. Then we characterize the E1E2 heterodimer by electron microscopy and model the structure by the coevolution based modeling strategy with Rosetta, revealing the potential interactions between E1 and E2. Moreover, the E1E2 heterodimer is applied to examine the interactions with the known HCV receptors, neutralizing antibodies as well as the inhibition of HCV infection, confirming the functionality of the E1E2 heterodimer and the binding profiles of E1E2 with the cellular receptors. Therefore, the expressed E1E2 heterodimer would be a valuable target for both viral studies and vaccination against HCV.

Hepatitis C virus (HCV) is an enveloped virus that infects millions of people worldwide and may lead to cirrhosis and hepatocellular carcinoma. HCV has two envelope proteins, E1 and E2, which form heterodimers on viral surface and are critical for HCV cell entry. However, current studies of HCV E1E2 are often limited by the poor quality of the in vitro expressed E1E2 heterodimers. Here we express the ectodomains of HCV E1E2 with different tags, and are able to isolate soluble E1E2 ectodomains suitable for structural and functional studies. Then we generate the 3D reconstruction of E1E2 heterodimer by electron microscopy and also model the E1E2 structure by the coevolution based strategy with Rosetta, showing the potential interactions between E1 and E2. Moreover, the E1E2 heterodimer is applied to examine the interactions with the HCV cellular receptors, neutralizing antibodies as well as the inhibition of HCV infection. These results suggest that the expressed E1E2 heterodimer would be a promising target for both viral studies and vaccination against HCV.

Function of the HCV E1 envelope glycoprotein in viral entry and assembly

HCV envelope glycoproteins, E1 and E2, are multifunctional proteins. Until recently, E2 glycoprotein was thought to be the fusion protein and was the focus of investigations. However, the recently obtained partial structures of E2 and E1 rather support a role for E1 alone or in association with E2 in HCV fusion. Moreover, they suggest that HCV harbors a new fusion mechanism, distinct from that of other members of the Flaviviridae family. In this context, E1 aroused a renewed interest. Recent functional characterizations of E1 revealed a more important role than previously thought in entry and assembly. Thus, E1 is involved in the viral genome encapsidation step and influences the association of the virus with lipoprotein components. Moreover, E1 modulates HCV–receptor interaction and participates in a late entry step potentially fusion. In this review, we outline our current knowledge on E1 functions in HCV assembly and entry.

Simultaneous Bayesian inference of phylogeny and molecular coevolution

Patterns of molecular coevolution can reveal structural and functional constraints within or among organic molecules. These patterns are better understood when considering the underlying evolutionary process, which enables us to disentangle the signal of the dependent evolution of sites (coevolution) from the effects of shared ancestry of genes. Conversely, disregarding the dependent evolution of sites when studying the history of genes negatively impacts the accuracy of the inferred phylogenetic trees. Although molecular coevolution and phylogenetic history are interdependent, analyses of the two processes are conducted separately, a choice dictated by computational convenience, but at the expense of accuracy. We present a Bayesian method and associated software to infer how many and which sites of an alignment evolve according to an independent or a pairwise dependent evolutionary process, and to simultaneously estimate the phylogenetic relationships among sequences. We validate our method on synthetic datasets and challenge our predictions of coevolution on the 16S rRNA molecule by comparing them with its known molecular structure. Finally, we assess the accuracy of phylogenetic trees inferred under the assumption of independence among sites using synthetic datasets, the 16S rRNA molecule and 10 additional alignments of protein-coding genes of eukaryotes. Our results demonstrate that inferring phylogenetic trees while accounting for dependent site evolution significantly impacts the estimates of the phylogeny and the evolutionary process.

Role of Hepatitis C Virus Envelope Glycoprotein E1 in Virus Entry and Assembly

Hepatitis C virus (HCV) glycoproteins E1 and E2 form a heterodimer to constitute viral envelope proteins, which play an essential role in virus entry. E1 does not directly interact with host receptors, and its functions in viral entry are exerted mostly through its interaction with E2 that directly binds the receptors. HCV enters the host cell via receptor-mediated endocytosis during which the fusion of viral and host endosomal membranes occurs to release viral genome to cytoplasm. A putative fusion peptide in E1 has been proposed to participate in membrane fusion, but its exact role and underlying molecular mechanisms remain to be deciphered. Recently solved crystal structures of the E2 ectodomains and N-terminal of E1 fail to reveal a classical fusion-like structure in HCV envelope glycoproteins. In addition, accumulating evidence suggests that E1 also plays an important role in virus assembly. In this mini-review, we summarize current knowledge on HCV E1 including its structure and biological functions in virus entry, fusion, and assembly, which may provide clues for developing HCV vaccines and more effective antivirals.