Insight Into the Binding Mechanism of p53/pDIQ-MDMX/MDM2 With the Interaction Entropy Method

Structural optimization of pyrrolopyrimidine BTK inhibitors based on molecular simulation

CONTEXT

BTK is a critical regulator involved in the proliferation, differentiation, and apoptosis of B cells. BTK inhibitors can effectively alleviate various diseases such as tumors, leukemia, and asthma. During this study, a range of novel BTK inhibitors were designed using 3D-QSAR, molecular docking, and molecular dynamics (MD) simulation.

METHODS

We selected 41 pyrrolopyrimidine derivatives as BTK inhibitors to structure a 3D-QSAR model. Comparative molecular field analysis (CoMFA) and comparative molecular similarity index analysis (CoMSIA) were adopted to research the connection between the pharmacological activities and chemical structures of the compounds. The CoMFA model (q2 = 0.519, R2 = 0.971), CoMSIA model (q2 = 0.512, R2 = 0.990), and external validation demonstrated excellent predictive performance and reliability of the 3D-QSAR model. We designed eight novel molecules with higher inhibitory activities according to the three-dimensional equipotential fields and explored the interactions between the compounds and BTK by molecular docking, which showed that the novel molecules had higher binding affinities with BTK than the template molecule 18. Then, the results of molecular docking were further verified by MD simulation, which showed that amino acid residues such as Leu528, Val416, and Met477 played vital parts in the interaction, and the binding free energy analysis showed that the novel molecules had higher stability with BTK. Finally, the ADME/T properties were predicted for all of the novel compounds, and the results showed that the majority of them had favorable pharmacokinetic properties. Therefore, this study provides strong support for the development of novel BTK inhibitors.

Hesperidin Inhibits the p53-MDMXInteraction-Induced Apoptosis of Non-Small-Cell Lung Cancer and Enhances the Antitumor Effect of Carboplatin

Objective

This study aimed to observe the effect of hesperidin on the apoptosis, proliferation, and invasion of non-small-cell lung cancer, as well as to explore the possible mechanism. The inhibitory effect of hesperidin combined with carboplatin on non-small-cell lung cancer was also investigated.

Methods

A549 and NCI-H460 cells were treated with different concentrations of hesperidin (10, 50, and 100 μM). The effect of siRNA knockdown on MDMX on the antitumor effect of hesperidin was observed. CCK-8 was used to detect cell activity. The apoptosis rate was determined by TUNEL. The transwell assay detects the ability of cell migration and invasion. The expression levels of the apoptosis-related proteins p53, MDM2, MDMX, p21, PUMA, Bcl-2, and Bax were detected by qRT-PCR. Cell-proliferation and transwell assays were used to detect the effects of the combined use of hesperidin and carboplatin on lung cancer cells.

Results

Hesperidin significantly inhibited the activity and invasion of A549 and NCI-H460 cells in a dose-dependent manner. Hesperidin also induced the apoptosis of A549 and NCI-H460 cells. Hesperidin further inhibited the interaction between p53 and MDMX, increased the expression of p53, and played an anticancer role. The combination of hesperidin and carboplatin showed the most obvious antitumor effect.

Conclusion

Hesperidin can inhibit lung cancer by inhibiting the interaction between p53 and MDMX. Moreover, the combination of hesperidin and carboplatin can inhibit the migration and invasion of lung cancer cell lines through p53 upregulation, thereby increasing the antitumor effect of carboplatin.

Molecular investigation of the dual inhibition mechanism for targeted P53 regulator MDM2/MDMX inhibitors

