My Cousin, My Enemy: quasispecies suppression of drug resistance

Molecular mechanism of rhinovirus escape from the Pyrazolo[3,4-d]pyrimidine capsid-binding inhibitor OBR-5-340 via mutations distant from the binding pocket: Derivatives that brake resistance

Rhinoviruses (RVs) cause the common cold. Attempts at discovering small molecule inhibitors have mainly concentrated on compounds supplanting the medium chain fatty acids residing in the sixty icosahedral symmetry-related hydrophobic pockets of the viral capsid of the Rhinovirus-A and -B species. High-affinity binding to these pockets stabilizes the capsid against structural changes necessary for the release of the ss(+) RNA genome into the cytosol of the host cell. However, single-point mutations may abolish this binding. RV-B5 is one of several RVs that are naturally resistant against the well-established antiviral agent pleconaril. However, RV-B5 is strongly inhibited by the pyrazolopyrimidine OBR-5-340. Here, we report on isolation and characterization of RV-B5 mutants escaping OBR-5-340 inhibition and show that substitution of amino acid residues not only within the binding pocket but also remote from the binding pocket hamper inhibition. Molecular dynamics network analysis revealed that strong inhibition occurs when an ensemble of several sequence stretches of the capsid proteins enveloping OBR-5-340 move together with OBR-5-340. Mutations abrogating this dynamic, regardless of whether being localized within the binding pocket or distant from it result in escape from inhibition. Pyrazolo [3,4-d]pyrimidine derivatives overcoming OBR-5-340 escape of various RV-B5 mutants were identified. Our work contributes to the understanding of the properties of capsid-binding inhibitors necessary for potent and broad-spectrum inhibition of RVs.

Differential inhibition of intra- and inter-molecular protease cleavages by antiviral compounds

IMPORTANCE

Most protease-targeted antiviral development evaluates the ability of small molecules to inhibit the cleavage of artificial substrates. However, before they can cleave any other substrates, viral proteases need to cleave themselves out of the viral polyprotein in which they have been translated. This can occur either intra- or inter-molecularly. Whether this process occurs intra- or inter-molecularly has implications for the potential for precursors to accumulate and for the effectiveness of antiviral drugs. We argue that evaluating candidate antivirals for their ability to block these cleavages is vital to drug development because the buildup of uncleaved precursors can be inhibitory to the virus and potentially suppress the selection of drug-resistant variants.

Puzzles, challenges, and information reservoir of SARS-CoV-2 quasispecies

Upon the emergence of SARS-CoV-2 in the human population, it was conjectured that for this coronavirus the dynamic intra-host heterogeneity typical of RNA viruses would be toned down. Nothing of this sort is observed. Here we review the main observations on the complexity and diverse composition of SARS-CoV-2 mutant spectra sampled from infected patients, within the framework of quasispecies dynamics. The analyses suggest that the information provided by myriads of genomic sequences within infected individuals may have a predictive value of the genomic sequences that acquire epidemiological relevance. Possibilities to reconcile the presence of broad mutant spectra in the large RNA coronavirus genome with its encoding a 3' to 5' exonuclease proofreading-repair activity are considered. Indeterminations in the behavior of individual viral genomes provide a benefit for the survival of the ensemble. We propose that this concept falls in the domain of "stochastic thinking," a notion that applies also to cellular processes, as a means for biological systems to face unexpected needs.

SPLASH: a statistical, reference-free genomic algorithm unifies biological discovery

Today's genomics workflows typically require alignment to a reference sequence, which limits discovery. We introduce a new unifying paradigm, SPLASH (Statistically Primary aLignment Agnostic Sequence Homing), an approach that directly analyzes raw sequencing data to detect a signature of regulation: sample-specific sequence variation. The approach, which includes a new statistical test, is computationally efficient and can be run at scale. SPLASH unifies detection of myriad forms of sequence variation. We demonstrate that SPLASH identifies complex mutation patterns in SARS-CoV-2 strains, discovers regulated RNA isoforms at the single cell level, documents the vast sequence diversity of adaptive immune receptors, and uncovers biology in non-model organisms undocumented in their reference genomes: geographic and seasonal variation and diatom association in eelgrass, an oceanic plant impacted by climate change, and tissue-specific transcripts in octopus. SPLASH is a new unifying approach to genomic analysis that enables an expansive scope of discovery without metadata or references.

