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 KD 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 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 KD 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.

Mutational analysis of the spike protein of SARS-COV-2 isolates revealed atomistic features responsible for higher binding and infectivity.

Introduction: The perpetual appearance of Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-COV-2), and its new variants devastated the public health and social fabric around the world. Understanding the genomic patterns and connecting them to phenotypic attributes is of great interest to devise a treatment strategy to control this pandemic. Materials and Methods: In this regard, computational methods to understand the evolution, dynamics and mutational spectrum of SARS-CoV-2 and its new variants are significantly important. Thus, herein, we used computational methods to screen the genomes of SARS-CoV-2 isolated from Pakistan and connect them to the phenotypic attributes of spike protein; we used stability-function correlation methods, protein-protein docking, and molecular dynamics simulation. Results: Using the Global initiative on sharing all influenza data (GISAID) a total of 21 unique mutations were identified, among which five were reported as stabilizing while 16 were destabilizing revealed through mCSM, DynaMut 2.0, and I-Mutant servers. Protein-protein docking with Angiotensin-converting enzyme 2 (ACE2) and monoclonal antibody (4A8) revealed that mutation G446V in the receptor-binding domain; R102S and G181V in the N-terminal domain (NTD) significantly affected the binding and thus increased the infectivity. The interaction pattern also revealed significant variations in the hydrogen bonding, salt bridges and non-bonded contact networks. The structural-dynamic features of these mutations revealed the global dynamic trend and the finding energy calculation further established that the G446V mutation increases the binding affinity towards ACE2 while R102S and G181V help in evading the host immune response. The other mutations reported supplement these processes indirectly. The binding free energy results revealed that wild type-RBD has a TBE of -60.55 kcal/mol while G446V-RBD reported a TBE of -73.49 kcal/mol. On the other hand, wild type-NTD reported -67.77 kcal/mol of TBE, R102S-NTD reported -51.25 kcal/mol of TBE while G181V-NTD reported a TBE of -63.68 kcal/mol. Conclusions: In conclusion, the current findings revealed basis for higher infectivity and immune evasion associated with the aforementioned mutations and structure-based drug discovery against such variants.

Mutational analysis of the spike protein of SARS-COV-2 isolates revealed atomistic features responsible for higher binding and infectivity

Introduction: The perpetual appearance of Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-COV-2), and its new variants devastated the public health and social fabric around the world. Understanding the genomic patterns and connecting them to phenotypic attributes is of great interest to devise a treatment strategy to control this pandemic. Materials and Methods: In this regard, computational methods to understand the evolution, dynamics and mutational spectrum of SARS-CoV-2 and its new variants are significantly important. Thus, herein, we used computational methods to screen the genomes of SARS-CoV-2 isolated from Pakistan and connect them to the phenotypic attributes of spike protein;we used stability-function correlation methods, protein-protein docking, and molecular dynamics simulation. Results: Using the Global initiative on sharing all influenza data (GISAID) a total of 21 unique mutations were identified, among which five were reported as stabilizing while 16 were destabilizing revealed through mCSM, DynaMut 2.0, and I-Mutant servers. Protein-protein docking with Angiotensin-converting enzyme 2 (ACE2) and monoclonal antibody (4A8) revealed that mutation G446V in the receptor-binding domain;R102S and G181V in the N-terminal domain (NTD) significantly affected the binding and thus increased the infectivity. The interaction pattern also revealed significant variations in the hydrogen bonding, salt bridges and non-bonded contact networks. The structural-dynamic features of these mutations revealed the global dynamic trend and the finding energy calculation further established that the G446V mutation increases the binding affinity towards ACE2 while R102S and G181V help in evading the host immune response. The other mutations reported supplement these processes indirectly. The binding free energy results revealed that wild type-RBD has a TBE of −60.55 kcal/mol while G446V-RBD reported a TBE of −73.49 kcal/mol. On the other hand, wild type-NTD reported −67.77 kcal/mol of TBE, R102S-NTD reported −51.25 kcal/mol of TBE while G181V-NTD reported a TBE of −63.68 kcal/mol. Conclusions: In conclusion, the current findings revealed basis for higher infectivity and immune evasion associated with the aforementioned mutations and structure-based drug discovery against such variants.