Inhibitors that competitively bind MDM2/MDMX can block the inhibition of P53 by MDM2/MDMX and restart its tumor-suppressive effect. Molecular studies targeting MDM2/MDMX inhibitors have always been a hot topic in anticancer drug design. Although numerous inhibitors have been designed previously against MDM2/MDMX, their dual inhibition efficacy has not been demonstrated, and few studies assessed the general causes affecting the dual inhibition of MDM2/MDMX by these inhibitors. Here, molecular dynamics simulations and alanine scanning combined with the interaction entropy method were employed to precisely investigate whether 16 inhibitors could dually inhibit MDM2/MDMX and the similarities and differences in the interaction modes. Thereby addressing the key residue sites affecting dual inhibition. Residues L54/M53, I61/60, M62/61, Y67/66, and V93/92 of MDM2/MDMX, which are in corresponding positions in both protein structures, provide significant conditions for these inhibitors to bind to MDM2/MDMX tightly. In addition, most of these inhibitors prefer to bind MDM2 than MDMX, and residues H96 and I99 in MDM2 are attractive targets for inhibitors, resulting in inhibitors binding to MDM2/MDMX with different affinity. These key residues should be considered in the development of dual inhibitors. For these 16 inhibitors, most have dual inhibitory potential for MDM2/MDMX based on the binding affinity of the complexes. Still, it is questionable whether they can exert excellent dual inhibition considering the assessment of the hot-spots. At least their binding affinity for MDMX is not superior to that for MDM2 due to the difference in energy of the van der Waals interactions at the key sites. Furthermore, based on the analysis of three representative inhibitors (TUZ/HRH and HRQ with different binding preferences for MDM2/MDMX), 3-chloropyridine in TUZ leads to the differential binding affinity between the inhibitor and MDM2/MDMX. It readily forms hydrophobic interactions with the surrounding residues H96 and I99. But this phenomenon does not occur in the TUZ-MDMX system, implying the critical role of residues H96/P95 and I99/L98. And the completely different binding mechanism of HRQ binding to MDM2/MDMX explains its inability to inhibit MDM2 well. Thus, we are cautious about its dual inhibitory ability. Besides, HRH is more prone to strong van der Waals interactions with MDM2 than MDMX whereas its 2-chlorofluorobenzene is detrimental to this. We hope that these findings will provide reliable molecular insights for the screening and optimization of targeting MDM2/MDMX dual inhibitors.

HISNAPI: a bioinformatic tool for dynamic hot spot analysis in nucleic acid-protein interface with a case study

Protein-nucleic acid interactions play essential roles in many biological processes, such as transcription, replication and translation. In protein-nucleic acid interfaces, hotspot residues contribute the majority of binding affinity toward molecular recognition. Hotspot residues are commonly regarded as potential binding sites for compound molecules in drug design projects. The dynamic property is a considerable factor that affects the binding of ligands. Computational approaches have been developed to expedite the prediction of hotspot residues on protein-nucleic acid interfaces. However, existing approaches overlook hotspot dynamics, despite their essential role in protein function. Here, we report a web server named Hotspots In silico Scanning on Nucleic Acid and Protein Interface (HISNAPI) to analyze hotspot residue dynamics by integrating molecular dynamics simulation and one-step free energy perturbation. HISNAPI is capable of not only predicting the hotspot residues in protein-nucleic acid interfaces but also providing insights into their intensity and correlation of dynamic motion. Protein dynamics have been recognized as a vital factor that has an effect on the interaction specificity and affinity of the binding partners. We applied HISNAPI to the case of SARS-CoV-2 RNA-dependent RNA polymerase, a vital target of the antiviral drug for the treatment of coronavirus disease 2019. We identified the hotspot residues and characterized their dynamic behaviors, which might provide insight into the target site for antiviral drug design. The web server is freely available via a user-friendly web interface at http://chemyang.ccnu.edu.cn/ccb/server/HISNAPI/ and http://agroda.gzu.edu.cn:9999/ccb/server/HISNAPI/.

Gonçalves Vianez JL, Santana da Costa K, de Molfetta FA, Nahum Alves C. Analysis of Kojic Acid Derivatives as Competitive Inhibitors of Tyrosinase: A Molecular Modeling Approach

Tyrosinases belong to the functional copper-containing proteins family, and their structure contains two copper atoms, in the active site, which are coordinated by three histidine residues. The biosynthesis of melanin in melanocytes has two stages depending on the actions of the natural substrates L-DOPA and L-tyrosine. The dysregulation of tyrosinase is involved in skin cancer initiation. In the present study, using molecular modeling tools, we analyzed the inhibition activity of tyrosinase activity using kojic acid (KA) derivatives designed from aromatic aldehydes and malononitrile. All derivatives showed conformational affinity to the enzyme active site, and a favorable distance to chelate the copper ion, which is essential for enzyme function. Molecular dynamics simulations revealed that the derivatives formed promising complexes, presenting stable conformations with deviations between 0.2 and 0.35 Å. In addition, the investigated KA derivatives showed favorable binding free energies. The most stable KA derivatives showed the following binding free energies: -17.65 kcal mol-1 (D6), -18.07 kcal mol-1 (D2), -18.13 (D5) kcal mol-1, and -10.31 kcal mol-1 (D4). Our results suggest that these derivatives could be potent competitive inhibitors of the natural substrates of L-DOPA (-12.84 kcal mol-1) and L-tyrosine (-9.04 kcal mol-1) in melanogenesis.

Structure-based designing efficient peptides based on p53 binding site residues to disrupt p53-MDM2/X interaction