Viral Fitness, Population Complexity, Host Interactions, and Resistance to Antiviral Agents

Fitness of viruses has become a standard parameter to quantify their adaptation to a biological environment. Fitness determinations for RNA viruses (and some highly variable DNA viruses) meet with several uncertainties. Of particular interest are those that arise from mutant spectrum complexity, absence of population equilibrium, and internal interactions among components of a mutant spectrum. Here, concepts, fitness measurements, limitations, and current views on experimental viral fitness landscapes are discussed. The effect of viral fitness on resistance to antiviral agents is covered in some detail since it constitutes a widespread problem in antiviral pharmacology, and a challenge for the design of effective antiviral treatments. Recent evidence with hepatitis C virus suggests the operation of mechanisms of antiviral resistance additional to the standard selection of drug-escape mutants. The possibility that high replicative fitness may be the driver of such alternative mechanisms is considered. New broad-spectrum antiviral designs that target viral fitness may curtail the impact of drug-escape mutants in treatment failures. We consider to what extent fitness-related concepts apply to coronaviruses and how they may affect strategies for COVID-19 prevention and treatment.

SARS-CoV-2 Mutant Spectra at Different Depth Levels Reveal an Overwhelming Abundance of Low Frequency Mutations

Populations of RNA viruses are composed of complex and dynamic mixtures of variant genomes that are termed mutant spectra or mutant clouds. This applies also to SARS-CoV-2, and mutations that are detected at low frequency in an infected individual can be dominant (represented in the consensus sequence) in subsequent variants of interest or variants of concern. Here we briefly review the main conclusions of our work on mutant spectrum characterization of hepatitis C virus (HCV) and SARS-CoV-2 at the nucleotide and amino acid levels and address the following two new questions derived from previous results: (i) how is the SARS-CoV-2 mutant and deletion spectrum composition in diagnostic samples, when examined at progressively lower cut-off mutant frequency values in ultra-deep sequencing; (ii) how the frequency distribution of minority amino acid substitutions in SARS-CoV-2 compares with that of HCV sampled also from infected patients. The main conclusions are the following: (i) the number of different mutations found at low frequency in SARS-CoV-2 mutant spectra increases dramatically (50- to 100-fold) as the cut-off frequency for mutation detection is lowered from 0.5% to 0.1%, and (ii) that, contrary to HCV, SARS-CoV-2 mutant spectra exhibit a deficit of intermediate frequency amino acid substitutions. The possible origin and implications of mutant spectrum differences among RNA viruses are discussed.

A Targeted Computational Screen of the SWEETLEAD Database Reveals FDA-Approved Compounds with Anti-Dengue Viral Activity

