Molecular Dynamic Behavior and Binding Affinity of Flavonoid Analogues to the Cyclin Dependent Kinase 6/cyclin D Complex

Kaempferol regulates the thermogenic function of adipocytes in high-fat-diet-induced obesity via the CDK6/RUNX1/UCP1 signaling pathway

Activation of adipose tissue thermogenesis is a promising strategy in the treatment of obesity and obesity-related metabolic disorders. Kaempferol (KPF) is a predominant dietary flavonoid with multiple pharmacological properties, such as anti-inflammatory and antioxidant activities. In this study, we sought to characterize the role of KPF in adipocyte thermogenesis. We demonstrated that KPF-treated mice were protected from diet-induced obesity, glucose tolerance, and insulin resistance, accompanied by markedly increased energy expenditure, ex vivo oxygen consumption of white fat, and increased expression of proteins related to adaptive thermogenesis. KPF-promoted beige cell formation is a cell-autonomous effect, since the overexpression of cyclin-dependent kinase 6 (CDK6) in preadipocytes partially reversed browning phenotypes observed in KPF-treated cells. Overall, these data implicate that KPF is involved in promoting beige cell formation by suppressing CDK6 protein expression. This study provides evidence that KPF is a promising natural product for obesity treatment by boosting energy expenditure.

Characterization of Promising Cytotoxic Metabolites from Tabebuia guayacan Hemsl.: Computational Prediction and In Vitro Testing

Genus Tabebuia is famous for its traditional uses and valuable phytoconstituents. Our previous investigation of Tabebuia species noted the promising anticancer activity of T. guayacan Hemsl. leaves extract, however, the mechanism underlying the observed anticancer activity is still unexplored. The current research was designed to explore the phytochemical content as well as to address the phytoconstituent(s) responsible for the recorded anticancer activity. Accordingly, sixteen compounds were isolated, and their structures were elucidated using different spectroscopic techniques. The drug-likeness of the isolated compounds, as well as their binding affinity with four anticancer drug target receptors: CDK-2/6, topoisomerase-1, and VEGFR-2, were evaluated. Additionally, the most promising compounds were in vitro evaluated for inhibitory activities against CDK-2/6 and VEGFR-2 enzymes using kinase assays method. Corosolic acid (3) and luteolin-7-O-β-glucoside (16) were the most active inhibitors against CDK-2 (-13.44 kcal/mol) and topoisomerase 1 (-13.83 kcal/mol), respectively. Meanwhile, quercetin 3-O-β-xyloside (10) scored the highest binding free energies against both CDK-6 (-16.23 kcal/mol) as well as against VEGFR-2 protein targets (-10.39 kcal/mol). Molecular dynamic simulation indicated that quercetin 3-O-β-xyloside (10) exhibited the least fluctuations and deviations from the starting binding pose with RMSD (2.6 Å). Interestingly, in vitro testing results confirmed the potent activity of 10 (IC50 = 0.154 µg/mL) compared to IC50 = 0.159 µg/mL of the reference drug ribociclib. These findings suggest the three noted compounds (3, 10, and 16) for further in vivo anticancer studies.

Navigating Calcium and Reactive Oxygen Species by Natural Flavones for the Treatment of Heart Failure

Heart failure (HF), the leading cause of death among men and women world-wide, causes great health and economic burdens. HF can be triggered by many factors, such as coronary artery disease, heart attack, cardiomyopathy, hypertension, obesity, etc., all of which have close relations with calcium signal and the level of reactive oxygen species (ROS). Calcium is an essential second messenger in signaling pathways, playing a pivotal role in regulating the life and death of cardiomyocytes via the calcium-apoptosis link mediated by the cellular level of calcium. Meanwhile, calcium can also control the rate of energy production in mitochondria that are the major resources of ROS whose overproduction can lead to cell death. More importantly, there are bidirectional interactions between calcium and ROS, and such interactions may have therapeutic implications in treating HF through finely tuning the balance between these two by certain drugs. Many naturally derived products, e.g., flavones and isoflavones, have been shown to possess activities in regulating calcium and ROS simultaneously, thereby leading to a balanced microenvironment in heart tissues to exert therapeutic efficacies in HF. In this mini review, we aimed to provide an updated knowledge of the interplay between calcium and ROS in the development of HF. In addition, we summarized the recent studies (in vitro, in vivo and in clinical trials) using natural isolated flavones and isoflavones in treating HF. Critical challenges are also discussed. The information collected may help to evoke multidisciplinary efforts in developing novel agents for the potential prevention and treatment of HF.

Understanding of the drug resistance mechanism of hepatitis C virus NS3/4A to paritaprevir due to D168N/Y mutations: A molecular dynamics simulation perspective

