Search

Since its introduction in the human population, SARS-CoV-2 has evolved into multiple clades, but the events in its intrahost diversification are not well understood. Here, we compare three-dimensional (3D) self-organized neural haplotype maps (SOMs) of SARS-CoV-2 from thirty individual nasopharyngeal diagnostic samples obtained within a 19-day interval in Madrid (Spain), at the time of transition between clades 19 and 20. SOMs have been trained with the haplotype repertoire present in the mutant spectra of the nsp12- and spike (S)-coding regions. Each SOM consisted of a dominant neuron (displaying the maximum frequency), surrounded by a low-frequency neuron cloud. The sequence of the master (dominant) neuron was either identical to that of the reference Wuhan-Hu-1 genome or differed from it at one nucleotide position. Six different deviant haplotype sequences were identified among the master neurons. Some of the substitutions in the neural clouds affected critical sites of the nsp12-nsp8-nsp7 polymerase complex and resulted in altered kinetics of RNA synthesis in an in vitro primer extension assay. Thus, the analysis has identified mutations that are relevant to modification of viral RNA synthesis, present in the mutant clouds of SARS-CoV-2 quasispecies. These mutations most likely occurred during intrahost diversification in several COVID-19 patients, during an initial stage of the pandemic, and within a brief time period.

There is a need for effective anti-COVID-19 treatments, mainly for individuals at risk of severe disease such as the elderly and the immunosuppressed. Drug repositioning has proved effective in identifying drugs that can find a new application for the control of coronavirus disease, in particular COVID-19. The purpose of the present study was to find synergistic antiviral combinations for COVID-19 based on lethal mutagenesis.

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.

Mutant spectra of RNA viruses are important to understand viral pathogenesis and response to selective pressures. There is a need to characterize the complexity of mutant spectra in coronaviruses sampled from infected patients. In particular, the possible relationship between SARS-CoV-2 mutant spectrum complexity and disease associations has not been established. In the present study, we report an ultradeep sequencing (UDS) analysis of the mutant spectrum of amplicons from the nsp12 (polymerase)- and spike (S)-coding regions of 30 nasopharyngeal isolates (diagnostic samples) of SARS-CoV-2 of the first COVID-19 pandemic wave (Madrid, Spain, April 2020) classified according to the severity of ensuing COVID-19. Low-frequency mutations and deletions, counted relative to the consensus sequence of the corresponding isolate, were overwhelmingly abundant. We show that the average number of different point mutations, mutations per haplotype, and several diversity indices was significantly higher in SARS-CoV-2 isolated from patients who developed mild disease than in those associated with moderate or severe disease (exitus). No such bias was observed with RNA deletions. Location of amino acid substitutions in the three-dimensional structures of nsp12 (polymerase) and S suggest significant structural or functional effects. Thus, patients who develop mild symptoms may be a richer source of genetic variants of SARS-CoV-2 than patients with moderate or severe COVID-19.

Mutant spectra of RNA viruses are important to understand viral pathogenesis, and response to selective pressures. There is a need to characterize the complexity of mutant spectra in coronaviruses sampled from infected patients. In particular, the possible relationship between SARS-CoV-2 mutant spectrum complexity and disease associations has not been established. In the present study, we report an ultra-deep sequencing (UDS) analysis of the mutant spectrum of amplicons from the nsp12 (polymerase)- and spike (S)-coding regions of thirty nasopharyngeal isolates (diagnostic samples) of SARS-CoV-2 of the first COVID-19 pandemic wave (Madrid, Spain, April 2020) classified according to the severity of ensuing COVID-19. Low frequency mutations and deletions, counted relative to the consensus sequence of the corresponding isolate, were overwhelmingly abundant. We show that the average number of different point mutations, mutations per haplotype and several diversity indices was significantly higher in SARS-CoV-2 isolated from patients who developed mild disease than in those associated with moderate or severe disease (exitus). No such bias was observed with RNA deletions. Location of amino acid substitutions in the three dimensional structures of nsp12 (polymerase) and S suggest significant structural or functional effects. Thus, patients who develop mild symptoms may be a richer source of genetic variants of SARS-CoV-2 than patients with moderate or severe COVID-19.

