Higher binding affinity of furin for SARS-CoV-2 spike (S) protein D614G mutant could be associated with higher SARS-CoV-2 infectivity

Impact of African-Specific ACE2 Polymorphisms on Omicron BA.4/5 RBD Binding and Allosteric Communication Within the ACE2-RBD Protein Complex

Severe acute respiratory symptom coronavirus 2 (SARS-CoV-2) infection occurs via the attachment of the spike (S) protein's receptor binding domain (RBD) to human ACE2 (hACE2). Natural polymorphisms in hACE2, particularly at the interface, may alter RBD-hACE2 interactions, potentially affecting viral infectivity across populations. This study identified the effects of six naturally occurring hACE2 polymorphisms with high allele frequency in the African population (S19P, K26R, M82I, K341R, N546D and D597Q) on the interaction with the S protein RBD of the BA.4/5 Omicron sub-lineage through post-molecular dynamics (MD), inter-protein interaction and dynamic residue network (DRN) analyses. Inter-protein interaction analysis suggested that the K26R variation, with the highest interactions, aligns with reports of enhanced RBD binding and increased SARS-CoV-2 susceptibility. Conversely, S19P, showing the fewest interactions and largest inter-protein distances, agrees with studies indicating it hinders RBD binding. The hACE2 M82I substitution destabilized RBD-hACE2 interactions, reducing contact frequency from 92 (WT) to 27. The K341R hACE2 variant, located distally, had allosteric effects that increased RBD-hACE2 contacts compared to WThACE2. This polymorphism has been linked to enhanced affinity for Alpha, Beta and Delta lineages. DRN analyses revealed that hACE2 polymorphisms may alter the interaction networks, especially in key residues involved in enzyme activity and RBD binding. Notably, S19P may weaken hACE2-RBD interactions, while M82I showed reduced centrality of zinc and chloride-coordinating residues, hinting at impaired communication pathways. Overall, our findings show that hACE2 polymorphisms affect S BA.4/5 RBD stability and modulate spike RBD-hACE2 interactions, potentially influencing SARS-CoV-2 infectivity-key insights for vaccine and therapeutic development.

Misinformation and COVID-19 vaccine uptake hesitancy among frontline workers in Tanzania: Do demographic variables matter?

Although COVID-19 vaccination has been widely considered as an important remedy to confront COVID-19, people remain hesitant to take it. The objective of this study was to assess the moderation effects of demographic characteristics on the relationship between forms of misinformation and COVID-19 vaccine uptake hesitancy among frontline workers in Dar es Salaam and Dodoma, Tanzania. Using a sample of 200 respondents, it assessed the differences in ratings on misinformation regarding COVID-19 vaccine based on respondents' demographics. The study used a Five-point Likert scale questionnaire distributed through snowball sampling to frontline workers from Dar es Salaam and Dodoma regions. Data was analyzed using binary logistic regression. It was found that the forms of misinformation revealed were manipulated imposters, satire, fabricated contents and false contents with their connection, which they influenced COVID-19 hesitancy significantly. With exception of age, that significantly moderated hesitancy, this study uncovers that, sex and education level moderated insignificantly in predicting those who are misinformed; misinformed individuals are not any less educated or not based on one's sex, different than individuals who are informed. The study informs policy makers on devising appropriate strategies to promote COVID-19 vaccination uptake among the different contextual demographic variables. Promotion of information, media and health literacy to the general public should be considered to deter spreading of vaccine-related misinformation.

Preventing SARS-CoV-2 infection using Fv-antibodies targeting the proprotein convertase (PPC) cleavage site

