Mutational landscape and in silico structure models of SARS-CoV-2 spike receptor binding domain reveal key molecular determinants for virus-host interaction

Computational screening combined with well-tempered metadynamics simulations identifies potential TMPRSS2 inhibitors

Type-II transmembrane serine proteases are effective pharmacological targets for host defence against viral entry and in certain cancer cell progressions. These serine proteases cleave viral spike proteins to expose the fusion peptide for cell entry, which is essential to the life cycle of the virus. TMPRSS2 inhibitors can also fight against respiratory viruses that employ them for cell entry. Our study combining virtual screening, all-atom molecular dynamics, and well-tempered metadynamics simulation identifies vicenin-2, neohesperidin, naringin, and rhoifolin as promising TMPRSS2 antagonists. The binding energies obtained are - 16.3, - 15.4, - 13.6, and - 13.8 kcal/mol for vicenin-2, neohesperidin, naringin, and rhoifolin respectively. The RMSD, RMSF, PCA, DCCM, and binding free energy profiles also correlate with the stable binding of these ligands at the active site of TMPRSS2. The study reveals that these molecules could be promising lead molecules for combating future outbreaks of coronavirus and other respiratory viruses.

Discovery of small molecule entry inhibitors targeting the linoleic acid binding pocket of SARS-CoV-2 spike protein

Spike glycoprotein has a significant role in the entry of SARS-CoV-2 to host cells, which makes it a potential drug target. Continued accumulation of non-synonymous mutations in the receptor binding domain of spike protein poses great challenges in identifying antiviral drugs targeting this protein. This study aims to identify potential entry inhibitors of SARS-CoV-2 using virtual screening and molecular dynamics (MD) simulations from three distinct chemical libraries including Pandemic Response Box, Drugbank and DrugCentral, comprising 6971 small molecules. The molecules were screened against a binding pocket identified in the receptor-binding domain (RBD) region of the spike protein which is known as the linoleic acid binding pocket, a highly conserved motif among several SARS-CoV-2 variants. Through virtual screening and binding free energy calculations, we identified four top-scoring compounds, MMV1579787 ([2-Oxo-2-[2-(3-phenoxyphenyl)ethylamino]ethyl]phosphonic acid), Tretinoin, MMV1633963 ((2E,4E)-5-[3-(3,5-dichlorophenoxy)phenyl]penta-2,4-dienoic acid) and Polydatin, which were previously reported to have antibacterial, antifungal or antiviral properties. These molecules showed stable binding on MD simulations over 100 ns and maintained stable interactions with TYR365, PHE338, PHE342, PHE377, TYR369, PHE374 and LEU368 of the spike protein RBD that are found to be conserved among SARS-CoV-2 variants. Our findings were further validated with free energy landscape, principal component analysis and dynamic cross-correlation analysis. Our in silico analysis of binding mode and MD simulation analyses suggest that the identified compounds may impede viral entrance by interacting with the linoleic acid binding site of the spike protein of SARS-CoV-2 regardless of its variants, and they thus demand for further in vitro and in vivo research.Communicated by Ramaswamy H. Sarma.

In silico study identifies peptide inhibitors that negate the effect of non-synonymous mutations in major drug targets of SARS-CoV-2 variants

Since its advent in December 2019, SARS-CoV-2 has diverged into multiple variants with differing levels of virulence owing to the accumulation of mutations in its genome. The structural changes induced by non-synonymous mutations in major drug targets of the virus are known to alter the binding of potential antagonistic inhibitors. Here, we analyzed the effects of non-synonymous mutations in major targets of SARS-CoV-2 in response to potential peptide inhibitors. We screened 12 peptides reported to have anti-viral properties against RBD and 5 peptides against Mpro of SARS-CoV-2 variants using molecular docking and simulation approaches. The mutational landscape of RBD among SARS-CoV-2 variants had 21 non-synonymous mutations across 18 distinct sites. Among these, 14 mutations were present in the RBM region directly interacting with the hACE2 receptor. However, Only 3 non-synonymous mutations were observed in Mpro. We found that LCB1 - a de novo-synthesized peptide has the highest binding affinity to RBD despite non-synonymous mutations in variants and engages key residues of RBD-hACE2 interaction such as K417, E484, N487, and N501. Similarly, an antimicrobial peptide; 2JOS, was identified against Mpro with high binding affinity as it interacts with key residues in dimerization sites such as E166 and F140 crucial for viral replication. MD simulations affirm the stability of RBD-LCB1 and Mpro-2JOS complexes with an average RMSD of 1.902 and 2.476 respectively. We ascertain that LCB1 and 2JOS peptides are promising inhibitors to combat emerging variants of SARS-CoV-2 and thus warrant further investigations using in-vitro and in-vivo analysis.Communicated by Ramaswamy H. Sarma.

Computational analysis of missense variant CYP4F2*3 (V433M) in association with human CYP4F2 dysfunction: a functional and structural impact

