The H163A Mutation Unravels an Oxidized Conformation of the SARS-CoV-2 Main Protease and Opens a New Avenue for Anti-Viral Therapeutic Design

The main protease of SARS-CoV-2 (Mpro) is an important target for developing COVID-19 therapeutics. Recent work has highlighted Mpros susceptibility to undergo redox-associated conformational changes in response to cellular and immune-system-induced oxidation. Despite structural evidence indicating large-scale rearrangements upon oxidation, the mechanisms of conformational change and its functional consequences are poorly understood. Here, we present the crystal structure of a new Mpro point mutant (H163A) that shows an oxidized conformation with the catalytic cysteine in a disulfide bond. We hypothesize that Mpro adopts this conformation under oxidative stress to protect against over-oxidation. Our metadynamics simulations illustrated a potential mechanism by which H163 modulates this transition and suggest that this equilibrium exists in the wild-type enzyme. We show that other point mutations can also significantly shift the equilibrium towards this state by altering conformational free energies. New therapeutic strategies against SARS-CoV-2 can be explored by understanding how H163 modulates this equilibrium.

Molecular dynamics and in silico mutagenesis on the reversible inhibitor-bound SARS-CoV-2 Main Protease complexes reveal the role of lateral pocket in enhancing the ligand affinity

The 2019 novel coronavirus pandemic caused by SARS-CoV-2 remains a serious health threat to humans and a number of countries are already in the middle of the second wave of infection. There is an urgent need to develop therapeutics against this deadly virus. Recent scientific evidences have suggested that the main protease (Mpro) enzyme in SARS-CoV-2 can be an ideal drug target due to its crucial role in the viral replication and transcription processes. Therefore, there are ongoing research efforts to identify drug candidates against SARS-CoV-2 Mpro that resulted in hundreds of X-ray crystal structures of ligand bound Mpro complexes in the protein data bank (PDB) that describe structural details of different chemotypes of fragments binding within different sites in Mpro. In this work, we perform rigorous molecular dynamics (MD) simulation of 62 reversible ligand-Mpro complexes in the PDB to gain mechanistic insights about their interactions at atomic level. Using a total of ~2.25 s long MD trajectories, we identified and characterized different pockets and their conformational dynamics in the apo Mpro structure. Later, using the published PDB structures, we analyzed the dynamic interactions and binding affinity of small ligands within those pockets. Our results identified the key residues that stabilize the ligands in the catalytic sites and other pockets in Mpro. Our analyses unraveled the role of a lateral pocket in the catalytic site in Mpro that is critical for enhancing the ligand binding to the enzyme. We also highlighted the important contribution from HIS163 in this lateral pocket towards ligand binding and affinity against Mpro through computational mutation analyses. Further, we revealed the effects of explicit water molecules and Mpro dimerization in the ligand association with the target. Thus, comprehensive molecular level insights gained from this work can be useful to identify or design potent small molecule inhibitors against SARS-CoV-2 Mpro.

A novel strategy for mechanism based computational drug discovery.

Glioma, the common brain tumor, which arises from the glial cells, offers worse prognosis and therapy than any other tumors. Despite the genetic and pathological diversities of malignant gliomas, common signaling pathways that drive cellular proliferation, survival, invasion and angiogenesis have been identified. Very often, various tyrosine kinase receptors are inappropriately activated in human brain tumors and contribute to tumor malignancy. During such tumourous states where multiple pathways are involved, a few of them are responsbile for cell differentiation, proliferation and anti-apoptosis. Computational simulation studies of normal EGFR signaling in glioma together with the mutant EGFR mediated signaling and the MAPK signaling in glioma were carried out. There were no significant cross talks observed between the mutant EGFR and the MAPK pathways and thus from the simulation results, we propose a novel concept of 'multiple-targeting' that combines EGFR and Ras targeted therapy thereby providing a better therapeutic value against glioma. Diallyl Disulfide (DADS) that has been commonly used for Ras inhibition in glioma was taken for analyses and the effect of inhibiting the EGFR downstream signaling protein with this DADS was analyzed using the simulation and docking studies.