Fv-antibodies targeting the proprotein convertase (PPC) region of the SARS-CoV-2 spike protein (SP) were screened from an Fv-antibody library to inhibit SARS-CoV-2 infection. Two selected Fv-antibodies were expressed as soluble recombinant proteins, and their binding affinities were assessed using a surface plasmon resonance biosensor. The binding regions of these Fv-antibodies corresponded to the cleavage sites of furin (S1/S2) and transmembrane serine protease 2 (TMPRSS2, S2'). The neutralizing activities of the two Fv-antibodies were demonstrated using a cell-based infection assay with pseudo-viruses carrying the SP of four different SARS-CoV-2 variants: wild-type (D614), delta (B.1.617.2), omicron (BA.2), and omicron (BA.4/5).

Temperature impacts SARS-CoV-2 spike fusogenicity and evolution

SARS-CoV-2 infects both the upper and lower respiratory tracts, which are characterized by different temperatures (33°C and 37°C, respectively). In addition, fever is a common COVID-19 symptom. SARS-CoV-2 has been shown to replicate more efficiently at low temperatures, but the effect of temperature on different viral proteins remains poorly understood. Here, we investigate how temperature affects the SARS-CoV-2 spike function and evolution. We first observed that increasing temperature from 33°C to 37°C or 39°C increased spike-mediated cell-cell fusion. We then experimentally evolved a recombinant vesicular stomatitis virus expressing the SARS-CoV-2 spike at these different temperatures. We found that spike-mediated cell-cell fusion was maintained during evolution at 39°C but was lost in a high proportion of viruses that evolved at 33°C or 37°C. Consistently, sequencing of the spikes evolved at 33°C or 37°C revealed the accumulation of mutations around the furin cleavage site, a region that determines cell-cell fusion, whereas this did not occur in spikes evolved at 39°C. Finally, using site-directed mutagenesis, we found that disruption of the furin cleavage site had a temperature-dependent effect on spike-induced cell-cell fusion and viral fitness. Our results suggest that variations in body temperature may affect the activity and diversification of the SARS-CoV-2 spike.

IMPORTANCE

When it infects humans, SARS-CoV-2 is exposed to different temperatures (e.g., replication site and fever). Temperature has been shown to strongly impact SARS-CoV-2 replication, but how it affects the activity and evolution of the spike protein remains poorly understood. Here, we first show that high temperatures increase the SARS-CoV-2 spike fusogenicity. Then, we demonstrate that the evolution of the spike activity and variants depends on temperature. Finally, we show that the functional effect of specific spike mutations is temperature-dependent. Overall, our results suggest that temperature may be a factor influencing the activity and adaptation of the SARS-CoV-2 spike in vivo, which will help understanding viral tropism, pathogenesis, and evolution.

Biosensing circulating MicroRNAs in autoinflammatory skin diseases: Focus on Hidradenitis suppurativa

MicroRNAs (miRNAs) play a crucial role in the early diagnosis of autoinflammatory diseases, with Hidradenitis Suppurativa (HS) being a notable example. HS, an autoinflammatory skin disease affecting the pilosebaceous unit, profoundly impacts patients' quality of life. Its hidden nature, with insidious initial symptoms and patient reluctance to seek medical consultation, often leads to a diagnostic delay of up to 7 years. Recognizing the urgency for early diagnostic tools, recent research identified significant differences in circulating miRNA expression, including miR-24-1-5p, miR-146a-5p, miR26a-5p, miR-206, miR338-3p, and miR-338-5p, between HS patients and healthy controls. These miRNAs serve as potential biomarkers for earlier disease detection. Traditional molecular biology techniques, like reverse transcription quantitative-polymerase chain reaction (RT-qPCR), are employed for their detection using specific primers and probes. Alternatively, short peptides offer a versatile and effective means for capturing miRNAs, providing specificity, ease of synthesis, stability, and multiplexing potential. In this context, we present a computational simulation pipeline designed for crafting peptide sequences that can capture circulating miRNAs in the blood of patients with autoinflammatory skin diseases, including HS. This innovative approach aims to expedite early diagnosis and enhance therapeutic follow-up, addressing the critical need for timely intervention in HS and similar conditions.

Elucidating the binding mechanism of SARS-CoV-2 NSP6-TBK1 and structure-based designing of phytocompounds inhibitors for instigating the host immune response

SARS-CoV-2, also referred to as severe acute respiratory syndrome coronavirus 2, is the virus responsible for causing COVID-19, an infectious disease that emerged in Wuhan, China, in December 2019. Among its crucial functions, NSP6 plays a vital role in evading the human immune system by directly interacting with a receptor called TANK-binding kinase (TBK1), leading to the suppression of IFNβ production. Consequently, in the present study we used the structural and biophysical approaches to analyze the effect of newly emerged mutations on the binding of NSP6 and TBK1. Among the identified mutations, four (F35G, L37F, L125F, and I162T) were found to significantly destabilize the structure of NSP6. Furthermore, the molecular docking analysis highlighted that the mutant NSP6 displayed its highest binding affinity with TBK1, exhibiting docking scores of -1436.2 for the wildtype and -1723.2, -1788.6, -1510.2, and -1551.7 for the F35G, L37F, L125F, and I162T mutants, respectively. This suggests the potential for an enhanced immune system evasion capability of NSP6. Particularly, the F35G mutation exhibited the strongest binding affinity, supported by a calculated binding free energy of -172.19 kcal/mol. To disrupt the binding between NSP6 and TBK1, we conducted virtual drug screening to develop a novel inhibitor derived from natural products. From this screening, we identified the top 5 hit compounds as the most promising candidates with a docking score of -6.59 kcal/mol, -6.52 kcal/mol, -6.32 kcal/mol, -6.22 kcal/mol, and -6.21 kcal/mol. The molecular dynamic simulation of top 3 hits further verified the dynamic stability of drugs-NSP6 complexes. In conclusion, this study provides valuable insight into the higher infectivity of the SARS-CoV-2 new variants and a strong rationale for the development of novel drugs against NSP6.

Tissue-type plasminogen activator (tPA) homozygous Tyr471His mutation associates with thromboembolic disease

Tissue-type plasminogen activator (tPA) encoded by PLAT is a major mediator that promotes fibrinolysis and prevents thrombosis. Pathogenetic mutations in PLAT associated with venous thromboembolism have rarely been reported. Here, we report the first case of a homozygous point mutation c.1411T>C (p.Y471H) in PLAT leading to thromboembolic events and conduct related functional studies. The corresponding tPA mutant protein (tPA-Y471H) and wild-type tPA (tPA-WT) were synthesized in vitro, and mutant mice (PLATH/H mice) were constructed. The molecular docking and surface plasmon resonance results indicated that the mutation impeded the hydrogen-bonding interactions between the protease domain of tPA and the kringle 4 domain of plasminogen, and the binding affinity of tPA and plasminogen was significantly reduced with a difference of one order of magnitude. mRNA half-life assay showed that the half-life of tPA-Y471H was shortened. The inferior vena cava thrombosis model showed that the rate of venous thrombosis in PLATH/H mice was 80% compared with 53% in wild-type mice. Our data suggested a novel role for the protease domain of tPA in efficient plasminogen activation, and demonstrated that this tPA mutation could reduce the fibrinolysis function of the body and lead to an increased propensity for thrombosis.

Targeting SARS-CoV-2 Macrodomain-1 to Restore the Innate Immune Response Using In Silico Screening of Medicinal Compounds and Free Energy Calculation Approaches

Among the different drug targets of SARS-CoV-2, a multi-domain protein known as NSP3 is a critical element of the translational and replication machinery. The macrodomain-I, in particular, has been reported to have an essential role in the viral attack on the innate immune response. In this study, we explore natural medicinal compounds and identify potential inhibitors to target the SARS-CoV-2-NSP3 macrodomain-I. Computational modeling and simulation tools were utilized to investigate the structural-dynamic properties using triplicates of 100 ns MD simulations. In addition, the MM/GBSA method was used to calculate the total binding free energy of each inhibitor bound to macrodomain-I. Two significant hits were identified: 3,5,7,4'-tetrahydroxyflavanone 3'-(4-hydroxybenzoic acid) and 2-hydroxy-3-O-beta-glucopyranosyl-benzoic acid. The structural-dynamic investigation of both compounds with macrodomain-I revealed stable dynamics and compact behavior. In addition, the total binding free energy for each complex demonstrated a robust binding affinity, of ΔG -61.98 ± 0.9 kcal/mol for Compound A, while for Compound B, the ΔG was -45.125 ± 2.8 kcal/mol, indicating the inhibitory potential of these compounds. In silico bioactivity and dissociation constant (KD) determination for both complexes further validated the inhibitory potency of each compound. In conclusion, the aforementioned natural products have the potential to inhibit NSP3, to directly rescue the host immune response. The current study provides the basis for novel drug development against SARS-CoV-2 and its variants.

Rationale for combined therapies in severe-to-critical COVID-19 patients

An unprecedented global social and economic impact as well as a significant number of fatalities have been brought on by the coronavirus disease 2019 (COVID-19), produced by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Acute SARS-CoV-2 infection can, in certain situations, cause immunological abnormalities, leading to an anomalous innate and adaptive immune response. While most patients only experience mild symptoms and recover without the need for mechanical ventilation, a substantial percentage of those who are affected develop severe respiratory illness, which can be fatal. The absence of effective therapies when disease progresses to a very severe condition coupled with the incomplete understanding of COVID-19's pathogenesis triggers the need to develop innovative therapeutic approaches for patients at high risk of mortality. As a result, we investigate the potential contribution of promising combinatorial cell therapy to prevent death in critical patients.

Insights on the interaction of SARS-CoV-2 variant B.1.617.2 with antibody CR3022 and analysis of antibody resistance

BACKGROUND

The existence of mutated Delta (B.1.617.2) variants of SARS-CoV-2 causes rapid transmissibility, increase in virulence, and decrease in the effectiveness of public health. Majority of mutations are seen in the surface spike, and they are considered as antigenicity and immunogenicity of the virus. Hence, finding suitable cross antibody or natural antibody and understanding its biomolecular recognition for neutralizing surface spike are crucial for developing many clinically approved COVID-19 vaccines. Here, we aim to design SARS-CoV-2 variant and hence, to understand its mechanism, binding affinity and neutralization potential with several antibodies.

RESULTS

In this study, we modelled six feasible spike protein (S1) configurations for Delta SARS-CoV-2 (B.1.617.2) and identified the best structure to interact with human antibodies. Initially, the impact of mutations at the receptor-binding domain (RBD) of B.1.617.2 was tested, and it is found that all mutations increase the stability of proteins (ΔΔG) and decrease the entropies. An exceptional case is noted for the mutation of G614D variant for which the vibration entropy change is found to be within the range of 0.133-0.004 kcal/mol/K. Temperature-dependent free energy change values (ΔG) for wild type is found to be - 0.1 kcal/mol, whereas all other cases exhibit values within the range of - 5.1 to - 5.5 kcal/mol. Mutation on spike increases the interaction with the glycoprotein antibody CR3022 and the binding affinity (CLUSpro energy =  - 99.7 kcal/mol). The docked Delta variant with the following antibodies, etesevimab, bebtelovimab, BD-368-2, imdevimab, bamlanivimab, and casirivimab, exhibit a substantially decreased docking score (- 61.7 to - 112.0 kcal/mol) and the disappearance of several hydrogen bond interactions.

CONCLUSION

Characterization of antibody resistance for Delta variant with respect to the wild type gives understanding regarding why Delta variant endures the resistance boosted through several trademark vaccines. Several interactions with CR3022 have appeared compared to Wild for Delta variant, and hence, it is suggested that modification on the CR3022 antibody could further improve for the prevention of viral spread. Antibody resistance decreased significantly due to numerous hydrogen bond interactions which clearly indicate that these marketed/launched vaccines (etesevimab) will be effective for Delta variants.

Structural insights into the effect of mutations in the spike protein of SARS-CoV-2 on the binding with human furin protein

The SARS COV-2 and its variants are spreading around the world at an alarming speed, due to its higher transmissibility and the conformational changes caused by mutations. The resulting COVID-19 pandemic has imposed severe health consequences on human health. Several countries of the world including Pakistan have studied its genome extensively and provided productive findings. In the current study, the mCSM, DynaMut2, and I-Mutant servers were used to analyze the effect of identified mutations on the structural stability of spike protein however, the molecular docking and simulations approaches were used to evaluate the dynamics of the bonding network between the wild-type and mutant spike proteins with furin. We addressed the mutational modifications that have occurred in the spike protein of SARS-COV-2 that were found in 215 Pakistani's isolates of COVID-19 patients to study the influence of mutations on the stability of the protein and its interaction with the host cell. We found 7 single amino acid substitute mutations in various domains that reside in spike protein. The H49Y, N74K, G181V, and G446V were found in the S1 domain while the D614A, V622F, and Q677H mutations were found in the central helices of the spike protein. Based on the observation, G181V, G446V, D614A, and V622F mutants were found highly destabilizing and responsible for structural perturbation. Protein-protein docking and molecular simulation analysis with that of furin have predicted that all the mutants enhanced the binding efficiency however, the V622F mutant has greatly altered the binding capacity which is further verified by the K_D value (7.1 E^-14) and therefore may enhance the spike protein cleavage by Furin and increase the rate of infectivity by SARS-CoV-2. On the other hand, the total binding energy for each complex was calculated which revealed -50.57 kcal/mol for the wild type, for G181V -52.69 kcal/mol, for G446V -56.44 kcal/mol, for D614A -59.78 kcal/mol while for V622F the TBE was calculated to be -85.84 kcal/mol. Overall, the current finding shows that these mutations have increased the binding of Furin for spike protein and shows that D614A and V622F have significant effects on the binding and infectivity.

Structural evolution of Delta lineage of SARS-CoV-2

One of the main obstacles in prevention and treatment of COVID-19 is the rapid evolution of the SARS-CoV-2 Spike protein. Given that Spike is the main target of common treatments of COVID-19, mutations occurring at this virulent factor can affect the effectiveness of treatments. The B.1.617.2 lineage of SARS-CoV-2, being characterized by many Spike mutations inside and outside of its receptor-binding domain (RBD), shows high infectivity and relative resistance to existing cures. Here, utilizing a wide range of computational biology approaches, such as immunoinformatics, molecular dynamics (MD), analysis of intrinsically disordered regions (IDRs), protein-protein interaction analyses, residue scanning, and free energy calculations, we examine the structural and biological attributes of the B.1.617.2 Spike protein. Furthermore, the antibody design protocol of Rosetta was implemented for evaluation the stability and affinity improvement of the Bamlanivimab (LY-CoV55) antibody, which is not capable of interactions with the B.1.617.2 Spike. We observed that the detected mutations in the Spike of the B1.617.2 variant of concern can cause extensive structural changes compatible with the described variation in immunogenicity, secondary and tertiary structure, oligomerization potency, Furin cleavability, and drug targetability. Compared to the Spike of Wuhan lineage, the B.1.617.2 Spike is more stable and binds to the Angiotensin-converting enzyme 2 (ACE2) with higher affinity.

Mutational profile confers increased stability of SARS-CoV-2 spike protein in Brazilian isolates

Spike (S) protein has been recognized as a promising molecular target for diagnostic, vaccines and antiviral drugs development for COVID-19. In this study, we analyzed the most predominant mutations in the S protein of Brazilian isolates and predicted the effect of these amino acid alterations to protein conformation. A total of 25,924 sequences were obtained from GISAID for five regions of Brazilian territory (Midwest, North, Northeast, South, and Southeast), according to exclusion criteria. Most of the SARS-CoV-2 isolates belongs to the G clade and showed a large occurrence of D614G, N501Y and L18F substitutions. Prediction effects of these amino acid substitutions on the structure dynamics of the spike protein indicated a positive ΔΔG values and negative ΔΔS_Vib in most cases which is associated to structural stabilization and flexibility reduction of the S protein. Mutations E484K, N501Y and K417N belong to several SARS-CoV-2 variants of concern such as Alpha, Beta, Gamma and Delta, and showed high incidence among Brazilian isolates. These mutations have been described to increase RBD affinity to ACE-2 host and abolishment of RBD affinity to potent neutralizing ant-RBD. The increase in rates of infection and reinfection requires continuous genomic surveillance studies in order to characterize emerging mutations and monitor vaccine efficacy, and thus consideration structural data and dynamics in the observed phenotypes.Communicated by Ramaswamy H. Sarma.

Molecular Dynamics Analysis of Fast-Spreading Severe Acute Respiratory Syndrome Coronavirus 2 Variants and Their Effects on the Interaction with Human Angiotensin-Converting Enzyme 2

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is evolving with mutations in the spike protein, especially in the receptor-binding domain (RBD). The failure of public health measures in some countries to contain the spread of the disease has given rise to novel viral variants with increased transmissibility. However, key questions about how quickly the variants can spread remain unclear. Herein, we performed a structural investigation using molecular dynamics simulations and determined dissociation constant (K _D) values using surface plasmon resonance assays of three fast-spreading SARS-CoV-2 variants, alpha, beta, and gamma, as well as genetic factors in host cells that may be related to the viral infection. Our results suggest that the SARS-CoV-2 variants facilitate their entry into the host cell by moderately increased binding affinities to the human ACE2 receptor, different torsions in hACE2 mediated by RBD variants, and an increased spike exposure time to proteolytic enzymes. We also found that other host cell aspects, such as gene and isoform expression of key genes for the infection (ACE2, FURIN, and TMPRSS2), may have few contributions to the SARS-CoV-2 variant infectivity. In summary, we concluded that a combination of viral and host cell factors allows SARS-CoV-2 variants to increase their abilities to spread faster than the wild type.

Modeling the change in European and US COVID-19 death rates

Motivated by several possible differences in Covid-19 virus strains, age demographics, and face mask wearing between continents and countries, we focussed on changes in Covid death rates in 2020. We have extended our Covid-19 multicompartment model (Khan et al., 2020) to fit cumulative case and death data for 49 European countries and 52 US states and territories during the recent pandemic, and found that the case mortality rate had decreased by at least 80% in most of the US and at least 90% in most of Europe. We found that death rate decreases do not have strong correlations to other model parameters (such as contact rate) or other standard state/national metrics such as population density, GDP, and median age. Almost all the decreases occurred between mid-April and mid-June 2020, which corresponds to the time when many state and national lockdowns were relaxed resulting in surges of new cases. We examine here several plausible causes for this drop—improvements in treatment, face mask wearing, new virus strains, testing, potentially changing demographics of infected patients, and changes in data collection and reporting—but none of their effects are as significant as the death rate changes suggest. In conclusion, this work shows that a two death rate model is effective in quantifying the reported drop in death rates.

Discovery of natural products to block SARS-CoV-2 S-protein interaction with Neuropilin-1 receptor: A molecular dynamics simulation approach

Neuropilin-1 (NRP1) is a widely expressed cell surface receptor protein characterized by its pleiotropic function. Recent reports highlighted NRP1 as an additional entry point of the SARS-CoV-2 virus, enhancing viral infectivity by interacting with the S-protein of SARS-CoV-2. The ubiquitous distribution and mechanism of action of NRP1 enable the SARS-CoV-2 virus to attack multiple organs in the body simultaneously. Therefore, blocking NRP1 is a potential therapeutic approach against SARS-CoV-2 infection. The current study screened the South African natural compounds database (SANCDB) for molecules that can disrupt the SARS-CoV-2 S protein-NRP1 interaction as a potential antiviral target for SARS-CoV-2 cellular entry. Following excessive screening and validation analysis 3-O-Methylquercetin and Esculetin were identified as potential compounds to disrupt the S-protein-NRP1 interaction. Furthermore, to understand the conformational stability and dynamic features between NRP1 interaction with the selected natural products, we performed 200 ns molecular dynamics (MD) simulations. In addition, molecular mechanics-generalized Born surface area (MM/GBSA) was utilized to calculate the free binding energies of the natural products interacting with NRP1. 3-O-methylquercetin showed an inhibitory effect with binding energies ΔG of -25.52 ± 0.04 kcal/mol to NRP1, indicating the possible disruption of the NRP1-S-protein interaction. Our analysis demonstrated that 3-O-methylquercetin presents a potential antiviral compound against SARS-CoV-2 infectivity. These results set the path for future functional in-vitro and in-vivo studies in SARS-CoV-2 research.

Sequence similarity of SARS-CoV-2 and humans: Implications for SARS-CoV-2 detection

Detecting severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) needs human samples, which inevitably contain trace human DNA and RNA. Sequence similarity may cause invalid detection results; however, there is still a lack of gene similarity analysis of SARS-CoV-2 and humans. All publicly reported complete genome assemblies in the Entrez genome database were collected for multiple sequence alignment, similarity and phylogenetic analysis. The complete genomes showed high similarity (>99.88% sequence identity). Phylogenetic analysis divided these viruses into three major clades with significant geographic group effects. Viruses from the United States showed considerable variability. Sequence similarity analysis revealed that SARS-CoV-2 has 612 similar sequences with the human genome and 100 similar sequences with the human transcriptome. The sequence characteristics and genome distribution of these similar sequences were confirmed. The sequence similarity and evolutionary mutations provide indispensable references for dynamic updates of SARS-CoV-2 detection primers and methods.

