The coronavirus pandemic is caused by intense acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Identifying the atomic structure of this virus can lead to the treatment of related diseases in medical cases. In the current computational study, the atomic evolution of the coronavirus in an aqueous environment using the Molecular Dynamics (MD) approach is explained. The virus behaviors by reporting the physical attributes such as total energy, temperature, potential energy, interaction energy, volume, entropy, and radius of gyration of the modeled virus are reported. The MD results indicated the atomic stability of the simulated virus significantly reduced after 25.33 ns. Furthermore, the volume of simulated virus changes from 182397 Å3 to 372589 Å3 after t = 30 ns. This result shows the atomic interaction between various atoms in coronavirus structure decreases in the vicinity of H2O molecules. Numerically, the interaction energy between virus and aqueous environment converges to -12387 eV and -251 eV values in the initial and final time steps of the MD study procedure, respectively.
The emergence outbreak caused by a novel severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has received significant attention on the global risks. Due to itscrucial role in viral replication, the main protease 3CLpro is an important target for drug discovery and development to combat COVID-19. In this work, the structural and dynamic behaviors as well as binding efficiency of the four peptidomimetic inhibitors (N3, 11a, 13b, and 14b) recently co-crystalized with SARS-CoV-2 3CLpro were studied and compared using all-atom molecular dynamics (MD) simulations and solvated interaction energy-based binding free energy calculations. The per-residue decomposition free energy results suggested that the key residues involved in inhibitors binding were H41, M49, L141-C145, H163-E166, P168, and Q189-T190 in the domains I and II. The van der Waals interaction yielded the main energy contribution stabilizing all the focused inhibitors. Besides, their hydrogen bond formations with F140, G143, C145, H164, E166, and Q189 residues in the substrate-binding pocket were also essential for strengthening the molecular complexation. The predicted binding affinity of the four peptidomimetic inhibitors agreed with the reported experimental data, and the 13b showed the most efficient binding to SARS-CoV-2 3CLpro. From rational drug design strategies based on 13b, the polar moieties (e.g., benzamide) and the bulky N-terminal protecting groups (e.g., thiazole) should be introduced to P1' and P4 sites in order to enhance H-bonds and hydrophobic interactions, respectively. We hope that the obtained structural and energetic information could be beneficial for developing novel SARS-CoV-2 3CLpro inhibitors with higher inhibitory potency to combat COVID-19.