Structural plasticity of omicron BA.5 and BA.2.75 for enhanced ACE-dependent entry into cells.

The current study investigated the binding variations among the wilt type, Omicron sub-variants BA.2.75 and BA.5, using protein-protein docking, protein structural graphs (P SG), and molecular simulation methods. HADDOCK predicted docking scores and dissociation constant (KD) revealed tighter binding of these sub-variants in contrast to the WT. Further investigation revealed variations in the hub residues, protein sub-networks, and GlobalMetapath in these variants as compared to the WT. A very unusual dynamic for BA.2.75 and BA.5 was observed, and secondary structure transition can also be witnessed in the loops (44-505). The results show that the flexibility of these three loops is increased by the mutations as an allosteric effect and thus enhances the chances of bonding with the nearby residues to connect and form a stable connection. Furthermore, the additional hydrogen bonding contacts steer the robust binding of these variants in contrast to the wild type. The total binding free energy for the wild type was calculated to be -61.38 kcal/mol, while for BA.2.75 and BA.5 variants the T BE was calculated to be -70.42 kcal/mol and 69.78 kcal/mol, respectively. We observed that the binding of BA.2.75 is steered by the electrostatic interactions while the BA.5 additional contacts are due to the vdW (Van der Waal) energy. From these findings, it can be observed the Spike (S) protein is undergoing structural adjustments to bind efficiently to the hACE2 (human angiotensin-converting enzyme 2) receptor and, in turn, increase entry to the host cells. The current study will aid the development of structure-based drugs against these variants.Communicated by Ramaswamy H. Sarma.

Ultra-Large-Scale Screening of Natural Compounds and Free Energy Calculations Revealed Potential Inhibitors for the Receptor-Binding Domain (RBD) of SARS-CoV-2.

The emergence of immune-evading variants of SARS-CoV-2 further aggravated the ongoing pandemic. Despite the deployments of various vaccines, the acquired mutations are capable of escaping both natural and vaccine-induced immune responses. Therefore, further investigation is needed to design a decisive pharmacological treatment that could efficiently block the entry of this virus into cells. Hence, the current study used structure-based methods to target the RBD of the recombinant variant (Deltacron) of SARS-CoV-2, which was used as a model variant. From the virtual drug screenings of various databases, a total of four hits were identified as potential lead molecules. Key residues were blocked by these molecules with favorable structural dynamic features. The binding free energies further validated the potentials of these molecules. The TBE for MNP was calculated to be -32.86 ± 0.10 kcal/mol, for SANC00222 the TBE was -23.41 ± 0.15 kcal/mol, for Liriodenine the TBE was -34.29 ± 0.07 kcal/mol, while for Carviolin the TBE was calculated to be -27.67 ± 0.12 kcal/mol. Moreover, each complex demonstrated distinct internal motion and a free energy profile, indicating a different strategy for the interaction with and inhibition of the RBD. In conclusion, the current study demands further in vivo and in vitro validation for the possible usage of these compounds as potential drugs against SARS-CoV-2 and its variants.

Blocking key mutated hotspot residues in the RBD of the omicron variant (B.1.1.529) with medicinal compounds to disrupt the RBD-hACE2 complex using molecular screening and simulation approaches† † Electronic supplementary information (ESI) available. See DOI: 10.1039/d2ra00277a

A new variant of SARS-CoV-2 known as the omicron variant (B.1.1.529) reported in South Africa with 30 mutations in the whole spike protein, among which 15 mutations are in the receptor-binding domain, is continuously spreading exponentially around the world. The omicron variant is reported to be highly contagious with antibody-escaping activity. The emergence of antibody-escaping variants is alarming, and thus the quick discovery of small molecule inhibitors is needed. Hence, the current study uses computational drug screening and molecular dynamics simulation approaches (replicated) to identify novel drugs that can inhibit the binding of the receptor-binding domain (RBD) with hACE2. Screening of the North African, East African and North-East African medicinal compound databases by employing a multi-step screening approach revealed four compounds, namely (−)-pipoxide (C1), 2-(p-hydroxybenzyl) benzofuran-6-ol (C2), 1-(4-hydroxy-3-methoxyphenyl)-2-{4-[(E)-3-hydroxy-1-propenyl]-2-methoxyphenoxy}-1,3-propanediol (C3), and Rhein (C4), with excellent anti-viral properties against the RBD of the omicron variant. Investigation of the dynamics demonstrates stable behavior, good residue flexibility profiles, and structural compactness. Validation of the top hits using computational bioactivity analysis, binding free energy calculations and dissociation constant (KD) analysis also indicated the anti-viral properties of these compounds. In conclusion, this study will help in the design and discovery of novel drug therapeutics, which may be used against the emerging omicron variant of SARS-CoV-2. A new variant of SARS-CoV-2 known as the omicron variant (B.1.1.529) reported in South Africa with 30 mutations in the whole spike protein, among which 15 mutations are in the receptor-binding domain, is continuously spreading exponentially around the world.

