Stabilization of microtubules by taxane diterpenoids: insight from docking and MD simulations

Machine Learning Integrating Protein Structure, Sequence, and Dynamics to Predict the Enzyme Activity of Bovine Enterokinase Variants

Despite recent advances in computational protein science, the dynamic behavior of proteins, which directly governs their biological activity, cannot be gleaned from sequence information alone. To overcome this challenge, we propose a framework that integrates the peptide sequence, protein structure, and protein dynamics descriptors into machine learning algorithms to enhance their predictive capabilities and achieve improved prediction of the protein variant function. The resulting machine learning pipeline integrates traditional sequence and structure information with molecular dynamics simulation data to predict the effects of multiple point mutations on the fold improvement of the activity of bovine enterokinase variants. This study highlights how the combination of structural and dynamic data can provide predictive insights into protein functionality and address protein engineering challenges in industrial contexts.

Identification of natural marine compounds as potential inhibitors of CDK2 using molecular docking and molecular dynamics simulation approach

The multifunctional enzyme cyclin-dependent kinase 2 (CDK2) protein is essential for cell proliferation, transcription and modulation of the cell cycle. There is a dysfunction that is connected to various diseases, such as cancer, making it an important treatment target in oncology and beyond. The goal of this study is to identify novel CDK2 ATP binding site inhibitors using in silico drug designing. To find competitive inhibitors for the ATP site, molecular docking, molecular dynamics (MD) simulation and free-binding energy calculations were used. Natural compounds retrieved from marine sources (fungi and algae) were docked against protein, and the best-binding compounds were further evaluated using MD simulations. LIG1, LIG2 and LIG3 (ΔGPB = -19.98, -15.82 and -12.98 kcal/mol, respectively) were placed in the top positions based on their overall binding energy calculated using MMPBSA approach. Stability of the complexes was confirmed by carefully analyzing the rmsd and rmsf patterns retrieved from the MD trajectories. Several residues and areas (Leu124, Val123, Phe80, Leu83, Glu81, Arg 126, Asn132, Leu134, Gln131, Lys88 and Glu195) appear to be critical for inhibitor retention across the active pocket, according to RMSD and RMSF. The dynamics of the ligands inside the active pocket were mapped using principle component analysis. It has been observed that LIG1-3 appear to be the best possible inhibitors due to their high binding energies, interaction pattern and retention inside the active pocket.Communicated by Ramaswamy H. Sarma.

Rational Engineering Docetaxel Prodrug Nanoassemblies: Response Modules Guiding Efficacy Enhancement and Toxicity Reduction

Prodrug-based nanoassemblies have been developed to solve the bottlenecks of chemotherapeutic drugs. The fabricated prodrugs usually consist of active drug modules, response modules, and modification modules. Among three modules, the response modules play a vital role in controlling the intelligent drug release at tumor sites. Herein, various locations of disulfide bond linkages were selected as response modules to construct three Docetaxel (DTX) prodrugs. Interestingly, the small structural difference caused by the length of response modules endowed corresponding prodrug nanoassemblies with unique characteristic. α-DTX-OD nanoparticles (NPs) possessed the advantages of high redox-responsiveness due to their shortest linkages. However, they were too sensitive to retain the intact structure in the blood circulation, leading to severe systematic toxicity. β-DTX-OD NPs significantly improved the pharmacokinetics of DTX but may induce damage to the liver. In comparison, γ-DTX-OD NPs with the longest linkages greatly ameliorated the delivery efficiency of DTX as well as improved DTX's tolerance dose.

Citrus nutraceutical eriocitrin and its metabolites are partial agonists of peroxisome proliferator-activated receptor gamma (PPARγ): a molecular docking and molecular dynamics study