Picornavirus genome replication takes place in specialized intracellular membrane compartments that concentrate viral RNA and proteins as well as a number of host factors that also participate in the process. The core enzyme in the replication machinery is the viral RNA-dependent RNA polymerase (RdRP) 3Dpol. Replication requires the primer protein 3B (or VPg) attached to two uridine molecules. 3B uridylylation is also catalysed by 3Dpol. Another critical interaction in picornavirus replication is that between 3Dpol and the precursor 3AB, a membrane-binding protein responsible for the localization of 3Dpol to the membranous compartments at which replication occurs. Unlike other picornaviruses, the animal pathogen foot-and-mouth disease virus (FMDV), encodes three non-identical copies of the 3B (3B1, 3B2, and 3B3) that could be specialized in different functions within the replication complex. Here, we have used a combination of biophysics, molecular and structural biology approaches to characterize the functional binding of FMDV 3B1 to the base of the palm of 3Dpol. The 1.7 Å resolution crystal structure of the FMDV 3Dpol -3B1 complex shows that 3B1 simultaneously links two 3Dpol molecules by binding at the bottom of their palm subdomains in an almost symmetric way. The two 3B1 contact surfaces involve a combination of hydrophobic and basic residues at the N- (G5-P6, R9; Region I) and C-terminus (R16, L19-P20; Region II) of this small protein. Enzyme-Linked Immunosorbent Assays (ELISA) show that the two 3B1 binding sites play a role in 3Dpol binding, with region II presenting the highest affinity. ELISA assays show that 3Dpol has higher binding affinity for 3B1 than for 3B2 or 3B3. Membrane-based pull-down assays show that 3B1 region II, and to a lesser extent also region I play essential roles in mediating the interaction of 3AB with the polymerase and its recruitment to intracellular membranes.

(A) FMDV 3B-3Dpol complex showing the primer peptide in green bound to active site cleft of the polymerase in yellow [20](PDB id. 2F8E), the CVB3 3B-3Dpol complex, showing 3B bound to the back side of the polymerase (in slate) [10] (PDB id. 3CDW) and the EV71 3B-3Dpol complex bound to the base of the polymerase palm (in sand) [11] (PDB id. IKA4). (B) Schematic drawing of the FMDV genome.

(A) 3Dpol- 3Dpol interactions in the AB plane. The reference molecule is shown in green cartoons and the contacting neighbours in grey. (B) Table showing the different contact surfaces calculated with the PISA software [41]. (C) Close up views showing the main interacting regions

Filoviruses depend on the nucleoprotein (NP) to accomplish multiple functions during the viral life cycle. NP is the most abundantly expressed viral protein in infected cells and the main component of the viral nucleocapsid. It can be structurally divided into amino- and carboxy- terminal domains (NTD and CTD). The NTD can homo-oligomerize to interact and protect the (−) ssRNA genome, forming long helical structures. The flexible CTD is responsible for the binding of other nucleocapsid proteins and is involved in the formation of inclusion bodies (IBs)—the cytoplasmic sites of nucleocapsid formation and genome replication. The CTD ends in a ~100-residue globular tail. Měnglà virus (MLAV) is the only member of the new Dianlovirus genus within the Filoviridae family. Their differential characteristics and the possibility of becoming a threat for human health justify the interest in better understanding of its structure and function. In this work, we present the structure of the globular tail of the MLAV NP CTD, showing an overall conformation closely related to that previously reported for the equivalent NP region in MARV. Moreover, analyses of the CTD-CTD interactions in the crystal asymmetric unit revealed a higher-order helicoidal structure. Mutational studies underscore the crucial role of a number of residues, located at the CTD-CTD contact interface, for IB formation. Site-directed mutagenesis of amino acids L653 and F687, involved in the formation of these helicoidal assemblies, abrogate the growth of IBs when the full-length MLAV NP was ectopically expressed in cells. Our findings confirm the role of NP CTD in IB formation, also for MLAV, and focus the attention on particular residues as starting point for further analysis.

