Unraveling the binding mechanism of the active form of Remdesivir to RdRp of SARS-CoV-2 and designing new potential analogues: Insights from molecular dynamics simulations.

The binding of the active form of Remdesivir (RTP) to RNA-dependent RNA Polymerase (RdRp) of SARS-CoV-2 was studied using molecular dynamics simulation. The RTP maintained the interactions observed in the experimental cryo-EM structure. Next, we designed new analogues of RTP, which not only binds to the RNA primer strand in a similar pose as that of RTP, but also binds more strongly than RTP does as predicted by MM-PBSA binding energy. This suggest that these analogues might be able to covalently link to the primer strand as RTP, but their 3' modification would terminate the primer strand growth.

The Search for Cyclooxygenase-2 (COX-2) Inhibitors for the Treatment of Inflammation Disease: An in-silico Study.

Purpose: Cyclooxygenase (COX-2) has been validated as a molecular target for treating inflammatory diseases. The present work was performed to identify potential COX-2 inhibitors by employing pharmacophore modeling. Methods: The pharmacophore features consisted of seven features, ie, three hydrophobic, one negative ion, and three hydrogen bond acceptors, which were developed based on the structure of COX-2 inhibitor, (R)-naproxen. Results: The pharmacophore model was validated with a Goodness of Hit (GH score) of 0.754 and the values of AUC100% 0.51. Screening against the ZINC database retrieved 1675 hits, while the molecular docking procedure identified four best hit molecules in term of binding orientation and binding energies, ie, Lig_1805/ZINC103584272 (E = -11.03 kcal/mol), Lig_553/ZINC408573132 (E = -10.92 kcal/mol), Lig_680/ZINC103584263 (E = -10.90 kcal/mol), and Lig_2006/ZINC19324645 (E = -10.62 kcal/mol). Conclusion: The interactions of the four hits occurred in the binding site as (R)-naproxen did, and interestingly, their binding affinities were stronger than (R)-naproxen, implying their potential as COX-2 inhibitors.

Novel inhibitors to ADP ribose phosphatase of SARS-CoV-2 identified by structure-based high throughput virtual screening and molecular dynamics simulations.

The outbreak of a new coronavirus (SARS-CoV-2) was first identified in Wuhan, People's Republic of China, in 2019, which has led to a severe, life-threatening form of pneumonia (COVID-19). Research scientists all around the world have been trying to find small molecule drugs to treat COVID-19. In the present study, a conserved macrodomain, ADP Ribose phosphatase (ADRP), of a critical non-structural protein (Nsp3) in all coronaviruses was probed using large-scale Molecular Dynamics (MD) simulations to identify novel inhibitors. In our virtual screening workflow, the recently-solved X-ray complex structure, 6W6Y, with a substrate-mimics was used to screen 17 million ZINC15 compounds using drug property filters and Glide docking scores. The top twenty output compounds each underwent 200 ns MD simulations (i.e. 20 × 200 ns) to validate their individual stability as potential inhibitors. Eight out of the twenty compounds showed stable binding modes in the MD simulations, as well as favorable drug properties from our predctions. Therefore, our computational data suggest that the resulting top eight out of twenty compounds could potentially be novel inhibitors to ADRP of SARS-CoV-2.

GC-MS, LC-MS/MS, Docking and Molecular Dynamics Approaches to Identify Potential SARS-CoV-2 3-Chymotrypsin-Like Protease Inhibitors from Zingiber officinale Roscoe.

This study aims to identify and isolate the secondary metabolites of Zingiber officinale using GC-MS, preparative TLC, and LC-MS/MS methods, to evaluate the inhibitory potency on SARS-CoV-2 3 chymotrypsin-like protease enzyme, as well as to study the molecular interaction and stability by using docking and molecular dynamics simulations. GC-MS analysis suggested for the isolation of terpenoids compounds as major compounds on methanol extract of pseudostems and rhizomes. Isolation and LC-MS/MS analysis identified 5-hydro-7, 8, 2'-trimethoxyflavanone (9), (E)-hexadecyl-ferulate (1), isocyperol (2), N-isobutyl-(2E,4E)-octadecadienamide (3), and nootkatone (4) from the rhizome extract, as well as from the leaves extract with the absence of 9. Three known steroid compounds, i.e., spinasterone (7), spinasterol (8), and 24-methylcholesta-7-en-3ß-on (6), were further identified from the pseudostem extract. Molecular docking showed that steroids compounds 7, 8, and 6 have lower predictive binding energies (MMGBSA) than other metabolites with binding energy of -87.91, -78.11, and -68.80 kcal/mole, respectively. Further characterization on the single isolated compound by NMR showed that 6 was identified and possessed 75% inhibitory activity on SARS-CoV-2 3CL protease enzyme that was slightly different with the positive control GC376 (77%). MD simulations showed the complex stability with compound 6 during 100 ns simulation time.

