Computational identification of potential lead molecules targeting rho receptor of Neisseria gonorrhoeae

Response regulator GacA and transcriptional activator RhlR proteins involved in biofilm formation of Pseudomonas aeruginosa are prospective targets for natural lead molecules: Computational modelling, molecular docking and dynamic simulation studies

This study aimed to screen putative drug targets associated with biofilm formation of multidrug-resistant Acinetobacter baumannii and Pseudomonas areugenosa and prioritize carbon nano-fullerene as potential lead molecule by structure-based virtual screening. Based on the functional role, 36 and 83 genes that are involved in biofilm formation of A. baumannii and P. areugenosa respectively were selected and metabolic network was computationally constructed. The genes that lack three-dimensional structures were predicted and validated. Carbon nano-fullerene selected as lead molecule and their drug-likeliness and pharmacokinetics properties were computationally predicted. The binding potential of carbon nano-fullerene toward selected drug targets was modeled and compared with the binding of conventional drugs, doripenem, and polymyxin-B with their usual targets. The stabilities of four best-docked complexes were confirmed by molecular dynamic (MD) simulation. This study suggested that selected genes demonstrated relevant interactions in the constructed metabolic pathways. Carbon fullerene exhibited significant binding abilities to most of the prioritized targets in comparison with the binding of last-resort antibiotics and their usual target. The four best ligand–receptor interactions predicted by molecular docking revealed that stability throughout MD simulation. Notably, carbon fullerene exhibited profound binding with outer membrane protein (OmpA) and ribonuclease-HII (rnhB) of A. baumannii and 2-heptyl-4(1H)-quinolone synthase (pqsBC) and chemotaxis protein (wspA) of P. aeruginosa. Thus, the current study suggested that carbon fullerene was probably used as potential lead molecules toward selected targets of A. baumannii and P. aeruginosa and the applied aspects probably scaled up to design promising lead molecules toward these pathogens. Communicated by Ramaswamy H. Sarma.

Lymphatic filariasis is a debilitating vector borne parasitic disease that infects human lymphatic system by nematode Brugia malayi. Currently available anti-filarial drugs are effective only on the larval stages of parasite. So far, no effective drugs are available for humans to treat filarial infections. In this regard, aspartate semialdehyde dehydrogenase (ASDase) in lysine biosynthetic pathway from Wolbachia endosymbiont Brugia malayi represents an attractive therapeutic target for the development of novel anti-filarial agents. In this present study, molecular modeling combined with molecular dynamics simulations and structure-based virtual screening were performed to identify potent lead molecules against ASDase. Based on Glide score, toxicity profile, binding affinity and mode of interactions with the ASDase, five potent lead molecules were selected. The molecular docking and dynamics results revealed that the amino acid residues Arg103, Asn133, Cys134, Gln161, Ser164, Lys218, Arg239, His246, and Asn321 plays a crucial role in effective binding of Top leads into the active site of ASDase. The stability of the ASDase-lead complexes was confirmed by running the 30 ns molecular dynamics simulations. The pharmacokinetic properties of the identified lead molecules are in the acceptable range. Furthermore, density functional theory and binding free energy calculations were performed to rank the lead molecules. Thus, the identified lead molecules can be used for the development of anti-filarial agents to combat the pathogenecity of Brugia malayi.

