The recent occurrence of COVID-19 pandemic causing fever, respiratory illness, and pneumonia, has drawn the attention across the globe and has been declared a global health emergency. Among humans, coronaviruses were not believed to be pathogenic until the severe acute respiratory syndrome (SARS) epidemic occurred in late 2002. Ten years later in 2012, a new coronavirus epidemic named as Middle East respiratory syndrome (MERS) emerged in Asian and Middle Eastern countries. Several studies have been published to understand the cell entry, replication, and protein structures of virus. The genomic organization of SARS and MERS is similar, comprising an enveloped, positive-sense, single-stranded RNA genome that translates into various structural and non-structural proteins (nsps). After viral entry into the host cell, virus releases its RNA into the cytoplasm where it is translated into polyproteins pp1a and pp1ab. These polyproteins are cleaved by proteases that produce 16 nsps. The viral genome is translated, replicated, and assembled in the cytoplasm which is then released via exocytosis.
To combat with the coronavirus diseases, it is important to block viral infection or replication. Therefore, two potential target proteins were selected for the drug development. The first target is an accessory protein 4a of MERS-CoV which has been proved through in vitro and in vivo studies that it is involved in inhibition of IFN production which influences the viral pathogenicity. 3D structure of protein 4a was modeled followed by structural validation and computational alanine mutagenesis. Residues N8, K27, W45, K63, and K67 were selected based on the dStability and dAffinity values. The mutant and wild type complexes were simulated for 200 ns each. The MD simulation results revealed that K27, W45, K63, and K67 are key residues that play crucial role in RNA binding. Thus, suggesting that blocking these residual interactions could destabilize the 4a-RNA complex. Furthermore, considering these residues a part of binding site, virtual screening of
antiviral as well as druglike and leadlike compounds from ZINC database is conducted using the most stable confirmation extracted from MD simulations. The ten compounds from both datasets were selected based on binding affinity scores and their binding pose. The interactions were analyzed in detail along with the free binding energy calculations. Most of the selected compounds possess antiviral properties. These compounds can be tested through in vitro and in vivo experiments.
The 3CL protease is another potential target due to its proteolytic activity that cleaves the polyprotein. Many drugs have been repositioned, and a few ligands are cocrystalized with SARS-CoV-2 3CL protease. The flexibility of the catalytic site leads to the conformational changes that influences the ligand binding. Thus, it is necessary to probe the flexibility and rigidity of protein’s active site. 100 ns MD simulations were performed for apo form of 3CL protease and its complexes with different cocrystalized ligands. The trajectory analysis and binding free energy calculation revealed that the ligand that can break the catalytic dyad bond in the active site would be a potential candidate to block the proteolytic activity. The alpha-keto group present in the inhibitor of 6LU7 make strong bonds with both residues Cys and His of catalytic dyad. Therefore, addition of such chemical groups that break the bond between Cys and His would be important to block the functional 3CL protease.
Conclusively, both target proteins 4a and 3CL protease are important for viral pathogenicity and replication, respectively. Therefore, understanding of structure dynamics and insights into the binding site is crucial for development of an efficient drug molecule. 3CL protease is structurally and functionally conserved protein among all three coronaviruses. Thus, targeting this protein would be a hope to combat all coronavirus diseases.