Peroxisome proliferator-activated receptor gamma (PPARγ) agonists are potent insulin sensitizers in treating type 2 diabetes. Despite being very effective in the fight against diabetes-mediated complications, PPARγ agonists are accompanied by severe side effects leading to complicated health problems, making the discovery of novel safe ligands highly pertinent. A significant intense research effort is in progress to explore the PPARγ activating potential of a wide range of natural compounds. Lemon (Citrus limon) contains various bioactive flavonoids, and eriocitrin is the major flavonoid. It possesses substantial antioxidant and anticancer, lipid-lowering activities and prevents obesity-associated metabolic diseases. Eriocitrin is metabolized to eriodictyol in the intestine, and the absorbed eriodictyol undergoes conversion to numerous metabolites in vivo. It is unclear if eriocitrin or its metabolites are responsible for their beneficial effects. We have used molecular docking, ADMET properties, drug-likeness score and molecular dynamics simulation study to find if eriocitrin and its metabolites are potent binders for PPARγ. Docking studies revealed that eriocitrin binds to PPARγ with the highest binding energy, but ADMET properties and in vivo studies show that the bioavailability of eriocitrin is very poor. Molecular dynamics studies were carried out to validate the docking results, and multiple parameters like RMSD, RMSF, Radius of gyration, SASA, hydrogen bond analysis, interaction energy, principal component analysis, Gibbs free energy and MM-PBSA were calculated. Based on our studies, eriodictyol, eriodictyol 7-O-glucuronide, eriodictyol 3'-O-glucuronide, homoeriodictyol and homoeriodictyol 7-O-glucuronide which are metabolites of eriocitrin appear to be potent partial agonists of PPARγ under physiological conditions.Communicated by Ramaswamy H. Sarma.

Purification of Two Taxanes from Taxus cuspidata by Preparative High-Performance Liquid Chromatography

In the present study, an effective method of preparative high-performance liquid chromatography (Prep-HPLC) was established to purify two taxanes in Taxus cuspidata. During the experimental operation, the effects of flow rate, injection volume, and column temperature on the purity of 10-deacetyltaxol (10-DAT) and paclitaxel (PTX) were investigated, and the optimized conditions were as follows: flow rate of 10 mL/min, injection volume of 0.5 mL, and column temperature of 30 °C. Under these conditions, the purity of 10-DAT and PTX reached 95.33% and 99.15%, respectively. The purified products were characterized by scanning electron microscopy (SEM), high-performance liquid chromatography (HPLC), and electrospray ionization-high resolution mass spectrometry (ESI-HRMS). The results demonstrated that preparative HPLC can effectively purify 10-DAT and PTX from Taxus cuspidata with a purity of >95%, which was suitable for the large-scale preparation of 10-DAT and PTX.

Separation and Purification of Taxanes from Crude Taxus cuspidata Extract by Antisolvent Recrystallization Method

Taxanes are natural compounds with strong antitumor activity. In this study, we first enriched taxanes by ultrasonic extraction and liquid–liquid extraction from Taxus cuspidata, then purified these taxanes by the antisolvent recrystallization method, and discussed the effects of four recrystallization conditions on the purity of eight target compounds. The most promising purification results were obtained using methanol as a solvent and water as an antisolvent. Response surface methodology (RSM) was used to further optimize the optimal purification conditions: when the crude extraction concentration was 555.28 mg/mL, an antisolvent to solvent volume ratio was 28.16 times, the deposition temperature was 22.91 °C, and the deposition time was 1.76 min, the purity of the taxanes reached its maximum. The scanning electron microscopy (SEM) results showed that recrystallization could effectively reduce the particle size of crude Taxus cuspidata and control the particle morphology. X-ray diffraction (XRD) and Raman spectrum experiments demonstrated that the amorphous state of the crude Taxus cuspidata did not change during the recrystallization process, and always remained amorphous. This recrystallization method can effectively improve the purity of taxanes in Taxus cuspidata, and is suitable for the preliminary purification of taxanes.

Anti-proliferative activity of labdane diterpenes isolated from Polyalthia cerasoides and their molecular interaction studies

BACKGROUND

Polyalthia cerasoides is well known for its therapeutic effects and is extensively used by the tribal people of South India and Africa to treat infertility, toothache, inflammation, rheumatism, fever, and to combat stress.

OBJECTIVE

In the present research, the anti-proliferative potential of two bioactive compounds isolated from the stem bark of P. cerasoides (Roxb.) Bedd. of the Annonaceae family were investigated.

METHOD