Mycolic acids are the key constituents of mycobacterial cell wall, which protect the bacteria from antibiotic susceptibility, helping to subvert and escape from the host immune system. Thus, the enzymes involved in regulating and biosynthesis of mycolic acids can be explored as potential drug targets to kill Mycobacterium tuberculosis (Mtb). Herein, Kyoto Encyclopedia of Genes and Genomes is used to understand the fatty acid metabolism signaling pathway and integrative computational approach to identify the novel lead molecules against the mtFabH (β-ketoacyl-acyl carrier protein synthase III), the key regulatory enzyme of the mycolic acid pathway. The structure-based virtual screening of antimycobacterial compounds from ChEMBL library against mtFabH results in the selection of 10 lead molecules. Molecular binding and drug-likeness properties of lead molecules compared with mtFabH inhibitor suggest that only two compounds, ChEMBL414848 (C1) and ChEMBL363794 (C2), may be explored as potential lead molecules. However, the spatial stability and binding free energy estimation of thiolactomycin (TLM) and compounds C1 and C2 with mtFabH using molecular dynamics simulation, followed by molecular mechanics Poisson–Boltzmann surface area (MM/PBSA) indicate the better activity of C2 (ΔG = −14.18 kcal/mol) as compared with TLM (ΔG = −9.21 kcal/mol) and C1 (ΔG = −13.50 kcal/mol). Thus, compound C1 may be explored as promising drug candidate for the structure-based drug designing of mtFabH inhibitors in the therapy of Mtb.

Hormone based birth control often causes various side effects. A recent study revealed that temporary infertility without changing hormone levels can be attained by inhibiting Katanin p60 ATPase-containing subunit A-like 1 protein (KATNAL1) which is critical for sperm maturation in the testes. This study aimed at attaining the most energetically stable three dimensional (3D) structure of KATNAL1 protein using comparative modeling followed by screening of a ligand library of known natural spermicidal compounds for their binding affinity with KATNAL1. This in turn may inhibit the development of mature sperm in the seminiferous epithelium. A series of computational techniques were used for building the 3D structure of KATNAL1 which was further optimized by molecular dynamics (MD) simulation. For revealing the ATP binding mode of KATNAL1, docking study was carried out using the optimized model obtained from the MD simulation. The docking study was also employed to test the binding efficiency of the ligand library. Molecular docking study confirmed the ATP binding of KATNAL1 with various hydrophobic and hydrogen bond interactions. Binding efficiency of the ligand library suggested that calotropin, a cardenolide of Calotropis procera showed the highest binding efficiency against the target protein without toxicity. MD simulation of the docked complex validated the results of the docking study. This study revealed the ATP binding mode of KATNAL1 and identified calotropin as a potential lead molecule against it showing high binding efficiency with good bioavailability and no mutagenicity. Further in vitro and in vivo bioassay of calotropin could facilitate the development of novel non-hormonal male-specific contraceptive in near future.

Background Triple-negative breast cancers are defined as tumors that lack the expression of the estrogen receptor (ER), progesterone receptor (PR), and human epidermal growth factor receptor 2 (HER2). It exhibits unique clinical and pathological features, is highly aggressive, and has a relatively poor prognosis and poor clinical outcome. Objective To identify a novel drug target protein against triple-negative breast cancer (TNBC) and potential phytochemical lead molecules against novel drug targets. Methods In this study, we retrieved TNBC samples from NGS and microarray datasets in the Gene Expression Omnibus database and employed a combination of differential gene expression studies, protein-protein interaction analysis, and network topology investigation to identify the target protein. Using molecular docking and molecular dynamics simulation studies, followed by Molecular Mechanics with Generalised Born Surface Area solvation, a potential lead molecule was identified. Result The androgen receptor (AR) was found to be the target protein, and 2-hydroxynaringenin was discovered to be a possible phytochemical lead molecule to combat TNBC. Upregulated genes with LogFC > 1.25 and P-value < 0.05 from the TNBC gene expression dataset were given to STRING tool to investigate the network topology, and androgen receptor (AR) was found to be an appropriate hub gene in the protein-protein interaction network. Phytochemicals that inhibit breast cancer were retrieved from the PubChem database and virtual screening was performed using PyRx against the AR protein. Based on Lipinski’s rule and ADMET properties, molecular interaction studies were analyzed using induced fit docking, wherein significant binding interactions were displayed by 2-hydroxynaringenin. Molecular dynamics studies and MM-GBSA of AR and the 2-hydroxynaringenin complex revealed strong and stable interactions. Conclusion AR was identified as a hub protein that is highly expressed in breast cancer and 2-hydroxynaringenin efficacy of counter TNBC needs to be investigated further in vitro and in vivo.