Subtractive proteomics assisted therapeutic targets mining and designing ensemble vaccine against Candida auris for immune response induction.

The emergence of variants and the reports of co-infection caused by Candida auris in COVID-19 patients adds a further complication to the global pandemic situation. To date, no effective therapy is available for C. auris infections. Thus, characterization of therapeutic targets and designing effective vaccine candidates using subtractive proteomics and immune-informatics approaches is useful tool in controlling the emerging infections associated with SARS-CoV-2. In the current study, subtractive proteomics-assisted annotation of the vaccine targets was performed, which revealed seven vaccine targets. An immunoinformatic-driven approach was then employed to map protein-specific and proteome-wide immunogenic peptides (CTL, B cell, and HTL) for the design of multi-epitope vaccine candidates (MEVCs). The results demonstrated that the vaccine candidates possess strong antigenic features (>0.4 threshold score) and are classified as non-allergenic. Validation of the designed MEVCs through molecular docking, in-silico cloning, and immune simulation further demonstrated the efficacy of the vaccines by producing immune factor titers (ranging from 2500 to 16000 au/mL) i.e., IgM, IgG, IL-6, and Interferon-α. In conclusion, the current study provides a strong impetus in designing anti-fungal strategies against Candida auris.

The Omicron (B.1.1.529) variant of SARS-CoV-2 binds to the hACE2 receptor more strongly and escapes the antibody response: Insights from structural and simulation data.

As SARS-CoV-2 (severe acute respiratory syndrome coronavirus 2) continues to inflict chaos globally, a new variant officially known as B.1.1.529 was reported in South Africa and was found to harbor 30 mutations in the spike protein. It is too early to speculate on transmission and hospitalizations. Hence, more analyses are required, particularly to connect the genomic patterns to the phenotypic attributes to reveal the binding differences and antibody response for this variant, which can then be used for therapeutic interventions. Given the urgency of the required analysis and data on the B.1.1.529 variant, we have performed a detailed investigation to provide an understanding of the impact of these novel mutations on the structure, function, and binding of RBD to hACE2 and mAb to the NTD of the spike protein. The differences in the binding pattern between the wild type and B.1.1.529 variant complexes revealed that the key substitutions Asn417, Ser446, Arg493, and Arg498 in the B.1.1.529 RBD caused additional interactions with hACE2 and the loss of key residues in the B.1.1.529 NTD resulted in decreased interactions with three CDR regions (1-3) in the mAb. Further investigation revealed that B.1.1.529 displayed a stable dynamic that follows a global stability trend. In addition, the dissociation constant (KD), hydrogen bonding analysis, and binding free energy calculations further validated the findings. Hydrogen bonding analysis demonstrated that significant hydrogen bonding reprogramming took place, which revealed key differences in the binding. The total binding free energy using MM/GBSA and MM/PBSA further validated the docking results and demonstrated significant variations in the binding. This study is the first to provide a basis for the higher infectivity of the new SARS-CoV-2 variants and provides a strong impetus for the development of novel drugs against them.

Computational modelling of potentially emerging SARS-CoV-2 spike protein RBDs mutations with higher binding affinity towards ACE2: A structural modelling study.

