Pharmacoinformatics-based identification of chemically active molecules against Ebola virus

Quantum biochemical analysis of the binding interactions between a potential inhibitory drug and the Ebola viral glycoprotein

Ebola virus disease (EVD) causes outbreaks and epidemics in West Africa that persist until today. The envelope glycoprotein of Ebola virus (GP) consists of two subunits, GP1 and GP2, and plays a key role in anchoring or fusing the virus to the host cell in its active form on the virion surface. Toremifene (TOR) is a ligand that mainly acts as an estrogen receptor antagonist; however, a recent study showed a strong and efficient interaction with GP. In this context, we aimed to evaluate the energetic affinity features involved in the interaction between GP and toremifene by computer simulation techniques using the Molecular Fractionation Method with Conjugate Caps (MFCC) scheme and quantum-mechanical (QM) calculations, as well as missense mutations to assess protein stability. We identified ASP522, GLU100, TYR517, THR519, LEU186, LEU515 as the most attractive residues in the EBOV glycoprotein structure that form the binding pocket. We divided toremifene into three regions and evaluated that region i was more important than region iii and region ii for the formation of the TOR-GP1/GP2 complex, which might control the molecular remodeling process of TOR. The mutations that caused more destabilization were ARG134, LEU515, TYR517 and ARG559, while those that caused stabilization were GLU523 and ASP522. TYR517 is a critical residue for the binding of TOR, and is highly conserved among EBOV species. Our results may help to elucidate the mechanism of drug action on the GP protein of the Ebola virus and subsequently develop new pharmacological approaches against EVD.Communicated by Ramaswamy H. Sarma.

Computational Study on Potential Novel Anti-Ebola Virus Protein VP35 Natural Compounds

Ebola virus (EBOV) is one of the most lethal pathogens that can infect humans. The Ebola viral protein VP35 (EBOV VP35) inhibits host IFN-α/β production by interfering with host immune responses to viral invasion and is thus considered as a plausible drug target. The aim of this study was to identify potential novel lead compounds against EBOV VP35 using computational techniques in drug discovery. The 3D structure of the EBOV VP35 with PDB ID: 3FKE was used for molecular docking studies. An integrated library of 7675 African natural product was pre-filtered using ADMET risk, with a threshold of 7 and, as a result, 1470 ligands were obtained for the downstream molecular docking using AutoDock Vina, after an energy minimization of the protein via GROMACS. Five known inhibitors, namely, amodiaquine, chloroquine, gossypetin, taxifolin and EGCG were used as standard control compounds for this study. The area under the curve (AUC) value, evaluating the docking protocol obtained from the receiver operating characteristic (ROC) curve, generated was 0.72, which was considered to be acceptable. The four identified potential lead compounds of NANPDB4048, NANPDB2412, ZINC000095486250 and NANPDB2476 had binding affinities of −8.2, −8.2, −8.1 and −8.0 kcal/mol, respectively, and were predicted to possess desirable antiviral activity including the inhibition of RNA synthesis and membrane permeability, with the probable activity (Pa) being greater than the probable inactivity (Pi) values. The predicted anti-EBOV inhibition efficiency values (IC₅₀), found using a random forest classifier, ranged from 3.35 to 11.99 μM, while the Ki values ranged from 0.97 to 1.37 μM. The compounds NANPDB4048 and NANPDB2412 had the lowest binding energy of −8.2 kcal/mol, implying a higher binding affinity to EBOV VP35 which was greater than those of the known inhibitors. The compounds were predicted to possess a low toxicity risk and to possess reasonably good pharmacological profiles. Molecular dynamics (MD) simulations of the protein–ligand complexes, lasting 50 ns, and molecular mechanisms Poisson-Boltzmann surface area (MM-PBSA) calculations corroborated the binding affinities of the identified compounds and identified novel critical interacting residues. The antiviral potential of the molecules could be confirmed experimentally, while the scaffolds could be optimized for the design of future novel anti-EBOV chemotherapeutics.

Identification of bioactive compounds from Glycyrrhiza glabra as possible inhibitor of SARS-CoV-2 spike glycoprotein and non-structural protein-15: a pharmacoinformatics study