Mechanistic insight on the remdesivir binding to RNA-Dependent RNA polymerase (RdRp) of SARS-cov-2.

The RNA-dependent RNA polymerase (RdRp) is a key enzyme which regulates the viral replication of SARS-CoV-2. Remdesivir (RDV) is clinically used drug which targets RdRp, however its mechanism of action remains elusive. This study aims to find out the binding dynamics of active Remdesivir-triphosphate (RDV-TP) to RdRp by means of molecular dynamics (MD) simulation. We built a homology model of RdRp along with RNA and manganese ion using RdRp hepatitis C virus and recent SARS-CoV-2 structures. We determined that the model was stable during the 500 ns MD simulations. We then employed the model to study the binding of RDV-TP to RdRp during three independent 500 ns MD simulations. It was revealed that the interactions of protein and template-primer RNA were dominated by salt bridge interactions with phosphate groups of RNA, while interactions with base pairs of template-primer RNA were minimal. The binding of RDV-TP showed that the position of phosphate groups was at the entry of the NTP channel and it was stabilized by the interactions with K551, R553, and K621, while the adenosine group on RDV-TP was pairing with U2 of the template strand. The manganese ion was located close to D618, D760, and D761, and helps in stabilization of the phosphate groups of RDV-TP. Further we identified three hits from the natural product database that pose similar to RDV-TP while having lower binding energies than that of RDV-TP, and that SN00359915 had binding free energy about three times lower than that of RDV-TP.

Pharmacophore-based virtual screening for identifying ß5 subunit inhibitor of 20S proteasome.

Proteasomal system plays an important role in maintaining cell homeostatis. Overexpression of proteasomes leads to several major diseases, such as cancer and autoimmune disorder. The ß5 subunit of proteasome is a crucial active site in proteolysis, and targeting proteasome ß5 subunit is essential for proteasome inhibition. In the present study, a pharmacophore-based virtual screening and molecular docking were employed to identify ligands as inhibitors of ß5 subunit of proteasome. The pharmacophore features were built with one hydrogen bond donor, two hydrogen bond acceptors, and one hydrophobic feature using native ligand of proteasome (HU10), which was then used to screen ZINC database using ZINCPharmer. The retrieved virtual hits were subjected to molecular docking analysis using iDock. The best six hits were subjected to molecular dynamics (MD) simulation and each complex was stable during 40 ns MD simulation as indicated by root-mean-square-deviation (RMSD) and root-mean-square-fluctuation (RMSF) values. The current study identifies 5 best hits having better binding potentials than HU10 as predicted by molecular mechanics Poisson-Boltzmann Surface Area (MM-PBSA) method, i.e. Lig1540/ZINC33356240, Lig1546/ZINC33356235, Lig1522/ZINC20854878, Lig980/ZINC12391945, and Lig1119/ZINC19865241, which can be used in the development of new proteasome inhibitors.

Molecular modeling on porphyrin derivatives as ß5 subunit inhibitor of 20S proteasome.

The ubiquitin-proteasome system plays an important role in protein quality control. Currently, inhibition of the proteasome has been validated as a promising approach in anticancer therapy. The 20S core particle of the proteasome harbors ß5 subunit which is a crucial active site in proteolysis. Targeting proteasome ß5 subunit which is responsible for the chymotrypsin-like activity of small molecules has been regarded as an important way for achieving therapeutics target. In the present study, a series of porphyrin derivatives bearing either pyridine or pyrazole rings as meso-substituents were designed and evaluated as an inhibitor for the ß5 subunit of the proteasome by employing molecular docking and dynamics simulations. The molecular docking was performed with the help of AutoDock 4.2, while molecular dynamics simulation was done using AMBER 14. All compounds bound to the proteasome with similar binding modes, and each porphyrin-proteasome complex was stable during 30 ns MD simulation as indicated by root-mean-square-deviation (RMSD) value. An analysis on protein residue fluctuation of porphyrin binding demonstrates that in all complexes, porphyrin binding produces minor fluctuation on amino acid residues. The molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) free energy calculation shows that the binding affinities of mono-H2PyP, bis-H2PzP, and tetra-H2PyP were comparable with that of the potential inhibitor, HU10. It is noted that the electrostatic interaction increases with the number of meso-substituents, which was favourable for porphyrin binding. The present study shows that both electrostatic and van der Waals interaction are the main force which controls the interaction of porphyrin compounds with the proteasome.