The dried stem bark was powdered and subjected to extraction using methanol and further partitioned using petroleum ether. Yellow viscous oil was isolated from the petroleum ether fraction using column and preparative thin-layer chromatography. The chromatographic fractions were characterized using GC-MS. The anti-proliferative effect of the isolated compounds was assessed against HepG2 Cells using MTT- Cytotoxicity test. Furthermore, comparative in-silico docking studies were performed to predict the binding pattern of isolated molecules individually, as well as simultaneously with α, β -tubulin, a critical protein involved in the molecular mechanism of microtubule formation.

RESULTS

GC-MS analysis of yellow viscous oil from petroleum fraction confirmed the presence of two labdane diterpenes and were identified as 12E-3,4-Seco-labda-4(18),8(17),12,14-tetraen-3-oic acid , and methyl harvadate C by mass fragmentation analysis. The MTT-Cytotoxicity assay showed the dose-dependent cytotoxic effect on HepG2 Cells. The comparative docking studies of the isolated compounds exhibited strong interactions with the α, β -tubulin protein.

CONCLUSION

The prominent anti-proliferative effect exhibited by the isolated compounds along with effective binding to α, β -tubulin protein encourages their future utilization as prominent anti-cancer molecules.

Discovery of microtubule stabilizers with novel scaffold structures based on virtual screening, biological evaluation, and molecular dynamics simulation

Disrupting the dynamics and structures of microtubules can perturb mitotic spindle formation, cause cell cycle arrest in G2/M phase, and subsequently lead to cellular death via apoptosis. In this investigation, the structure-based virtual screening methods, including molecular docking and rescoring, and similarity analysis of interaction molecular fingerprints, were developed to discover novel tubulin inhibitors from ChemDiv database with 1,601,806 compounds. The screened compounds were further filtered by PAINS, ADME/T, Toxscore, SAscore, and Drug-likeness analysis. Finally, 17 hit compounds were selected, and then submitted to the biologic evaluation. Among these hits, the P2 exhibited the strongest antiproliferative activity against four tumor cells including HeLa, HepG2, MCF-7, and A549. The in vitro tubulin polymerization assay revealed P2 could promote tubulin polymerization in a dose dependent manner. Finally, in order to analyze the interaction modes of complexes, the molecular dynamics simulation was performed to investigate the interactions between P2 and tubulin. The molecular dynamics simulation analysis showed that P2 could stably bind to taxane site, induced H6-H7, B9-B10, and M-loop regions changes. The ΔG_bind energies of tubulin-P2 and tubulin-paclitaxel were -68.25 ± 12.98 and -146.05 ± 16.17 kJ mol^-1, respectively, which were in line with the results of the experimental test. Therefore, P2 has been well characterized as lead compounds for developing new tubulin inhibitors with potential anticancer activity.

Mechanisms of influence of the microtubule over-stabilizing ligands on the structure and intrinsic dynamics of α,β-Tubulin

The intervention into the cell cycle progression by administering microtubule over-stabilizing ligands that arrest the mitotic cell division by preventing spindle dissociation, is a promising strategy to fight against cancers. The building blocks of the microtubules and the spindles, i.e. the α,β-tubulin dimer, upon binding of such ligands, stay more comfortably in the microtubular multimeric form; the phenomenon of which is the key to the said over-stabilization. Using two such over-stabilizing ligands, Taxol and Taxotere, the present work reports the collective changes that these ligands induce on the structure and dynamics of the α,β-tubulin dimer which could be reconciled as the molecular basis of the over-stabilization of the microtubules; the trends have been found to be statistically significant across all independent calculations on them. The ligand binding increases the coherence between the residue communities of the two opposite faces of the β-subunit, which in a periodic arrangement in microtubule are knwon to form intermolecular contact with each other. This is likely to create an indirect cooperativity between those structural regions and this is a consequence of the reshuffling of the internal network of interactions upon ligand binding. Such reorganizations are also complemented by the increased contributions of the softer modes of the intrinsic dynamics more, which is likely to increase the plasticity of the system favourable for making structural adjustments in a multimer. Further, the ligands are able to compensate the drawback of lacking one phosphate group in protein-GDP interactions compared to the same for protein-GTP and this is in agreement with the hints form the earlier reports. The findings form a mechanistic basis of the enhanced capacity of the α,β-tubulin dimer to get more favourably accommodated into the microtubule superstructure upon binding either of Taxol and Taxotere.