Affordable and effective antiviral therapies are needed worldwide, especially against agents such as dengue virus that are endemic in underserved regions. Many antiviral compounds have been studied in cultured cells but are unsuitable for clinical applications due to pharmacokinetic profiles, side effects, or inconsistent efficacy across dengue serotypes. Such tool compounds can, however, aid in identifying clinically useful treatments. Here, computational screening (Rapid Overlay of Chemical Structures) was used to identify entries in an in silico database of safe-in-human compounds (SWEETLEAD) that display high chemical similarities to known inhibitors of dengue virus. Inhibitors of the dengue proteinase NS2B/3, the dengue capsid, and the host autophagy pathway were used as query compounds. Three FDA-approved compounds that resemble the tool molecules structurally, cause little toxicity, and display strong antiviral activity in cultured cells were selected for further analysis. Pyrimethamine (50% inhibitory concentration [IC₅₀] = 1.2 μM), like the dengue proteinase inhibitor ARDP0006 to which it shows structural similarity, inhibited intramolecular NS2B/3 cleavage. Lack of toxicity early in infection allowed testing in mice, in which pyrimethamine also reduced viral loads. Niclosamide (IC₅₀ = 0.28 μM), like dengue core inhibitor ST-148, affected structural components of the virion and inhibited early processes during infection. Vandetanib (IC₅₀ = 1.6 μM), like cellular autophagy inhibitor spautin-1, blocked viral exit from cells and could be shown to extend survival in vivo Thus, three FDA-approved compounds with promising utility for repurposing to treat dengue virus infections and their potential mechanisms were identified using computational tools and minimal phenotypic screening.IMPORTANCE No antiviral therapeutics are currently available for dengue virus infections. By computationally overlaying the three-dimensional (3D) chemical structures of compounds known to inhibit dengue virus over those of compounds known to be safe in humans, we identified three FDA-approved compounds that are attractive candidates for repurposing as antivirals. We identified targets for two previously identified antiviral compounds and revealed a previously unknown potential anti-dengue compound, vandetanib. This computational approach to analyze a highly curated library of structures has the benefits of speed and cost efficiency. It also leverages mechanistic work with query compounds used in biomedical research to provide strong hypotheses for the antiviral mechanisms of the safer hit compounds. This workflow to identify compounds with known safety profiles can be expanded to any biological activity for which a small-molecule query compound has been identified, potentially expediting the translation of basic research to clinical interventions.

Full-length three-dimensional structure of the influenza A virus M1 protein and its organization into a matrix layer

Matrix proteins are encoded by many enveloped viruses, including influenza viruses, herpes viruses, and coronaviruses. Underneath the viral envelope of influenza virus, matrix protein 1 (M1) forms an oligomeric layer critical for particle stability and pH-dependent RNA genome release. However, high-resolution structures of full-length monomeric M1 and the matrix layer have not been available, impeding antiviral targeting and understanding of the pH-dependent transitions involved in cell entry. Here, purification and extensive mutagenesis revealed protein–protein interfaces required for the formation of multilayered helical M1 oligomers similar to those observed in virions exposed to the low pH of cell entry. However, single-layered helical oligomers with biochemical and ultrastructural similarity to those found in infectious virions before cell entry were observed upon mutation of a single amino acid. The highly ordered structure of the single-layered oligomers and their likeness to the matrix layer of intact virions prompted structural analysis by cryo-electron microscopy (cryo-EM). The resulting 3.4-Å–resolution structure revealed the molecular details of M1 folding and its organization within the single-shelled matrix. The solution of the full-length M1 structure, the identification of critical assembly interfaces, and the development of M1 assembly assays with purified proteins are crucial advances for antiviral targeting of influenza viruses.

Multi-subunit shells of matrix proteins line the interior of infectious influenza virus particles. In this study, biochemical purification of wild-type and mutant influenza M1 proteins allows the structural determination of an oligomer whose shape corresponds to that of infectious virions and suggests mechanisms for its formation and dismantling during infection.

Viral fitness: history and relevance for viral pathogenesis and antiviral interventions

The quasispecies dynamics of viral populations (continuous generation of variant genomes and competition among them) has as one of its frequent consequences variations in overall multiplication capacity, a major component of viral fitness. This parameter has multiple implications for viral pathogenesis and viral disease control, some of them unveiled thanks to deep sequencing of viral populations. Darwinian fitness is an old concept whose quantification dates back to the early developments of population genetics. It was later applied to viruses (mainly to RNA viruses) to quantify relative multiplication capacities of individual mutant clones or complex populations. The present article reviews the fitness concept and its relevance for the understanding of the adaptive dynamics of viruses in constant and changing environments. Many studies have addressed the fitness cost of escape mutations (to antibodies, cytotoxic T cells or inhibitors) as an influence on the efficacy of antiviral interventions. Here, we summarize the evidence that the basal fitness level can be a determinant of inhibitor resistance.

