Global prevalence of adaptive and prolonged infections’ muta-tions in the receptor-binding domain of the SARS-CoV-2 spike protein

Several vaccines with varying efficacies have been developed and are currently administered globally to minimize the spread of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Despite having an RNA-dependent RNA polymerase with a proofreading activity, new variants of SARS-CoV-2 are on the rise periodically. Some of the mutations in these variants, especially mutations on the spike protein, aids the virus in transmission, infectivity and host immune evasion. Further, these mutations also reduce the effectiveness of some of the current vaccines and monoclonal antibodies (mAbs). In the present study, using the available 984,769 SARS-CoV-2 nucleotide sequences on the NCBI database from the end of 2019 till 28 July 2021, we have estimated the global prevalence of so-called ‘adaptive mutations’ and ‘mutations identified in the prolonged infections’, in the receptor-binding domain (RBD) of the spike (S) protein. Irrespective of the geographical regions, in the case of the adaptive mutations, N501Y (48.38 %) was found to be the dominant mutation followed by L452R (17.52 %), T478K (14.31 %), E484K (4.69 %), S477N (3.29 %), K417T (1.64 %), N439K (0.7 %) and S494P (0.7 %). Other mutations were found to be less prevalent (less than 0.7 %). Since the last two months, there has been a massive increase of L452R and T478K mutations (delta variant) in certain areas. In the case of prolonged infections’ mutations (long-term SARS-CoV-2 infections), V483A (0.009 %) was found to be dominant followed by Q493R (0.009 %), while other mutations were found in less than 0.007 % of the studied sequences. The data obtained in this study will aid in the development of better infection control policies thereby curbing the spread of this virus

Biophysical Studies on HCV 1a NS3/4A Protease and Its Catalytic Triad in Wild Type and Mutants by the In Silico Approach

The hepatitis C virus (HCV), of the family flaviviridae, is one of the major causes of chronic liver diseases. Until the year 2012, HCV infections were treated using PEG-interferon and ribavirin combinations, which have a low cure rate and severe side effects. Currently, many direct-acting antivirals (DAAs) are available, e.g. protease inhibitors, NS5A and polymerase inhibitors. These drugs have proven to be efficient in interferon-free treatment combinations and capable of enhancing the cure rate to above 90 %. Unlike PEG-interferon and ribavirin combinations, DAAs select for resistance in HCV. The R155K mutation in the HCV was found to resist all the currently available protease inhibitors. Here, we studied biophysical parameters like pocket (cavity) geometries and stabilizing residues of HCV 1a NS3/4A protease in wild type and mutants. We also studied HCV 1a NS3/4A protease’s catalytic residues: their accessibility, energy, flexibility and binding to Phase II oral protease inhibitor vedroprevir (GS-9451), and compared these parameters between wild type and mutant(s). All these studies were performed using various bioinformatics tools (e.g. Swiss-PdbViewer and Schrödinger’s Maestro) and web servers (e.g. DoGSiteScorer, SRide, ASA-View, WHAT IF, elNémo, CABS-flex, PatchDock and PLIP). From our study, we found that introduction of R155K, A156T or D168A mutation to wild-type NS3/4A protease increases the pocket’s volume, surface (in the R155K mutant, surface decreases), lipo surface and depth and decreases the number of stabilizing residues. Additionally, differences in catalytic residues’ solvent accessibility, energy, root-mean-square deviation (RMSD) and flexibility between wild type and mutants might explain changes in the protease activity and the resistance to protease inhibitors.N/

SARS-CoV-2 Omicron variants and their susceptibility towards Monoclonal Antibodies in a Swedish cohort during 2022-23, studied by next-generation sequencing

Monoclonal antibodies (mAbs) are an important treatment option for COVID-19 caused by SARS-CoV-2, especially in immunosuppressed patients. However, this treatment option can become ineffective due to mutations in the SARS-CoV-2 genome, mainly in the receptor binding domain (RBD) of the spike (S) protein. In the present study 7950 SARS-CoV-2 positive samples from the Uppsala and Örebro regions of central Sweden collected between March 2022 and May 2023 were whole-genome sequenced using next-generation sequencing, mainly with the Nanopore sequencing method. Pango lineages were determined and all single nucleotide polymorphism (SNP) mutations that occurred in these samples were identified. The dominant sublineages changed over time and mutations conferring resistance to currently available mAbs became common. Notable ones are R346T and K444T mutations in the RBD that confer significant resistance against tixagevimab and cilgavimab mAbs. Further, mutations conferring a high-fold resistance to bebtelovimab, such as the K444T and V445P mutations, were also observed in the samples. This study highlights that resistance mutations have over time rendered currently available mAbs ineffective against SARS-CoV-2 in most patients. Therefore, there is a need for continued surveillance of resistance mutations and the development of new mAbs that target more conserved regions of the RBD.N/