Evaluation of binding of potential ADMET/tox screened saquinavir analogues for inhibition of HIV-protease via molecular dynamics and binding free energy calculations

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.

A computational prospective on isoform-selective CB2 inhibitors

The strategy of enlarging the molecular size, increasing hydrogen bond interactions with the key residues of CB2, and furthering geometric distribution of hydrophobic groups would improve the selectivity towards CB2 receptor.

Computational and In Vitro Investigation of (-)-Epicatechin and Proanthocyanidin B2 as Inhibitors of Human Matrix Metalloproteinase 1

Matrix metalloproteinases 1 (MMP-1) energetically triggers the enzymatic proteolysis of extracellular matrix collagenase (ECM), resulting in progressive skin aging. Natural flavonoids are well known for their antioxidant properties and have been evaluated for inhibition of matrix metalloproteins in human. Recently, (-)-epicatechin and proanthocyanidin B2 were reported as essential flavanols from various natural reservoirs as potential anti-inflammatory and free radical scavengers. However, their molecular interactions and inhibitory potential against MMP-1 are not yet well studied. In this study, sequential absorption, distribution, metabolism, and excretion (ADME) profiling, quantum mechanics calculations, and molecular docking simulations by extra precision Glide protocol predicted the drug-likeness of (-)-epicatechin (−7.862 kcal/mol) and proanthocyanidin B2 (−8.145 kcal/mol) with the least reactivity and substantial binding affinity in the catalytic pocket of human MMP-1 by comparison to reference bioactive compound epigallocatechin gallate (−6.488 kcal/mol). These flavanols in docked complexes with MMP-1 were further studied by 500 ns molecular dynamics simulations that revealed substantial stability and intermolecular interactions, viz. hydrogen and ionic interactions, with essential residues, i.e., His218, Glu219, His222, and His228, in the active pocket of MMP-1. In addition, binding free energy calculations using the Molecular Mechanics Generalized Born Surface Area (MM/GBSA) method suggested the significant role of Coulomb interactions and van der Waals forces in the stability of respective docked MMP-1-flavonol complexes by comparison to MMP-1-epigallocatechin gallate; these observations were further supported by MMP-1 inhibition assay using zymography. Altogether with computational and MMP-1–zymography results, our findings support (-)-epicatechin as a comparatively strong inhibitor of human MMP-1 with considerable drug-likeness against proanthocyanidin B2 in reference to epigallocatechin gallate.

Identification of Phosphoinositide-3 Kinases Delta and Gamma Dual Inhibitors Based on the p110δ/γ Crystal Structure

Background:

Phosphoinositide-3 kinases (PI3Ks) are key signaling molecules that affect a diverse array of biological processes in cells, including proliferation, differentiation, survival, and metabolism. The abnormal activity of PI3K signals is closely related to the occurrence of many diseases, which has become a very promising drug target, especially for the treatment of cancer. PI3Kδ/γ inhibitors can reduce toxicity concerns for chronic indications such as asthma and rheumatoid arthritis compared with pan PI3Ks inhibitors.

Methods:

With the aim of finding more effective PI3Kδ/γ dual inhibitors, virtual screening, ADMET prediction Molecular Dynamics (MD) simulations and MM-GBSA were executed based on the known p110δ/γ crystal structure. Compound ZINC28564067 with high docking score and low toxicity was obtained.

Results:

By MD simulations and MM-GBSA, we could observe that ZINC28564067 had more favorable conformation binding to the PI3Kδ/γ than the original ligands.

Conclusion:

The results provided a rapid approach for the discovery of novel PI3Kδ/γ dual inhibitors which might be a potential anti-tumor lead compound.

Rational drug repurposing for cancer by inclusion of the unbiased molecular dynamics simulation in the structure-based virtual screening approach: Challenges and breakthroughs

