Computational Investigation of Increased Virulence and Pathogenesis of SARS-CoV-2 Lineage B.1.1.7

New variants of SARS-CoV-2 are being reported worldwide. More specifically, the variants reported in South Africa (501Y.V2) and United Kingdom (B.1.1.7) were found to be more contagious than the wild type. There are also speculations that the variants might evade the host immune responses induced by currently available vaccines and develop resistance to drugs under consideration. The first step of viral infection in COVID-19, occurs through the interaction of receptor binding domain (RBD) of the spike protein with peptidase domain of the human ACE-2 (hACE-2) receptor. So, possibly the mutations in the RBD domain of spike protein in the new variants could modulate the protein-protein interaction with hACE-2 receptor leading to the increased virulence. In this study, we aim to get molecular level understanding into the mechanism behind the increased infection rate due to such mutations in these variants. We have computationally studied the interaction of the spike protein in both wild-type and B.1.1.7 variant with hACE-2 receptor using combined molecular dynamics and binding free energy calculations using molecular mechanics-Generalized Born surface area (MM-GBSA) approach. The binding free energies computed using configurations from minimization run and low temperature simulation show that mutant variant of spike protein has increased binding affinity for hACE-2 receptor (i.e. {Delta}{Delta}G(N501Y,A570D) is in the range -20.4 to -21.4 kcal/mol)The residue-wise decomposition analysis and intermolecular hydrogen bond analysis evidenced that the N501Y mutation has increased interaction between RBD of spike protein with ACE-2 receptor. We have also carried out calculations using density functional theory and the results evidenced the increased interaction between three pairs of residues (TYR449 (spike)-ASP38 (ACE-2), TYR453-HIE34 and TYR501-LYS353) in the variant that could be attributed to its increased virulence. The free energies of wild-type and mutant variants of the spike protein computed from MM-GBSA approach suggests that latter variant is stable by about -10.4 kcal/mol when compared to wild type suggesting that it will be retained in the evolution due to increased stability. We demonstrate that with the use of the state-of-the art of computational approaches, we can in advance predict the more virulent nature of variants of SARS-CoV-2 and alert the world health-care system.

Structural and functional insights of STAT2-NS5 interaction for the identification of NS5 antagonist - An approach for restoring interferon signaling.

Dengue is a mosquito-borne viral infection caused by Dengue virus (DENV) and is an emerging concern in public health affecting billions of people worldwide annually with no effective drugs available till now. Immunogenic and highly conserved properties of Non-Structural Protein 5(NS5) in DENV makes it a potent marker to identify DENV infection. DENV interfere in the innate immune signaling and thereby decreases antiviral responses and favors viral replication. Viral recognition by host pathogen recognition receptors facilitates binding of interferon (IFN) to the interferon receptors that further activates both the Signal Transducer and Activator of Transcription-2 (STAT-2) a factor producing an antiviral response. The most debilitating factor of DENV infection is emaciation of human immune system by DENV- NS5. NS5 counters the antiviral response by STAT2 degradation impeding the transcriptional activation of interferon stimulated genes through interferon stimulated response elements. The present study aims to identify inhibitors for NS5 Methyl Transferase (MTase) domain and to provide an insight into the mechanism of STAT2 degradation in the host infected with DENV. Virtual screening and molecular docking studies identified five potential inhibitors ZINC84154300, ZINC08762321, ZINC08762323, ZINC12659408 and ZINC12285470 with docking scores of -10.55, -10.53, -10.78, -11.28 and -10.78 kcal/mol respectively. To further investigate the stability of the complexes, we have used Molecular Dynamics Simulations (MD). Besides, the binding free energy of top 5 docked ligands were estimated through Molecular Mechanics Generalized Born and Surface Area Solvation (MM/GBSA) methods. This study also provides an insight on the mechanism of immunological processes involved in alleviating the antiviral immune response and computational identification of potent inhibitors for viral NS5 protein.