Picornavirus genome replication takes place in specialized intracellular membrane compartments that concentrate viral RNA and proteins as well as a number of host factors that also participate in the process. The core enzyme in the replication machinery is the viral RNA-dependent RNA polymerase (RdRP) 3Dpol. Replication requires the primer protein 3B (or VPg) attached to two uridine molecules. 3B uridylylation is also catalysed by 3Dpol. Another critical interaction in picornavirus replication is that between 3Dpol and the precursor 3AB, a membrane-binding protein responsible for the localization of 3Dpol to the membranous compartments at which replication occurs. Unlike other picornaviruses, the animal pathogen foot-and-mouth disease virus (FMDV), encodes three non-identical copies of the 3B (3B1, 3B2, and 3B3) that could be specialized in different functions within the replication complex. Here, we have used a combination of biophysics, molecular and structural biology approaches to characterize the functional binding of FMDV 3B1 to the base of the palm of 3Dpol. The 1.7 Å resolution crystal structure of the FMDV 3Dpol -3B1 complex shows that 3B1 simultaneously links two 3Dpol molecules by binding at the bottom of their palm subdomains in an almost symmetric way. The two 3B1 contact surfaces involve a combination of hydrophobic and basic residues at the N- (G5-P6, R9; Region I) and C-terminus (R16, L19-P20; Region II) of this small protein. Enzyme-Linked Immunosorbent Assays (ELISA) show that the two 3B1 binding sites play a role in 3Dpol binding, with region II presenting the highest affinity. ELISA assays show that 3Dpol has higher binding affinity for 3B1 than for 3B2 or 3B3. Membrane-based pull-down assays show that 3B1 region II, and to a lesser extent also region I play essential roles in mediating the interaction of 3AB with the polymerase and its recruitment to intracellular membranes.

The strictly conserved residues are in red blocks and similar residues in red characters. The FMDV 3B1 residues interacting with FMDV 3Dpol are marked by green asterisks. Residues of EV71 3B previously shown to contact the bottom of the palm of EV71 3Dpol in the X-ray structure of the complex [11](PDB:4IKA) are highlighted in yellow boxes.

E.coli membranes containing, wild-type 3AB1 or 3AB1(P6S/R9A), 3AB1(R16A/L19S) or 3AB1(P6S/R9A/R16A/L19S) mutants, disrupting the 3B1-3Dpol binding site were used in the analysis. (A) Western blot experiment that shows the 3Dpol protein recruited to the membrane. (B) Bars diagram indicating the percentage of 3Dpol bound by the membranes. To calculate the fraction of 3Dpol bound to 3AB, the volume of the 3Dpol band from a control membrane was subtracted from each of the other 3Dpol bands, and these were then normalized to the 3Dpol band pulled down by wild-type 3AB. Assays were performed in triplicate, and the mean value and the standard deviation are reported.

(A) Structure-based sequence alignment of the picornavirus 3Dpol residues located the base of the palm that would participate in interactions with 3B, (B) Structural superimposition of the two quasi-equivalent 3B1 binding sites in FMDV 3Dpol. The polymerase residues and the bound 3B1 regions are shown in sticks, coloured as in Fig 1, but with molecule I shown in semi-transparent. (C) The 3B binding site in EV71 3Dpol, as seen in the X-ray structure of the EV71 3Dpol -3B complex [11] (PDB:4IKA). (D-G) Structural comparisons of the putative 3B binding region in 3Dpol of other representative picornaviruses whose structure is known: the enteroviruses PV [5] (PDB: 1RA7; light blue) (D) and HRV1B [7] (PDB: 1XR6; salmon) (E), the cardiovirus EMCV [14] (PDB: 4NYZ; cyan) (F), and the kobusvirus porcine aichi virus [13](PDB: 6R1I; green).

† Rwork = ∑hkl ||Fobs(hkl)|—|Fcalc(hkl)|| / ∑hkl |Fobs(hkl)|, where Fobs and Fcalc are the structure factors, deduced from measured intensities and calculated from the model, respectively. ‡ Rfree = as for Rwork but for 5% of the total reflections chosen at random and omitted from refinement.

The two 3Dpol molecules (green and yellow cartoons) linked by 3B1 (grey ribbon) also interact with each other in the crystal packing through contacts between the finger subdomains, forming long fibres along the a-b diagonal. (B) Close up of the interactions involving the direct 3Dpol-3Dpol contacts that facilitate fibre formation.

(A) Scheme representing the fluorescence microscopy experiments comparing the distribution of 3Dpol in the presence of 3AB1, wild type and mutant. The left panel shows 3Dpol interacting with the 3AB1 precursor bound to the membrane. The right panel mimics the scenario in the presence of the 3AB1 mutant 3AB1(P6S/R9A/R16A/L19S), unable to bind 3Dpol. (B) Fluorescence images of HeLa cells showing the different distribution of 3Dpol bound to 3AB1 wild type or the 3AB1(P6A/R9A/R16A/L19S) mutant. Upper panels show the control cells transfected with the polymerase only (green), which appears distributed throughout the cell. Middle panels show cells transfected with 3AB1 wild type (red) and 3Dpol (green), where 3Dpol mostly co-localizes with the 3AB1 protein in a continuous compartment in the cytoplasm. The lower panels show cells transfected with the 3Dpol and the 3AB1 mutant (P6A/R9A/R16A/L19S), where 3Dpol recovers its localization throughout the cell. The images shown are representative of the total number of images obtained. (C) Fluorescence Intensity plots comparing the relative distribution of 3Dpol (green) in presence of wild type and mutant 3AB1 proteins (red), left and right panels, respectively. The nucleus is labelled in blue (DAPI).