Characteristic analysis of Omicron-included SARS-CoV-2 variants of concern

In view of the rapid development of the COVID-19 pandemic and SARS-CoV-2 mutation, we characterized the emerging SARS-CoV-2 variants of concern (VOCs) by both bioinformatics methods and experiments. The representative genomic sequences of SARS-CoV-2 VOCs were first downloaded from NCBI, including the prototypic strain, Alpha (B.1.1.7) strain, Beta (B.1.351) strain, Delta (B.1.617.2), and Omicron (B1.1.529) strain. Bioinformatics analysis revealed that the D614G mutation led to formation of a protruding spike (S) in the tertiary structure of spike protein, which could be responsible for the enhanced binding to angiotensin-converting enzyme 2 (ACE2) receptor. The epitope analysis further showed that the S protein antigenicity of the Omicron variant changed dramatically, which was possibly associated with its enhanced ability of immune escape. To verify the bioinformatics results, we performed experiments of pseudovirus infection and protein affinity assay. Notably, we found that the spike protein of Omicron variant showed the weakest infectivity and binding ability among all tested strains. Finally, we also proved this through virus infection experiments, and found that the cytotoxicity of Omicron seems to be not strong enough. The results in this study provide guidelines for prevention and control of COVID-19.

Molecular variants of SARS-CoV-2: antigenic properties and current vaccine efficacy