Quasispecies dynamics in disease prevention and control

Medical interventions to prevent and treat viral disease constitute evolutionary forces that may modify the genetic composition of viral populations that replicate in an infected host and influence the genomic composition of those viruses that are transmitted and progress at the epidemiological level. Given the adaptive potential of viruses in general and the RNA viruses in particular, the selection of viral mutants that display some degree of resistance to inhibitors or vaccines is a tangible challenge. Mutant selection may jeopardize control of the viral disease. Strategies intended to minimize vaccination and treatment failures are proposed and justified based on fundamental features of viral dynamics explained in the preceding chapters. The recommended use of complex, multiepitopic vaccines, and combination therapies as early as possible after initiation of infection falls under the general concept that complexity cannot be combated with simplicity. It also follows that sociopolitical action to interrupt virus replication and spread as soon as possible is as important as scientifically sound treatment designs to control viral disease on a global scale.

Viral quasispecies

Viral quasispecies refers to a population structure that consists of extremely large numbers of variant genomes, termed mutant spectra, mutant swarms or mutant clouds. Fueled by high mutation rates, mutants arise continually, and they change in relative frequency as viral replication proceeds. The term quasispecies was adopted from a theory of the origin of life in which primitive replicons) consisted of mutant distributions, as found experimentally with present day RNA viruses. The theory provided a new definition of wild type, and a conceptual framework for the interpretation of the adaptive potential of RNA viruses that contrasted with classical studies based on consensus sequences. Standard clonal analyses and deep sequencing methodologies have confirmed the presence of myriads of mutant genomes in viral populations, and their participation in adaptive processes. The quasispecies concept applies to any biological entity, but its impact is more evident when the genome size is limited and the mutation rate is high. This is the case of the RNA viruses, ubiquitous in our biosphere, and that comprise many important pathogens. In virology, quasispecies are defined as complex distributions of closely related variant genomes subjected to genetic variation, competition and selection, and that may act as a unit of selection. Despite being an integral part of their replication, high mutation rates have an upper limit compatible with inheritable information. Crossing such a limit leads to RNA virus extinction, a transition that is the basis of an antiviral design termed lethal mutagenesis.

Enterovirus 71 targets the cardiopulmonary system in a robust oral infection mouse model

Severe infection with the re-emerging enterovirus 71 (EV71 or EV-A71) can cause cardiopulmonary failure. However, in patients' heart and lung, viral protein has not been detected. In mouse models, heart disease has not been reported. EV71-infected brainstem is generally believed to be responsible for the cardiopulmonary collapse. One major limitation in EV71 research is the lack of an efficient oral infection system using non-mouse-adapted clinical isolates. In a robust oral infection NOD/SCID mouse model, we detected EV71 protein at multiple organs, including heart and lung, in 100% of moribund mice with limb paralysis. Infiltrating leukocytes were always detected in heart and muscle, and VP1-positive M2 macrophages were abundant in the lung. Functional dissection on the pathogenesis mechanism revealed severe apoptosis, inflammatory cytokines, and abnormal electrocardiogram (EKG) in orally infected hearts. Therefore, cardiopulmonary disease could be one plausible cause of death in this mouse model. Inoculation of EV71 through an oral route resulted in viral infection in the intestine, viremia, and EV71 appeared to spread to peripheral tissues via blood circulation. Infectious virus was no longer detected in the blood on day 5 post-infection by the plaque formation assay. We demonstrated that both EV71 clinical isolate and cloned virus can target the cardiopulmonary system via a natural infection-like oral route.

Transmission genetics of drug-resistant hepatitis C virus