Managing cancer is now one of the biggest concerns of health organizations. Many strategies have been developed in drug discovery pipelines to help rectify this problem and two of the best ones are drug repurposing and computational methods. The combination of these approaches can have immense impact on the course of drug discovery. In silico drug repurposing can significantly reduce the time, the cost and the effort of drug development. Computational methods such as structure-based drug design (SBDD) and virtual screening can predict the potentials of small molecule binders, such as drugs, for having favorable effect on a particular molecular target. However, the demand for accuracy and efficiency of SBDD requires more sophisticated and complicated approaches such as unbiased molecular dynamics (UMD) simulation that has been recently introduced. As a complementary strategy, the knowledge acquired from UMD simulations can increase the chance of finding the right candidates and the pipeline of its administration is introduced and discussed in this review. An elaboration of this pipeline is also made by detailing an example, the binding and unbinding pathways of dasatinib-c-Src kinase complex, which shows that how influential this method can be in rational drug repurposing in cancer treatment.

Structure-based screening and validation of bioactive compounds as Zika virus methyltransferase (MTase) inhibitors through first-principle density functional theory, classical molecular simulation and QM/MM affinity estimation

Recent Zika virus (ZIKV) outbreak and association with human diseases such as neurological disorders have raised global health concerns. However, in the absence of an approved anti-ZIKV drug has generated urgency for the drug development against ZIKV infection. Here, structure-based virtual screening of 8589 bioactive compounds, screened at the substrate-binding site of ZIKV nonstructural 5 (NS5)-based structure N-terminal methyltransferase (MTase) domain followed by ADMET (absorption, distribution, metabolism, excretion and toxicity) profiling concluded the four potential lead inhibitors, i.e. (4-acetylamino-benzenesulfonylamino)-acetic acid (F3342-0450), 3-(5-methylfuran-2-yl)-N-(4-sulfamoylphenyl)propanamide (F1736-0142), 8-(2-hydroxy-ethylamino)-1,3-dimethyl-7-(3-methyl-benzyl)-3,7-dihydro-purine-2,6-dione (F0886-0080) and N-[4-(aminosulfonyl)phenyl]-2,3-dihydro-1,4-benzodioxine-2-carboxamide (F0451-2187). Collectively, extra precision docking and Density Functional Theory(DFT) calculations studies identified the F3342-0450 molecule, having strong interactions on the active site of MTase, further supported by molecular dynamics simulation, binding affinity and hybrid QM/MM calculations, suggest a new drug molecule for the antiviral drug development against ZIKV infection. Communicated by Ramaswamy H. Sarma.

Tubulin-Binding 3,5-Bis(styryl)pyrazoles as Lead Compounds for the Treatment of Castration-Resistant Prostate Cancer

The microtubule-binding taxanes, docetaxel and cabazitaxel, are administered intravenously for the treatment of castration-resistant prostate cancer (CRPC) as the oral administration of these drugs is largely hampered by their low and highly variable bioavailabilities. Using a simple, rapid, and environmentally friendly microwave-assisted protocol, we have synthesized a number of 3,5-bis(styryl)pyrazoles 2a-l, thus allowing for their screening for antiproliferative activity in the androgen-independent PC3 prostate cancer cell line. Surprisingly, two of these structurally simple 3,5-bis(styryl)pyrazoles (2a and 2l) had concentrations which gave 50% of the maximal inhibition of cell proliferation (GI₅₀) in the low micromolar range in the PC3 cell line and were thus selected for extensive further biologic evaluation (apoptosis and cell cycle analysis, and effects on tubulin and microtubules). Our findings from these studies show that 3,5-bis[(1E)-2(2,6-dichlorophenyl)ethenyl]-1H-pyrazole 2l 1) caused significant effects on the cell cycle in PC3 cells, with the vast majority of treated cells in the G₂/M phase (89%); 2) induces cell death in PC3 cells even after the removal of the compound; 3) binds to tubulin [dissociation constant (K_d) 0.4 ± 0.1 μM] and inhibits tubulin polymerization in vitro; 4) had no effect upon the polymerization of the bacterial cell division protein FtsZ (a homolog of tubulin); 5) is competitive with paclitaxel for binding to tubulin but not with vinblastine, crocin, or colchicine; and 6) leads to microtubule depolymerization in PC3 cells. Taken together, these results suggest that 3,5-bis(styryl)pyrazoles warrant further investigation as lead compounds for the treatment of CRPC. SIGNIFICANCE STATEMENT: The taxanes are important components of prostate cancer chemotherapy regimens, but their oral administration is hampered by very low and highly variable oral bioavailabilities resulting from their poor absorption, poor solubility, high first-pass metabolism, and efficient efflux by P-glycoprotein. New chemical entities for the treatment of prostate cancer are thus required, and we report here the synthesis and investigation of the mechanism of action of some bis(styryl)pyrazoles, demonstrating their potential as lead compounds for the treatment of prostate cancer.