Hepatitis C virus (HCV) NS3/4A protease is an attractive target for the development of antiviral therapy. However, the evolution of antiviral drug resistance is a major problem for treatment of HCV infected patients. Understanding of drug-resistance mechanisms at molecular level is therefore very important for the guidance of further design of antiviral drugs with high efficiency and specificity. Paritaprevir is a potent inhibitor against HCV NS3/4A protease genotype 1a. Unfortunately, this compound is highly susceptible to the substitution at D168 in the protease. In this work, molecular dynamics simulations of paritaprevir complexed with wild-type (WT) and two mutated strains (D168 N and D168Y) were carried out. Due to such mutations, the inhibitor-protein hydrogen bonding between them was weakened and the salt-bridge network among residues R123, R155 and D168 responsible for inhibitor binding was disrupted. Moreover, the per-residue free energy decomposition suggested that the main contributions from key residues such as Q80, V132, K136, G137 and R155 were lost in the D168 N/Y mutations. These lead to a lower binding affinity of paritaprevir for D168 N/Y variants of the HCV NS3/4A protease, consistent with the experimental data. This detailed information could be useful for further design of high potency anti-HCV NS3/4A inhibitors.

Exploring PfDHFR reaction surface: A combined molecular dynamics and QM/MM analysis

The substrate to the enzyme PfDHFR (Plasmodium falciparum Dihydrofolate Reductase) is a small molecule dihydrofolate (DHF), it gets converted to tetrahydrofolate (THF) in the active site of the enzyme. The PfDHFR reaction surface involves the protonation of DHF to DHFP as an initial step before the catalytic conversion. The binding affinities of all these species (DHF, DHFP and THF) contribute to the mechanism of DHFR catalytic action. Molecular dynamics (MD) simulations and Quantum Mechanics/Molecular Mechanics (QM/MM) analysis were performed to evaluate the binding affinity and molecular recognition interactions of the substrate DHF/DHFP and the product THF, in the active site of wild-type PfDHFR (wtPfDHFR). The binding affinities of the cofactor NADPH/NADP⁺ were also estimated in all the three complexes. The molecular dynamics (MD) simulations of the substrate, product and cofactor in the cavities of wtPfDHFR revealed the variation of the atomic level interactions during the course of the catalytic conversion. It was found that the DHFP binds very strongly to the PfDHFR active site and pulls the cofactor NADPH closer to itself. The QM/MM analysis revealed that the binding energy of DHFP (-59.82 kcal/mol) and NADPH (-100.24 kcal/mol) in DHFP-wtPfDHFR complex, is higher in comparison to the binding energy of DHF (-38.67 kcal/mol) and NADPH (-77.53 kcal/mol) in DHF-wtPfDHFR complex and the binding energy of THF (-30.72 kcal/mol) and NADP⁺ (-73.72 kcal/mol) in THF-wtPfDHFR complex. The hydride ion donor-acceptor distance (DAD) analysis was also carried out. This combined MD and QM/MM analysis revealed that the protonation of DHF increases the proximity between the substrate and the cofactor, thus facilitates the reaction profile of PfDHFR.

Mycobacterium tuberculosis serine/threonine protein kinases: structural information for the design of their specific ATP-competitive inhibitors

In the last decades, human protein kinases (PKs) have been relevant as targets in the development of novel therapies against many diseases, but the study of Mycobacterium tuberculosis PKs (MTPKs) involved in tuberculosis pathogenesis began much later and has not yet reached an advanced stage of development. To increase knowledge of these enzymes, in this work we studied the structural features of MTPKs, with focus on their ATP-binding sites and their interactions with inhibitors. PknA, PknB, and PknG are the most studied MTPKs, which were previously crystallized; ATP-competitive inhibitors have been designed against them in the last decade. In the current work, reported PknA, PknB, and PknG inhibitors were extracted from literature and their orientations inside the ATP-binding site were proposed by using docking method. With this information, interaction fingerprints were elaborated, which reveal the more relevant residues for establishing chemical interactions with inhibitors. The non-crystallized MTPKs PknD, PknF, PknH, PknJ, PknK, and PknL were also studied; their three-dimensional structural models were developed by using homology modeling. The main characteristics of MTPK ATP-binding sites (the non-crystallized and crystallized MTPKs, including PknE and PknI) were accounted; schemes of the main polar and nonpolar groups inside their ATP-binding sites were constructed, which are suitable for a major understanding of these proteins as antituberculotic targets. These schemes could be used for establishing comparisons between MTPKs and human PKs in order to increase selectivity of MTPK inhibitors. As a key tool for guiding medicinal chemists interested in the design of novel MTPK inhibitors, our work provides a map of the structural elements relevant for the design of more selective ATP-competitive MTPK inhibitors.

C-glycosyl flavone from Urginea indica inhibits proliferation & angiogenesis & induces apoptosis via cyclin-dependent kinase 6 in human breast, hepatic & colon cancer cell lines

Background & objectives:

Search for novel compounds beneficial to the treatment of cancer attracts a great deal of attention. We earlier demonstrated the isolation of 5,7-dihydroxy-2-[4’-hydroxy-3’-(methoxymethyl)phenyl]-6-C-β-glucopyranosyl flavone, a novel C-glycosyl flavone from Urginea indica bulb. The present study was undertaken to investigate the effect of this novel compound on human normal epithelial and breast, hepatic and colon cancer cell lines.

Methods:

