Insight into binding mechanisms of inhibitors MKP56, MKP73, MKP86, and MKP97 to HIV-1 protease by using molecular dynamics simulation

HIV-1 protease (PR) has been a significant target for design of potent inhibitors curing acquired immunodeficiency syndrome. Molecular dynamics simulations coupled with molecular mechanics Poisson–Boltzmann surface area method were performed to study interaction modes of four inhibitors MKP56, MKP73, MKP86, and MKP97 with PR. The results suggest that the main force controlling interactions of inhibitors with PR should be contributed by van der Waals interactions between inhibitors and PR. The cross-correlation analyses based on MD trajectories show that inhibitor binding produces significant effect on the flap dynamics of PR. Hydrogen bond analyses indicate that inhibitors can form stable hydrogen bonding interactions with the residues from the catalytic strands of PR. The contributions of separate residues to inhibitor bindings are evaluated by using residue-based free energy decomposition method and the results demonstrate that the CH–π and CH–CH interactions between the hydrophobic groups of inhibitors with residues drive the associations of inhibitors with PR. We expect that this study can provide a significant theoretical aid for design of potent inhibitors targeting PR.     

Computational study on the selective inhibition mechanism of MS402 to the first and second bromodomains of BRD4

As a member of the bromodomain and extraterminal domain (BET) family, BRD4 is considered as a potential target for cancer treatment. However, because of the highly conservation of its two homologous bromodomains (BD1/BD2), selective inhibition of each bromodomain remains a challenge. MS402 is a domain-selective inhibitor of BRD4-BD1 over BRD4-BD2 reported recently. Understanding the selectivity mechanism would be very useful for the further design of more potent BD1-selectivity inhibitors. Molecular dynamics simulation, adaptive biasing force and multiple-walker adaptive biasing force were performed to study the inhibition and domain-selective mechanism of MS402 toward BRD4-BD1 over BRD4-BD2 here. Results demonstrate BRD4-BD1 binds to MS402 with lower binding free energy than BRD4-BD2. Residues Gln85, Pro86, Asn140, and Ile146 are crucial for MS402's selectively binding to BRD4-BD1. MS402 needs to overcome more energy barrier to dissociate from BD1 than from BD2 pocket. These findings will be helpful for rational structural modification of existing inhibitors to increase their BD1-selectivity.     

Binding pocket-based design,synthesis and biological evaluation of novel selective BRD4-BD1 inhibitors

Bromodomain-containing protein 4 (BRD4), consisting of two tandem bromodomains (BD1 and BD2), is key epigenetic regulator in fibrosis and cancer, which has been reported that BD1 and BD2 have distinct roles in post-translational modification. But there are few selective inhibitors toward those two domains. Herein, this study designed and synthesized a series of novel selective BRD4-BD1 inhibitors, using computer-aided drug design (CADD) approach focused on exploring the difference of the binding pockets of BD1 and BD2, and finding the His437 a crucial way to achieve BRD4-BD1 selectivity. Our results revealed that the compound 3u is a potent selective BRD4-BD1 inhibitor with IC₅₀ values of 0.56?μM for BD1 but >100?μM for BD2. The compound exhibited a broad spectrum of anti-proliferative activity against several human cancer and fibroblastic cell lines, which might be related to its capability of reducing the expression of c-Myc and collagen I. Furthermore, it could induce apoptosis in A375 cells. To the contrary, the selective BD2 inhibitor, RVX-208, did not indicate any of these activities. Our findings highlight that the function of BRD4-BD1 might be predominant in fibrosis and cancer. And it is rational to further develop novel selective BRD4-BD1 inhibitors.     

Structure investigation,enrichment analysis and structure-based repurposing of FDA-approved drugs as inhibitors of BET-BRD4

We report herein detailed structural insights into the ligand recognition modes guiding bromodomain selectivity, enrichment analysis and docking-based database screening for the identification of the FDA-approved drugs that have potential to be the human BRD4 inhibitors. Analysis of multiple X-ray structures prevailed that the lysine-recognition sites are highly conserved, and apparently, the dynamic ZA loop guides the specific ligand-recognition. The protein–ligand interaction profiling revealed that both BRD2 and BRD4 shared hydrophobic interaction of bound ligands with PRO-98/PRO-82, PHE-99/PHE-83, LEU-108/LEU-92 and direct H-bonding with ASN-156/ASN-140 (BRD2/BRD4), while on the other hand the water-mediated H-bonding of bound ligands with PRO-82, GLN-85, PRO-86, VAL-87, ASP-88, LEU-92, TYR-97 and MET-132, and aromatic π–π stacking with TRP-81 prevailed as unique interaction in BRD4, and were not observed in BRD2. Subsequently, through ROC curve analysis, the best enrichment was found with PDB-ID 4QZS of BRD4 structures. Finally, through docking-based database screening study, we found that several drugs have better binding affinity than the control candidate lead (+)-JQ1 (Binding affinity?=?-7.9?kcal/mol), a well-known BRD4 inhibitor. Among the top-ranked drugs, azelastine, a selective histamine H1 receptor antagonist, showed the best binding affinity of –9.3?kcal/mol and showed interactions with several key residues of the acetyl lysine binding pocket. Azelastine may serve as a promising template for further medicinal chemistry. These insights may serve as basis for structure-based drug design, drug repurposing and the discovery of novel BRD4 inhibitors.