An ongoing pandemic of newly emerged SARS-CoV-2 has puzzled many scientists and health care policymakers around the globe. The appearance of the virus was accompanied by several distinct antigenic changes, specifically spike protein which is a key element for host cell entry of virus and major target of currently developing vaccines. Some of these mutations enable the virus to attach to receptors more firmly and easily. Moreover, a growing number of trials are demonstrating higher transmissibility and, in some of them, potentially more serious forms of illness related to novel variants. Some of these lineages, especially the Beta variant of concern, were reported to diminish the neutralizing activity of monoclonal and polyclonal antibodies present in both convalescent and vaccine sera. This could imply that these independently emerged variants could make antiviral strategies prone to serious threats. The rapid changes in the mutational profile of new clades, especially escape mutations, suggest the convergent evolution of the virus due to immune pressure. Nevertheless, great international efforts have been dedicated to producing efficacious vaccines with cutting-edge technologies. Despite the partial decrease in vaccines efficacy against worrisome clades, most current vaccines are still effective at preventing mild to severe forms of disease and hospital admission or death due to coronavirus disease 2019 (COVID-19). Here, we summarize existing evidence about newly emerged variants of SARS-CoV-2 and, notably, how well vaccines work against targeting new variants and modifications of highly flexible mRNA vaccines that might be required in the future.