Antiviral development is plagued by drug resistance and genetic barriers to resistance are needed. For HIV and hepatitis C virus (HCV), combination therapy has proved life-saving. The targets of direct-acting antivirals for HCV infection are NS3/4A protease, NS5A phosphoprotein and NS5B polymerase. Differential visualization of drug-resistant and -susceptible RNA genomes within cells revealed that resistant variants of NS3/4A protease and NS5A phosphoprotein are cis-dominant, ensuring their direct selection from complex environments. Confocal microscopy revealed that RNA replication complexes are genome-specific, rationalizing the non-interaction of wild-type and variant products. No HCV antivirals yet display the dominance of drug susceptibility shown for capsid proteins of other viruses. However, effective inhibitors of HCV polymerase exact such high fitness costs for drug resistance that stable genome selection is not observed. Barriers to drug resistance vary with target biochemistry and detailed analysis of these barriers should lead to the use of fewer drugs.

Viruses are simple organisms that consist of genetic information and a few types of proteins. They cannot replicate on their own, and instead hijack the molecular machinery of a host cell to produce more of themselves.

Inside an infected cell, the genetic information of the virus is replicated and ‘read’ to create viral proteins. These components are then assembled to form a new generation of viruses. During this process, genetic errors may occur that lead to modifications in the viral proteins, and help the virus become resistant to treatment. For instance, a viral protein that used to be targeted by a drug can change slightly and not be recognized anymore.

Currently, the most efficient way to fight drug resistance is to use combination therapy, where several drugs are given at the same time. This strategy is successful, for example to treat infections with the hepatitis C virus, but it is also expensive, especially for developing countries.

An alternative approach is dominant-drug targeting, which exploits the fact that both drug-resistant and drug-susceptible viruses are ‘born’ in the same cell. There, the susceptible viruses can overwhelm and ‘mask’ the benefits of the resistant ones. For example, proteins from resistant strains, which are no longer detected by a treatment, can bind to proteins from susceptible viruses; drugs will still be able to recognize these resulting viral structures. The proteins that operate in such ways are potential dominant-drug targets. However, resistant and susceptible strains can also cohabit without any contacts if their proteins do not interact with each other.

Now, van Buuren et al. screen several viral proteins, including one called NS5A, to test whether a dominant drug target exists for the hepatitis C virus. Only a few molecules of a drug that targets NS5A can stop the virus from growing. In theory, drug-bound NS5A proteins could block their non-drug-bound neighbors, but when these drugs have been used on their own, resistance quickly emerged.

Experiments showed that NS5A is not a dominant drug target because the drug-resistant and drug-susceptible proteins do not mix. Unless ‘forced’ in the laboratory, NS5A proteins only bind to the ones produced by the same strain of virus. This explains why resistant viruses quickly take over when NS5A drugs are the sole treatment. However, other hepatitis C proteins, such as the HCV core protein, are known to mix during the assembly of the virus, and thus are likely be dominant drug targets.

Quasispecies and virus

Quasispecies theory has been instrumental in the understanding of RNA virus population dynamics because it considered for the first time mutation as an integral part of the replication process. The key influences of quasispecies theory on experimental virology have been: (1) to disclose the mutant spectrum nature of viral populations and to evaluate its consequences; (2) to unveil collective properties of genome ensembles that can render a mutant spectrum a unit of selection; and (3) to identify new vulnerability points of pathogenic RNA viruses on three fronts: the need to apply multiple selective constraints (in the form of drug combinations) to minimize selection of treatment-escape variants, to translate the error threshold concept into antiviral designs, and to construct attenuated vaccine viruses through alterations of viral polymerase copying fidelity or through displacements of viral genomes towards unfavorable regions of sequence space. These three major influences on the understanding of viral pathogens preceded extensions of quasispecies to non-viral systems such as bacterial and tumor cell collectivities and prions. These developments are summarized here.

Discovery and Mechanistic Study of Benzamide Derivatives That Modulate Hepatitis B Virus Capsid Assembly