The maximum non-toxic concentration (MNTC) and cytotoxicity of C-glycosyl flavone were assayed by 3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide (MTT). Cell cycle was analyzed by flow cytometry. Docking studies were performed to predict possible targets. Levels of cyclin-dependent kinase 1 (CDK1) and CDK6, Bcl2 and BAX and cytochrome c were quantified by specific ELISA. Mitochondrial membrane potential was determined using JC-1 dye. Apoptosis was quantified by Annexin V ELISA method.

Results:

Flow cytometry analysis demonstrated G0/G1 arrest. In silico docking studies predicted CDK1 and CDK6 as a possible target of C-glycosyl flavone. In vitro study confirmed CDK6 as the main target in C-glycosyl flavone-treated cancer cell lines. C-glycosyl flavone treatment also induced membrane blebbing, chromatin fragmentation and nucleosome formation. C-glycosyl flavone treatment caused marked loss of mitochondrial membrane potential, decrease in Bcl2/BAX ratio and activation of caspase-3 and release of caspase-9 and cytochrome c. In addition, C-glycosyl flavone inhibited the tumour-induced angiogenesis and reduced the vascular endothelial growth factor levels. Similarly, CDK6 inhibitor significantly inhibited proliferation and angiogenesis and induced apoptosis in tested cell lines.

Interpretation & conclusions:

The results indicate that C-glycosyl flavone may exert induction of apoptosis, cell cycle arrest and inhibition of angiogenesis via CDK6. Thus, targeting CDK6 using C-glycosyl flavone may serve as a novel therapeutic approach for the treatment of breast, hepatic and colon cancers.

Deciphering the binding behavior of flavonoids to the cyclin dependent kinase 6/cyclin D complex

Flavonoids, a class of natural compounds with variable phenolic structures, have been found to possess anti-cancer activities by modulating different enzymes and receptors like CDK6. To understand the binding behavior of flavonoids that inhibit the active CDK6, molecular dynamics (MD) simulations were performed on six inhibitors, chrysin (M01), fisetin (M03), galangin (M04), genistein (M05), quercetin (M06) and kaempferol (M07), complexed with CDK6/cyclin D. For all six flavonoids, the 3'-OH and 4'-OH of B-ring were found to be favorable for hydrogen bond formation, but the 3-OH on the C-ring and 5-OH on the A-ring were unfavorable, which were confirmed by the MD simulation results of the test molecule, 3', 4', 7-trihydroxyflavone (M15). The binding efficiencies of flavonoids against the CDK6/cyclin D complex were mainly through the electrostatic (especially the H-bond force) and vdW interactions with residues ILE19, VAL27, ALA41, GLU61, PHE98, GLN103, ASP163 and LEU152. The order of binding affinities of these flavonoids toward the CDK6/cyclin D was M03 > M01 > M07 > M15 > M06 > M05 > M04. It is anticipated that the binding features of flavonoid inhibitors studied in the present work may provide valuable insights for the development of CDK6 inhibitors.

Molecular insights into inclusion complexes of mansonone E and H enantiomers with various β-cyclodextrins

The structural dynamics and stability of inclusion complexes of mansonone E (ME) and H (MH) including their stereoisomers with various βCDs (methylated- and hydroxypropylated-βCDs) were investigated by classical molecular dynamics (MD) simulations and binding free energy calculations. The simulation results revealed that mansonones are able to form inclusion complexes with βCDs. The guest molecules are not completely inserted into the host cavity, their preferably positions are nearby the secondary rim with the oxane ring dipping into the hydrophobic inner cavity. The encapsulation process leads to a higher rigidity of the βCDs enhancing the intramolecular hydrogen bond formation ability and decreasing the chance of glucopyranose rotation. According to the MM-PBSA binding free energy calculation, all considered inclusion complexes are stable and the binding energies are mainly caused by van der Waals interactions. Moreover, the free energy calculations showed significant differences in the complexation energies for the stereoisomers, which could enable the separation of the isomers by analytical techniques for further pharmaceutical applications.

Effect of D168V mutation in NS3/4A HCV protease on susceptibilities of faldaprevir and danoprevir

Hepatitis C virus (HCV) is a serious cause of liver inflammation, cirrhosis and the development of hepatocellular carcinoma. Its NS3/4A serine protease functions to cleave a specific peptide bond, which is an important step in HCV replication. Thus the NS3/4A protease has become one of the main drug-targets in the design and development of anti-HCV agents. Unfortunately, high mutation rates in HCV have been reported due to the lack of RNA proofreading activity resulting in drug resistance. Herein, all-atom molecular dynamics simulations were employed to understand and illustrate the effects of the NS3/4A D168V mutation on faldaprevir (FDV) and danoprevir (DNV) binding efficiency. The D168V mutation was shown to interrupt the hydrogen bonding network of Q80R155D168R123 embedded in the extended S2 and partial S4 subsites of the NS3 protein and as a result the R123 side chain was displaced and moved out from the binding pocket. By means of MM/PBSA and MM/GBSA binding free energy calculations, the FDV and DNV binding affinities were shown to be significantly reduced by ∼10-15 kcal mol^-1 and ∼4-9 kcal mol^-1 relative to the wild-type complexes, respectively, which somewhat agrees with the experimental resistance folds.