SYBR Green II Dye-Based Real-Time Assay for Measuring Inhibitor Activity Against HIV-1 Reverse Transcriptase

There are arrays of in vitro assays to quantify the activity of HIV-1 reverse transcriptase (HIV-1 RT). These assays utilize either chemically customized/labelled nucleotides, or TaqMan probes, or radiolabeled nucleotides/primers. Although several real-time PCR assays exist commercially for measuring the RT activity, which are usually used for quantifying the viral titres, these assays are not optimized for measuring the inhibitory concentrations (IC50) of HIV-1 RT inhibitors. Moreover, a recently established inorganic pyrophosphate-coupled enzyme assay cannot be employed for studying nonphosphorylated nucleoside reverse transcriptase inhibitors (NRTIs). In the present study, we have developed a novel one-step assay with native nucleotide substrates and SYBR Green II dye to determine IC50 values of triphosphorylated NRTIs against HIV-1 RT. Using exact batches of wild-type and mutant RT, and triphosphorylated NRTIs, we showed that our method gave IC50 values for inhibitors similar to that of an earlier published colorimetric assay with BrdUTP substrate (CABS). Our assay should be suitable for high-throughput screening of antiretroviral drugs and could also be suitable for studying drug resistance profiles. Additionally, we also used our assay to study inhibition by AZT in its nonphosphorylated form by supplementing the reaction mixture with necessary kinases and ATP.Swedish Physicians against AIDS Research Fun

Transmission of hepatitis C virus among intravenous drug users in the Uppsala region of Sweden

BACKGROUND: Epidemiology and transmission patterns of hepatitis C virus (HCV) are important subjects as we enter a new era of treatment with directly acting antivirals (DAAs). The highest prevalence of HCV in developed countries is found among intravenous drug users (IDUs), where unsafe needle sharing practices provide the main route of infection. Efforts to prohibit the continuous spread of HCV among these groups have been initiated by the community services and health care providers. Our goal was to understand how HCV was transmitted among IDUs within a limited population group. We provide a retrospective study (2005-2007) of the HCV transmission patterns in a population of IDUs in the Uppsala region of Sweden. METHOD: Eighty-two serum samples were collected from IDUs in Uppsala County. Our reverse transcription nested polymerase chain reaction (RT-nested PCR) and sequencing method enabled a comprehensive genetic analysis for a broad spectrum of genotypes of two relatively conserved regions, NS5B and NS3, that encodes for the viral polymerase and protease, respectively. HCV RNA in serum samples was amplified and sequenced with in-house primers. Sequence similarities between individuals and subgroups were analyzed with maximum likelihood (ML) phylogenetic trees. Published HCV reference sequences from other geographic regions and countries were also included for clarity. RESULTS: Phylogenetic analysis was possible for 59 NS5B (72%) and 29 NS3 (35%) sequences from Uppsala patients. Additionally, we also included 15 NS3 sequences from Örebro patients, making a total of 44 NS3 sequences for the analysis. By analyzing the NS3 sequences, two transmission sets were found between the IDUs (>98% sequence identity), with one set consisting of two individuals and another set consisting of three individuals. In addition, the phylogenetic analysis done with our serum samples displayed clusters that distinguished them from the reference sequences. CONCLUSION: Our method seems to enable us to trace the HCV transmission between IDUs. Furthermore, the method is fairly independent of the time of infection because the method uses relatively conserved HCV sequence regions (i.e. NS5B and NS3).N/

Prevalence of SARS-CoV-2 Omicron Sublineages and Spike Protein Mutations Conferring Resistance against Monoclonal Antibodies in a Swedish Cohort during 2022–2023