BACKGROUND

Cytochrome P450 4F2 (CYP4F2) enzyme is a member of the CYP4 family responsible for the metabolism of fatty acids, therapeutic drugs, and signaling molecules such as arachidonic acid, tocopherols, and vitamin K. Several reports have demonstrated that the missense variant CYP4F2*3 (V433M) causes decreased activity of CYP4F2 and inter-individual variations in warfarin dose in different ethnic groups. However, the molecular pathogenicity mechanism of missense V433M in CYP4F2 at the atomic level has not yet been completely elucidated.

METHODS AND RESULTS

In the current study, we evaluated the effect of the V433M substitution on CYP4F2 using 14 different bioinformatics tools. Further molecular dynamics (MD) simulations were performed to assess the impact of the V433M mutation on the CYP4F2 protein structure, stability, and dynamics. In addition, molecular docking was used to illustrate the effect of V433M on its interaction with vitamin K1. Based on our results, the CYP4F2*3 variant was a damaging amino acid substitution with a destabilizing nature. The simulation results showed that missense V433M affects the dynamics and stability of CYP4F2 by reducing its compactness and stability, which means that it tends to change the overall structural conformation and flexibility of CYP4F2. The docking results showed that the CYP4F2*3 variant decreased the binding affinity between vitamin K1 and CYP4F2, which reduced the activity of CYP4F2*3 compared to native CYP4F2.

CONCLUSIONS

This study determined the molecular pathogenicity mechanism of the CYP4F2*3 variant on the human CYP4F2 protein and provided new information for understanding the structure-function relationship of CYP4F2 and other CYP4 enzymes. These findings will aid in the development of effective drugs and treatment options.

Understanding Mutations in Human SARS-CoV-2 Spike Glycoprotein: A Systematic Review & Meta-Analysis

Genetic variant(s) of concern (VoC) of SARS-CoV-2 have been emerging worldwide due to mutations in the gene encoding spike glycoprotein. We performed comprehensive analyses of spike protein mutations in the significant variant clade of SARS-CoV-2, using the data available on the Nextstrain server. We selected various mutations, namely, A222V, N439K, N501Y, L452R, Y453F, E484K, K417N, T478K, L981F, L212I, N856K, T547K, G496S, and Y369C for this study. These mutations were chosen based on their global entropic score, emergence, spread, transmission, and their location in the spike receptor binding domain (RBD). The relative abundance of these mutations was mapped with global mutation D614G as a reference. Our analyses suggest the rapid emergence of newer global mutations alongside D614G, as reported during the recent waves of COVID-19 in various parts of the world. These mutations could be instrumentally imperative for the transmission, infectivity, virulence, and host immune system's evasion of SARS-CoV-2. The probable impact of these mutations on vaccine effectiveness, antigenic diversity, antibody interactions, protein stability, RBD flexibility, and accessibility to human cell receptor ACE2 was studied in silico. Overall, the present study can help researchers to design the next generation of vaccines and biotherapeutics to combat COVID-19 infection.

In Silico Study of Coumarins: Wedelolactone as a Potential Inhibitor of the Spike Protein of the SARS-CoV-2 Variants

Despite the rigorous global efforts to control SARS-CoV-2 transmission, it continues to pose a serious threat to humans with the frequent emergence of new variants. Thus, robust therapeutics to combat the virus are a desperate need. The SARS-CoV-2 spike (S) protein is an important target protein as it mediates the entry of the virus inside the host cells, which is initiated by the binding of the receptor-binding domain (RBD) to its cognate receptor, angiotensin-converting enzyme 2 (ACE-2). Herein, the inhibition potential of several naturally occurring coumarins was investigated against the spike proteins of SARS-CoV-2 variants using computational approaches. Molecular docking studies revealed 26 coumarins with better binding energies than the reference ligands, molnupiravir and ceftazidime, against the S-RBD of the omicron variant. The top 10 best-docked coumarins were further analyzed to understand their binding interactions against the spike proteins of other variants (wild-type, Alpha, Beta, Gamma, and Delta), and these studies also demonstrated decent binding energies. Physicochemical, QSAR, and pharmacokinetics analyses of the coumarins revealed wedelolactone as the best inhibitor of the spike protein with ideal Lipinski’s drug-likeness and optimal ADMET properties. Furthermore, coarse-grained molecular dynamics (MD) simulation studies of spike protein-wedelolactone complexes validated the stable binding of wedelolactone in the respective binding pockets. As an outcome, wedelolactone could be utilized to develop a potent drug candidate against COVID-19 by blocking the viral entry into the host cell.

Advanced Plasmonic Nanoparticle-Based Techniques for the Prevention, Detection, and Treatment of Current COVID-19