Communicated by Ramaswamy H. Sarma     

Design,synthesis and biological evaluation of novel 6-phenyl-1,3a,4,10b-tetrahydro-2H-benzo[c]thiazolo[4,5-e]azepin-2-one derivatives as potential BRD4 inhibitors

Bromodomain-containing protein 4 (BRD4) is a key epigenetic regulator in cancer, and inhibitors targeting BRD4 exhibit great anticancer activity. By replacing the methyltriazole ring of the BRD4 inhibitor I-BET-762 with an N-methylthiazolidone heterocyclic ring, fifteen novel BRD4 inhibitors were designed and synthesized. Compound 13f had a hydrophobic acetylcyclopentanyl side chain, showing the most potent BRD4 inhibitory activity in the BRD4-BD1 inhibition assay (IC₅₀ value of 110 nM), it also significantly suppressed the proliferation of MV-4-11 cells with high BRD4 level (IC₅₀ value of 0.42 μM). Furthermore, the potent apoptosis-promoting and G0/G1 cycle-arresting activity of compound 13f were indicated by flow cytometry. As the downstream-protein of BRD4, c-Myc was in significantly low expression by compound 13f treatment in a dose-dependent manner. All the findings supported that this novel compound 13f provided a perspective for developing effective BRD4 inhibitors.     

Probing Difference in Binding Modes of Inhibitors to MDMX by Molecular Dynamics Simulations and Different Free Energy Methods

The p53-MDMX interaction has attracted extensive attention of anti-cancer drug development in recent years. This current work adopted molecular dynamics (MD) simulations and cross-correlation analysis to investigate conformation changes of MDMX caused by inhibitor bindings. The obtained information indicates that the binding cleft of MDMX undergoes a large conformational change and the dynamic behavior of residues obviously change by the presence of different structural inhibitors. Two different methods of binding free energy predictions were employed to carry out a comparable insight into binding mechanisms of four inhibitors PMI, pDI, WK23 and WW8 to MDMX. The data show that the main factor controlling the inhibitor bindings to MDMX arises from van der Waals interactions. The binding free energies were further divided into contribution of each residue and the derived information gives a conclusion that the hydrophobic interactions, such as CH-CH, CH-π and π-π interactions, are responsible for the inhibitor associations with MDMX.     

Binding modes and conformational changes of FK506-binding protein 51 induced by inhibitor bindings: insight into molecular mechanisms based on multiple simulation technologies

Abstract

The FK506-binding protein 51 (FKBP51) is a cochaperone that modulates the signal transduction of steroid hormone receptors and has been involved in prostate cancer, indicating that FKBP51 is an attractive target of drug design curing the related cancers. In this work, multiple short molecular dynamics (MSMD) simulations are combined with MM-GBSA method to investigate binding modes of inhibitors 3JP, 3JR and 3JQ to FKBP51. The results show that the substitutions of diols (R)-19 and (S)-19 at the R position of 3JP strengthen binding of 3JR and 3JQ to FKBP51. Principal component (PC) analysis performed on the equilibrated MSMD trajectories suggests that three inhibitor bindings produce significant effect on dynamics behavior and conformational changes of the loops L1, L2 and the domain β-L-α-L-β in FKBP51. The calculations of residue-based free energy decomposition not only recognize the hot interaction spot of inhibitors with FKBP51, but also display that the substitutions of diols (R)-19 and (S)-19 at the R position of 3JP play significant role in stronger binding of 3JR and 3JQ to FKBP51 than 3JP. This work is expected to provide theoretical hints and molecular mechanism for design of highly efficient inhibitors toward FKBP51.     

In silico design and molecular basis for the selectivity of Olinone toward the first over the second bromodomain of BRD4