MDM2 and MDMX are known as overexpressed oncoproteins in several wild-type p53 cancer cells. The development of potent and dual antagonist peptides for p53-MDM2/X is a continuous challenge. In this study, we intended to investigate the pivotal structural points respecting the development of potent and dual inhibitors of MDM2/X. Correspondingly, MD simulation was performed on the experimentally confirmed peptides, comprising p53, pDI, pDIQ, PMI, and computationally screened mutant pDI and pDIQ. A follow-up secondary structure analysis showed the last three C-terminal residues provide the helicity reservation of peptides bound to MDM2/X. Furthermore, a delicate residue-residue examination displayed Met 11 and Ser12 in the modified peptides contribute significantly to dual inhibition of MDM2/X. Additionally, the peptides_MDM2/X complexes' ΔG_binding extracted by the umbrella sampling method were in agreement with the pattern of their experimental affinity values. It was concluded the screened pDI mutants were considered as suitable anti-MDM2/X peptides, and the data obtained could be exploited as the theoretical structure-based guide for rational peptide design. Taking account of results, the suitable C-terminal residues of p53-based peptides especially Met11, and Ser12, as well as higher umbrella sampling, generated ΔG_binding to MDM2/X would be considered as the positive structural markers of a promising anti-cancer agent.

An accurate free energy estimator: based on MM/PBSA combined with interaction entropy for protein-ligand binding affinity

The molecular mechanics/Poisson-Boltzmann surface area (MM/PBSA) method is constantly used to calculate the binding free energy of protein-ligand complexes, and has been shown to effectively balance computational cost against accuracy. The relative binding affinities obtained by the MM/PBSA approach are acceptable, while it usually overestimates the absolute binding free energy. This paper proposes four free energy estimators based on the MM/PBSA for enthalpy change combined with interaction entropy (IE) for entropy change using different weights for individual energy terms. The ΔG_{PBSA_IE} method is determined to be an optimal estimator based on its performance in terms of the correlation between experimental and theoretical values and error estimations. This approach is optimized using high-quality experimental values from a training set containing 84 protein-ligand systems, and the coefficients for the sum of electrostatic energy and polar solvation free energy, van der Waals (vdW) energy, non-polar solvation energy and entropy change are obtained by multivariate linear fitting to the corresponding experimental values. A comparison between the traditional MM/PBSA method and this method shows that the correlation coefficient is improved from 0.46 to 0.72 and the slope of the regression line increases from 0.10 to 1.00. More importantly, the mean absolute error (MAE) is significantly reduced from 22.52 to 1.59 kcal mol^-1. Furthermore, the numerical stability of this method is validated on a test set with a similar correlation coefficient, slope and MAE to those of the training set. Based on the above advantages, the ΔG_{PBSA_IE} method can be a powerful tool for a reliable and accurate estimation of binding free energy and plays a significant role in a detailed energetic investigation of protein-ligand interaction.

Computational Investigation on the p53-MDM2 Interaction Using the Potential of Mean Force Study

Murine double minute 2 (MDM2) proteins are found to be overproduced by many human tumors in order to inhibit the functioning of p53 molecules, a tumor suppressor protein. Thus, reactivating p53 functioning in cancer cells by disrupting p53-MDM2 interactions may offer a significant approach in cancer treatment. However, the structural characterization of the p53-MDM2 complex at the atomistic level and the mechanism of binding/unbinding of the p53-MDM2 complex still remain unclear. Therefore, we demonstrate here the probable binding (unbinding) pathway of transactivation domain 1 of p53 during the formation (dissociation) of the p53-MDM2 complex in terms of free energy as a function of reaction coordinate from the potential of mean force (PMF) study using two different force fields: ff99SB and ff99SB-ILDN. From the PMF plot, we noticed the PMF to have a minimum value at a p53-MDM2 separation of 12 Å, with a dissociation energy of 30 kcal mol^-1. We also analyzed the conformational dynamics and stability of p53 as a function of its distance of separation from MDM2. The secondary structure content (helix and turns) in p53 was found to vary with its distance of separation from MDM2. The p53-MDM2 complex structure with lowest potential energy was isolated from the ensemble at the reaction coordinate corresponding to the minimum PMF value and subjected to molecular dynamics simulation to identify the interface surface area, interacting residues at the interface, and the stability of the complex. The simulation results highlight the importance of hydrogen bonds and the salt bridge between Lys94 of MDM2 and Glu17 of p53 in the stability of the p53-MDM2 complex. We also carried out the binding free energy calculations and the per residue energy decomposition analyses of the interface residues of the p53-MDM2 complex. We found that the binding affinity between MDM2 and p53 is indeed high [ΔG _bind = -7.29 kcal mol^-1 from molecular mechanics/Poisson-Boltzmann surface area (MM/PBSA) and ΔG _bind = -53.29 kcal mol^-1 from molecular mechanics/generalized borne surface area]. The total binding energy obtained using the MM/PBSA method was noticed to be closer to the experimental values (-6.4 to -9.0 kcal mol^-1). The p53-MDM2 complex binding profile was observed to follow the same trend even in the duplicate simulation run and also in the simulation carried out with different force fields. We found that Lys51, Leu54, Tyr100, and Tyr104 from MDM2 and the residues Phe19, Trp23, and Leu26 from p53 provide the highest energy contributions for the p53-MDM2 interaction. Our findings highlight the prominent structural and binding characteristics of the p53-MDM2 complex that may be useful in designing potential inhibitors to disrupt the p53-MDM2 interactions.