Discovery of Ganoderma lucidum triterpenoids as potential inhibitors against Dengue virus NS2B-NS3 protease

Dengue virus (DENV) infection causes serious health problems in humans for which no drug is currently available. Recently, DENV NS2B-NS3 protease has been proposed as a primary target for anti-dengue drug discovery due to its important role in new virus particle formation by conducting DENV polyprotein cleavage. Triterpenoids from the medicinal fungus Ganoderma lucidum have been suggested as pharmacologically bioactive compounds and tested as anti-viral agents against various viral pathogens including human immunodeficiency virus. However, no reports are available concerning the anti-viral activity of triterpenoids from Ganoderma lucidum against DENV. Therefore, we employed a virtual screening approach to predict the functional triterpenoids from Ganoderma lucidum as potential inhibitors of DENV NS2B-NS3 protease, followed by an in vitro assay. From in silico analysis of twenty-two triterpenoids of Ganoderma lucidum, four triterpenoids, viz. Ganodermanontriol (−6.291 kcal/mol), Lucidumol A (−5.993 kcal/mol), Ganoderic acid C2 (−5.948 kcal/mol) and Ganosporeric acid A (−5.983 kcal/mol) were predicted to be viral protease inhibitors by comparison to reference inhibitor 1,8-Dihydroxy-4,5-dinitroanthraquinone (−5.377 kcal/mol). These results were further studied for binding affinity and stability using the molecular mechanics/generalized Born surface area method and Molecular Dynamics simulations, respectively. Also, in vitro viral infection inhibition suggested that Ganodermanontriol is a potent bioactive triterpenoid.

Computational aided mechanistic understanding of Camellia sinensis bioactive compounds against co-chaperone p23 as potential anticancer agent

Co-chaperon p23 has been well established as molecular chaperon for the heat shock protein 90 (Hsp90) that further leads to immorality in cancer cells by providing defense against Hsp90 inhibitors, and as stimulating agent for generating overexpressed antiapoptotic proteins, that is, Hsp70 and Hsp27. The natural compounds such as catechins from Camellia sinensis (green tea) are also well known for inhibition activity against various cancer. However, molecular interaction profile and potential lead bioactive compounds against co-chaperon p23 from green tea are not yet reported. To this context, we study the various secondary metabolites of green tea against co-chaperon p23 using structure-based virtual screening from Traditional Chinese Medicine (TCM) database. Following 26 compounds were obtained from TCM database and further studied for extra precision molecular docking that showed binding score between -10.221 and -2.276 kcal/mol with co-chaperon p23. However, relative docking score to known inhibitors, that is, ailanthone (-4.54 kcal/mol) and gedunin ( 3.60 kcal/mol) along with ADME profile analysis concluded epicatechin (-7.013 kcal/mol) and cis-theaspirone (-4.495 kcal/mol) as potential lead inhibitors from green tea against co-chaperone p23. Furthermore, molecular dynamics simulation and molecular mechanics generalized born surface area calculations validated that epicatechin and cis-theaspirone have significantly occupied the active region of co-chaperone p23 by hydrogen and hydrophobic interactions with various residues including most substantial amino acids, that is, Thr90, Ala94, and Lys95. Hence, these results supported the fact that green tea contained potential compounds with an ability to inhibit the cancer by disrupting the co-chaperon p23 activity.

The inhibitory mechanism of aurintricarboxylic acid targeting serine/threonine phosphatase Stp1 in Staphylococcus aureus: insights from molecular dynamics simulations