At present, the world is facing a pandemic named as COVID-19, caused by SARS-CoV-2. Traditional Chinese medicine has recommended the use of liquorice (Glycyrrhiza species) in the treatment of infections caused by SARS-CoV-2. Therefore, the present investigation was carried out to identify the active molecule from the liquorice against different protein targets of COVID-19 using an in-silico approach. The molecular docking simulation study of 20 compounds along with two standard antiviral drugs (Lopinavir and Rivabirin) was carried out with the help of Autodock vina software using two protein targets from COVID-19 i.e. spike glycoprotein (PDB ID: 6VSB) and Non-structural Protein-15 (Nsp15) endoribonuclease (PDB ID: 6W01). From the observed binding energy and the binding interactions, glyasperin A showed high affinity towards Nsp15 endoribonuclease with uridine specificity, while glycyrrhizic acid was found to be best suited for the binding pocket of spike glycoprotein and also prohibited the entry of the virus into the host cell. Further, the dynamic behavior of the best-docked molecules inside the spike glycoprotein and Nsp15 endoribonuclease were explored through all-atoms molecular dynamics (MD) simulation study. Several parameters from the MD simulation have substantiated the stability of protein-ligand stability. The binding free energy of both glyasperin A and glycyrrhizic acid was calculated from the entire MD simulation trajectory through the MM-PBSA approach and found to high binding affinity towards the respective protein receptor cavity. Thus, glyasperin A and glycyrrhizic acid could be considered as the best molecule from liquorice, which could find useful against COVID-19. Communicated by Ramaswamy H. Sarma.

Design, synthesis and evaluation of novel enzalutamide analogues as potential anticancer agents

The androgen receptor inhibitor, Enzalutamide, proved effective against castration resistance prostate cancer, has demonstrated clinical benefits and increased survival rate in men. However, AR mutation (F876L) converts Enzalutamide from antagonist to agonist indicating a rapid evolution of resistance. Hence, our goal is to overcome this resistance mechanism by designing and developing novel Enzalutamide analogues. We designed a dataset of Enzalutamide derivatives using Enzalutamide's shape and electrostatic features to match with pharmacophoric features essential for tight binding with the androgen receptor. Based on this design strategy ten novel derivatives were selected including 5,5-dimethyl-3-(6-substituted benzo[d]thia/oxazol-2-yl)-2-thioxo-1-(4-(trifluoromethyl)pyridin-2-yl)imidazolidin-4-one (6a-j) for synthesis. All the compounds were evaluated in-vitro on prostate cancer cell lines DU-145, LNCaP and PC3. Interestingly, two compounds 3-(6-hydroxybenzo[d]thiazol-2-yl)-5,5-dimethyl-2-thioxo-1-(4-(trifluoromethyl)pyridin-2-yl) imidazolidin-4-one (6c, IC50 - 18.26 to 20.31μM) and 3-(6-hydroxybenzo[d]oxazol-2-yl)-5,5-dimethyl -2-thioxo- 1- (4-(trifluoromethyl) pyridin-2-yl)imidazolidin-4-one (6h, IC50 - 18.26 to 20.31μM) were successful with promising in-vitro antiproliferative activity against prostate cancer cell lines. The binding mechanism of potential androgen receptor inhibitors was further studied by molecular docking, molecular dynamics simulations and MM-GBSA binding free energy calculations and found in agreement with the in vitro studies. It provided strong theoretical support to our hypothesis.

Molecular dynamics and combined docking studies for the identification of Zaire ebola virus inhibitors

Ebola virus (EBOV) is a lethal human pathogen with a risk of global spread of its zoonotic infections, and Ebolavirus Zaire specifically has the highest fatality rate amongst other species. There is a need for continuous effort towards having therapies, as a single licensed treatment to neutralize the EBOV is yet to come into reality. This present study virtually screened the MCULE database containing almost 36 million compounds against the structure of a Zaire Ebola viral protein (VP) 35 and a consensus scoring of both MCULE and CLCDDW docking programs remarked five compounds as potential hits. These compounds, with binding energies ranging from -7.9 to -8.9 kcal/mol, were assessed for predictions of their physicochemical and bioactivity properties, as well as absorption, distribution, metabolism, excretion, and toxicity (ADMET) criteria. The results of the 50 ns molecular dynamics simulations showed the presence of dynamic stability between ligand and protein complexes, and the structures remained significantly unchanged at the ligand-binding site throughout the simulation period. Both docking analysis and molecular dynamics simulation studies suggested strong binding affinity towards the receptor cavity and these selected compounds as potential inhibitors against the Zaire Ebola VP 35. With respect to inhibition constant values, bioavailability radar and other physicochemical properties, compound A (MCULE-1018045960-0-1) appeared to be the most promising hit compound. However, the ligand efficiency and ligand efficiency scale need improvement during optimization, and also validation via in vitro and in vivo studies are necessary to finally make a lead compound in treating Ebola virus diseases. Communicated by Ramaswamy H. Sarma.