Molecular dynamics (MD) simulations can provide not only plentiful dynamical structural information on biomacromolecules but also a wealth of energetic information about protein and ligand interactions. Such information is very important to understanding the structure-function relationship of the target and the essence of protein-ligand interactions and to guiding the drug discovery and design process. Thus, MD simulations have been applied widely and successfully in each step of modern drug discovery. Areas covered: In this review, the authors review the applications of MD simulations in novel drug discovery, including the pathogenic mechanisms of amyloidosis diseases, virtual screening and the interaction mechanisms between drugs and targets. Expert opinion: MD simulations have been used widely in investigating the pathogenic mechanisms of diseases caused by protein misfolding, in virtual screening, and in investigating drug resistance mechanisms caused by mutations of the target. These issues are very difficult to solve by experimental methods alone. Thus, in the future, MD simulations will have wider application with the further improvement of computational capacity and the development of better sampling methods and more accurate force fields together with more efficient analysis methods.

This study was aimed to design the first dual-target small-molecule inhibitor co-targeting poly (ADP-ribose) polymerase-1 (PARP1) and bromodomain containing protein 4 (BRD4), which had important cross relation in the global network of breast cancer, reflecting the synthetic lethal effect. A series of new BRD4 and PARP1 dual-target inhibitors were discovered and synthesized by fragment-based combinatorial screening and activity assays that together led to the chemical optimization. Among these compounds, 19d was selected and exhibited micromole enzymatic potencies against BRD4 and PARP1, respectively. Compound 19d was further shown to efficiently modulate the expression of BRD4 and PARP1. Subsequently, compound 19d was found to induce breast cancer cell apoptosis and stimulate cell cycle arrest at G1 phase. Following pharmacokinetic studies, compound 19d showed its antitumor activity in breast cancer susceptibility gene 1/2 ( BRCA1/2 ) wild-type MDA-MB-468 and MCF-7 xenograft models without apparent toxicity and loss of body weight. These results together demonstrated that a highly potent dual-targeted inhibitor was successfully synthesized and indicated that co-targeting of BRD4 and PARP1 based on the concept of synthetic lethality would be a promising therapeutic strategy for breast cancer.

BackgroundChinese traditional herbal medicine Fuzhengkangai (FZKA) formulation combination with gefitinib can overcome drug resistance and improve the prognosis of lung adenocarcinoma patients. However, the pharmacological and molecular mechanisms underlying the active ingredients, potential targets, and overcome drug resistance of the drug are still unclear. Therefore, it is necessary to explore the molecular mechanism of FZKA.MethodsA systems pharmacology and bioinformatics-based approach was employed to investigate the molecular pathogenesis of EGFR-TKI resistance with clinically effective herb formula. The differential gene expressions between EGFR-TKI sensitive and resistance cell lines were calculated and used to find overlap from targets as core targets. The prognosis of core targets was validated from the cancer genome atlas (TCGA) database by Cox regression. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment is applied to analysis core targets for revealing mechanism in biology.ResultsThe results showed that 35 active compounds of FZKA can interact with eight core targets proteins (ADRB2, BCL2, CDKN1A, HTR2C, KCNMA1, PLA2G4A, PRKCA and LYZ). The risk score of them were associated with overall survival and relapse free time (HR = 6.604, 95% CI: 2.314–18.850; HR = 5.132, 95% CI: 1.531–17.220). The pathway enrichment suggested that they involved in EGFR-TKI resistance and non-small cell lung cancer pathways, which directly affect EGFR-TKI resistance. The molecular docking showed that licochalcone a and beta-sitosterol can closely bind two targets (BCL2 and PRKCA) that involved in EGFR-TKI resistance pathway.ConclusionsThis study provided a workflow for understanding mechanism of CHM for against drug resistance.Electronic supplementary materialThe online version of this article (10.1186/s12906-018-2347-x) contains supplementary material, which is available to authorized users.

The reported crystal structures of β2 adrenergic receptor (β2AR) reveal that the open and closed states of the water channel are correlated with the inactive and active conformations of β2AR. However, more details about the process by which the water channel states are affected by the active to inactive conformational change of β2AR remain illusive. In this work, molecular dynamics simulations are performed to study the dynamical inactive and active conformational change of β2AR induced by inverse agonist ICI 118,551. Markov state model analysis and free energy calculation are employed to explore the open and close states of the water channel. The simulation results show that inverse agonist ICI 118,551 can induce water channel opening during the conformational transition of β2AR. Markov state model (MSM) analysis proves that the energy contour can be divided into seven states. States S1, S2 and S5, which represent the active conformation of β2AR, show that the water channel is in the closed state, while states S4 and S6, which correspond to the intermediate state conformation of β2AR, indicate the water channel opens gradually. State S7, which represents the inactive structure of β2AR, corresponds to the full open state of the water channel. The opening mechanism of the water channel is involved in the ligand-induced conformational change of β2AR. These results can provide useful information for understanding the opening mechanism of the water channel and will be useful for the rational design of potent inverse agonists of β2AR.