Comparative mutational analysis of SARS-CoV-2 isolates from Pakistan and structural-functional implications using computational modelling and simulation approaches

SARS-CoV-2, an RNA virus, has been prone to high mutations since its first emergence in Wuhan, China, and throughout its spread. Its genome has been sequenced continuously by many countries, including Pakistan, but the results vary. Understanding its genomic patterns and connecting them with phenotypic features will help in devising therapeutic strategies. Thus, in this study, we explored the mutation landscape of 250 Pakistani isolates of SARS-CoV-2 genomes to check the genome diversity and examine the impact of these mutations on protein stability and viral pathogenesis in comparison with a reference sequence (Wuhan NC 045512.2). Our results revealed that structural proteins mainly exhibit more mutations than others in the Pakistani isolates; in particular, the nucleocapsid protein is highly mutated. In comparison, the spike protein is the most mutated protein globally. Furthermore, nsp12 was found to be the most mutated NSP in the Pakistani isolates and worldwide. Regarding accessory proteins, ORF3A is the most mutated in the Pakistani isolates, whereas ORF8 is highly mutated in world isolates. These mutations decrease the structural stability of their proteins and alter different biological pathways. Molecular docking, the dissociation constant (KD), and MM/GBSA analysis showed that mutations in the S protein alter its binding with ACE2. The spike protein mutations D614G-S943T-V622F (-75.17 kcal/mol), D614G-Q677H (-75.78 kcal/mol), and N74K-D614G (-73.84 kcal/mol) exhibit stronger binding energy than the wild type (-66.34 kcal/mol), thus increasing infectivity. Furthermore, the simulation results strongly corroborated the predicted protein servers. Our analysis findings also showed that E, M, ORF6, ORF7A, ORF7B, and ORF10 are the most stable coding genes; they may be suitable targets for vaccine and drug development.

GlycoGrip: Cell Surface-Inspired Universal Sensor for Betacoronaviruses