The spike protein of SARS-CoV-2 and the host ACE2 receptor plays a vital role in the entry to the cell. Among which the hotspot residue 501 is continuously subjected to positive selection pressure and induces unusual virulence. Keeping in view the importance of the hot spot residue 501, we predicted the potentially emerging structural variants of 501 residue. We analyzed the binding pattern of wild type and mutants (Spike RBD) to the ACE2 receptor by deciphering variations in the amino acids' interaction networks by graph kernels along with evolutionary, network metrics, and energetic information. Our analysis revealed that N501I, N501T, and N501V increase the binding affinity and alter the intra and inter-residue bonding networks. The N501T has shown strong positive selection and fitness in other animals. Docking results and repeated simulations (three times) confirmed the structural stability and tighter binding of these three variants, correlated with the previous results following the global stability trend. Consequently, we reported three variants N501I, N501T, and N501V could worsen the situation further if they emerged. The relations between the viral fitness and binding affinity is a complicated game thus the emergence of high affinity mutations in the SARS-CoV-2 RBD brings up the question of whether or not positive selection favours these mutations or not?

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.

Annotation of Potential Vaccine Targets and Design of a Multi-Epitope Subunit Vaccine against Yersinia pestis through Reverse Vaccinology and Validation through an Agent-Based Modeling Approach.

Yersinia pestis is responsible for plague and major pandemics in Asia and Europe. This bacterium has shown resistance to an array of drugs commonly used for the treatment of plague. Therefore, effective therapeutics measurements, such as designing a vaccine that can effectively and safely prevent Y. pestis infection, are of high interest. To fast-track vaccine development against Yersinia pestis, herein, proteome-wide vaccine target annotation was performed, and structural vaccinology-assisted epitopes were predicted. Among the total 3909 proteins, only 5 (rstB, YPO2385, hmuR, flaA1a, and psaB) were shortlisted as essential vaccine targets. These targets were then subjected to multi-epitope vaccine design using different linkers. EAAK, AAY, and GPGPG as linkers were used to link CTL, HTL, and B-cell epitopes, and an adjuvant (beta defensin) was also added at the N-terminal of the MEVC. Physiochemical characterization, such as determination of the instability index, theoretical pI, half-life, aliphatic index, stability profiling, antigenicity, allergenicity, and hydropathy of the ensemble, showed that the vaccine is highly stable, antigenic, and non-allergenic and produces multiple interactions with immune receptors upon docking. In addition, molecular dynamics simulation confirmed the stable binding and good dynamic properties of the vaccine-TLR complex. Furthermore, in silico and immune simulation of the developed MEVC for Y. pestis showed that the vaccine triggered strong immune response after several doses at different intervals. Neutralization of the antigen was observed at the third day of injection. Conclusively, the vaccine designed here for Y. pestis produces an immune response; however, further immunological testing is needed to unveil its real efficacy.

Combined drug repurposing and virtual screening strategies with molecular dynamics simulation identified potent inhibitors for SARS-CoV-2 main protease (3CLpro).

The current coronavirus (SARS-COV-2) pandemic and phenomenal spread to every nook and cranny of the world has raised major apprehensions about the modern public health care system. So far as a result of this epidemic, 4,434,653 confirmed cases and 302,169 deaths are reported. The growing infection rate and death toll demand the use of all possible approaches to design novel drugs and vaccines to curb this disease. In this study, we combined drugs repurposing and virtual drug screening strategies to target 3CLpro, which has an essential role in viral maturation and replication. A total of 31 FDA approved anti-HIV drugs, and Traditional Chinese medicines (TCM) database were screened to find potential inhibitors. As a result, Saquinavir, and five drugs (TCM5280805, TCM5280445, TCM5280343, TCM5280863, and TCM5458190) from the TCM database were found as promising hits. Furthermore, results from molecular dynamics simulation and total binding free energy revealed that Saquinavir and TCM5280805 target the catalytic dyad (His41 and Cys145) and possess stable dynamics behavior. Thus, we suggest that these compounds should be tested experimentally against the SARS-COV-2 as Saquinavir has been reported to inhibit HIV protease experimentally. Considering the intensity of coronavirus dissemination, the present research is in line with the idea of discovering the latest inhibitors against the coronavirus essential pathways to accelerate the drug development cycle.Communicated by Ramaswamy H. Sarma.

Effect of temperature and humidity on coronavirus infection in Pakistan.

