Bromodomain结构域蛋白4(bromodomain-containing protein 4,BRD4)已成为治疗多种疾病药物设计的重要靶标. 最近在实验上发现了几种有效的靶向BRD4的抑制剂,但具体的抑制机理尚不清楚.此工作采用分子动力学模拟,动态相关性分析和结合自由能计算研究抑制剂8Q9和8QC与BRD4(1)的结合模式.分子动力学分析表明抑制剂结合对BRD4(1)的结构柔性产生重大影响.同时动态相关性分析进一步表明抑制剂结合极大地改变了BRD4(1)的运动模式.结合自由能计算结果表明范德华相互作用是抑制剂与BRD4(1)结合的主要驱动力.采用基于残基的自由能分解方法评估了分离残基对抑制剂结合的贡献,数据表明氢键相互作用和疏水相互作用是影响抑制剂与BRD4(1)结合的关键因素.本研究有望为设计和开发靶向BRD4的抑制剂提供有意义的理论指导. 相似文献
Inhibition of the p53–MDM2 interaction is a new therapeutic strategy to activate the wild-type function of p53 in tumors. Molecular dynamics (MD) simulations and calculations of binding free energies were performed to investigate the binding mechanisms of p53 and two inhibitors PMI and VZV to MDM2. The results show that van der Waals interaction is the main force to control the bindings of ligands to the hydrophobic cleft of MDM2, which basically agrees with the previous calculated and experimental studies. The results from the RMSF calculation, cross-correlation analysis and principal component (PC) analysis prove that the ligand bindings produce a significant effect on the conformation of the binding cleft of MDM2. In addition, the calculations of residue-based free energy decomposition suggest that the CH–CH, CH–π, and π–π interactions dominate the bindings of p53 and inhibitors to MDM2. This study can provide significant help for the design of potent inhibitors targeting the p53–MDM2 interaction. 相似文献
The C-C chemokine receptor 2 (CCR2) was proved as a multidrug target in many diseases like diabetes, inflammation and AIDS, but rational drug design on this target is still lagging behind as the information on the exact binding site and the crystal structure is not yet available. Therefore, for a successful structure-based drug design, an accurate receptor model in ligand-bound state is necessary. In this study, binding-site residues of CCR2 was determined using in silico alanine scanning mutagenesis and the interactions between TAK-779 and the developed homology model of CCR2. Molecular dynamic simulation and Molecular Mechanics-Generalized Born Solvent Area method was applied to calculate binding free energy difference between the template and mutated protein. Upon mutating 29 amino acids of template protein and comparison of binding free energy with wild type, six residues were identified as putative hot spots of CCR2. 相似文献
Despite the efficacy of imatinib therapy in chronic myelogenous leukemia, the development of drug-resistant Abl mutants, especially the most difficult overcoming T315I mutant, makes the search for new Abl T315I inhibitors a very interesting challenge in medicinal chemistry. In this work, a multistep computational framework combining the three dimensional quantitative structure-activity relationship (3D-QSAR), molecular docking, molecular dynamics (MD) simulation and binding free energy calculation, was performed to explore the structural requirements for the Abl T315I activities of benzimidazole/benzothiazole derivatives and the binding mechanism between the inhibitors and Abl T315I. The established 3D-QSAR models exhibited satisfactory internal and external predictability. Docking study elucidated the comformations of compounds and the key amino acid residues at the binding pocket, which were confirmed by MD simulation. The binding free energies correlated well with the experimental activities. The MM-GBSA energy decomposition revealed that the van der Waals interaction was the major driving force for the interaction between the ligands and Abl T315I. The hydrogen bond interactions between the inhibitors and Met318 also played an important role in stablizing the binding of compounds to Abl T315I. Finally, four new compounds with rather high Abl T315I activities were designed and presented to experimenters for reference. 相似文献
The protonation states of catalytic Asp25/25′ residues remarkably affect the binding mechanism of the HIV-1 protease–inhibitor complex. Here we report a molecular dynamics simulation study, which includes electrostatic polarisation effect, to investigate the influence of Asp25/25′ protonation states upon the binding free energy of the HIV-1 protease and a C2-symmetric inhibitor. Good agreements are obtained on inhibitor structure, hydrogen bond network, and binding free energy between our theoretical calculations and the experimental data. The calculations show that the Asp25 residue is deprotonated, and the Asp25′ residue is protonated. Our results reveal that the Asp25/25′ residues can have different protonation states when binding to different inhibitors although the protease and the inhibitors have the same symmetry. This study offers some insights into understanding the protonation state of HIV-1 protease–inhibitor complex, which could be helpful in designing new inhibitor molecules. 相似文献
Research on anticancer properties of natural compounds, as effective materials that are available while causing minimal side effects, is growing. Ellagic acid (EA) is a well-known polyphenolic compound, which has been found in both free and complex modes in several medicinal plants such as pomegranate, walnut, and berries. Although many articles have reported anticancer properties for this compound, its mechanism of action has not been fully elucidated. In this study, we used several online and offline bioinformatics tools and databases to identify the mechanism of action of EA on various types of human malignancies including bladder, blood, breast, cervical, colorectal, liver, pancreas, and prostate cancers. In this context, after identifying and extracting EA-affected human genes/proteins that have been reported in various references, we built the related gene networks and determined functional hub genes. In addition, docking was performed to recognize target proteins that react directly with EA and are in fact most affected by this compound. Our findings revealed that EA exerts its anticancer effects by influencing specific hub genes in various types of cancers. Moreover, different cellular signaling pathways are affected by this natural compound. Generally, it turned out that EA probably exerts most of its anticancer activities, through induction of apoptosis, as well as P53 and WNT signaling pathways, and also by affecting the expression of several hub genes such as CDKN1A, CDK4, CDK2, CDK6, TP53, JUN, CCNA2, MAPK14, CDK1, and CCNB1 and especially interactions with some related proteins including P53, CDK6, and MAPK14.
The formation mechanism for the icosahedral central structure of the He1 13 cluster is proposed and its total energy curve is calculated by the method of a Modified Arrangement Channel Quantum Mechanics. The energy is the function of separation R between two nuclei at the center and an apex of the icosahedral central structure. The result of the calculation has shown that the curve has a minimal energy -37.5765 (a.u.) at R = 2.70ao. The binding energy of He 13 with respect to He 12He was calculated to be 1.4046 a.u. This means that the cluster of He 13 may be formed in an icosahedral central structure with strong binding energy. 相似文献
The c-Met tyrosine kinase plays an important role in human cancers. Preclinical studies demonstrated that c-Met is over-expressed, mutated and amplified in a variety of human tumor types and design of more potent c-Met inhibitors is a priority. In this study, 14 molecular dynamics simulations of potent type II c-Met inhibitors were run to resolve the critical interactions responsible for high affinity of ligands towards c-Met considering the essential flexibility of protein–ligand interactions. Residues Phe1223 and Tyr1159, involved in pi-pi interactions were recognized as the most effective residues in the ligand binding in terms of binding free energies. Hydrogen bond interaction with Met1160 was also found necessary for effective type II ligand binding to c-Met.