Bromodomains (BrDs), a conserved structural module in chromatin-associated proteins, are well known for recognizing ε-N-acetyl lysine residues on histones. One of the most relevant BrDs is BRD4, a tandem BrD containing protein (BrD1 and BrD2) that plays a critical role in numerous diseases including cancer. Growing evidence shows that the two BrDs of BRD4 have different biological functions; hence selective ligands that can be used to study their functions are of great interest. Here, as a follow-up of our previous work, we first provide a detailed characterization study of the in silico rational design of Olinone as part of a series of five tetrahydropyrido indole-based compounds as BRD4 BrD1 inhibitors. Additionally, we investigated the molecular basis for Olinone's selective recognition by BrD1 over BrD2. Molecular dynamics simulations, free energy calculations, and conformational analyses of the apo-BRD4-BrD1|2 and BRD4-BrD1|2/Olinone complexes showed that Olinone's selectivity is facilitated by five key residues: Leu92 in BrD1|385 in BrD2 of ZA loop, Asn140|433, Asp144|His437 and Asp145|Glu438 of BC loop, and Ile146|Val49 of helix C. Furthermore, the difference in hydrogen bonds number and in mobility of the ZA and BC loops of the acetyl-lysine binding site between BRD4 BrD1/Olinone and BrD2/Olinone complexes also contribute to the difference in Olinone's binding affinity and selectivity toward BrD1 over BrD2. Altogether, our computer-aided molecular design techniques can effectively guide the development of small-molecule BRD4 BrD1 inhibitors, explain their selectivity origin, and further open doors to the design of new therapeutically improved derivatives.     

Computational revelation of binding mechanisms of inhibitors to endocellular protein tyrosine phosphatase 1B using molecular dynamics simulations

Endocellular protein tyrosine phosphatase 1B (PTP1B) is one of the most promising target for designing and developing drugs to cure type-II diabetes and obesity. Molecular dynamics (MD) simulations combined with molecular mechanics generalized Born surface area (MM-GBSA) and solvated interaction energy methods were applied to study binding differences of three inhibitors (ID: 901, 941, and 968) to PTP1B, the calculated results show that the inhibitor 901 has the strongest binding ability to PTP1B among the current inhibitors. Principal component (PC) analysis was also carried out to investigate the conformational change of PTP1B, and the results indicate that the associations of inhibitors with PTP1B generate a significant effect on the motion of the WPD-loop. Free energy decomposition method was applied to study the contributions of individual residues to inhibitor bindings, it is found that three inhibitors can generate hydrogen bonding interactions and hydrophobic interactions with different residues of PTP1B, which provide important forces for associations of inhibitors with PTP1B. This research is expected to give a meaningfully theoretical guidance to design and develop of effective drugs curing type-II diabetes and obesity.     

Molecular mechanism study of several inhibitors binding to BRD9 bromodomain based on molecular simulations

Bromodomain-containing protein 9 (BRD9) has been employed as a potential target for anticancer drugs in recent years. In this work, molecular docking, molecular dynamics (MD) simulations, binding free energy calculations, and per residue energy decomposition approaches were performed to elucidate the different binding modes between four pyridinone-like scaffold inhibitors and BRD9 bromodomain. Analysis results indicate that non-polar contribution mainly deriving from van der Waals energy is a critical impact on binding affinity of inhibitors against BRD9. Some key residues Phe44, Phe47, Val49, and Ile53 (at ZA loop) enhance the binding energy of inhibitors in BRD9 by means of providing hydrophobic interactions. Moreover, it is observed that BRD9 is anchored by the formation of a stable hydrogen bond between the carbonyl of the inhibitors and the residue Asn100 (at BC loop), and a strong π–π stacking interaction formed between the residue Tyr106 (at BC loop) and the inhibitors. The existence of dimethoxyphenyl structure and the aromatic ring merged to pyridinone scaffold are useful to enhance the BRD9 binding affinity. These findings should guide the rational design of more prospective inhibitors targeting BRD9.

Communicated by Ramaswamy H. Sarma     

An insight into the inhibitory selectivity of 4-(Pyrazol- 4-yl)-pyrimidines to CDK4 over CDK2