Palbociclib can overcome mutations in cyclin dependent kinase 6 that break hydrogen bonds between the drug and the protein

Inhibition of cyclin dependent kinases (CDKs) 4 and 6 prevent cells from entering the synthesis phase of the cell cycle. CDK4 and 6 are therefore important drug targets in various cancers. The selective CDK4/6 inhibitor palbociclib is approved for the treatment of breast cancer and has shown activity in a cellular model of mixed lineage leukaemia (MLL)‐rearranged acute myeloid leukaemia (AML). We studied the interactions of palbociclib and CDK6 using molecular dynamics simulations. Analysis of the simulations suggested several interactions that stabilized the drug in its binding site and that were not observed in the crystal structure of the protein‐drug complex. These included a hydrogen bond to His 100 that was hitherto not reported and several hydrophobic contacts. Evolutionary‐based bioinformatic analysis was used to suggest two mutants, D163G and H100L that would potentially yield drug resistance, as they lead to loss of important protein–drug interactions without hindering the viability of the protein. One of the mutants involved a change in the glycine of the well‐conserved DFG motif of the kinase. Interestingly, CDK6‐dependent human AML cells stably expressing either mutant retained sensitivity to palbociclib, indicating that the protein‐drug interactions are not affected by these. Furthermore, the cells were proliferative in the absence of palbociclib, indicating that the Asp to Gly mutation in the DFG motif did not interfere with the catalytic activity of the protein.

PDB Code(s): 2EUF

Prediction of luciferase inhibitors by the high-performance MIEC-GBDT approach based on interaction energetic patterns

The MIEC-GBDT model can be used as a powerful tool to identify potential interference compounds in luciferase-based high-throughput screening.

Role of R292K mutation in influenza H7N9 neuraminidase toward oseltamivir susceptibility: MD and MM/PB(GB)SA study

The H7N9 avian influenza virus is a novel re-assortment from at least four different strains of virus. Neuraminidase, which is a glycoprotein on the surface membrane, has been the target for drug treatment. However, some H7N9 strains that have been isolated from patient after drug treatment have a R292K mutation in neuraminidase. This substitution was found to facilitate drug resistance using protein- and virus- assays, in particular it gave a high resistance to the most commonly used drug, oseltamivir. The aim of this research is to understand the source of oseltamivir resistance using MD simulations and the MM/PB(GB)SA binding free energy approaches. Both methods can predict the reduced susceptibility of oseltamivir in good agreement to the IC ₅₀ binding energy, although MM/GBSA underestimates this prediction compared to the MM/PBSA calculation. Electrostatic interaction is the main contribution for oseltamivir binding in terms of both interaction and solvation. We found that the source of the drug resistance is a decrease in the binding interaction combined with the reduction of the dehydration penalty. The smaller K292 mutated residue has a larger binding pocket cavity compared to the wild-type resulting in the loss of drug carboxylate-K292 hydrogen bonding and an increased accessibility for water molecules around the K292 mutated residue. In addition, oseltamivir does not bind well to the R292K mutant complex as shown by the high degree of fluctuation in ligand RMSD during the simulation and the change in angular distribution of bulky side chain groups.

Two novel antimicrobial defensins from rice identified by gene coexpression network analyses

Defensins form an antimicrobial peptides (AMP) family, and have been widely studied in various plants because of their considerable inhibitory functions. However, their roles in rice (Oryza sativa L.) have not been characterized, even though rice is one of the most important staple crops that is susceptible to damaging infections. Additionally, a previous study identified 598 rice genes encoding cysteine-rich peptides, suggesting there are several uncharacterized AMPs in rice. We performed in silico gene expression and coexpression network analyses of all genes encoding defensin and defensin-like peptides, and determined that OsDEF7 and OsDEF8 are coexpressed with pathogen-responsive genes. Recombinant OsDEF7 and OsDEF8 could form homodimers. They inhibited the growth of the bacteria Xanthomonas oryzae pv. oryzae, X. oryzae pv. oryzicola, and Erwinia carotovora subsp. atroseptica with minimum inhibitory concentration (MIC) ranging from 0.6 to 63μg/mL. However, these OsDEFs are weakly active against the phytopathogenic fungi Helminthosporium oryzae and Fusarium oxysporum f.sp. cubense. This study describes a useful method for identifying potential plant AMPs with biological activities.

Screening of hepatitis C NS5B polymerase inhibitors containing benzothiadiazine core: a steered molecular dynamics

Hepatic C virus (HCV) is a global health problem, resulting in liver cirrhosis and inflammation that can develop to hepatocellular carcinoma and fatality. The NS5B polymerase of HCV plays an important role in viral RNA replication process, making it an attractive therapeutic target for design and development of anti-HCV drugs. To search new potent compounds against the HCV NS5B polymerase, the molecular docking and the steered molecular dynamics (SMD) simulation techniques were performed. The potential potent inhibitors of the NS5B polymerase were screened out from the ZINC database using structural similarity search and molecular docking technique. Five top-hit compounds (the ZINC compounds 49888724, 49054741, 49777239, 49793673, and 49780355) were then studied by the SMD simulations based on the hypothesis that a high rupture force relates to a high binding efficiency. The results demonstrated that the ZINC compound 49888724 had a greater maximum rupture force, reflecting a good binding strength and inhibitory potency than known inhibitors and the rest four ZINC compounds. Therefore, our finding indicated that the ZINC compound 49888724 is a potential candidate to be a novel NS5B inhibitor for further design. Besides, the van der Waals interaction could be considered as the main contribution for stabilizing the NS5B-ligand complex.