Inspired by the role of cell-surface glycoproteins as coreceptors for pathogens, we report the development of GlycoGrip: a glycopolymer-based lateral flow assay for detecting SARS-CoV-2 and its variants. GlycoGrip utilizes glycopolymers for primary capture and antispike antibodies labeled with gold nanoparticles for signal-generating detection. A lock-step integration between experiment and computation has enabled efficient optimization of GlycoGrip test strips which can selectively, sensitively, and rapidly detect SARS-CoV-2 and its variants in biofluids. Employing the power of the glycocalyx in a diagnostic assay has distinct advantages over conventional immunoassays as glycopolymers can bind to antigens in a multivalent capacity and are highly adaptable for mutated strains. As new variants of SARS-CoV-2 are identified, GlycoGrip will serve as a highly reconfigurable biosensor for their detection. Additionally, via extensive ensemble-based docking simulations which incorporate protein and glycan motion, we have elucidated important clues as to how heparan sulfate and other glycocalyx components may bind the spike glycoprotein during SARS-CoV-2 host-cell infection. GlycoGrip is a promising and generalizable alternative to costly, labor-intensive RT-PCR, and we envision it will be broadly useful, including for rural or low-income populations that are historically undertested and under-reported in infection statistics.

SARS-CoV-2 new variants: Characteristic features and impact on the efficacy of different vaccines

The emergence of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) and its new variants reported in different countries have posed a serious threat to human health and social fabrics worldwide. In addition, these new variants hindered the efforts of vaccines and other therapeutic developments. In this review article, we explained the emergence of new variants of SARS-CoV-2, their transmission risk, mortality rate, and, more importantly, the impact of each new variant on the efficacy of the developed vaccines reported in different literature and findings. The literature reported that with the emergence of new variants, the efficacy of different vaccines is declined, hospitalization and the risk of reinfection is increased. The reports concluded that the emergence of a variant that entirely evades the immune response triggered by the vaccine is improbable. The emergence of new variants and reports of re-infections are creating a more distressing situation and therefore demands further investigation to formulate an effective therapeutic strategy.

Computational prediction of the effect of amino acid changes on the binding affinity between SARS-CoV-2 spike RBD and human ACE2

The association of the receptor binding domain (RBD) of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) spike protein with human angiotensin-converting enzyme 2 (hACE2) represents the first required step for cellular entry. SARS-CoV-2 has continued to evolve with the emergence of several novel variants, and amino acid changes in the RBD have been implicated with increased fitness and potential for immune evasion. Reliably predicting the effect of amino acid changes on the ability of the RBD to interact more strongly with the hACE2 can help assess the implications for public health and the potential for spillover and adaptation into other animals. Here, we introduce a two-step framework that first relies on 48 independent 4-ns molecular dynamics (MD) trajectories of RBD-hACE2 variants to collect binding energy terms decomposed into Coulombic, covalent, van der Waals, lipophilic, generalized Born solvation, hydrogen bonding, π-π packing, and self-contact correction terms. The second step implements a neural network to classify and quantitatively predict binding affinity changes using the decomposed energy terms as descriptors. The computational base achieves a validation accuracy of 82.8% for classifying single-amino acid substitution variants of the RBD as worsening or improving binding affinity for hACE2 and a correlation coefficient of 0.73 between predicted and experimentally calculated changes in binding affinities. Both metrics are calculated using a fivefold cross-validation test. Our method thus sets up a framework for screening binding affinity changes caused by unknown single- and multiple-amino acid changes offering a valuable tool to predict host adaptation of SARS-CoV-2 variants toward tighter hACE2 binding.

Tomatidine and Patchouli Alcohol as Inhibitors of SARS-CoV-2 Enzymes (3CLpro, PLpro and NSP15) by Molecular Docking and Molecular Dynamics Simulations

Considering the current dramatic and fatal situation due to the high spreading of SARS-CoV-2 infection, there is an urgent unmet medical need to identify novel and effective approaches for prevention and treatment of Coronavirus disease (COVID 19) by re-evaluating and repurposing of known drugs. For this, tomatidine and patchouli alcohol have been selected as potential drugs for combating the virus. The hit compounds were subsequently docked into the active site and molecular docking analyses revealed that both drugs can bind the active site of SARS-CoV-2 3CLpro, PLpro, NSP15, COX-2 and PLA2 targets with a number of important binding interactions. To further validate the interactions of promising compound tomatidine, Molecular dynamics study of 100 ns was carried out towards 3CLpro, NSP15 and COX-2. This indicated that the protein-ligand complex was stable throughout the simulation period, and minimal backbone fluctuations have ensued in the system. Post dynamic MM-GBSA analysis of molecular dynamics data showed promising mean binding free energy 47.4633 ± 9.28, 51.8064 ± 8.91 and 54.8918 ± 7.55 kcal/mol, respectively. Likewise, in silico ADMET studies of the selected ligands showed excellent pharmacokinetic properties with good absorption, bioavailability and devoid of toxicity. Therefore, patchouli alcohol and especially, tomatidine may provide prospect treatment options against SARS-CoV-2 infection by potentially inhibiting virus duplication though more research is guaranteed and secured.

Emergence and Spread of a B.1.1.28-Derived P.6 Lineage with Q675H and Q677H Spike Mutations in Uruguay

Uruguay controlled the viral dissemination during the first nine months of the SARS-CoV-2 pandemic. Unfortunately, towards the end of 2020, the number of daily new cases exponentially increased. Herein, we analyzed the country-wide genetic diversity of SARS-CoV-2 between November 2020 and April 2021. We identified that the most prevalent viral variant during the first epidemic wave in Uruguay (December 2020-February 2021) was a B.1.1.28 sublineage carrying Spike mutations Q675H + Q677H, now designated as P.6, followed by lineages P.2 and P.7. P.6 probably arose around November 2020, in Montevideo, Uruguay's capital department, and rapidly spread to other departments, with evidence of further local transmission clusters; it also spread sporadically to the USA and Spain. The more efficient dissemination of lineage P.6 with respect to P.2 and P.7 and the presence of mutations (Q675H and Q677H) in the proximity of the key cleavage site at the S1/S2 boundary suggest that P.6 may be more transmissible than other lineages co-circulating in Uruguay. Although P.6 was replaced by the variant of concern (VOC) P.1 as the predominant lineage in Uruguay since April 2021, the monitoring of the concurrent emergence of Q675H + Q677H in VOCs should be of worldwide interest.