Molecular modeling of cationic porphyrin-anthraquinone hybrids as DNA topoisomerase IIß inhibitors.

Human DNA Topoisomerase II has been regarded as a promising target in anticancer drug discovery. In the present study, we designed six porphyrin-anthraquinone hybrids bearing pyrazole or pyridine group as meso substituents and evaluated their potentials as DNA Topoisomerase IIß inhibitor. First, we investigated the binding orientation of porphyrin hybrids into DNA topoisomerase IIß employing AutoDock 4.2 and then performed 20-ns molecular dynamics simulations to see the dynamic stability of each porphyrin-Topo IIß complex using Amber 14. We found that the binding of porphyrin hybrids occured through intercalation and groove binding mode in addition interaction with the amino acid residues constituting the active cavity of Topo IIß. Each porphyrin-Topo IIß complex was stabilized during 20-ns dynamics simulations. The MM-PBSA free energy calculation shows that the binding affinities of porphyrin hybrids were modified with the number of meso substituent. Interestingly, the affinity of all porphyrin hybrids to Topo IIß was stronger than that of native ligand (EVP), indicating the potential of the designed porphyrin to be considered in experimental research.

In silico study of porphyrin-anthraquinone hybrids as CDK2 inhibitor.

Cyclin-Dependent Kinases (CDKs) are known to play crucial roles in controlling cell cycle progression of eukaryotic cell and inhibition of their activity has long been considered as potential strategy in anti-cancer drug research. In the present work, a series of porphyrin-anthraquinone hybrids bearing meso-substituents, i.e. either pyridine or pyrazole rings were designed and computationally evaluated for their Cyclin Dependent Kinase-2 (CDK2) inhibitory activity using molecular docking, molecular dynamics simulation, and binding free energy calculation. The molecular docking simulation revealed that all six porphyrin hybrids were able to bind to ATP-binding site of CDK2 and interacted with key residues constituted the active cavity of CDK2, while molecular dynamics simulation indicated that all porphyrins bound to CDK2 were stable for 6ns. The binding free energies predicted by MM-PBSA method showed that most compounds exhibited higher affinity than that of native ligand (4-anilinoquinazoline, DTQ) and the affinity of mono-H2PyP-AQ was about three times better than that of DTQ, indicating its potential to be advanced as a new CDK2 inhibitor.

Molecular docking and dynamics simulations on the interaction of cationic porphyrin-anthraquinone hybrids with DNA G-quadruplexes.

A series of cationic porphyrin-anthraquinone hybrids bearing either pyridine, imidazole, or pyrazole rings at the meso-positions have been investigated for their interaction with DNA G-quadruplexes by employing molecular docking and molecular dynamics simulations. Three types of DNA G-quadruplexes were utilized, which comprise parallel, antiparallel, and mixed hybrid topologies. The porphyrin hybrids have a preference to bind with parallel and mixed hybrid structures compared to the antiparallel structure. This preference arises from the end stacking of porphyrin moiety following G-stem and loop binding of anthraquinone tail, which is not found in the antiparallel due to the presence of diagonal and lateral loops that crowd the G-quartet. The binding to the antiparallel, instead, occurred with poorer affinity through both the loop and wide groove. All sites of porphyrin binding were confirmed by 6 ns molecular dynamics simulation, as well as by the negative value of the total binding free energies that were calculated using the MMPBSA method. Free energy analysis shows that the favorable contribution came from the electrostatic term, which supposedly originated from the interaction of either cationic pyridinium, pyrazole, or imidazole groups and the anionic phosphate backbone, and also from the van der Waals energy, which primarily contributed through end stacking interaction.

The binding modes of cationic porphyrin-anthraquinone hybrids to DNA duplexes: in silico study.

Cationic porphyrin-anthraquinone hybrids bearing peripheral substituents, either pyridine, imidazole, or pyrazole rings have been investigated for their binding mode to DNA duplexes. The four kinds of DNA duplexes were used, which represent intercalation and groove binding modes. AutoDock 4.2 was used to dock nine hybrid compounds to four DNA duplexes, while monitoring of conformational changes of four best hybrid-DNA complexes during 2 ns was performed by Amber9 molecular dynamics package. The binding energy calculation of best four complexes was then carried out using MMPBSA method. The hybrid compounds interacted to DNA duplexes through intercalation and groove binding modes. The minor groove binding of DNA was energetically preferred by cationic porphyrin hybrids due to favorable electrostatic and van der Waals interactions. Both electrostatic and van der Waals contributions were able to distinguish the binding mode of porphyrin hybrid to DNA duplexes.