Improving the performance of the MM/PBSA and MM/GBSA methods in recognizing the native structure of the Bcl-2 family using the interaction entropy method

Correct discrimination of native structure plays an important role in drug design. IE method significantly improves the performance of MM/PB(GB)SA method in discriminating native and decoy structures in protein–ligand/protein systems of Bcl-2 family.

Combined use of murine double minute-2 promoter methylation and serum AFP improves diagnostic efficiency in hepatitis B virus-related hepatocellular carcinoma

Objective: Hepatocellular carcinoma (HCC) accounts for approximately 85% of all cases of liver cancer. In China, chronic hepatitis B virus-related HCC (HBV-related HCC) is the most common type of HCC. However, the majority of HBV-related HCC patients are asymptomatic, and the best opportunities for treating these patients are missed. The precise diagnosis of HBV-related HCC is crucial. The main purpose of this study was to evaluate the diagnostic value of murine double minute-2 (MDM2) promoter methylation in HBV-related HCC patients. Methods: The methylation status of the MDM2 promoter was detected by methylation-specific PCR. The MDM2 expression levels were validated by quantitative real-time PCR. Enzyme-linked immunosorbent assay was used to determine the levels of interleukin-6 (IL-6) and tumor-necrosis factor-α (TNF-α) in plasma. Results: The methylation frequency of the MDM2 promoter was decreased in HBV-related HCC patients. The MDM2 mRNA levels of patients with HBV-related HCC were higher than those of patients with liver cirrhosis and chronic hepatitis B. The plasma levels of IL-6 and TNF-α were significantly higher in HBV-related HCC patients than that in liver cirrhosis and chronic hepatitis B patients. The TNF-α levels were higher in the unmethylated MDM2 promoter group than in the methylated MDM2 promoter group in HBV-related HCC patients. Moreover, the combination of MDM2 promoter methylation and alpha-fetoprotein (AFP) improved the diagnosis of HBV-related HCC. Conclusions: Our study indicates, for the first time, that MDM2 promoter hypomethylation is present in HBV-related HCC patients. The combination of MDM2 promoter methylation and AFP can greatly improve diagnostic efficiency in HBV-related HCC, which might provide a new method for HBV-related HCC diagnosis.

Revealing binding selectivity of ligands toward murine double minute 2 and murine double minute X based on molecular dynamics simulations and binding free energy calculations

It is well known that the interactions of p53 with murine double minute 2 and murine double minute X, namely MDM2 and MDMX, have been significant targets of efficient anti-cancer drug design. In this study, molecular dynamics (MD) simulations, principal component (PC) analysis and binding free energy calculations are combined to recognize binding selectivity of three ligands to MDM2 and MDMX. The binding free energies were estimated by using molecular mechanics generalized Born surface area (MM-GBSA) method and the obtained results display that the increase in the binding enthalpy of three ligands to MDM2 relative to MDMX mainly drives the binding selectivity of them toward MDM2 and MDMX. The information obtained from PC analysis shows that the associations of ligands exert important impacts on internal dynamics of MDM2 and MDMX. Meanwhile, the calculations of residue-based free energy decomposition not only identify the hot interaction spots of ligands with MDM2 and MDMX, but also show the residues (L54, M53), (Y67, Y66), (V93, V92), (H96, P95), (I99, I98) and (Y100, Y99) in (MDM2, MDMX) are responsible for most contributions to the binding selectivity of three ligands toward MDM2 and MDMX. It is believed that this work can provide useful information for design of highly selective and dual inhibitors targeting MDM2 and MDMX.Communicated by Ramaswamy H. Sarma.

Drug-resistance mechanisms of three mutations in anaplastic lymphoma kinase against two inhibitors based on MM/PBSA combined with interaction entropy

As a promising drug target in the treatment of lung cancer, anaplastic lymphoma kinase (ALK) and its mutations have been studied widely. This work explored the origin of the resistance mechanism of the ALK mutants again two inhibitors.