Development of a Novel Pharmacophore Model Guided by the Ensemble of Waters and Small Molecule Fragments Bound to SARS-CoV-2 Main Protease

Structural and thermodynamic properties of conserved water molecules in Mpro native: A combined approach by MD simulation and Grid Inhomogeneous Solvation Theory

The new viral strains of severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) are continuously rising, becoming more virulent, and transmissible. Therefore, the development of new antiviral drugs is essential. Due to its significant role in the viral life cycle of SARS-CoV-2, the main protease (Mpro) enzyme is a leading target for antiviral drug design. The Mpro monomer consists of domain DI, DII, and DI-DII interface. Twenty-one conserved water molecules (W4-W24) are occupied at these domains according to multiple crystal structure analyses. The crystal and MD structures reveal the presence of eight conserved water sites in domain DI, DII and remaining in the DI-DII interface. Grid-based inhomogeneous fluid solvation theory (GIST) was employed on MD structures of Mpro native to predict structural and thermodynamic properties of each conserved water site for focusing to identify the specific conserved water molecules that can easily be displaced by proposed ligands. Finally, MD water W13 is emerged as a promising candidate for water mimic drug design due to its low mean interaction energy, loose binding character with the protein, and its involvement in a water-mediated H-bond with catalytic His41 via the interaction Thr25(OG)---W13---W---His41(NE2). In this context, water occupancy, relative interaction energy, entropy, and topologies of W13 are thermodynamically acceptable for the water displacement method. Therefore, the strategic use of W13's geometrical position in the DI domain may be implemented for drug discovery against COVID disease by designing new ligands with appropriately oriented chemical groups to mimic its structural, electronic, and thermodynamic properties.

In silico identification and molecular dynamic simulations of derivatives of 6,6-dimethyl-3-azabicyclo[3.1.0]hexane-2-carboxamide against main protease 3CLpro of SARS-CoV-2 viral infection

CONTEXT

The unavailability of target-specific antiviral drugs for SARS-CoV-2 viral infection kindled the motivation to virtually design derivatives of 6,6-dimethyl-3-azabicyclo[3.1.0]hexane-2-carboxamide as potential antiviral inhibitors against the concerned virus. The molecular docking and molecular dynamic results revealed that the reported derivatives have a potential to act as antiviral drug against SARS-CoV-2. The reported hit compounds can be considered for in vitro and in vivo analyses.

METHODS

Fragment-based drug designing was used to model the derivatives. Furthermore, DFT simulations were carried out using B3LYP/6-311G** basis set. Docking simulations were performed by using a combination of empirical free energy force field with a Lamarckian genetic algorithm under AutoDock 4.2. By the application of AMBER14 force field and SPCE water model, molecular dynamic simulations and MM-PBSA were calculated for 100 ns.