Exploring the molecular targets of dietary flavonoid fisetin in cancer

The last few decades have seen a resurgence of interest among the scientific community in exploring the efficacy of natural compounds against various human cancers. Compounds of plant origin belonging to different groups such as alkaloids, flavonoids and polyphenols evaluated for their cancer preventive effects have yielded promising data, thereby offering a potential therapeutic alternative against this deadly disease. The flavonol fisetin (3,3',4',7-tetrahydroxyflavone), present in fruits and vegetables such as strawberries, apple, cucumber, persimmon, grape and onion, was shown to possess anti-microbial, anti-inflammatory, anti-oxidant and more significantly anti-carcinogenic activity when assessed in diverse cell culture and animal model systems. The purpose of this review is to update and discuss key findings obtained till date from in vitro and in vivo studies on fisetin, with special focus on its anti-cancer role. The molecular mechanism(s) described in the observed growth inhibitory effects of fisetin in different cancer cell types is also summarized. Moreover, an attempt is made to delineate the direction of future studies that could lead to the development of fisetin as a potent chemopreventive/chemotherapeutic agent against cancer.

Physical properties and biological activities of hesperetin and naringenin in complex with methylated β-cyclodextrin

The aim of this work is to improve physical properties and biological activities of the two flavanones hesperetin and naringenin by complexation with β-cyclodextrin (β-CD) and its methylated derivatives (2,6-di-O-methyl-β-cyclodextrin, DM-β-CD and randomly methylated-β-CD, RAMEB). The free energies of inclusion complexes between hesperetin with cyclodextrins (β-CD and DM-β-CD) were theoretically investigated by molecular dynamics simulation. The free energy values obtained suggested a more stable inclusion complex with DM-β-CD. The vdW force is the main guest–host interaction when hesperetin binds with CDs. The phase solubility diagram showed the formation of a soluble complex of A_L type, with higher increase in solubility and stability when hesperetin and naringenin were complexed with RAMEB. Solid complexes were prepared by freeze-drying, and the data from differential scanning calorimetry (DSC) confirmed the formation of inclusion complexes. The data obtained by the dissolution method showed that complexation with RAMEB resulted in a better release of both flavanones to aqueous solution. The flavanones-β-CD/DM-β-CD complexes demonstrated a similar or a slight increase in anti-inflammatory activity and cytotoxicity towards three different cancer cell lines. The overall results suggested that solubilities and bioactivities of both flavanones were increased by complexation with methylated β-CDs.

Co-solvation effect on the binding mode of the α-mangostin/β-cyclodextrin inclusion complex

Cyclodextrins (CDs) have been extensively utilized as host molecules to enhance the solubility, stability and bioavailability of hydrophobic drug molecules through the formation of inclusion complexes. It was previously reported that the use of co-solvents in such studies may result in ternary (host:guest:co-solvent) complex formation. The objective of this work was to investigate the effect of ethanol as a co-solvent on the inclusion complex formation between α-mangostin (α-MGS) and β-CD, using both experimental and theoretical studies. Experimental phase-solubility studies were carried out in order to assess complex formation, with the mechanism of association being probed using a mathematical model. It was found that α-MGS was poorly soluble at low ethanol concentrations (0–10% v/v), but higher concentrations (10–40% v/v) resulted in better α-MGS solubility at all β-CD concentrations studied (0–10 mM). From the equilibrium constant calculation, the inclusion complex is still a binary complex (1:1), even in the presence of ethanol. The results from our theoretical study confirm that the binding mode is binary complex and the presence of ethanol as co-solvent enhances the solubility of α-MGS with some effects on the binding affinity with β-CD, depending on the concentration employed.

Inclusion complexation of pinostrobin with various cyclodextrin derivatives

Pinostrobin (PNS) is one of the important flavonoids and can be abundantly found in the rhizomes of fingerroot (Boesenbergia rotrunda) and galangal (Alpinia galangal and Alpinia officinarum), the herbal basis of Southeast Asian cooking. Similar to other flavonoids, PNS exhibits anti-oxidative, anti-inflammatory and anti-cancer properties. However, this compound has an extremely low water solubility that limits its use in pharmaceutical applications. Beta-cyclodextrin (βCD) and its derivatives, 2,6-dimethyl-βCD (2,6-DMβCD) and the three hydroxypropyl-βCDs (2-HPβCD, 6-HPβCD and 2,6-DHPβCD), have unique properties that enhance the stability and solubility of such low-soluble guest molecules. In the present study, molecular dynamics simulations were applied to investigate the dynamics and stability of PNS inclusion complexes with βCD and its derivatives (2,6-DMβCD, 2,6-DHPβCD, 2-HPβCD and 6-HPβCD). PNS was able to form complexes with βCD and all four of its derivatives by either the chromone (C-PNS) or phenyl (P-PNS) ring dipping toward the cavity. According to the molecular mechanics-generalized Born surface area binding free energy values, the stability of the different PNS/βCD complexes was ranked as 2,6-DHPβCD>2,6-DMβCD>2-HPβCD>6-HPβCD>βCD. These theoretical results were in good agreement with the stability constants that had been determined by the solubility method.

