Carbon fullerene acts as potential lead molecule against prospective molecular targets of biofilm-producing multidrug-resistant Acinetobacter baumanni and Pseudomonas aerugenosa: computational modeling and MD simulation studies

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

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.

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

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.

The bacterium Bacillus coagulans has attracted interest because of its ability to produce spores and advantageous probiotic traits, such as facilitating food digestion in the intestine, managing some disorders, and controlling the symbiotic microbiota. Spore-forming probiotic bacteria are especially important in the probiotic industry compared to non-spore-forming bacteria due to their stability during production and high resistance to adverse factors such as stomach acid. When spore-forming bacteria are exposed to environmental stresses, they enter the sporulation pathway to survive. This pathway is activated by the final phosphorylation of the master regulator of spore response, Spo0A, and upon achieving the phosphorylation threshold. Spo0A is indirectly inhibited by some enzymes of the aspartate response regulator phosphatase (Rap) family, such as RapJ. RapJ is one of the most important Rap enzymes in the sporogenesis pathway, which is naturally inhibited by the pentapeptides. This study used structure-based virtual screening and molecular dynamics (MD) simulation studies to find potential RapJ hits that could induce the sporulation pathway. The crystal structures of RapJ complexed with pentapeptide clearly elucidated their interactions with the enzyme active site. Based on the binding compartment, through molecular docking, MD simulation, hydrogen bonds, and binding-free energy calculations, a series of novel hits against RapJ named tandutinib, infigratinib, sitravatinib, linifanib, epertinib, surufatinib, and acarbose were identified. Among these compounds, acarbose obtained the highest score, especially in terms of the number of hydrogen bonds, which plays a major role in stabilizing RapJ-ligand complexes, and also according to the occupancy percentages of hydrogen bonds, its hydrogen bonds were more stable during the simulation time. Consequently, acarbose is probably the most suitable hit for RapJ enzyme. Notably, experimental validation is crucial to confirm the effectiveness of the selected ligands.

Herein, we report the synthesis, and characterization of a new series of 1,3,4–oxadiazole and 1,2,4–triazole derivatives based on azaindole acetamides and assigned as potential antibacterial and antitubercular substances. The structures of these compounds were established by 1H NMR, 13C NMR, and HRMS spectral analysis. In preliminary antibacterial studies, analogues 6b, 6d, and 6e were found to be most effective against S. aureus with MIC of 12.5, 6.25, and 12.5 μg/mL, whereas 8d displayed excellent activity against S. aureus, B. subtilis, E. coli bacterial strains with zones of inhibition 12.5, 25, and 12.5 μg/mL respectively. Particularly, the prepared scaffolds 8c, 8d, and 8e showed remarkable antifungal activity with MIC value 12.5, 12.5, and 6.25 μg/mL against A. flavus and 6d, 6c producing an increase in the activity against C. Albicans with zones of inhibition 12.5 and 12.5 μg/mL respectively. Also, through the antitubercular studies, we found that compounds 6e and 8b have a strong activity with M. tuberculosis H37Rv with MICs 3.26, and 6.48 μg/mL, respectively. The protein stability, fluctuations of APO-Protein, and protein-ligand complexes were investigated through Molecular Dynamics (MD) simulations studies using Desmond Maestro 11.3, and potential lead molecules were identified. Our findings were further confirmed using molecular docking, revealing that azaindole based ligand 6e, 6f, and 8a has strong hydrophobic Tyr179, Trp183, Ile177, Ile445, and H–bondings interactions Arg151 and Arg454 through molecular dynamics simulation studies, making it potential biological compound. These compounds were further evaluated for their ADMET and physicochemical properties by using SwissADME. Communicated by Ramaswamy H. Sarma