Summary of FMDV 3Dpol-3B complexes prepared and the number of crystals analysed that allowed obtaining X-ray diffraction data of sufficient quality to solve the structures.

We report that ribavirin exerts an inhibitory and mutagenic activity on SARS-CoV-2-infecting Vero cells, with a therapeutic index higher than 10. Deep sequencing analysis of the mutant spectrum of SARS-CoV-2 replicating in the absence or presence of ribavirin indicated an increase in the number of mutations, but not in deletions, and modification of diversity indices, expected from a mutagenic activity. Notably, the major mutation types enhanced by replication in the presence of ribavirin were A→G and U→C transitions, a pattern which is opposite to the dominance of G→A and C→U transitions previously described for most RNA viruses. Implications of the inhibitory activity of ribavirin, and the atypical mutational bias produced on SARS-CoV-2, for the search for synergistic anti-COVID-19 lethal mutagen combinations are discussed.

Usutu virus (USUV), is a mosquito-borne flavivirus currently spreading outside the African continent producing substantial avian mortality. In contrast, infected humans could exhibit mild neurological symptoms or remain asymptomatic. As in other flaviviruses, the capped USUV genome encodes three structural and seven non-structural (NS) proteins. Among the NS proteins, NS5 plays crucial roles in virus replication, harbouring the capping and methyltransferase (MTase) activities in its N-terminal domain and the RNA-dependent RNA polymerase (RdRP) activity at the C-terminus. In this work, we present the first structural and functional characterization of the USUV MTase domain. The first structure of the USUV MTase has been determined in complex with its natural ligands (S-adenosyl-L-methionine [SAM]) and S-adenosyl-L-homocysteine [SAH]) at 2.2 Å resolution, showing a molecular dimer in the crystal asymmetric unit. One molecule is bound to the methyl donor SAM while the second is bound to the reaction by-product SAH. Both molecules are almost identical and also show a high structural similarity to the MTase domains of other flaviviruses. The structure of the USUV MTase bound to the inhibitor sinefungin at 1.8 Å resolution is also described. Careful comparisons of the interactions in the SAM-binding cavity prompt us to hypothesize about the strength and weakness of the structure-based design of antivirals directed to the SAM/SAH binding site that could be effective to deal with this threat.

Despite the high sustained virological response rates achieved with current directly-acting antiviral agents (DAAs) against hepatitis C virus (HCV), around 2% to 5% of patients do not achieve such a response. Identification of amino acid substitutions associated with treatment failure requires analytical designs, such as subtype-specific ultra-deep sequencing (UDS) methods for HCV characterization and patient management. By deep sequencing analysis of 220 subtyped HCV samples from infected patients who failed therapy, collected from 39 Spanish hospitals, we determined amino acid sequences of the DAA-target proteins NS3, NS5A and NS5B, by UDS of HCV patient samples, in search of resistanceassociated substitutions (RAS). Using this procedure, we have identified six highly represented amino acid substitutions (HRSs) in NS5A and NS5B of HCV, which are not bona fide RAS. They were present frequently in basal and post-treatment virus of patients who failed therapy to different DAA-based therapies. Contrary to several RAS, HRSs belong to the acceptable subset of substitutions according to the PAM250 replacement matrix. Coherently, their mutant frequency, measured by the number of deep sequencing reads within the HCV quasispecies that encode the relevant substitutions, ranged between 90% and 100% in most cases. Also, they have limited predicted disruptive effects on the threedimensional structures of the proteins harboring them. The information on HRSs that will be gathered during sequencing should be relevant not only to help predict treatment outcomes and disease progression but also to further understand HCV population dynamics, which appears much more complex than thought prior to the introduction of deep sequencing.