Binding Activity Prediction of Cyclin-Dependent Inhibitors

The Cyclin-Dependent Kinases (CDKs) are the core components coordinating eukaryotic cell division cycle. Generally the crystal structure of CDKs provides information on possible molecular mechanisms of ligand binding. However, reliable and robust estimation of ligand binding activity has been a challenging task in drug design. In this regard, various machine learning techniques, such as Support Vector Machine, Naive Bayesian classifier, Decision Tree, and K-Nearest Neighbor classifier, have been used. The performance of these heterogeneous classification techniques depends on proper selection of features from the data set. This fact motivated us to propose an integrated classification technique using Genetic Algorithm (GA), Rotational Feature Selection (RFS) scheme, and Ensemble of Machine Learning methods, named as the Genetic Algorithm integrated Rotational Ensemble based classification technique, for the prediction of ligand binding activity of CDKs. This technique can automatically find the important features and the ensemble size. For this purpose, GA encodes the features and ensemble size in a chromosome as a binary string. Such encoded features are then used to create diverse sets of training points using RFS in order to train the machine learning method multiple times. The RFS scheme works on Principal Component Analysis (PCA) to preserve the variability information of the rotational nonoverlapping subsets of original data. Thereafter, the testing points are fed to the different instances of trained machine learning method in order to produce the ensemble result. Here accuracy is computed as a final result after 10-fold cross validation, which also used as an objective function for GA to maximize. The effectiveness of the proposed classification technique has been demonstrated quantitatively and visually in comparison with different machine learning methods for 16 ligand binding CDK docking and rescoring data sets. In addition, the best possible features have been reported for CDK docking and rescoring data sets separately. Finally, the Friedman test has been conducted to judge the statistical significance of the results produced by the proposed technique. The results indicate that the integrated classification technique has high relevance in predicting of protein-ligand binding activity.

Structural insight of DNA topoisomerases I from camptothecin-producing plants revealed by molecular dynamics simulations

DNA topoisomerase I (Top1) catalyzes changes in DNA topology by cleaving and rejoining one strand of the double stranded (ds)DNA. Eukaryotic Top1s are the cellular target of the plant-derived anticancer indole alkaloid camptothecin (CPT), which reversibly stabilizes the Top1-dsDNA complex. However, CPT-producing plants, including Camptotheca acuminata, Ophiorrhiza pumila and Ophiorrhiza liukiuensis, are highly resistant to CPT because they possess point-mutated Top1. Here, the adaptive convergent evolution is reported between CPT production ability and mutations in their Top1, as a universal resistance mechanism found in all tested CPT-producing plants. This includes Nothapodytes nimmoniana, one of the major sources of CPT. To obtain a structural insight of the resistance mechanism, molecular dynamics simulations of CPT- resistant and -sensitive plant Top1s complexed with dsDNA and topotecan (a CPT derivative) were performed, these being compared to that for the CPT-sensitive human Top1. As a result, two mutations, Val617Gly and Asp710Gly, were identified in O. pumila Top1 and C. acuminata Top1, respectively. The substitutions at these two positions, surprisingly, are the same as those found in a CPT derivative-resistant human colon adenocarcinoma cell line. The results also demonstrated an increased linker flexibility of the CPT-resistant Top1, providing an additional explanation for the resistance mechanism found in CPT-producing plants. These mutations could reflect the long evolutionary adaptation of CPT-producing plant Top1s to confer a higher degree of resistance.

Binding mode and free energy prediction of fisetin/β-cyclodextrin inclusion complexes

In the present study, our aim is to investigate the preferential binding mode and encapsulation of the flavonoid fisetin in the nano-pore of β-cyclodextrin (β-CD) at the molecular level using various theoretical approaches: molecular docking, molecular dynamics (MD) simulations and binding free energy calculations. The molecular docking suggested four possible fisetin orientations in the cavity through its chromone or phenyl ring with two different geometries of fisetin due to the rotatable bond between the two rings. From the multiple MD results, the phenyl ring of fisetin favours its inclusion into the β-CD cavity, whilst less binding or even unbinding preference was observed in the complexes where the larger chromone ring is located in the cavity. All MM- and QM-PBSA/GBSA free energy predictions supported the more stable fisetin/β-CD complex of the bound phenyl ring. Van der Waals interaction is the key force in forming the complexes. In addition, the quantum mechanics calculations with M06-2X/6-31G(d,p) clearly showed that both solvation effect and BSSE correction cannot be neglected for the energy determination of the chosen system.

Investigations on the mechanisms of interactions between matrix metalloproteinase 9 and its flavonoid inhibitors using a combination of molecular docking, hybrid quantum mechanical/molecular mechanical calculations, and molecular dynamics simulations