Comparative study of predicted miRNA between Indonesia and China (Wuhan) SARS-CoV-2: a bioinformatics analysis

BACKGROUND

Several reports on the discovery of SARS-CoV-2 mutations and variations in Indonesia COVID-19 cases led to genomic dysregulation with the first pandemic cases in Wuhan, China. MicroRNA (miRNA) plays an important role in this genetic regulation and contributes to the enhancement of viral RNA binding through the host mRNA.

OBJECTIVE

This research is aimed to detect miRNA targets of SARS-CoV-2 and examines their role in Indonesia cases against Wuhan cases.

METHODS

SARS-CoV-2 sequences were obtained from GISAID ( https://www.gisaid.org/ ), NCBI ( https://ncbi.nlm.nih.gov ), and National Genomics Data Center ( https://bigd.big.ac.cn/gwh/ ) databases. MiRDB ( https://github.com/gbnegrini/mirdb-custom-target-search ) was used to annotate and predict target human mature miRNAs. For statistical analysis, we utilized a series chi-square test to obtain significant miRNA. DIANA-miRPath v3.0 ( http://www.microrna.gr/miRPathv3 ) analyzed the Gene Ontology of mature miRNAs.

RESULT

The statistical results detected five significant miRNAs. Two miRNAs: hsa-miR-4778-5p and hsa-miR-4531 were consistently found in the majority of Wuhan samples, while they were only found in less than half of the Indonesia samples. The other three miRNA, hsa-miR-6844, hsa-miR-627-5p, and hsa-miR-3674, were discovered in most samples in both groups but with a significant difference ratio. Among these five significant miRNA targets, hsa-miR-6844 is the only miRNA that has an association with the ORF1ab gene of SARS-CoV-2.

CONCLUSION

The Gene Ontology analysis of five significant miRNA targets indicates a significant role in inflammation and the immune system. The specific detection of host miRNAs in this study shows that there are differences in the characteristics of SARS-CoV-2 between Indonesia and Wuhan.

Recent progress on the mutations of SARS-CoV-2 spike protein and suggestions for prevention and controlling of the pandemic

Severe acute respiratory syndrome coronavirus type 2 (SARS-CoV-2) infection has caused a global pandemic in the past year, which poses continuing threat to human beings. To date, more than 3561 mutations in the viral spike protein were identified, including 2434 mutations that cause amino acid changes with 343 amino acids located in the viral receptor-binding domain (RBD). Among these mutations, the most representative ones are substitution mutations such as D614G, N501Y, Y453F, N439K/R, P681H, K417N/T, and E484K, and deletion mutations of ΔH69/V70 and Δ242-244, which confer the virus with enhanced infectivity, transmissibility, and resistance to neutralization. In this review, we discussed the recent findings of SARS-CoV-2 for highlighting mutations and variants on virus transmissibility and pathogenicity. Moreover, several suggestions for prevention and controlling the pandemic are also proposed.

Accelerating COVID-19 Research Using Molecular Dynamics Simulation

The COVID-19 pandemic has emerged as a global medico-socio-economic disaster. Given the lack of effective therapeutics against SARS-CoV-2, scientists are racing to disseminate suggestions for rapidly deployable therapeutic options, including drug repurposing and repositioning strategies. Molecular dynamics (MD) simulations have provided the opportunity to make rational scientific breakthroughs in a time of crisis. Advancements in these technologies in recent years have become an indispensable tool for scientists studying protein structure, function, dynamics, interactions, and drug discovery. Integrating the structural data obtained from high-resolution methods with MD simulations has helped in comprehending the process of infection and pathogenesis, as well as the SARS-CoV-2 maturation in host cells, in a short duration of time. It has also guided us to identify and prioritize drug targets and new chemical entities, and to repurpose drugs. Here, we discuss how MD simulation has been explored by the scientific community to accelerate and guide translational research on SARS-CoV-2 in the past year. We have also considered future research directions for researchers, where MD simulations can help fill the existing gaps in COVID-19 research.

D614G Substitution of SARS-CoV-2 Spike Protein Increases Syncytium Formation and Virus Titer via Enhanced Furin-Mediated Spike Cleavage

Since the D614G substitution in the spike (S) protein of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) emerged, the variant strain has undergone a rapid expansion to become the most abundant strain worldwide. Therefore, this substitution may provide an advantage for viral spreading. To explore the mechanism, we analyzed 18 viral isolates containing S proteins with either G614 or D614 (S-G614 and S-D614, respectively). The plaque assay showed a significantly higher virus titer in S-G614 than in S-D614 isolates. We further found increased cleavage of the S protein at the furin substrate site, a key event that promotes syncytium formation, in S-G614 isolates. The enhancement of the D614G substitution in the cleavage of the S protein and in syncytium formation has been validated in cells expressing S protein. The effect on the syncytium was abolished by furin inhibitor treatment and mutation of the furin cleavage site, suggesting its dependence on cleavage by furin. Our study pointed to the impact of the D614G substitution on syncytium formation through enhanced furin-mediated S cleavage, which might increase the transmissibility and infectivity of SARS-CoV-2 strains containing S-G614.

The virological impacts of SARS-CoV-2 D614G mutation