RNA virus populations are composed of complex mixtures of genomes that are termed mutant spectra. SARS-CoV-2 replicates as a viral quasispecies, and mutations that are detected at low frequencies in a host can be dominant in subsequent variants. We have studied mutant spectrum complexities of SARS-CoV-2 populations derived from thirty nasopharyngeal swabs of patients infected during the first wave (April 2020) in the Hospital Universitario Fundación Jiménez Díaz. The patients were classified according to the COVID-19 severity in mild (non-hospitalized), moderate (hospitalized) and exitus (hospitalized with ICU admission and who passed away due to COVID-19). Using ultra-deep sequencing technologies (MiSeq, Illumina), we have examined four amplicons of the nsp12 (polymerase)-coding region and two amplicons of the spike-coding region. Ultra-deep sequencing data were analyzed with different cut-off frequency for mutation detection. Average number of different point mutations, mutations per haplotype and several diversity indices were significantly higher in SARS-CoV-2 isolated from patients who developed mild disease. A feature that we noted in the SARS-CoV-2 mutant spectra from diagnostic samples is the remarkable absence of mutations at intermediate frequencies, and an overwhelming abundance of mutations at frequencies lower than 10%. Thus, the decrease of the cut-off frequency for mutation detection from 0.5% to 0.1% revealed an increasement (50- to 100 fold) in the number of different mutations. The significantly higher frequency of mutations in virus from patients displaying mild than moderate or severe disease was maintained with the 0.1% cut- off frequency. To evaluate whether the frequency repertoire of amino acid substitutions differed between SARS-CoV-2 and the well characterized hepatitis C virus (HCV), we performed a comparative study of mutant spectra from infected patients using the same bioinformatics pipelines. HCV did not show the deficit of intermediat

Vaults are ubiquitous ribonucleoprotein particles involved in a diversity of cellular processes, with promising applications as nanodevices for delivery of multiple cargos. The vault shell is assembled by the symmetrical association of multiple copies of the major vault protein that, initially, generates half vaults. The pairwise, anti-parallel association of two half vaults produces whole vaults. Here, using a combination of vault recombinant reconstitution and structural techniques, we characterized the molecular determinants for the vault opening process. This process commences with a relaxation of the vault waist, causing the expansion of the inner cavity. Then, local disengagement of amino-terminal domains at the vault midsection seeds a conformational change that leads to the aperture, facilitating access to the inner cavity where cargo is hosted. These results inform a hitherto uncharacterized step of the vault cycle and will aid current engineering efforts leveraging vault for tailored cargo delivery.

Viruses are extraordinary biological entities that can only thrive as obligate intracellular parasites, exploiting other living cellular components in order to reproduce. Despite the diversity in the virosphere, a central component of the viral replication machinery is the nucleic acid polymerases that are responsible for copying the viral genome and for generating genetic diversity. In the RNA virus world, with the exception of retroviruses, the RNA-dependent RNA polymerases (RdRPs) catalyze RNA synthesis, functioning either as single polypeptides, or in a complex with other viral or host components to transcribe and replicate viral RNA genomes. In this Special Issue, we seek to highlight recent progress in our understanding of the structure and function of different viral replication complexes. Four original studies and seven reviews contributed to the issue to increase our general knowledge of RNA virus replication.

RNA viruses typically encode their own RNA-dependent RNA polymerase (RdRP) to ensure genome replication and transcription. The closed “right hand” architecture of RdRPs encircles seven conserved structural motifs (A to G) that regulate the polymerization activity. The four palm motifs, arranged in the sequential order A to D, are common to all known template dependent polynucleotide polymerases, with motifs A and C containing the catalytic aspartic acid residues. Exceptions to this design have been reported in members of the Permutotetraviridae and Birnaviridae families of positive single stranded (+ss) and double-stranded (ds) RNA viruses, respectively. In these enzymes, motif C is located upstream of motif A, displaying a permuted C–A–B–D connectivity. Here we study the details of the replication elongation process in the non-canonical RdRP of the Thosea asigna virus (TaV), an insect virus from the Permutatetraviridae family. We report the X-ray structures of three replicative complexes of the TaV polymerase obtained with an RNA template-primer in the absence and in the presence of incoming rNTPs. The structures captured different replication events and allowed to define the critical interactions involved in: (i) the positioning of the acceptor base of the template strand, (ii) the positioning of the 3’-OH group of the primer nucleotide during RNA replication and (iii) the recognition and positioning of the incoming nucleotide. Structural comparisons unveiled a closure of the active site on the RNA template-primer binding, before rNTP entry. This conformational rearrangement that also includes the repositioning of the motif A aspartate for the catalytic reaction to take place is maintained on rNTP and metal ion binding and after nucleotide incorporation, before translocation.