Structural Differences In 3C-like protease (Mpro) From SARS-CoV and SARS-CoV-2: Molecular Insights For Drug Repurposing Against COVID-19 Revealed by Molecular Dynamics Simulations

Abstract

A recent fatal outbreak of novel coronavirus SARS-CoV-2, identified preliminary as a causative agent for series of unusual pneumonia cases in Wuhan city, China has infected more than 20 million individuals with more than 4 million mortalities. Since, the infection crossed geographical barriers, the WHO permanently named the causing disease as COVID-2019 by declaring it a pandemic situation. SARS-CoV-2 is an enveloped single-stranded RNA virus causing a wide range of pathological conditions from common cold symptoms to pneumonia and fatal severe respiratory syndrome. Genome sequencing of SARS-CoV-2 has revealed 96% identity to the bat coronavirus and 79.6% sequence identity to the previous SARS-CoV. The main protease (known as 3C-like proteinase/ Mpro) plays a vital role during the infection with the processing of replicase polyprotein thus offering an attractive target for therapeutic interventions. SARS-CoV and SARS-CoV-2 Mpro shares 97% sequence identity, with 12 variable residues but none of them are present in the catalytic and substrate binding site. With the high level of sequence and structural similarity and absence of any drug/vaccine against SARS-CoV-2, drug repurposing against Mpro is an effective strategy to combat COVID-19. Here, we report a detailed comparison of SARS-CoV-2 Mpro with SARS-CoV Mpro using molecular dynamics simulations to assess the impact of 12 divergent residues on the molecular microenvironment of Mpro. Structural comparison and analysis are made on how these variable residues affect the intra-molecular interactions between key residues in the monomer and biologically active dimer form of Mpro. The present MD simulations study concluded the change in microenvironment of active-site residues at the entrance (T25, T26, M49 and Q189), near the catalytic region (F140, H163, H164, M165 and H172) and other residues in substrate binding site (V35T, N65S, K88R and N180K) due to 12 mutation incorporated in the SARS-CoV-2 Mpro. It is also evident that SARS-CoV-2 dimer is more stable and less flexible state compared to monomer which may be due to these variable residues, mainly F140, E166 and H172 which are involved in dimerization. This also warrants a need for inhibitor design considering the more stable dimer form. The mutation accumulated in SARS-CoV-2 Mpro indirectly reconfigures the key molecular networks around the active site conferring a potential change in SARS-CoV-2, thus posing a challenge in drug repurposing SARS drugs for COVID-19. The new networks and changes in the microenvironment identified by our work might guide attempts needed for repurposing and identification of new Mpro inhibitors.

Development of Masitinib Derivatives with Enhanced Mpro Ligand Efficiency and Reduced Cytotoxicity

“…It was recently shown that HCoV-OC43 is a good biosafety-level 2 (BSL-2) model system to study antiviral drugs against SARS-CoV-2 infection [ 6 ]. The two homologous proteins share roughly 50% of sequence identity [ 16 ], with the catalytic dyad of Cys-His conserved throughout the two sequences, as they are essential for enzyme activity. Figure S1 shows a comparison between both structures, where the high structural identity is evident (RMSD for the active site main chain is lower than 1.7 Å).…”