Coronavirus is an ongoing global pandemic caused by severe acute respiratory syndrome coronavirus 2. Coronavirus disease 2019 known as COVID-19 is the worst pandemic since World War II. The outbreak of COVID-19 had a significant repercussion on the health, economy, politics, and environment, making coronavirus-related issues more complicated and becoming one of the most challenging pandemics of the last century with deadly outcomes and a high rate of the reproduction number. There are thousands of different types - or variants - of COVID circulating across the world. Viruses mutate all the time; it emphasizes the critical need for the designing of efficient vaccines to prevent virus infection, early and fast diagnosis, and effective antiviral and protective therapeutics. In this regard, the use of nanotechnology offers new opportunities for the development of novel strategies in terms of prevention, diagnosis, and treatment of COVID-19. This review presents an outline of the platforms developed using plasmonic nanoparticles in the detection, treatment, and prevention of SARS-CoV-2. We select the best strategies in each of these approaches. The properties of metallic plasmon NPs and their relevance in the development of novel point-of-care diagnosis approaches for COVID-19 are highlighted. Also, we discuss the current challenges and the future perspectives looking towards the clinical translation and the commercial aspects of nanotechnology and plasmonic NP-based diagnostic tools and therapy to fight COVID-19 pandemic. The article could be of significance for researchers dedicated to developing suitable plasmonic detection tools and therapy approaches for COVID-19 viruses and future pandemics.

Evolution of Anti-SARS-CoV-2 Therapeutic Antibodies

Since the first COVID-19 reports back in December of 2019, this viral infection caused by SARS-CoV-2 has claimed millions of lives. To control the COVID-19 pandemic, the Food and Drug Administration (FDA) and/or European Agency of Medicines (EMA) have granted Emergency Use Authorization (EUA) to nine therapeutic antibodies. Nonetheless, the natural evolution of SARS-CoV-2 has generated numerous variants of concern (VOCs) that have challenged the efficacy of the EUA antibodies. Here, we review the most relevant characteristics of these therapeutic antibodies, including timeline of approval, neutralization profile against the VOCs, selection methods of their variable regions, somatic mutations, HCDR3 and LCDR3 features, isotype, Fc modifications used in the therapeutic format, and epitope recognized on the receptor-binding domain (RBD) of SARS-CoV-2. One of the conclusions of the review is that the EUA therapeutic antibodies that still retain efficacy against new VOCs bind an epitope formed by conserved residues that seem to be evolutionarily conserved as thus, critical for the RBD:hACE-2 interaction. The information reviewed here should help to design new and more efficacious antibodies to prevent and/or treat COVID-19, as well as other infectious diseases.

SARS-CoV-2 Genetic Diversity and Lineage Dynamics in Egypt during the First 18 Months of the Pandemic

COVID-19 was first diagnosed in Egypt on 14 February 2020. By the end of November 2021, over 333,840 cases and 18,832 deaths had been reported. As part of the national genomic surveillance, 1027 SARS-CoV-2 near whole-genomes were generated and published by the end of July 2021. Here we describe the genomic epidemiology of SARS-CoV-2 in Egypt over this period using a subset of 976 high-quality Egyptian genomes analyzed together with a representative set of global sequences within a phylogenetic framework. A single lineage, C.36, introduced early in the pandemic was responsible for most of the cases in Egypt. Furthermore, to remain dominant in the face of mounting immunity from previous infections and vaccinations, this lineage acquired several mutations known to confer an adaptive advantage. These results highlight the value of continuous genomic surveillance in regions where VOCs are not predominant and the need for enforcement of public health measures to prevent expansion of the existing lineages.

Exploring the role of framework mutations in enabling breadth of a cross-reactive antibody (CR3022) against the SARS-CoV-2 RBD and its variants of concern

Cross-reactive and broadly neutralizing antibodies against surface proteins of diverse strains of rapidly evolving viral pathogens like SARS-CoV-2 can prevent infection and therefore are crucial for the development of effective universal vaccines. While antibodies typically incorporate mutations in their complementarity determining regions during affinity maturation, mutations in the framework regions have been reported as players in determining properties of broadly neutralizing antibodies against HIV and the Influenza virus. We propose an increase in the cross-reactive potential of CR3022 against the emerging SARS- CoV-2 variants of concern through enhanced conformational flexibility. In this study, we use molecular dynamics simulations, in silico mutagenesis, structural modeling, and docking to explore the role of light chain FWR mutations in CR3022, a SARS-CoV anti-spike (S)-protein antibody cross-reactive to the S-protein receptor binding domain of SARS-CoV-2. Our study shows that single substitutions in the light chain framework region of CR3022 with conserved epitopes across SARS-CoV strains allow targeting of diverse antibody epitope footprints that align with the epitopes of recently-categorized neutralizing antibody classes while enabling binding to more than one strain of SARS-CoV-2. Our study has implications for rapid and evolution-based engineering of broadly neutralizing antibodies and reaffirms the role of framework mutations in effective change of antibody orientation and conformation via improved flexibility.Communicated by Ramaswamy H. Sarma.