Chronic hepatitis B virus (HBV) infection is a global public health problem. Although the currently approved medications can reliably reduce the viral load and prevent the progression of liver diseases, they fail to cure the viral infection. In an effort toward discovery of novel antiviral agents against HBV, a group of benzamide (BA) derivatives that significantly reduced the amount of cytoplasmic HBV DNA were discovered. The initial lead optimization efforts identified two BA derivatives with improved antiviral activity for further mechanistic studies. Interestingly, similar to our previously reported sulfamoylbenzamides (SBAs), the BAs promote the formation of empty capsids through specific interaction with HBV core protein but not other viral and host cellular components. Genetic evidence suggested that both SBAs and BAs inhibited HBV nucleocapsid assembly by binding to the heteroaryldihydropyrimidine (HAP) pocket between core protein dimer-dimer interfaces. However, unlike SBAs, BA compounds uniquely induced the formation of empty capsids that migrated more slowly in native agarose gel electrophoresis from A36V mutant than from the wild-type core protein. Moreover, we showed that the assembly of chimeric capsids from wild-type and drug-resistant core proteins was susceptible to multiple capsid assembly modulators. Hence, HBV core protein is a dominant antiviral target that may suppress the selection of drug-resistant viruses during core protein-targeting antiviral therapy. Our studies thus indicate that BAs are a chemically and mechanistically unique type of HBV capsid assembly modulators and warranted for further development as antiviral agents against HBV.IMPORTANCE HBV core protein plays essential roles in many steps of the viral replication cycle. In addition to packaging viral pregenomic RNA (pgRNA) and DNA polymerase complex into nucleocapsids for reverse transcriptional DNA replication to take place, the core protein dimers, existing in several different quaternary structures in infected hepatocytes, participate in and regulate HBV virion assembly, capsid uncoating, and covalently closed circular DNA (cccDNA) formation. It is anticipated that small molecular core protein assembly modulators may disrupt one or multiple steps of HBV replication, depending on their interaction with the distinct quaternary structures of core protein. The discovery of novel core protein-targeting antivirals, such as benzamide derivatives reported here, and investigation of their antiviral mechanism may lead to the identification of antiviral therapeutics for the cure of chronic hepatitis B.

Internal Disequilibria and Phenotypic Diversification during Replication of Hepatitis C Virus in a Noncoevolving Cellular Environment

Viral quasispecies evolution upon long-term virus replication in a noncoevolving cellular environment raises relevant general issues, such as the attainment of population equilibrium, compliance with the molecular-clock hypothesis, or stability of the phenotypic profile. Here, we evaluate the adaptation, mutant spectrum dynamics, and phenotypic diversification of hepatitis C virus (HCV) in the course of 200 passages in human hepatoma cells in an experimental design that precluded coevolution of the cells with the virus. Adaptation to the cells was evidenced by increase in progeny production. The rate of accumulation of mutations in the genomic consensus sequence deviated slightly from linearity, and mutant spectrum analyses revealed a complex dynamic of mutational waves, which was sustained beyond passage 100. The virus underwent several phenotypic changes, some of which impacted the virus-host relationship, such as enhanced cell killing, a shift toward higher virion density, and increased shutoff of host cell protein synthesis. Fluctuations in progeny production and failure to reach population equilibrium at the genomic level suggest internal instabilities that anticipate an unpredictable HCV evolution in the complex liver environment.IMPORTANCE Long-term virus evolution in an unperturbed cellular environment can reveal features of virus evolution that cannot be explained by comparing natural viral isolates. In the present study, we investigate genetic and phenotypic changes that occur upon prolonged passage of hepatitis C virus (HCV) in human hepatoma cells in an experimental design in which host cell evolutionary change is prevented. Despite replication in a noncoevolving cellular environment, the virus exhibited internal population disequilibria that did not decline with increased adaptation to the host cells. The diversification of phenotypic traits suggests that disequilibria inherent to viral populations may provide a selective advantage to viruses that can be fully exploited in changing environments.