Potential inhibitors for peroxiredoxin 6 of W. bancrofti: A combined study of modelling, structure-based drug design and MD simulation

Developing novel drug molecules against HIV is a scientific quest necessitated by development of drug resistance against used drugs. We report comparative results of molecular dynamics simulation studies on 11 structural analogues of Saquinavir (SQV) – against HIV-protease that were earlier examined for pharmacodynamic and pharmacokinetic properties. We reported analogues S1, S5 and S8 to qualify the ADMET criterion and may be considered as potential lead molecules. In this study the designed molecules were successively docked with native HIV-protease at AutoDock. Docking scores established relative goodness of the 11 analogues against the benchmark for Saquinavir. The docked complexes were subjected to molecular dynamics simulation studies using GROMACS 4.6.2. Four parameters viz. H-bonding, RMSD, Binding energy and Protein–Ligand Distance were used for comparative analyses of the analogues relative to Saquinavir. The comparison and analysis of the results are indicative that analogues S8, S9 and S1 are promising candidates among all the analogues studied. From our earlier work and present study it is evident that among the three S8 and S1 qualify the ADMET criterion and between S1 and S8, the analogue S8 shows more target efficacy and specificity over S1 and have better molecular dynamics simulation results. Thus, of the 11 de novo Saquinavir analogues, the S8 appears to be the most promising candidate as lead molecule for HIV-protease inhibitor and is best suited for testing under biological system. Further validation of the proposed lead molecules through wet lab studies involving antiviral assays however is required. Communicated by Ramaswamy H. Sarma

The prevalence of multi-drug resistance S. aureus is one of the most challenging tasks for the treatment of nosocomial infections. Proteins and enzymes of peptidoglycan biosynthesis pathway are one among the well-studied targets, but many of the enzymes are unexplored as targets. MurE is one such enzyme featured to be a promising target. As MurE plays an important role in ligating the L-lys to stem peptide at third position that is crucial for peptidoglycan synthesis. To screen the potential MurE inhibitor by in silico approach and evaluate the best potential lead molecule by in vitro methods. In the current study, we have employed structure based virtual screening targeting the active site of MurE, followed by Molecular dynamics and in vitro studies. Virtual screening resulted in successful screening of potential lead molecule ((2R)-2-[[1-[(2R)- 2-(benzyloxycarbonylamino) propanoyl] piperidine-4-carbonyl]amino]-5-guanidino-pentan). The molecular dynamics of the MurE and Lead molecule complex emphasizes that lead molecule has shown stable interactions with active site residues Asp 406 and with Glu 460. In vitro studies demonstrate that the lead molecule shows antibacterial activity close to standard antibiotic Vancomycin and higher than that of Ampicillin, Streptomycin and Rifampicin. The MIC of lead molecule at 50μg/mL was observed to be 3.75 μg/mL, MBC being bactericidal with value of 6.25 μg/mL, cytotoxicity showing 34.44% and IC50 of 40.06μg/mL. These results suggest ((2R)-2-[[1-[(2R)-2-(benzyloxycarbonylamino) propanoyl] piperidine-4-carbonyl]amino]-5-guanidino-pentan) as a promising lead molecule for developing a MurE inhibitor against treatment of S. aureus infections.

Identification of potential inhibitors for oncogenic target of dihydroorotate dehydrogenase using in silico approaches