Ongoing Coronavirus epidemic (COVID-19) identified first in Wuhan, China posed huge impact on public health and economy around the globe. Both cough and sneeze based droplets or aerosols encapsulated COVID-19 particles are responsible for airborne transmission of this virus and caused an unexpected escalation and high mortality worldwide. Current study intends to investigate the correlation of COVID-19 epidemic with meteorological parameters, particularly temperature and humidity. A data set of Epidemiological data of COVID-19 for highly infected provinces of Pakistan was collected from the official website of (https://www.covid.gov.pk/) and weather data was collected from (https://www.timeanddate.com/) during the time period of 1st March to 30th September 2020. The GrapPad prism 5 Software was used to calculate the mean and standard error of mean (SEM). In the current study the incident of daily covid cases is recorded higher in the month of June while the less number of case were reported in the month of May as compared to the other months (April, May, June, July, September and August) in the four province of Pakistan. We also find out that the incident of Covid19 were high at higher temperature (like the average temperature in the month of June 37 °C) while less cases were reported in May the average temperature was 29.5 °C. Furthermore the incident of covid cases were less reported at low humidity while more intendant with high humidity. Pearson's (r) determine the strength of the relationship between the variables. Pearson's correlation coefficient test employed for data analysis revealed that temperature average (TA) and average humidity is not a significant correlated with COVID-19 pandemic. The results obtained from the current analysis for selected parameters indirect correlation of COVID-19 transmission with temperature variation, and humidity. In the present study association of parameters is not correlated with COVID-19 pandemic, suggested need of more strict actions and control measures for highly populated cities. These findings will be helpful for health regulatory authorities and policy makers to take specific measures to combat COVID-19 epidemic in Pakistan.

Structural-Dynamics and Binding Analysis of RBD from SARS-CoV-2 Variants of Concern (VOCs) and GRP78 Receptor Revealed Basis for Higher Infectivity.

Glucose-regulated protein 78 (GRP78) might be a receptor for SARS-CoV-2 to bind and enter the host cell. Recently reported mutations in the spike glycoprotein unique to the receptor-binding domain (RBD) of different variants might increase the binding and pathogenesis. However, it is still not known how these mutations affect the binding of RBD to GRP78. The current study provides a structural basis for the binding of GRP78 to the different variants, i.e., B.1.1.7, B.1.351, B.1.617, and P.1 (spike RBD), of SARS-CoV-2 using a biomolecular simulation approach. Docking results showed that the new variants bound stronger than the wild-type, which was further confirmed through the free energy calculation results. All-atom simulation confirmed structural stability, which was consistent with previous results by following the global stability trend. We concluded that the increased binding affinity of the B.1.1.7, B.1.351, and P.1 variants was due to a variation in the bonding network that helped the virus induce a higher infectivity and disease severity. Consequently, we reported that the aforementioned new variants use GRP78 as an alternate receptor to enhance their seriousness.

Bioinformatics analysis of the differences in the binding profile of the wild-type and mutants of the SARS-CoV-2 spike protein variants with the ACE2 receptor.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is the causative agent of coronavirus disease 2019 (COVID-19). Reports of new variants that potentially increase virulence and viral transmission, as well as reduce the efficacy of available vaccines, have recently emerged. In this study, we computationally analyzed the N439K, S477 N, and T478K variants for their ability to bind Angiotensin-converting enzyme 2 (ACE2). We used the protein-protein docking approach to explore whether the three variants displayed a higher binding affinity to the ACE2 receptor than the wild type. We found that these variants alter the hydrogen bonding network and the cluster of interactions. Additional salt bridges, hydrogen bonds, and a high number of non-bonded contacts (i.e., non-bonded interactions between atoms in the same molecule and those in other molecules) were observed only in the mutant complexes, allowing efficient binding to the ACE2 receptor. Furthermore, we used a 2.0-µs all-atoms simulation approach to detect differences in the structural dynamic features of the resulting protein complexes. Our findings revealed that the mutant complexes possessed stable dynamics, consistent with the global trend of mutations yielding variants with improved stability and enhanced affinity. Binding energy calculations based on molecular mechanics/generalized Born surface area (MM/GBSA) further revealed that electrostatic interactions principally increased net binding energies. The stability and binding energies of N439K, S477 N, and T478K variants were enhanced compared to the wild-type-ACE2 complex. The net binding energy of the systems was -31.86 kcal/mol for the wild-type-ACE2 complex, -67.85 kcal/mol for N439K, -69.82 kcal/mol for S477 N, and -69.64 kcal/mol for T478K. The current study provides a basis for exploring the enhanced binding abilities and structural features of SARS-CoV-2 variants to design novel therapeutics against the virus.

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.

Immunogenomics guided design of immunomodulatory multi-epitope subunit vaccine against the SARS-CoV-2 new variants, and its validation through in silico cloning and immune simulation.

Reports of the novel and more contagious strains of SARS-CoV-2 originating in different countries have further aggravated the pandemic situation. The recent substitutions in spike protein may be critical for the virus to evade the host's immune system and therapeutics that have already been developed. Thus, this study has employed an immunoinformatics pipeline to target the spike protein of this novel strain to construct an immunogenic epitope (CTL, HTL, and B cell) vaccine against the new variant. Our investigation revealed that 12 different epitopes imparted a critical role in immune response induction. This was validated by an exploration of physiochemical properties and experimental feasibility. In silico and host immune simulation confirmed the expression and induction of both primary and secondary immune factors such as IL, cytokines, and antibodies. The current study warrants further lab experiments to demonstrate its efficacy and safety.

Higher infectivity of the SARS-CoV-2 new variants is associated with K417N/T, E484K, and N501Y mutants: An insight from structural data.

The evolution of the SARS-CoV-2 new variants reported to be 70% more contagious than the earlier one is now spreading fast worldwide. There is an instant need to discover how the new variants interact with the host receptor (ACE2). Among the reported mutations in the Spike glycoprotein of the new variants, three are specific to the receptor-binding domain (RBD) and required insightful scrutiny for new therapeutic options. These structural evolutions in the RBD domain may impart a critical role to the unique pathogenicity of the SARS-CoV-2 new variants. Herein, using structural and biophysical approaches, we explored that the specific mutations in the UK (N501Y), South African (K417N-E484K-N501Y), Brazilian (K417T-E484K-N501Y), and hypothetical (N501Y-E484K) variants alter the binding affinity, create new inter-protein contacts and changes the internal structural dynamics thereby increases the binding and eventually the infectivity. Our investigation highlighted that the South African (K417N-E484K-N501Y), Brazilian (K417T-E484K-N501Y) variants are more lethal than the UK variant (N501Y). The behavior of the wild type and N501Y is comparable. Free energy calculations further confirmed that increased binding of the spike RBD to the ACE2 is mainly due to the electrostatic contribution. Further, we find that the unusual virulence of this virus is potentially the consequence of Darwinian selection-driven epistasis in protein evolution. The triple mutants (South African and Brazilian) may pose a serious threat to the efficacy of the already developed vaccine. Our analysis would help to understand the binding and structural dynamics of the new mutations in the RBD domain of the Spike protein and demand further investigation in in vitro and in vivo models to design potential therapeutics against the new variants.

In Silico Mutagenesis-Based Remodelling of SARS-CoV-1 Peptide (ATLQAIAS) to Inhibit SARS-CoV-2: Structural-Dynamics and Free Energy Calculations.

The prolific spread of COVID-19 caused by a novel coronavirus (SARS-CoV-2) from its epicenter in Wuhan, China, to every nook and cranny of the world after December 2019, jeopardize the prevailing health system in the world and has raised serious concerns about human safety. Multi-directional efforts are made to design small molecule inhibitors, and vaccines and many other therapeutic options are practiced, but their final therapeutic potential is still to be tested. Using the old drug or vaccine or peptides could aid this process to avoid such long experimental procedures. Hence, here, we have repurposed a small peptide (ATLQAIAS) from the previous study, which reported the inhibitory effects of this peptide. We used in silico mutagenesis approach to design more peptides from the native wild peptide, which revealed that substitutions (T2W, T2Y, L3R, and A5W) could increase the binding affinity of the peptide towards the 3CLpro. Furthermore, using MD simulation and free energy calculation confirmed its dynamics stability and stronger binding affinities. Per-residue energy decomposition analysis revealed that the specified substitution significantly increased the binding affinity at the residue level. Our wide-ranging analyses of binding affinities disclosed that our designed peptide owns the potential to hinder the SARS-CoV-2 and will reduce the progression of SARS-CoV-2-borne pneumonia. Our research strongly suggests the experimental and clinical validation of these peptides to curtail the recent corona outbreak.