Snake venom metalloproteinase (SVMP) (Echis coloratus (Carpet viper) is a multifunctional enzyme that is involved in producing several symptoms that follow a snakebite, such as severe local hemorrhage, nervous system effects and tissue necrosis. Because the three-dimensional (3D) structure of SVMP is not known, models were constructed, and the best model was selected based on its stereo-chemical quality. The stability of the modeled protein was analyzed through molecular dynamics (MD) simulation studies. Structure-based virtual screening was performed, and 15 potential molecules with the highest binding energies were selected. Further analysis was carried out with induced fit docking, Prime/MM–GBSA (ΔGBind calculations), quantum-polarized ligand docking, and density functional theory calculations. Further, the stability of the lead molecules in the SVMP-active site was examined using MD simulation. The results showed that the selected lead molecules were highly stable in the active site of SVMP. Hence, these molecules could potentially be selective inhibitors of SVMP. These lead molecules can be experimentally validated, and their backbone structural scaffold could serve as building blocks in designing drug-like molecules for snake antivenom.

Gonorrhea, one of the sexually transmitted disease caused by a gram negative diplococcus bacteria Neisseria gonorrhoeae. Rho protein is indispensable for bacterial viability due to its versatile functions in physiology apart from RNA dependent transcription termination. Based on conserved function and wider role in several cellular processes, inhibitors specifically targeting Rho proteins are largely in use these days to treat various bacterial infections. In this study, three dimensional structure of Rho protein was modeled using the template protein from E. coli and further the optimized model was simulated for 100 ns to understand the structural stability and compactness. Owing to the therapeutic potential of Rho, traditional structure-based virtual screening was applied to identify potential inhibitors for the selected target. Based on empirical glide scoring functions two potent lead molecules (ChemBridge_6121956 and ChemBridge_5232688) were selected from ChemBridge database. The pharmacokinetic properties of these lead molecules are within the permissible range. DFT descriptor revealed that the lead molecules are more reactive, which also supports the molecular docking studies. The stability of Rho and Rho-inhibitor complexes was studied using molecular dynamics simulation. Parameters include binding free energy calculation, RMSD, RMSF and hydrogen bond analysis depicts the stability of Rho and Rho-inhibitors throughout the simulation. Altogether, the identified lead molecules require further optimization towards the design and development of new antibiotics against N. gonorrhoeae. Communicated by Ramaswamy H. Sarma

Carboxymuconolactone decarboxylase is a prospective molecular target for multi-drug resistant Acinetobacter baumannii-computational modeling, molecular docking and dynamic simulation studies

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

Carbon fullerene and nanotube are probable binders to multiple targets of SARS-CoV-2: Insights from computational modeling and molecular dynamic simulation studies.

Alzheimer is a severe memory and cognitive impairment neurodegenerative disease that is the most common cause of dementia worldwide and characterized by the pathological accumulation of tau protein and amyloid-beta peptides. In this study, we have developed E-pharmacophore modeling to screen the eMolecules database with the help of a reported co-crystal structure bound with Beta-Site Amyloid Precursor Protein Cleaving Enzyme 1 (BACE-1). Flumemetamol, florbetaben, and florbetapir are currently approved drugs for use in the clinical diagnosis of Alzheimer’s disease. Despite the benefits of commercially approved drugs, there is still a need for novel diagnostic agents with enhanced physicochemical and pharmacokinetic properties compared to those currently used in clinical practice and research. In the E-pharmacophore modeling results, it is revealed that two aromatic rings (R19, R20), one donor (D12), and one acceptor (A8) are obtained, and also that similar pharmacophoric features of compounds are identified from pharmacophore-based virtual screening. The identified screened hits were filtered for further analyses using structure-based virtual screening and MM/GBSA. From the analyses, top hits such as ZINC39592220 and en1003sfl.46293 are selected based on their top docking scores (−8.182 and −7.184 Kcal/mol, respectively) and binding free energy (−58.803 and −56.951 Kcal/mol, respectively). Furthermore, a molecular dynamics simulation and MMPBSA study were performed, which revealed admirable stability and good binding free energy throughout the simulation period. Moreover, Qikprop results revealed that the selected, screened hits have good drug-likeness and pharmacokinetic properties. The screened hits ZINC39592220 and en1003sfl.46293 could be used to develop drug molecules against Alzheimer’s disease.

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.

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.

The study aimed to screen prospective molecular targets of BCC and potential natural lead candidates as effective binders by computational modeling, molecular docking, and dynamic (MD) simulation studies. Based on the virulent functions, tRNA 5-methylaminomethyl-2-thiouridine biosynthesis bifunctional protein (mnmC) and pyrimidine/purine nucleoside phosphorylase (ppnP) were selected as the prospective molecular targets. In the absence of experimental data, the three-dimensional (3D) structures of these targets were computationally predicted. After a thorough literature survey and database search, the drug-likeness, and pharmacokinetic properties of 70 natural molecules were computationally predicted and the effectual binding of the best lead molecules against both the targets was predicted by molecular docking. The stabilities of the best-docked complexes were validated by MD simulation and the binding energy calculations were carried out by MM-GBSA approaches. The present study revealed that the hypothetical models of mnmC and ppnP showed stereochemical accuracy. The study also showed that among 70 natural compounds subjected to computational screening, Honokiol (3',5-Di(prop-2-en-1-yl) [1,1'-biphenyl]-2,4'-diol) present in Magnolia showed ideal drug-likeness, pharmacokinetic features and showed effectual binding with mnmC and ppnP (binding energies -7.3kcal/mol and -6.6kcal/mol, respectively). The MD simulation and GBSA calculation studies showed that the ligand-protein complexes stabilized throughout tMD simulation. The present study suggests that Honokiol can be used as a potential lead molecule against mnmC and ppnP targets of BCC and this study provides insight into further experimental validation for alternative lead development against drug resistant BCC.