Identification of repurposing therapeutics toward SARS-CoV-2 main protease by virtual screening.

SARS-CoV-2 causes the current global pandemic coronavirus disease 2019. Widely-available effective drugs could be a critical factor in halting the pandemic. The main protease (3CLpro) plays a vital role in viral replication; therefore, it is of great interest to find inhibitors for this enzyme. We applied the combination of virtual screening based on molecular docking derived from the crystal structure of the peptidomimetic inhibitors (N3, 13b, and 11a), and experimental verification revealed FDA-approved drugs that could inhibit the 3CLpro of SARS-CoV-2. Three drugs were selected using the binding energy criteria and subsequently performed the 3CLpro inhibition by enzyme-based assay. In addition, six common drugs were also chosen to study the 3CLpro inhibition. Among these compounds, lapatinib showed high efficiency of 3CLpro inhibition (IC50 value of 35 ± 1 µM and Ki of 23 ± 1 µM). The binding behavior of lapatinib against 3CLpro was elucidated by molecular dynamics simulations. This drug could well bind with 3CLpro residues in the five subsites S1', S1, S2, S3, and S4. Moreover, lapatinib's key chemical pharmacophore features toward SAR-CoV-2 3CLpro shared important HBD and HBA with potent peptidomimetic inhibitors. The rational design of lapatinib was subsequently carried out using the obtained results. Our discovery provides an effective repurposed drug and its newly designed analogs to inhibit SARS-CoV-2 3CLpro.

Computational study on peptidomimetic inhibitors against SARS-CoV-2 main protease.

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.

In Silico Screening of DNA Gyrase B Potent Flavonoids for the Treatment of Clostridium difficile Infection from PhytoHub Database

Abstract Clostridium difficile infection (CDI) is the most common hospital acquired diarrheal disease with its increasing incidence and mortality rate globally. DNA Gyrase B (GyrB) is a key component of DNA replication process across all bacterial genera; thus, this offers a potential target for the treatment of CDI. In the present study, several virtual screening approaches were employed to identify a novel C. difficile GyrB inhibitor. The 139 known metabolites were screened out from the 480 flavonoids in PhytoHub database. Molinspiration and PROTOX II servers were used to calculate the ADME properties and oral toxicity of the metabolites, whereas mutagenicity, tumorigenicity, irritant, and reproductive effect were predicted using DataWarrior program. The binding mode and the binding efficiency of the screened flavonoids against the GyrB were studied using FlexX docking program. From virtual screening of 139 metabolites, we found 25 flavonoids with no mutagenicity, tumorigenicity, irritant, and reproductive effect. Docking study suggested that flavonoids 1030 ((-)-epicatechin 3'-O-sulfate), 1032 ((-)-epicatechin 4'-O-sulfate), 1049 (3'-O-methyl-(-)-epicatechin 4-O-sulfate), 1051 (3'-O-methyl-(-)-epicatechin 7-O-sulfate), 1055 (4'-O-methyl-(-)-epicatechin 7-O-sulfate) and 1317 (quercetin sulfate) have significantly higher binding affinity than the known GyrB inhibitor novobiocin. The results from molecular dynamics simulation and free energy calculations based on solvated interaction energy suggested that (-)-epicatechin 3'-O-sulfate could be a potential drug candidate in the management of CDI.

Torsional flexibility of undecorated catechol diether compound as potent NNRTI targeting HIV-1 reverse transcriptase.

Conformational adaptation of non-nucleoside reverse transcriptase inhibitor (NNRTI) via torsional flexibility is found to be very significant for targeting human immunodeficiency virus type 1 reverse transcriptase (HIV-1 RT) mutants. Catechol diether derivative including flexible torsions is new potent NNRTI with picomolar activity. Moreover, this derivative also reveals the good solubility, low toxicity and potent inhibition for HIV-1 mutants. In this study, torsional flexibility of an undecorated catechol diether compound in the binding pocket of wild type and mutants (Y181C and K103N/Y181C) HIV-1 RT is investigated by using QM/MM calculations. From the results, the uracil ring is found to exhibit more flexibility in the NNIBP. On the contrary, potential energy surfaces show that high energy is encountered by changing of the corresponding torsion of the cyanovinyl aryl ring indicating the limitation for torsional flexibility. For pointing out the key interaction for the binding, the residual interaction energies are performed by means of QM calculations. Important attractive interactions through hydrogen bonds between the inhibitor and K102, K/N103, V106, and Y188 are observed. The catechol ring is proposed to be modified in order to strengthen interactions with surrounding amino acids. The results may help for the designing of new potent NNRTIs.

Mechanistic insights into the catalytic reaction of plant allene oxide synthase (pAOS) via QM and QM/MM calculations.

QM cluster and QM/MM protein models have been employed to understand aspects of the reaction mechanism of plant allene oxide synthase (pAOS). In this study we have investigated two reaction mechanisms for pAOS. The standard pAOS mechanism was contrasted with an alternative involving an additional active site molecule which has been shown to facilitate proton coupled electron transfer (PCET) in related systems. Firstly, we found that the results from QM/MM protein model are comparable with those from the QM cluster model, presumably due to the large active site used. Furthermore, the results from the QM cluster model show that the Fe(III) and Fe(IV) pathways for the standard mechanism have similar energetic and structural properties, indicating that the reaction mechanism may well proceed via both pathways. However, while the PCET process is facilitated by an additional active site bound water in other related families, in pAOS it is not, suggesting this type of process is not general to all closely related family members.

Insight into the reaction mechanism of cis,cis-muconate lactonizing enzymes: a DFT QM/MM study.

MLEs derived from mycobacterium smegmatis and seudomonas fluorescens share ∼76% identity and have a very similar arrangement of catalytic residues in their active site configuration. However, while they catalyze the conversion of cis,cis-muconate to the same achiral product, muconolactone, studies in deuterated solvent surprisingly show that the cyclo-isomerization proceeds with the formation of a chiral product. In this paper we discuss the application of DFT QM/MM calculations on both MLEs, to our knowledge the first reported in the literature on this protein. We investigate the proposal that the base involved in the catalytic reaction is the lysine residue found at the end of the 2(nd) strand given: (a) that the lysine residue at the end of the 6(th) strand is in an apparently equally effective position to catalyze reaction and (b) that the structural related epimerase in-fact achieve their stereo-specific outcomes by relying on either the base from the 2(nd) or 6(th) strand.