The coronavirus diseases 2019 (COVID-19) caused by the infection of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) emerged in December 2019 has caused more than 140 million infections worldwide by the end of April 2021. As an enveloped single-stranded positive-sense RNA virus, SARS-CoV-2 underwent constant evolution that produced novel variants carrying mutation conferring fitness advantages. The current prevalent D614G variant, with glycine substituted for aspartic acid at position 614 in the spike glycoprotein, is one of such variants that became the main circulating strain worldwide in a short period of time. Over the past year, intensive studies from all over the world had defined the epidemiological characteristics of this highly contagious variant and revealed the underlying mechanisms. This review aims at presenting an overall picture of the impacts of D614G mutation on virus transmission, elucidating the underlying mechanisms of D614G in virus pathogenicity, and providing insights into the development of effective therapeutics.

SARS-CoV-2: Possible recombination and emergence of potentially more virulent strains

COVID-19 is challenging healthcare preparedness, world economies, and livelihoods. The infection and death rates associated with this pandemic are strikingly variable in different countries. To elucidate this discrepancy, we analyzed 2431 early spread SARS-CoV-2 sequences from GISAID. We estimated continental-wise admixture proportions, assessed haplotype block estimation, and tested for the presence or absence of strains' recombination. Herein, we identified 1010 unique missense mutations and seven different SARS-CoV-2 clusters. In samples from Asia, a small haplotype block was identified, whereas samples from Europe and North America harbored large and different haplotype blocks with nonsynonymous variants. Variant frequency and linkage disequilibrium varied among continents, especially in North America. Recombination between different strains was only observed in North American and European sequences. In addition, we structurally modelled the two most common mutations, Spike_D614G and Nsp12_P314L, which suggested that these linked mutations may enhance viral entry and replication, respectively. Overall, we propose that genomic recombination between different strains may contribute to SARS-CoV-2 virulence and COVID-19 severity and may present additional challenges for current treatment regimens and countermeasures. Furthermore, our study provides a possible explanation for the substantial second wave of COVID-19 presented with higher infection and death rates in many countries.

Cellular and Humoral Immune Responses in Covid-19 and Immunotherapeutic Approaches

Coronavirus disease 2019 (Covid-19), caused by the novel coronavirus severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), can range in severity from asymptomatic to severe/critical disease. SARS-CoV-2 uses angiotensin-converting enzyme 2 to infect cells leading to a strong inflammatory response, which is most profound in patients who progress to severe Covid-19. Recent studies have begun to unravel some of the differences in the innate and adaptive immune response to SARS-CoV-2 in patients with different degrees of disease severity. These studies have attributed the severe form of Covid-19 to a dysfunctional innate immune response, such as a delayed and/or deficient type I interferon response, coupled with an exaggerated and/or a dysfunctional adaptive immunity. Differences in T-cell (including CD4+ T-cells, CD8+ T-cells, T follicular helper cells, γδ-T-cells, and regulatory T-cells) and B-cell (transitional cells, double-negative 2 cells, antibody-secreting cells) responses have been identified in patients with severe disease compared to mild cases. Moreover, differences in the kinetic/titer of neutralizing antibody responses have been described in severe disease, which may be confounded by antibody-dependent enhancement. Importantly, the presence of preexisting autoantibodies against type I interferon has been described as a major cause of severe/critical disease. Additionally, priorVaccine and multiple vaccine exposure, trained innate immunity, cross-reactive immunity, and serological immune imprinting may all contribute towards disease severity and outcome. Several therapeutic and preventative approaches have been under intense investigations; these include vaccines (three of which have passed Phase 3 clinical trials), therapeutic antibodies, and immunosuppressants.

Free SARS-CoV-2 Spike Protein S1 Particles May Play a Role in the Pathogenesis of COVID-19 Infection

The imbalance of the renin-angiotensin system is currently considered as a potentially important factor of the pathogenesis of COVID-19 disease. It has been shown previously in the murine model, that the expression of angiotensin-converting enzyme 2 (ACE2) on the cell surface is downregulated in response to the infection by SARS-CoV virus or recombinant spike protein (S protein) alone. In the case of natural infection, circulation of the S protein in a soluble form is unlikely. However, in SARS-CoV-2, a large fraction of S protein trimers is pre-processed during virion morphogenesis due to the presence of furin protease cleavage site between the S1 and S2 subunits. Therefore, S protein transition into the fusion conformation may be accompanied by the separation of the S1 subunits carrying the receptor-binding domains from the membrane-bound S2 subunits. The fate of the S1 particles shed due to the spontaneous "firing" of some S protein trimers exposed on the virions and on the surface of infected cells has been never investigated. We hypothesize that the soluble S1 subunits of the SARS-CoV-2 S protein shed from the infected cells and from the virions in vivo may bind to the ACE2 and downregulate cell surface expression of this protein. The decrease in the ACE2 activity on the background of constant or increased ACE activity in the lungs may lead to the prevalence of angiotensin II effects over those of angiotensin (1-7), thus promoting thrombosis, inflammation, and pulmonary damage. This hypothesis also suggests the association between less pronounced shedding of the S1 particles reported for the S protein carrying the D614G mutation (vs. the wild type D614 protein), and lack of increased severity of the COVID-19 infection caused by the mutant (D614G) SARS-CoV-2 strain, despite its higher infectivity and higher in vivo viral load.

Functional importance of the D614G mutation in the SARS-CoV-2 spike protein

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is an enveloped virus which binds its cellular receptor angiotensin-converting enzyme 2 (ACE2) and enters hosts cells through the action of its spike (S) glycoprotein displayed on the surface of the virion. Compared to the reference strain of SARS-CoV-2, the majority of currently circulating isolates possess an S protein variant characterized by an aspartic acid-to-glycine substitution at amino acid position 614 (D614G). Residue 614 lies outside the receptor binding domain (RBD) and the mutation does not alter the affinity of monomeric S protein for ACE2. However, S(G614), compared to S(D614), mediates more efficient ACE2-mediated transduction of cells by S-pseudotyped vectors and more efficient infection of cells and animals by live SARS-CoV-2. This review summarizes and synthesizes the epidemiological and functional observations of the D614G spike mutation, with focus on the biochemical and cell-biological impact of this mutation and its consequences for S protein function. We further discuss the significance of these recent findings in the context of the current global pandemic.