Serine/threonine phosphatase (Stp1) is a member of the bacterial Mg2+- or Mn2+- dependent protein phosphatase/protein phosphatase 2C family, which is involved in the regulation of Staphylococcus aureus virulence. Aurintricarboxylic acid (ATA) is a known Stp1 inhibitor with an IC50 of 1.03 μM, but its inhibitory mechanism has not been elucidated in detail because the Stp1-ATA cocrystal structure has not been determined thus far. In this study, we performed 400 ns molecular dynamics (MD) simulations of the apo-Stp1 and Stp1-ATA complex models. During MD simulations, the flap subdomain of the Stp1-ATA complex experienced a clear conformational transition from an open state to a closed state, whereas the flap domain of apo-Stp1 changed from an open state to a semi-open state. In the Stp1-ATA complex model, the hydrogen bond (H-bond) between D137 and N142 disappeared, whereas critical H-bond interactions were formed between Q160 and H13, Q160/R161 and ATA, as well as N162 and D198. Finally, four residues (D137, N142, Q160, and R161) in Stp1 were mutated to alanine and the mutant enzymes were assessed using phosphate enzyme activity assays, which confirmed their important roles in maintaining Stp1 activity. This study indicated the inhibitory mechanism of ATA targeting Stp1 using MD simulations and sheds light on the future design of allosteric Stp1 inhibitors.

Evaluation of caffeine as inhibitor against collagenase, elastase and tyrosinase using in silico and in vitro approach

Skin ageing results from enhanced activation of intracellular enzymes such as collagenases, elastases and tyrosinase, stimulated by intrinsic ageing and photoageing factors. Recently, caffeine-based cosmetics are introduced that demonstrates to slow down skin photoageing process. However, no attempts have been done so for to understand caffeine functional inhibitory activity against photoageing related enzymes. Hence, this study established the caffeine molecular interaction and inhibition activity profiles against respective enzymes using in silico and in vitro methods, respectively. Results from in silico study indicates that caffeine has comparatively good affinity with collagenase (−4.6 kcal/mol), elastase (−3.36 kcal/mol) and tyrosinase (−2.86 kcal/mol) and formed the stable protein-ligand complex as validated by molecular dynamics simulation (protein-ligand contacts, RMSD, RMSF and secondary structure changes analysis). Moreover, in vitro data showed that caffeine (1000 µg/mL) has statistically significant maximum inhibition activity of 41.86, 36.44 and 13.72% for collagenase, elastase and tyrosinase, respectively.

Computational and In-Vitro Validation of Natural Molecules as Potential Acetylcholinesterase Inhibitors and Neuroprotective Agents

Background:

Cholinesterase inhibitors are the first line of therapy for the management of Alzheimer’s disease (AD), however, it is now established that they provide only temporary and symptomatic relief, besides, having several inherited side-effects. Therefore, an alternative drug discovery method is used to identify new and safer ‘disease-modifying drugs’.

Methods:

Herein, we screened 646 small molecules of natural origin having reported pharmacological and functional values through in-silico docking studies to predict safer neuromodulatory molecules with potential to modulate acetylcholine metabolism. Further, the potential of the predicted molecules to inhibit acetylcholinesterase (AChE) activity and their ability to protect neurons from degeneration was determined through in-vitro assays.

Results:

Based on in-silico AChE interaction studies, we predicted quercetin, caffeine, ascorbic acid and gallic acid to be potential AChE inhibitors. We confirmed the AChE inhibitory potential of these molecules through in-vitro AChE inhibition assay and compared results with donepezil and begacestat. Herbal molecules significantly inhibited enzyme activity and inhibition for quercetin and caffeine did not show any significant difference from donepezil. Further, the tested molecules did not show any neurotoxicity against primary (E18) hippocampal neurons. We observed that quercetin and caffeine significantly improved neuronal survival and efficiently protected hippocampal neurons from HgCl2 induced neurodegeneration, which other molecules, including donepezil and begacestat, failed to do.

Conclusion:

Quercetin and caffeine have the potential as “disease-modifying drugs” and may find application in the management of neurological disorders such as AD.

Identification of isobavachalcone as a potential drug for rice blast disease caused by the fungus Magnaporthe grisea