Many experimental studies have found that flavonoids can inhibit the activities of matrix metalloproteinases (MMPs), but the relevant mechanisms are still unclear. In this paper, the interaction mechanisms of MMP-9 with its five flavonoid inhibitors are investigated using a combination of molecular docking, hybrid quantum mechanical and molecular mechanical (QM/MM) calculations, and molecular dynamics simulations. The molecular dynamics simulation results show a good linear correlation between the calculated binding free energies of QM/MM−Poisson–Boltzmann surface area (PBSA) and the experimental −log(EC50) regarding the studied five flavonoids on MMP-9 inhibition in explicit solvent. It is found that compared with the MM−PBSA method, the QM/MM−PBSA method can obviously improve the accuracy for the calculated binding free energies. The predicted binding modes of the five flavonoid−MMP-9 complexes reveal that the different hydrogen bond networks can form besides producing the Zn−O coordination bonds, which can reasonably explain previous experimental results. The agreement between our calculated results and the previous experimental facts indicates that the force field parameters used here are effective and reliable for investigating the systems of flavonoid−MMP-9 interactions, and thus, these simulations and analyses could be reproduced for the other related systems involving protein−ligand interactions. This paper may be helpful for designing the new MMP-9 inhibitors having higher biological activities by carrying out the structural modifications of flavonoid molecules.

Molecular dynamics of interactions between rigid and flexible antifolates and dihydrofolate reductase from pyrimethamine-sensitive and pyrimethamine-resistant Plasmodium falciparum

Currently, the usefulness of antimalarials such as pyrimethamine (PYR) is drastically reduced due to the emergence of resistant Plasmodium falciparum (Pf) caused by its dihydrofolate reductase (PfDHFR) mutations, especially the quadruple N51I/C59R/S108N/I164L mutations. The resistance was due to the steric conflict of PYR with S108N. WR99210 (WR), a dihydrotriazine antifolate with a flexible side chain that can avoid such conflict, can overcome this resistance through tight binding with the mutant. To understand factors contributing to different binding affinities of PYR/WR to the wild type (WT) and quadruple mutant (QM), we performed simulations on WR-WT, WR-QM, PYR-WT, and PYR-QM complexes and found that Ile14 and Asp54 were crucial for PYR/WR binding to PfDHFR due to strong hydrogen bonds. The quadruple mutations cause PYR to form, on average, fewer hydrogen bonds with Ile14 and Leu164, and to be displaced from its optimal orientation for Asp54 interaction. The predicted binding affinity ranking (WR-QM ≈ WR-WT ≈ PYR-WT >> PYR-QM) reasonably agrees with the inhibition constant (K(i)) ranking. Our results reveal important residues for tight binding of PYR/WR to WT/QM, which may be used to evaluate the inhibition effectiveness of antimalarials and to provide fundamental information for designing new drugs effective against drug-resistant P. falciparum.

Key binding and susceptibility of NS3/4A serine protease inhibitors against hepatitis C virus

Hepatitis C virus (HCV) causes an infectious disease that manifests itself as liver inflammation, cirrhosis, and can lead to the development of liver cancer. Its NS3/4A serine protease is a potent target for drug design and development since it is responsible for cleavage of the scissile peptide bonds in the polyprotein important for the HCV life cycle. Herein, the ligand-target interactions and the binding free energy of the four current NS3/4A inhibitors (boceprevir, telaprevir, danoprevir, and BI201335) were investigated by all-atom molecular dynamics simulations with three different initial atomic velocities. The per-residue free energy decomposition suggests that the key residues involved in inhibitor binding were residues 41-43, 57, 81, 136-139, 155-159, and 168 in the NS3 domain. The van der Waals interactions yielded the main driving force for inhibitor binding at the protease active site for the cleavage reaction. In addition, the highest number of hydrogen bonds was formed at the reactive P1 site of the four studied inhibitors. Although the hydrogen bond patterns of these inhibitors were different, their P3 site was most likely to be recognized by the A157 backbone. Both molecular mechanic (MM)/Poisson-Boltzmann surface area and MM/generalized Born surface area approaches predicted the relative binding affinities of the four inhibitors in a somewhat similar trend to their experimentally derived biological activities.

Enhanced stability of a naringenin/2,6-dimethyl β-cyclodextrin inclusion complex: molecular dynamics and free energy calculations based on MM- and QM-PBSA/GBSA

The structure, dynamic behavior and binding affinity of the inclusion complexes between naringenin and the two cyclodextrins (CDs), β-CD and its 2,6-dimethyl derivative (DM-β-CD), were theoretically studied by multiple molecular dynamics simulations and free energy calculations. Naringenin most likely prefers to bind with CDs through the phenyl ring. Although a lower hydrogen bond formation of naringenin with the 3-hydroxyl group of DM-β-CD (relative to β-CD) was observed, the higher cavity could encapsulate almost the whole naringenin molecule. In contrast for the naringenin/β-CD complex, the phenyl ring feasibly passed through the primary rim resulting in the chromone ring binding inside instead. MM-PBSA/GBSA and QM-PBSA/GBSA binding free energies strongly suggested a greater stability of the naringenin/DM-β-CD inclusion complex. Van der Waals force played an important role as the key guest-host interaction for the complexation between naringenin and each cyclodextrin.