Revealing selectivity mechanism of cyclin-dependent kinases (CDK) and their inhibitors is an important issue to develop potential anticancer drugs. The substituted 4-(Pyrazol-4-yl)-pyrimidines are potent inhibitors of CDK4 but not of the highly homologous CDK2. In order to reveal the inhibitory selectivity of these inhibitors to CDK4 over CDK2, we select one of substituted 4-(Pyrazol-4-yl)-pyrimidines as a representative (marked as A1 hereunder) and perform molecular docking, molecular dynamics simulations and binding free energy analysis for CDK4/A1 and CDK2/A1, respectively. The electrostatic and van der Waals (vdW) interactions of the A1 inhibitor with CDK4/CDK2 are discussed. The computed binding free energies based on the MM-PBSA method are consistent with experimental bioactivity ranking of A1 inhibitor to CDK4/CDK2. On the other hand, the conformational characteristics of CDK2 and CDK4 induced by A1 inhibitor are analysed and revealed. Results demonstrate that the vdW interactions considerably contribute to binding of CDK4/CDK2 with A1 inhibitor and are similar in size. The hydrogen bonding between A1 inhibitor and CDK4/CDK2 is considerably favourable to the binding, in which the hydrogen bond between the NH group of the pyrazole group of A1 and the residue Asp158 of CDK4 plays a crucial role in inhibitory selectivity of A1 inhibitor to CDK4 over CDK2. The electrostatic interaction energy differences between the corresponding residues of CDK4/A1 and CDK2/A1 confirm the above inference. The conformational changes of CDK2 and CDK4 induced by A1 inhibitor influence the selectivity of A1 inhibitor to CDK4/CDK2.     

Small-Molecule Dual PLK1 and BRD4 Inhibitors are Active Against Preclinical Models of Pediatric Solid Tumors

Simultaneous inhibition of multiple molecular targets is an established strategy to improve the continuance of clinical response to therapy. Here, we screened 49 molecules with dual nanomolar inhibitory activity against BRD4 and PLK1, best classified as dual kinase-bromodomain inhibitors, in pediatric tumor cell lines for their antitumor activity. We identified two candidate dual kinase-bromodomain inhibitors with strong and tumor-specific activity against neuroblastoma, medulloblastoma, and rhabdomyosarcoma tumor cells. Dual PLK1 and BRD4 inhibitor treatment suppressed proliferation and induced apoptosis in pediatric tumor cell lines at low nanomolar concentrations. This was associated with reduced MYCN-driven gene expression as assessed by RNA sequencing. Treatment of patient-derived xenografts with dual inhibitor UMB103 led to significant tumor regression. We demonstrate that concurrent inhibition of two central regulators of MYC protein family of protooncogenes, BRD4, and PLK1, with single small molecules has strong and specific antitumor effects in preclinical pediatric cancer models.     

Unravelling novel congeners from acetyllysine mimicking ligand targeting a lysine acetyltransferase PCAF bromodomain

p300/CBP Associated Factor (PCAF) bromodomain (BRD), a lysine acetyltransferases, has emerged as a promising drug target as its dysfunction is linked to onset and progression of several diseases like cancer, diabetes, AIDS, etc. In this study, a three featured E-Pharmacophore (ARR) was generated based on acetyllysine mimicking inhibitor of PCAF BRD which is available as co-crystal structure (PDB ID: 5FDZ). It was used for filtering small molecule databases followed by molecular docking and consequently validated using enrichment calculation. The resulted hits were found to be congeners which show the predictive power of E-Pharmacophore hypothesis. Further, Induced Fit Docking method, Binding energy calculation, ADME prediction, Single Point Energy calculation and Molecular Dynamics simulation were performed to find better hits against PCAF BRD. Based on the results, it was concluded that Asn803, Tyr809 and Tyr802 along with a water molecule (HOH1001) plays crucial role in binding with inhibitor. It is also proposed that four hits from Life Chemicals database namely, F2276-0099, F2276-0008, F2276-0104 and F2276-0106 could act as potent drug molecules for PCAF BRD. Thus, the present study is strongly believed to have bright impact on rational drug design of potent and novel congeners of PCAF BRD inhibitors.     

Design,synthesis and biological evaluation of novel indole derivatives as potential HDAC/BRD4 dual inhibitors and anti-leukemia agents

HDAC inhibitors and BRD4 inhibitors were considered to be potent anti-cancer agents. Recent studies have demonstrated that HDAC and BRD4 participate in the regulation of some signal paths like PI3K-AKT. In this work, a series of indole derivatives that combine the inhibitory activities of BRD4 and HDAC into one molecule were designed and synthesized through the structure-based design method. Most compounds showed potent HDAC inhibitory activity and moderate BRD4 inhibitory activity. In vitro anti-proliferation activities of the synthesized compounds were also evaluated. Among them, 19f was the most potent inhibitor against HDAC3 with IC₅₀ value of 5 nM and BRD4 inhibition rate of 88% at 10 μM. It was confirmed that 19f could up-regulate the expression of Ac-H3 and reduce the expression of c-Myc by western blot analysis. These results indicated that 19f was a potent dual HDAC/BRD4 inhibitor and deserved further investigation.     