Triple-negative breast cancers (TNBC) are defined as tumors that lack the expression of the estrogen receptor (ER), progesterone receptor (PR), and human epidermal growth factor receptor 2 (HER2). It exhibits unique clinical and pathological features, demonstrates high aggressiveness, and has a relatively poor prognosis and clinical outcome. To identify a novel drug target protein against TNBC and potential phytochemical lead molecules against the identified target. In this study, we retrieved TNBC samples from NGS and microarray datasets in the Gene Expression Omnibus database. We employed a combination of differential gene expression studies, protein-protein interaction analysis, and network topology investigation to identify the target protein. Additionally, the molecular docking and molecular dynamics (MD) simulation studies followed by Molecular Mechanics with Generalised Born Surface Area salvation was used to identify potential lead molecule. The upregulated genes with LogFC > 1.25 and P-value < 0.05 from the TNBC gene expression dataset were identified. Androgen receptor (AR) was found to be an appropriate hub target in the protein-protein interaction network. Phytochemicals that inhibit breast cancer target were retrieved from the PubChem database and virtual screening was performed using PyRx against the AR protein. Thereby, the AR was found to be the target protein and 2-hydroxynaringenin was discovered to be a possible phytochemical lead molecule for combating TNBC. Moreover, the AR and the 2-hydroxynaringenin complex showed structural stability and higher binding affinity through molecular dynamics and MM-GBSA studies. AR was identified as a hub protein that is highly expressed in breast cancer and 2-hydroxynaringenin efficacy of counter TNBC requires further investigation both in vitro and in vivo.

Alzheimer’s disease (AD) is the most common neurodegenerative disease followed by various other disorders such as Parkinson’s disease, amyloid lateral sclerosis, Huntington’s, and glioblastoma. The target that has been a key focus in this study is N-methyl-D-aspartate (NMDA) receptor that is responsible for the transmission, integration, and plasticity of excitatory signals needed for the proper functioning of the central nervous system. NMDA is also responsible for the Ca2+ influx. Any damage to the NMDA receptor will lead to neuronal cell death which, in turn, leads to AD. Damage to the receptor also affects Ca2+ influx by developing a lead molecule that can partially bind to the receptor and when it is damaged can prevent neural death. An alkaloid compound called moupinamide (N-trans-feruloyltyramine) which was determined as a lead molecule was further optimized utilizing various in silico methods. Using moupinamide as the basic framework, initial and secondary drug designing is done and then it is followed by docking studies, absorption, distribution, metabolism, and excretion analysis, and molecular simulations. In this study, it was observed that optimized lead molecule high modulus polyethylene could be a potential lead molecule as it showed a great potential.

Ligand-based pharmacophore filtering, atom based 3D-QSAR, virtual screening and ADME studies for the discovery of potential ck2 inhibitors

Computational insights into potential marine natural products as selective inhibitors of Mycobacterium tuberculosis InhA: A structure-based virtual screening study

One of the significant therapeutic targets for Alzheimer’s disease (AD) is Glycogen Synthase Kinase-3β (GSK-3β). Inhibition of GSK-3β can prevent hyperphosphorylation of tau, and thus prevent formation and accumulation of neurofibrillary tangles and neuropil threads that block intracellular transport, trigger unfolded protein response, and increase oxidative stress, cumulatively leading to neurodegeneration. In this study, we have performed structure-based virtual screening of two small-molecule libraries from ChemDiv CNS databases using AutoDock Vina to identify hit molecules based on their binding affinities compared to that of an established GSK-3β inhibitor, indirubin-3’-monoxime (IMO). Pharmacoinformatic screening on SwissADME and pkCSM servers enabled identification of lead molecules with favorable pharmacoinformatic properties for drug likeliness, including blood brain barrier (BBB) permeability. Further, molecular dynamic simulations identified six candidate lead molecules that show stable complex formation with GSK-3β in dynamic state under physiological conditions. Principal component analysis of the dynamic state was used to plot Free Energy Landscapes (FELs) of GSK-3β-ligand complexes. STRIDE secondary structure analysis of the lowest energy conformations identified from FEL plots, and binding free energy calculations by Molecular Mechanics Poisson-Boltzmann Surface Area ((ΔGbind )MM-PBSA) of the simulation trajectories led to the identification of two ligands as potential lead molecules having favorable free energy landscape profiles as well as significantly lower (ΔGbind )MM-PBSA in dynamic state compared to that of reference inhibitor IMO. Hence, this study identifies two novel brain penetrant GSK-3β inhibitors that are likely to have therapeutic potential against Alzheimer’s disease.