Predicting the mixture effects of three pesticides by integrating molecular simulation with concentration addition modeling

Molecular simulation techniques are used to identify the mode of inhibition of chemicals at the ligand–receptor level.

Related flavonoids cause cooperative inhibition of the sarcoplasmic reticulum Ca²⁺ ATPase by multimode mechanisms

Flavonoids are group of plant-derived hydroxylated polycyclic molecules found in fruit and vegetables. They are known to bio-accumulate within humans and are considered to have beneficial health effects, including cancer chemoprotection. One mechanism proposed to explain this is that they are able to induce apoptosis in cancer cells by inhibiting a variety of kinases and also the Ca²⁺ ATPase. An investigation was undertaken with respect to the mechanism of inhibition for three flavonoids: quercetin, galangin and 3,6 dihydroxyflavone (3,6-DHF). Each inhibited the Ca²⁺ ATPase with K(i) values of 8.7, 10.3 and 5.4 μM, respectively, showing cooperative inhibition with n ~ 2. Given their similar structures, the flavonoids showed several differences in their mechanisms of inhibition. All three flavonoids stabilized the ATPase in the E₁ conformation and reduced [³²P]-ATP binding. However, both galangin and 3,6-DHF increased the affinity of Ca²⁺ for the ATPase by decreasing the Ca²⁺-dissociation rate constant, whereas quercetin had little effect. Ca²⁺-induced changes in tryptophan fluorescence levels were reduced in the presence of 3,6-DHF and galangin (but not with quercetin), indicating that Ca²⁺-associated changes within the transmembrane helices are altered. Both galangin and quercetin reduced the rates of ATP-dependent phosphorylation and dephosphorylation, whereas 3,6-DHF did not. Modelling studies suggest that flavonoids could potentially bind to two sites: one directly where nucleotides bind within ATP binding site and the other at a site close by. We hypothesize that interactions of these two neighbouring sites may account for both the cooperative inhibition and the multimode mechanisms of action seen with related flavonoids.

Computational study of the effects of protein tyrosine nitrations on the catalytic activity of human thymidylate synthase

Tyrosine nitration is a widespread post-translational modification capable of affecting both the function and structure of the host protein molecule. Enzyme thymidylate synthase (TS), a homodimer, is a molecular target for anticancer therapy. Recently purified TS preparations, isolated from mammalian tissues, were found to be nitrated, suggesting this modification to appear endogenously in normal and tumor tissues. Moreover, human TS (hTS) nitration in vitro led to a by twofold lowered catalytic activity following nitration in average of 1 tyrosine residue per monomer (Dąbrowska-Maś et al. in Org Biomol Chem 10:323-331, 2012), with the modification identified by mass spectrometry at seven different sites (Y33, Y65, Y135, Y213, Y230, Y258 and Y301). In the present paper, combined computational approach, including molecular and essential dynamics and free energy computations, was used to predict the influence on the activity of hTS of nitration of each of the seven tyrosine residues. The simulations were based on the crystal structure of hTS ternary complex with dUMP and Tomudex (PDB code: 1I00), with the Tomudex molecule replaced by the molecule of TS cofactor analogue, tetrahydrofolate. The present results indicate that while with nitration of five out of seven residues (Y33, Y135, Y230, Y258 and Y301), single residue modification appears to have a strong reducing effect on the activity, with the remaining two, Y65 and Y213, no or a weaker influence is apparent. Taken together, these results demonstrate that tyrosine nitrations in the hTS enzyme show clear tendency to influence the structure and dynamics and, in turn, catalytic properties of the host enzyme. These effects are overall distance-dependent.

Apigenin and its impact on gastrointestinal cancers

Apigenin (4',5,7-trihydroxyflavone, 5,7-dihydroxy-2-(4-hydroxyphenyl)-4H-1-benzopyran-4-one) is a flavonoid found in many fruits, vegetables, and herbs, the most abundant sources being the leafy herb parsley and dried flowers of chamomile. Present in dietary sources as a glycoside, it is cleaved in the gastrointestinal lumen to be absorbed and distributed as apigenin itself. For this reason, the epithelium of the gastrointestinal tract is exposed to higher concentrations of apigenin than tissues at other locations. This would also be true for epithelial cancers of the gastrointestinal tract. We consider the evidence for actions of apigenin that might hinder the ability of gastrointestinal cancers to progress and spread. Apigenin has been shown to inhibit cell growth, sensitize cancer cells to elimination by apoptosis, and hinder the development of blood vessels to serve the growing tumor. It also has actions that alter the relationship of the cancer cells with their microenvironment. Apigenin is able to reduce cancer cell glucose uptake, inhibit remodeling of the extracellular matrix, inhibit cell adhesion molecules that participate in cancer progression, and oppose chemokine signaling pathways that direct the course of metastasis into other locations. As such, apigenin may provide some additional benefit beyond existing drugs in slowing the emergence of metastatic disease.