Monoclonal antibodies (mAbs) are an important treatment option for COVID-19 caused by SARS-CoV-2, especially in immunosuppressed patients. However, this treatment option can become ineffective due to mutations in the SARS-CoV-2 genome, mainly in the receptor binding domain (RBD) of the spike (S) protein. In the present study, 7950 SARS-CoV-2 positive samples from the Uppsala and Örebro regions of central Sweden, collected between March 2022 and May 2023, were whole-genome sequenced using amplicon-based sequencing methods on Oxford Nanopore GridION, Illumina MiSeq, Illumina HiSeq, or MGI DNBSEQ-G400 instruments. Pango lineages were determined and all single nucleotide polymorphism (SNP) mutations that occurred in these samples were identified. We found that the dominant sublineages changed over time, and mutations conferring resistance to currently available mAbs became common. Notable ones are R346T and K444T mutations in the RBD that confer significant resistance against tixagevimab and cilgavimab mAbs. Further, mutations conferring a high-fold resistance to bebtelovimab, such as the K444T and V445P mutations, were also observed in the samples. This study highlights that resistance mutations have over time rendered currently available mAbs ineffective against SARS-CoV-2 in most patients. Therefore, there is a need for continued surveillance of resistance mutations and the development of new mAbs that target more conserved regions of the RBD.Acknowledgements: We would like to thank the personnel at Clinical Microbiology and Hospital Hygiene, Uppsala University Hospital for their help with handling samples, RNA extraction, PCR, and whole-genome sequencing.Funder: Regional Research Council Mid Sweden; Grant(s): RFR-98011

Identification of a C2-symmetric diol based human immunodeficiency virus protease inhibitor targeting Zika virus NS2B-NS3 protease

Zika virus (ZIKV) is an emerging mosquito-borne flavivirus and infection by ZIKV Asian lineage is known to cause fetal brain anomalies and Guillain-Barrés syndrome. The WHO declared ZIKV a global public health emergency in 2016. However, currently neither vaccines nor antiviral prophylaxis/treatments are available. In this study, we report the identification of a C2-symmetric diol-based Human immunodeficiency virus type-1 (HIV) protease inhibitor active against ZIKV NS2B-NS3 protease. The compound, referred to as 9b, was identified by in silico screening of a library of 6265 protease inhibitors. Molecular dynamics (MD) simulation studies revealed that compound 9b formed a stable complex with ZIKV protease. Interaction analysis of compound 9b's binding pose from the cluster analysis of MD simulations trajectories predicted that 9b mostly interacted with ZIKV NS3. Although designed as an aspartyl protease inhibitor, compound 9b was found to inhibit ZIKV serine protease in vitro with IC50 = 143.25 ± 5.45 µM, in line with the in silico results. Additionally, linear interaction energy method (LIE) was used to estimate binding affinities of compounds 9b and 86 (a known panflavivirus peptide hybrid with IC50 = 1.64 ± 0.015 µM against ZIKV protease). The LIE method correctly predicted the binding affinity of compound 86 to be lower than that of 9b, proving to be superior to the molecular docking methods in scoring and ranking compounds. Since most of the reported ZIKV protease inhibitors are positively charged peptide-hybrids, with our without electrophilic warheads, compound 9b represents a less polar and more drug-like non-peptide hit compound useful for further optimization.Communicated by Ramaswamy Sarma.Scandinavian Society for Antimicrobial Chemotherap

Baseline dasabuvir resistance in Hepatitis C virus from the genotypes 1, 2 and 3 and modeling of the NS5B-dasabuvir complex by the in silico approach

Background: Current combination treatments with direct-acting antiviral agents (DAAs) can cure more than 95% of hepatitis C virus (HCV) infections. However, resistance-associated substitutions (RASs) may emerge and can also be present in treatment-naïve patients. Methods, results and discussion: In this study, a semi-pan-genotypic population sequencing method was developed and used to assess all NS5B amino acid variants between residue positions 310 and 564. Our method successfully sequenced more than 90% of genotype (GT) 1a, 1b, 2b and 3a samples. By using the population sequencing method with a cut-off of 20%, we found the dasabuvir RASs A553V and C445F to be a baseline polymorphism of GT 2b (8 out of 8) and GT 3a (18 out of 18) sequences, respectively. In GT 1a and 1b treatment-naïve subjects (n=25), no high-fold resistance polymorphism/RASs were identified. We further predicted dasabuvir's binding pose with the NS5B polymerase using the in silico methods to elucidate the reasons associated with the resistance of clinically relevant RASs. Dasabuvir was docked at the palm-I site and was found to form hydrogen bonds with the residues S288, I447, Y448, N291 and D318. The RAS positions 316, 414, 448, 553 and 556 were found to constitute the dasabuvir binding pocket.Uppsala-Örebro Regional Research Council, the Erik, Karin and Gösta Selander Foundation and the Scandinavian Society for Antimicrobial Chemotherap