The rice blast disease caused by the fungus Magnaporthe grisea is one of the most devastating rice diseases, but there is no effective fungicide toward chitinase which is a key enzyme of M. grisea. In this study, we observed that distortion and cell-wall damage of M. grisea hyphae were significantly under the scanning electron micrograph after a 24-h treatment with 10 mg/L isobavachalcone (IBC) extracted from Psoralea corylifolia L. To further explore the effect of IBC on the cell wall of M. grisea, we examined changes in enzymes associated with cell wall degradation by enzyme activity experiments, treated liquid culture mycelia with 10 mg/L IBC for 1 h. Results displayed that chitinase was obviously more active than control group. To illustrate the interactions between IBC and chitinase, the studies of homology modeling and molecular docking were carried out successively. The results revealed that IBC had hydrogen bonds with residues ASP267 and ARG276 of chitinase. Besides, these nonpolar residues TYR270, PRO271, VAL272, LEU310, PRO311, TYR316, and LEU317 were able to form strong hydrophobic interactions. Binding energies of the chitinase-IBC complexes were calculated by MM-GBSA showed that the ΔG_bind score of molecular dynamics had lower binding energy and more stable than docking complexes. All above, IBC owns significant agonistic activity in chitinase and would be a potent fungicide to inhibit the growth of M. grisea. We hope the above information provides an important insight for understanding the interactions between IBC and chitinase, which may be useful in the discovery of a novel potent agonist. Communicated by Ramaswamy H. Sarma.

Syntheses, spectroscopic and crystallographic characterizations of cis- and trans-dispirocyclic ferrocenylphosphazenes: molecular dockings, cytotoxic and antimicrobial activities

The crystallographic and molecular docking investigations and the determinations of the cytotoxic activities of the dispirocyclic ferrocenylphosphazenes were performed.

Design, sythesis and evaluation of a series of 3- or 4-alkoxy substituted phenoxy derivatives as PPARs agonists

Peroxisome proliferators-activated receptors (PPARα, γ and δ) are potentially effective targets for Type 2 diabetes mellitus therapy. The severe effects of known glitazones and the successfully approved agents (saroglitazar and lobeglitazone) motivated us to study novelly potent PPARs drugs with improved safety profile. In this work, we received 15 carboxylic acids based on the combination principle to integrate the polar head of bezafibrate with the hydrophobic tail of pioglitazone. Another 12 tetrazoles based on the bioisosterism principle were obtained accordingly. Furthermore, in vitro PPARs transactivation assays on these 3- or 4-alkoxy substituted phenoxy derivatives afforded six compounds. Interactions and binding stability from the docking analysis and 20 ns molecular dynamic simulations confirmed the representative compounds to be suitable and plausible for PPARs pockets. The above-mentioned results demonstrated that the compounds may be used as reference for further optimization for enhanced PPARs activities and wide safety range.

Novel anti-tubulin agents from plant and marine origins: insight from a molecular modeling and dynamics study

The screening of a variety of botanical species and marine organisms provided satisfactory novel tubulin binding agents (TBAs).

Identification of dual ligands targeting angiotensin II type 1 receptor and peroxisome proliferator-activated receptor-γ by core hopping of telmisartan

It has been reported previously that some angiotensin II receptor blockers not only antagonize angiotensin II type 1 receptor (AT₁R), but also exert stimulation in peroxisome proliferator-activated receptor γ (PPARγ) partial activation, among which telmisartan displays the best. Telmisartan has been tested as a bifunctional ligand with antihypertensive and hypoglycemic activity. Aiming at more potent leads with selective AT₁R antagonism and PPARγ partial agonism, the three parts of telmisartan including the distal benzimidazole ring, the biphenyl moiety, and the carboxylic acid group experienced modification by core hopping method in our study. The central benzimidazole ring, however, remained intact considering its great affinity toward AT₁R and PPARγ. We utilized computational techniques for the sake of details on the binding interactions and conformational stability. Standard precision docking analysis and absorption, distribution, metabolism, excretion, and toxicity prediction received 10 molecules with higher Glide scores, similar interactions, and improved pharmacokinetic profiles compared to telmisartan. Comp#91 with highest scores for AT₁R (-11.92 kcal/mol) and PPARγ (-13.88 kcal/mol) exhibited excellent binding modes and pharmacokinetic parameters. Molecular dynamics trajectories on best docking pose of comp#91 confirmed the docking results and verified the conformational stability with both receptors throughout the course of 20-ns simulations. Thus, comp#91 could be identified as a promising lead in the development of dual AT₁R antagonist and PPARγ partial agonist against hypertension and type 2 diabetes.