Recently, small-molecule compounds have been reported to block the PD-1/PD-L1 interaction by inducing the dimerization of PD-L1. All these inhibitors had a common scaffold and interacted with the cavity formed by two PD-L1 monomers. This special interactive mode provided clues for the structure-based drug design, however, also showed limitations for the discovery of small-molecule inhibitors with new scaffolds. In this study, we revealed the structure-activity relationship of the current small-molecule inhibitors targeting dimerization of PD-L1 by predicting their binding and unbinding mechanism via conventional molecular dynamics and metadynamics simulation. During the binding process, the representative inhibitors (BMS-8 and BMS-1166) tended to have a more stable binding mode with one PD-L1 monomer than the other and the small-molecule inducing PD-L1 dimerization was further stabilized by the non-polar interaction of Ile54, Tyr56, Met115, Ala121, and Tyr123 on both monomers and the water bridges involved in ALys124. The unbinding process prediction showed that the PD-L1 dimerization kept stable upon the dissociation of ligands. It's indicated that the formation and stability of the small-molecule inducing PD-L1 dimerization was the key factor for the inhibitory activities of these ligands. The contact analysis, R-group based quantitative structure-activity relationship (QSAR) analysis and molecular docking further suggested that each attachment point on the core scaffold of ligands had a specific preference for pharmacophore elements when improving the inhibitory activities by structural modifications. Taken together, the results in this study could guide the structural optimization and the further discovery of novel small-molecule inhibitors targeting PD-L1.

is an NAD + -dependent protein deacetylase localized primarily in the mitochondria with many links to different types of human cancers. Autophagy, which is a highly conserved lysosomal degradation process in eukaryotic cells, has been recently reported to be positively regulated by SIRT3 in cancer; therefore, activating SIRT3-modulated autophagy may be a promising strategy for drug discovery. In this study, we discovered a small-molecule activator of SIRT3 compound 33c (ADTL-SA1215) with specific SIRT3 deacetylase activity by structure-guided design and high-throughput screening. Subsequently, compound 33c inhibited the proliferation and migration of human breast carcinoma MDA-MB-231 cells by SIRT3-driven autophagy/mitophagy signaling pathways in vitro and in vivo. Collectively, these results demonstrate that pharmacological activation of SIRT3 is a potential therapeutic approach of triple negative breast cancer (TNBC). More importantly, compound 33c may be a first-in-class specific small-molecule activator of SIRT3 that would be utilized for future cancer drug development.

As co-chaperones of the 90-kDa heat shock protein(HSP90), FK506 binding protein 51 (FKBP51) and FK506 binding protein 52 (FKBP52) modulate the maturation of steroid hormone receptor through their specific FK1 domains (FKBP12-like domain 1). The inhibitors targeting FK1 domains are potential therapies for endocrine-related physiological disorders. However, the structural conservation of the FK1 domains between FKBP51 and FKBP52 make it difficult to obtain satisfactory selectivity in FK506-based drug design. Fortunately, a series of iFit ligands synthesized by Hausch et al exhibited excellent selectivity for FKBP51, providing new opportunity for design selective inhibitors. We performed molecular dynamics simulation, binding free energy calculation and unbinding pathway analysis to reveal selective mechanism for the inhibitor iFit4 binding with FKBP51 and FKBP52. The conformational stability evaluation of the "Phe67-in" and "Phe67-out" states implies that FKBP51 and FKBP52 have different preferences for "Phe67-in" and "Phe67-out" states, which we suggest as the determinant factor for the selectivity for FKBP51. The binding free energy calculations demonstrate that nonpolar interaction is favorable for the inhibitors binding, while the polar interaction and entropy contribution are adverse for the inhibitors binding.According to the results from binding free energy decomposition, the electrostatic difference of residue 85 causes the most significant thermodynamics effects on the binding of iFit4 to FKBP51 and FKBP52. Furthermore, the importance of substructure units on iFit4 were further evaluated by unbinding pathway analysis and residue-residue contact analysis between iFit4 and the proteins. The results will provide new clues for the design of selective inhibitors for FKBP51.

Xanthatin (XT), a naturally occurring sesquiterpene lactone presented in cocklebur ( Xanthium strumarium L.), is under development as a potential anticancer agent. Despite the promising anticancer effect of XT, the molecular mechanism underlying its cellular action has not been well elucidated. The mammalian thioredoxin reductase (TrxR) enzymes, the essential seleno-flavoproteins containing a penultimate selenocysteine (Sec) residue at the C-terminus, represent a promising target for cancer chemotherapeutic agents. In this study, XT inhibits both the purified TrxR and the enzyme in cells. The possible binding mode of XT with the TrxR protein is predicted by the covalent docking method. Mechanism studies reveal that XT targets the Sec residue of TrxR and inhibits the enzyme activity irreversibly. Simultaneously, the inhibition of TrxR by XT promotes the oxidative stress-mediated apoptosis of HeLa cells. Importantly, the knockdown of the enzyme sensitizes the cells to XT treatment. Targeting TrxR thus discloses a novel molecular mechanism in accounting for the cellular action of XT and provides insights into the development of XT as an anticancer agent.