Discovery of a new class of PROTAC BRD4 degraders based on a dihydroquinazolinone derivative and lenalidomide/pomalidomide

BRD4 has emerged as an attractive target for anticancer therapy. However, BRD4 inhibitors treatment leads to BRD4 protein accumulation, together with the reversible nature of inhibitors binding to BRD4, which may limit the efficacy of BRD4 inhibitors. To address these problems, a protein degradation strategy based on the proteolysis targeting chimera (PROTAC) technology has been developed to target BRD4 recently. Herein, we present our design, synthesis and biological evaluation of a new class of PROTAC BRD4 degraders, which were based on a potent dihydroquinazolinone-based BRD4 inhibitor compound 6 and lenalidomide/pomalidomide as ligand for E3 ligase cereblon. Gratifyingly, several compounds showed excellent inhibitory activity against BRD4, and high anti-proliferative potency against human monocyte lymphoma cell line THP-1. Especially, compound 21 (BRD4 BD1, IC₅₀ = 41.8 nM) achieved a submicromolar IC₅₀ value of 0.81 μM in inhibiting the growth of THP-1 cell line, and was 4 times more potent than compound 6. Moreover, the mechanism study established that 21 could effectively induce the degradation of BRD4 protein and suppression of c-Myc. All of these results suggested that 21 was an efficacious BRD4 degrader for further investigation.     

Comparative density functional theory based study of the reactivity of Cu,Ag, and Au nanoparticles and of (111) surfaces toward CO oxidation and NO2 reduction

Development of multi-target drugs is becoming increasingly attractive in the repertoire of protein kinase inhibitors discovery. In this study, we carried out molecular docking, molecular dynamics simulations, molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) binding free energy calculations, principal component analysis (PCA), and dynamical cross-correlation matrices (DCCM) to dissect the molecular mechanism for the valmerin-19 acting as a dual inhibitor for glycogen synthase kinase 3β (GSK3β) and cyclin-dependent kinase 5 (CDK5). Detailed MM-PBSA calculations revealed that the binding free energies of the valmerin-19 to GSK3β/CDK5 were calculated to be ?12.60?±?2.28 kcal mol^-1 and ?11.85?±?2.54 kcal mol^-1, respectively, indicating that valmerin-19 has the potential to act as a dual inhibitor of GSK3β/CDK5. The analyses of PCA and DCCM results unraveled that binding of the valmerin-19 reduced the conformational dynamics of GSK3β/CDK5 and the valmerin-19 bound to GSK3β/CDK5 might occur mostly through a conformational selection mechanism. This study may be helpful for the future design of novel and potent dual GSK3β/CDK5 inhibitors.     

Molecular dynamics simulation,binding free energy calculation and unbinding pathway analysis on selectivity difference between FKBP51 and FKBP52: Insight into the molecular mechanism of isoform selectivity

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.     

Insights into the crystal structure of BRD2-BD2 – phenanthridinone complex and theoretical studies on phenanthridinone analogs

Bromodomain and extra-terminal family proteins recognize the acetylated histone code on chromatin and participate in downstream processes like DNA replication, modification, and repair. As part of epigenetic approaches, BRD2 and BRD4 were identified as putative targets, for the management of chronic diseases. We have recently reported the discovery of a new scaffold of the phenanthridinone-based inhibitor (L10) of the second bromodomain of BRD2 (BRD2-BD2). Here, we present the crystal structure of the BRD2-BD2, refined to 1.4 Å resolution, in complex with β-mercaptoethanol (a component of the protein buffer). The β-mercaptoethanol covalently links to C425 of BD2 in the acetyl-lysine binding pocket, to form a modified cysteine mercaptoethanol (CME). The CME modification significantly hinders the entry of ligands into the BD2 binding pocket, suggesting that β-mercaptoethanol should be removed during protein production process. Next, to confirm whether phenanthridionone scaffold is a new inhibitor family of BRD2-BD2, we have determined the crystal structure of BD2 in complex with 6(5H)-Phenanthridinone (a core moiety of L10), refined to 1.28 Å resolution. It confirmed that the phenanthridinone molecule, unambiguously, binds to BD2. Moreover, we performed molecular docking and molecular dynamic studies on selected phenanthridinone analogs. The predicted L10 analogs are stable with essential hydrophobic and hydrophilic interactions with BD2 during molecular dynamic simulations. We propose that the predicted phenanthridinone analogs may be potential molecules for inhibiting the BD2 function of acetylated histone recognition.