Insight into mechanism of small molecule inhibitors of the MDM2–p53 interaction: Molecular dynamics simulation and free energy analysis

Insights into the Interaction Mechanisms of Peptide and Non-Peptide Inhibitors with MDM2 Using Gaussian-Accelerated Molecular Dynamics Simulations and Deep Learning

Inhibiting MDM2-p53 interaction is considered an efficient mode of cancer treatment. In our current study, Gaussian-accelerated molecular dynamics (GaMD), deep learning (DL), and binding free energy calculations were combined together to probe the binding mechanism of non-peptide inhibitors K23 and 0Y7 and peptide ones PDI6W and PDI to MDM2. The GaMD trajectory-based DL approach successfully identified significant functional domains, predominantly located at the helixes α2 and α2', as well as the β-strands and loops between α2 and α2'. The post-processing analysis of the GaMD simulations indicated that inhibitor binding highly influences the structural flexibility and collective motions of MDM2. Calculations of molecular mechanics-generalized Born surface area (MM-GBSA) and solvated interaction energy (SIE) not only suggest that the ranking of the calculated binding free energies is in agreement with that of the experimental results, but also verify that van der Walls interactions are the primary forces responsible for inhibitor-MDM2 binding. Our findings also indicate that peptide inhibitors yield more interaction contacts with MDM2 compared to non-peptide inhibitors. Principal component analysis (PCA) and free energy landscape (FEL) analysis indicated that the piperidinone inhibitor 0Y7 shows the most pronounced impact on the free energy profiles of MDM2, with the piperidinone inhibitor demonstrating higher fluctuation amplitudes along primary eigenvectors. The hot spots of MDM2 revealed by residue-based free energy estimation provide target sites for drug design toward MDM2. This study is expected to provide useful theoretical aid for the development of selective inhibitors of MDM2 family members.

Promising disruptors of p53-MDM2 dimerization from some medicinal plant phytochemicals: a molecular modeling study

Cancer is a major global health issue that has a high mortality rate. p53, which functions as a tumor suppressor, is critical in preventing tumor development by regulating the cell cycle and inducing apoptosis in damaged cells. However, the tumor suppressor function of p53 is effectively inhibited by its direct interaction with the hydrophobic cleft of MDM2 protein via multiple mechanisms As a result, restoring p53 activity by blocking the p53-MDM2 protein-protein interaction has been proposed as a compelling therapeutic strategy for cancer treatment. The use of molecular docking and phytochemical screening procedures are appraised to inhibit MDM2's hydrophobic cleft and disrupt the p53-MDM2 interaction. For this purpose, a library of 51 bioactive compounds from 10 medicinal plants was compiled and subjected to structure-based virtual screening. Out of these, only 3 compounds (Atalantoflavone, Cudraxanthone 1, and Ursolic acid) emerged as promising inhibitors of MDM2-p53 based on their binding affinities (-9.1 kcal/mol, -8.8 kcal/mol, and -8.8 kcal/mol respectively) when compared to the standard (-8.8 kcal/mol). Moreover, these compounds showed better pharmacokinetic and drug-like profiling than the standard inhibitor (Chromonotriazolopyrimidine 1). Finally, the 100 ns MD simulation analysis confirmed no significant perturbation in the conformational dynamics of the simulated binary complexes when compared to the standard. In particular, Ursolic acid was found to satisfy the molecular enumeration the most compared to the other inhibitors. Our overall molecular modeling finding shows why these compounds may emerge as potent arsenals for cancer therapeutics. Nonetheless, extensive experimental and clinical research is needed to augment their use in clinics.Communicated by Ramaswamy H. Sarma.

A systematic drug repurposing approach to identify promising inhibitors from FDA-approved drugs against Nsp4 protein of SARS-CoV-2

COVID-19 is caused by SARS-CoV-2 and responsible for the ongoing global pandemic in the world. After more than a year, we are still in lurch to combat and control the situation. Therefore, new therapeutic options to control the ongoing COVID-19 are urgently in need. In our study, we found that nonstructural protein 4 (Nsp4) of SARS-CoV-2 could be a potential target for drug repurposing. Due to availability of only the crystal structure of C-terminal domain of Nsp4 (Ct-Nsp4) and its crucial participation in viral RNA synthesis, we have chosen Ct-Nsp4 as a target for screening the 1600 FDA-approved drugs using molecular docking. Top 102 drugs were found to have the binding energy equal or less than -7.0 kcal/mol. Eribulin and Suvorexant were identified as the two most promising drug molecules based on the docking score. The dynamics of Ct-Nsp4-drug binding was monitored using 100 ns molecular dynamics simulations. From binding free energy calculation over the simulation, both the drugs were found to have considerable binding energy. The present study has identified Eribulin and Suvorexant as promising drug candidates. This finding will be helpful to accelerate the drug discovery process against COVID-19 disease.Communicated by Ramaswamy H. Sarma.

Structural Analysis of the Spike Protein of SARS-CoV-2 Variants and Other Betacoronaviruses Using Molecular Dynamics

A structural analysis over various spike proteins from three highly pathogenic Betacoronavirus was done to understand their structural differences. The proteins were modeled using crystal structures from SARS-CoV, MERS-CoV, and other Betacoronavirus that infect bats and pangolins. The group was split in two sets; the first set corresponds to the non-mutated spike proteins, while the second set corresponds to mutated spike variants alpha, beta, gamma, delta, omicron and mu; five of them classified as variants of concern and the last one as variant of interest. A conformational space exploration was carried out for every protein by using molecular dynamic simulations. Root mean square fluctuations, principal component and cross-correlation analysis were carried out over the dynamics to analyze the flexibility and rigidity of every protein in comparison to the wild type Spike protein from the SARS-CoV-2. The obtained results indicate that the proteins, which are not spread among humans, have smooth movements compared to those of SARS-CoV-2 and its variants. In addition, a relationship between the speed of the virulence and the movement of the protein can explain the behavior of delta and omicron variants.

Integrated virtual screening and molecular dynamics simulation revealed promising drug candidates of p53-MDM2 interaction

In the vast majority of malignancies, the p53 tumor suppressor pathway is compromised. In some cancer cells, high levels of MDM2 polyubiquitinate p53 and mark it for destruction, thereby leading to a corresponding downregulation of the protein. MDM2 interacts with p53 via its hydrophobic pocket, and chemical entities that block the dimerization of the protein-protein complex can restore p53 activity. Thus far, only a few chemical compounds have been reported as potent arsenals against p53-MDM2. The Protein Data Bank has crystallogaphic structures of MDM2 in complex with certain compounds. Herein, we have exploited one of the complexes in the identification of new p53-MDM2 antagonists using a hierarchical virtual screening technique. The initial stage was to compile a targeted library of structurally appropriate compounds related to a known effective inhibitor, Nutlin 2, from the PubChem database. The identified 57 compounds were subjected to virtual screening using molecular docking to discover inhibitors with high binding affinity for MDM2. Consequently, five compounds with higher binding affinity than the standard emerged as the most promising therapeutic candidates. When compared to Nutlin 2, four of the drug candidates (CID_140017825, CID_69844501, CID_22721108, and CID_22720965) demonstrated satisfactory pharmacokinetic and pharmacodynamic profiles. Finally, MD simulation of the dynamic behavior of lead-protein complexes reveals the stability of the complexes after a 100,000 ps simulation period. In particular, when compared to the other three leads, overall computational modeling found CID_140017825 to be the best pharmacological candidate. Following thorough experimental trials, it may emerge as a promising chemical entity for cancer therapy.

Targeting p53-MDM2 interactions to identify small molecule inhibitors for cancer therapy: beyond "Failure to rescue"

At present, disrupting p53-MDM2 interactions through small molecule ligands is a promising approach to safe treatment and management of human cancer. Tumor cells unlike the normal cells, are rapidly evolving affecting the efficacy of many approved anti-cancer agents due to drug resistance. Therefore, identifying a potential anticancer compound is crucial. Pharmacophore based virtual screening, followed by molecular docking, ADMET evaluation, and molecular dynamics studies against MDM2 protein was investigated to identify potential ligands that may act as inhibitors. The model (AHRR_1) with survival score (4.176) was selected among the top ranked generated Pharmacophore hypothesis. Validation of the model hypothesis by an external dataset of actives and inactive compounds produced significant validation attributes including; AUC = 0.85, BEDROC = 0.56 at α = 20.0, RIE = 8.18, AUAC = 0.88, and EF of 6.2 at the top 2% of the dataset. The model was use for screening the ZINC database, and the top 1375 hits satisfying the model hypothesis were subjected to molecular docking studies to understand the molecular and structural basis of selectivity of compounds for MDM2 protein. A sub-set of 25 compounds with binding energy lower than the reference inhibitors were evaluated for pharmacokinetic properties. Four compounds (ZINC02639178, ZINC06752762, ZINC38933175, and ZINC77969611) showed the most desired pharmacokinetic profile. Lastly, investigation of the dynamic behaviour of leads-protein complexes through MD simulation showed similar RMSD, RMSF, and H-bond occupancy profile compared to a reference inhibitor, suggesting stability throughout the simulation time. However, ZINC02639178 was found to satisfy the molecular enumeration the most compared to the other three leads. It may emerge as potential treatment option after extensive experimental studies. Communicated by Ramaswamy H. Sarma.

Effect of Stapling on the Thermodynamics of mdm2-p53 Binding

Protein-protein interaction (PPI) is one of the key regulatory features driving biomolecular processes and hence is targeted for designing therapeutics against diseases. Small peptides are a new and emerging class of therapeutics owing to their high specificity and low toxicity. For achieving efficient targeting of the PPI, amino acid side chains are often stapled together, resulting in the rigidification of these peptides. Exploring the scope of these peptides demands a comprehensive understanding of their working principle. In this work, two stapled p53 peptides have been considered to delineate their binding mechanism with mdm2 using computational approaches. The addition of stapling agent protects the secondary structure of the peptides even in the case of thermal and chemical denaturation. Although the introduction of a stapling agent increases the hydrophobicity of the peptide, the enthalpic stabilization decreases. This is overcome by the lowering of the entropic penalty, and the overall binding affinity improves. The mechanistic insights into the benefit of peptide stapling can be adopted for further improvement of peptide therapeutics.

Design and activity detection of BH3-mimetic peptide based on Bcl-XL binding potency

Communicated by Ramaswamy H. Sarma.

Molecular Dynamics Studies on the Enzalutamide Resistance Mechanisms Induced by Androgen Receptor Mutations

The second-generation antiandrogen enzalutamide, targeting androgen receptor (AR), was approved to treat castration resistant prostate cancer (CRPC) in 2012. Its resistance was observed when it was in the clinical research stage. AR mutation is the main factor of enzalutamide resistance. AR F876L and F876L_T877A mutations were reported to switch enzalutamide from AR antagonist to agonist, but W741C cannot. There are various mutations in the ligand binding domain of AR LBD, such as L701H, W741L, H874Y, T877A, and M895T, if these mutations can lead to drug resistance problem or not is not known. In this work, molecular dynamics (MD) simulations and molecular mechanics Generalized Born (GB) surface area (MM-GBSA) calculations were employed to explore the interaction mechanisms between enzalutamide and wild-type (WT)/mutant ARs. The simulation results indicate that helix 12 (H12), which lies on the top of the AR LBD like a cover, plays a vital role for the function of enzalutamide. When C-ring of enzalutamide locates near to H12, the distance between enzalutamide and H12 is reduced, which prevents H12 from closing and distort the coactivator binding site, resulting in the inactivation of transcription. In this case, enzalutamide acts as an AR antagonist. However, when the C-ring of enzalutamide is near to helix H11 or the Loop 11-12, H12 tends to close to form a coactivator binding site to facilitate transcription, enzalutamide acts as an AR agonist. Moreover, per-residue free energy decomposition analysis indicates that M895 and I899 are key residues in the antagonist mechanism of enzalutamide. J. Cell. Biochem. 118: 2792-2801, 2017. © 2017 Wiley Periodicals, Inc.

Association between positive murine double minute 2 expression and clinicopathological characteristics of esophageal squamous cell carcinoma: a meta-analysis

The correlations of murine double minute 2 (MDM2) T309G and esophageal cancer were elucidated because the association between MDM2 expression states and clinicopathological parameters of esophageal squamous cell carcinoma (ESCC) is controversial. We conducted a meta-analysis on studies screened from PubMed, Web of Science, Embase, the Cochrane Library, and the Chinese Biomedical Literature Databases that were published before October 2014. All studies describing the association between MDM2 and ESCC were traced. Meta-analysis was performed using the STATA software (Stata Corp., College Station, TX, USA). A total of 9 studies with 707 cases and 324 controls were included. MDM2 expression was higher in ESCC than in normal esophageal epithelium (odds ratio [OR] 10.38, 95% confidence interval [CI] 6.42-16.78, P < 0.001). High MDM2 expression was associated with early primary tumor stage (T1/T2 vs. T3/T4, OR 0.59, 95% CI 0.38-0.92, P = 0.018) and increased risk of regional lymph node metastasis (N0 vs. N1, OR 1.66, 95% CI 1.03-2.67, P = 0.039). However, no relationship was observed between MDM2 expression and the risk of distant metastasis (OR = 2.09, 95% CI 1.00-4.36, P = 0.050), and MDM2 was not significantly correlated with TP53 expression (OR 1.22, 95% CI 0.53-2.77, P = 0.643). Our analysis suggests that MDM2 acts as a potent marker of early primary tumor stage but higher risk of regional lymph node metastasis in ESCC. However, because of the limited number of studies included, the result should be further clarified by well-designed prospective studies.

Energetic Landscape of MDM2-p53 Interactions by Computational Mutagenesis of the MDM2-p53 Interaction

The ubiquitin ligase MDM2, a principle regulator of the tumor suppressor p53, plays an integral role in regulating cellular levels of p53 and thus a prominent role in current cancer research. Computational analysis used MUMBO to rotamerize the MDM2-p53 crystal structure 1YCR to obtain an exhaustive search of point mutations, resulting in the calculation of the ΔΔG comprehensive energy landscape for the p53-bound regulator. The results herein have revealed a set of residues R65-E69 on MDM2 proximal to the p53 hydrophobic binding pocket that exhibited an energetic profile deviating significantly from similar residues elsewhere in the protein. In light of the continued search for novel competitive inhibitors for MDM2, we discuss possible implications of our findings on the drug discovery field.

Probing Origin of Binding Difference of inhibitors to MDM2 and MDMX by Polarizable Molecular Dynamics Simulation and QM/MM-GBSA Calculation

Binding abilities of current inhibitors to MDMX are weaker than to MDM2. Polarizable molecular dynamics simulations (MD) followed by Quantum mechanics/molecular mechanics generalized Born surface area (QM//MM-GBSA) calculations were performed to investigate the binding difference of inhibitors to MDM2 and MDMX. The predicted binding free energies not only agree well with the experimental results, but also show that the decrease in van der Walls interactions of inhibitors with MDMX relative to MDM2 is a main factor of weaker bindings of inhibitors to MDMX. The analyses of dihedral angles based on MD trajectories suggest that the closed conformation formed by the residues M53 and Y99 in MDMX leads to a potential steric clash with inhibitors and prevents inhibitors from arriving in the deep of MDMX binding cleft, which reduces the van der Waals contacts of inhibitors with M53, V92, P95 and L98. The calculated results using the residue-based free energy decomposition method further prove that the interaction strength of inhibitors with M53, V92, P95 and L98 from MDMX are obviously reduced compared to MDM2. We expect that this study can provide significant theoretical guidance for designs of potent dual inhibitors to block the p53-MDM2/MDMX interactions.

Effect of introducing aib in a designed helical inhibitor of hdm2-p53 interaction: A molecular dynamics study

Although p53 is an intrinsically disordered protein, upon binding to Hdm2, a short stretch (residues 19-25) comprising the binding epitope assumes a helical backbone. Because the allowed conformational space of α-aminoisobutyric acid (Aib) is restricted to only the helical basin, Aib-containing helical mimics of p53 (binding epitope) are expected to inhibit interaction between p53 and Hdm2 with a much stronger affinity than the wild type p53 peptide (binding epitope), due to the entropic advantage associated with Aib. However, the IC50 values for the disruption of p53-Hdm2 interaction by Aib-p53 peptides and wild type p53 peptide were found to be comparable (J. Peptide Res. 2002, 60:88-94). To understand why incorporation of Aib didn't substantially increase Hdm2 affinity of Aib-p53 peptides, a series of molecular dynamics simulations were performed. It was found that despite stabilizing a helical backbone in the unbound state, the Aib residues in Aib-p53 peptide arrested two functionally important side-chains (F19 and W23) in non-productive conformations, resulting in relative side-chain orientations of the binding triad F19-W23-L26 incompatible with the bound conformation. Therefore, although a Aib-induced pre-formed helical peptide backbone in the unbound state is expected to favor binding, the locked side-chain orientations of the binding triad in non-productive modes would disfavor binding. This study shows that when using Aib to design functionally important helical peptides, care must be taken to consider potential interactions between side-chains of neighboring residues and Aib in the unbound state.

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.

Elucidation of Ligand-Dependent Modulation of Disorder-Order Transitions in the Oncoprotein MDM2

Numerous biomolecular interactions involve unstructured protein regions, but how to exploit such interactions to enhance the affinity of a lead molecule in the context of rational drug design remains uncertain. Here clarification was sought for cases where interactions of different ligands with the same disordered protein region yield qualitatively different results. Specifically, conformational ensembles for the disordered lid region of the N-terminal domain of the oncoprotein MDM2 in the presence of different ligands were computed by means of a novel combination of accelerated molecular dynamics, umbrella sampling, and variational free energy profile methodologies. The resulting conformational ensembles for MDM2, free and bound to p53 TAD (17-29) peptide identify lid states compatible with previous NMR measurements. Remarkably, the MDM2 lid region is shown to adopt distinct conformational states in the presence of different small-molecule ligands. Detailed analyses of small-molecule bound ensembles reveal that the ca. 25-fold affinity improvement of the piperidinone family of inhibitors for MDM2 constructs that include the full lid correlates with interactions between ligand hydrophobic groups and the C-terminal lid region that is already partially ordered in apo MDM2. By contrast, Nutlin or benzodiazepinedione inhibitors, that bind with similar affinity to full lid and lid-truncated MDM2 constructs, interact additionally through their solubilizing groups with N-terminal lid residues that are more disordered in apo MDM2.

Life as we know it depends on interactions between proteins. There is substantial evidence that many interactions between proteins involve very flexible protein regions. These disordered regions may undergo disorder/order transitions upon forming an interaction with another protein. Many successful approaches to medicinal chemistry are based on mimicking the interactions of biological molecules with man-made small molecules. However how drug-like small-molecules may modulate protein disorder is currently poorly understood, largely because it is difficult to measure in details this type of interaction with experimental methods. Here we have used computer simulations to resolve with great details the process by which different small-molecules modulate the flexibility of a disordered region of the protein MDM2. This protein is overexpressed in many cancers and small-molecules that recognize MDM2 have been developed over the last decade as possible novel anti-cancer agents. We show that the flexible MDM2 “lid” region adopts different conformational states in the presence of different small-molecules. Our results suggest why some classes of small-molecules form favorable interactions with the lid region, whereas others do not. These findings may prove crucial to develop new and more effective MDM2 inhibitors, and more generally to help drug designers target disordered proteins regions with small-molecules.

CH-π hydrogen bonds in biological macromolecules

This is a sequel to the previous Perspective "The CH-π hydrogen bond in chemistry. Conformation, supramolecules, optical resolution and interactions involving carbohydrates", which featured in a PCCP themed issue on "Weak Hydrogen Bonds - Strong Effects?": Phys. Chem. Chem. Phys., 2011, 13, 13873-13900. Evidence that weak hydrogen bonds play an enormously important role in chemistry and biochemistry has now accumulated to an extent that the rigid classical concept of hydrogen bonds formulated by Pauling needs to be seriously revised and extended. The concept of a more generalized hydrogen bond definition is indispensable for understanding the folding mechanisms of proteins. The CH-π hydrogen bond, a weak molecular force occurring between a soft acid CH and a soft base π-electron system, among all is one of the most important and plays a functional role in defining the conformation and stability of 3D structures as well as in many molecular recognition events. This concept is also valuable in structure-based drug design efforts. Despite their frequent occurrence in organic molecules and bio-molecules, the importance of CH-π hydrogen bonds is still largely unknown to many chemists and biochemists. Here we present a review that deals with the evidence, nature, characteristics and consequences of the CH-π hydrogen bond in biological macromolecules (proteins, nucleic acids, lipids and polysaccharides). It is hoped that the present Perspective will show the importance of CH-π hydrogen bonds and stimulate interest in the interactions of biological macromolecules, one of the most fascinating fields in bioorganic chemistry. Implication of this concept is enormous and valuable in the scientific community.

Binding modes of three inhibitors 8CA, F8A and I4A to A-FABP studied based on molecular dynamics simulation

Adipocyte fatty-acid binding protein (A-FABP) is an important target of drug designs treating some diseases related to lipid-mediated biology. Molecular dynamics (MD) simulations coupled with solvated interaction energy method (SIE) were carried out to study the binding modes of three inhibitors 8CA, F8A and I4A to A-FABP. The rank of our predicted binding affinities is in accordance with experimental data. The results show that the substitution in the position 5 of N-benzyl and the seven-membered ring of N-benzyl-indole carboxylic acids strengthen the I4A binding, while the substitution in the position 2 of N-benzyl weakens the F8A binding. Computational alanine scanning and dynamics analyses were performed and the results suggest that the polar interactions of the positively charged residue R126 with the three inhibitors provide a significant contribution to inhibitor bindings. This polar interaction induces the disappearance of the correlated motion of the C terminus of A-FABP relative to the N terminus and favors the stability of the binding complex. This study is helpful for the rational design of potent inhibitors within the fields of metabolic disease, inflammation and atherosclerosis.

Investigating the impact of Asp181 point mutations on interactions between PTP1B and phosphotyrosine substrate

Protein tyrosine phosphatase 1B (PTP1B) is a key negative regulator of insulin and leptin signaling, which suggests that it is an attractive therapeutic target in type II diabetes and obesity. The aim of this research is to explore residues which interact with phosphotyrosine substrate can be affected by D181 point mutations and lead to increased substrate binding. To achieve this goal, molecular dynamics simulations were performed on wild type (WT) and two mutated PTP1B/substrate complexes. The cross-correlation and principal component analyses show that point mutations can affect the motions of some residues in the active site of PTP1B. Moreover, the hydrogen bond and energy decomposition analyses indicate that apart from residue 181, point mutations have influence on the interactions of substrate with several residues in the active site of PTP1B.

A computational analysis of binding modes and conformation changes of MDM2 induced by p53 and inhibitor bindings

Molecular dynamics (MD) simulations followed by principal component analysis were performed to study the conformational change of MDM2 induced by p53 and two inhibitor (P4 and MI63a) bindings. The results show that the hydrophobic cleft of MDM2 is very flexible and adaptive to different structural binding partners. The cleft tends to become wider and more stable as MDM2 binds to the three binding partners, while unbound MDM2 shows a narrower and pretty flexible cleft, which agrees with recent experimental data and theoretical studies. It was also found that the binding of P4 and p53 stabilizes the motion of the loop L2 linking the helix α2 and β strand (β3), but the presence of MI63a makes the motion of L2 disordered. In addition, the binding free energies of the three partners to MDM2 were calculated using molecular mechanics generalized Born surface area to explain the binding modes of these three partners to MDM2. This study will be helpful not only for better understanding the functional, concerted motion of MDM2, but also for the rational design of potent anticancer drugs targeting the p53-MDM2 interaction.

Molecular dynamic simulation insights into the normal state and restoration of p53 function

As a tumor suppressor protein, p53 plays a crucial role in the cell cycle and in cancer prevention. Almost 50 percent of all human malignant tumors are closely related to a deletion or mutation in p53. The activity of p53 is inhibited by over-active celluar antagonists, especially by the over-expression of the negative regulators MDM2 and MDMX. Protein-protein interactions, or post-translational modifications of the C-terminal negative regulatory domain of p53, also regulate its tumor suppressor activity. Restoration of p53 function through peptide and small molecular inhibitors has become a promising strategy for novel anti-cancer drug design and development. Molecular dynamics simulations have been extensively applied to investigate the conformation changes of p53 induced by protein-protein interactions and protein-ligand interactions, including peptide and small molecular inhibitors. This review focuses on the latest MD simulation research, to provide an overview of the current understanding of interactions between p53 and its partners at an atomic level.

Simulating molecular mechanisms of the MDM2-mediated regulatory interactions: a conformational selection model of the MDM2 lid dynamics

Diversity and complexity of MDM2 mechanisms govern its principal function as the cellular antagonist of the p53 tumor suppressor. Structural and biophysical studies have demonstrated that MDM2 binding could be regulated by the dynamics of a pseudo-substrate lid motif. However, these experiments and subsequent computational studies have produced conflicting mechanistic models of MDM2 function and dynamics. We propose a unifying conformational selection model that can reconcile experimental findings and reveal a fundamental role of the lid as a dynamic regulator of MDM2-mediated binding. In this work, structure, dynamics and energetics of apo-MDM2 are studied as a function of posttranslational modifications and length of the lid. We found that the dynamic equilibrium between "closed" and "semi-closed" lid forms may be a fundamental characteristic of MDM2 regulatory interactions, which can be modulated by phosphorylation, phosphomimetic mutation as well as by the lid size. Our results revealed that these factors may regulate p53-MDM2 binding by fine-tuning the thermodynamic equilibrium between preexisting conformational states of apo-MDM2. In agreement with NMR studies, the effect of phosphorylation on MDM2 interactions was more pronounced with the truncated lid variant that favored the thermodynamically dominant closed form. The phosphomimetic mutation S17D may alter the lid dynamics by shifting the thermodynamic equilibrium towards the ensemble of "semi-closed" conformations. The dominant "semi-closed" lid form and weakened dependence on the phosphorylation seen in simulations with the complete lid can provide a rationale for binding of small p53-based mimetics and inhibitors without a direct competition with the lid dynamics. The results suggested that a conformational selection model of preexisting MDM2 states may provide a robust theoretical framework for understanding MDM2 dynamics. Probing biological functions and mechanisms of MDM2 regulation would require further integration of computational and experimental studies and may help to guide drug design of novel anti-cancer therapeutics.

Computational studies of difference in binding modes of peptide and non-peptide inhibitors to MDM2/MDMX based on molecular dynamics simulations

Inhibition of p53-MDM2/MDMX interaction is considered to be a promising strategy for anticancer drug design to activate wild-type p53 in tumors. We carry out molecular dynamics (MD) simulations to study the binding mechanisms of peptide and non-peptide inhibitors to MDM2/MDMX. The rank of binding free energies calculated by molecular mechanics generalized Born surface area (MM-GBSA) method agrees with one of the experimental values. The results suggest that van der Waals energy drives two kinds of inhibitors to MDM2/MDMX. We also find that the peptide inhibitors can produce more interaction contacts with MDM2/MDMX than the non-peptide inhibitors. Binding mode predictions based on the inhibitor-residue interactions show that the π–π, CH–π and CH–CH interactions dominated by shape complimentarity, govern the binding of the inhibitors in the hydrophobic cleft of MDM2/MDMX. Our studies confirm the residue Tyr99 in MDMX can generate a steric clash with the inhibitors due to energy and structure. This finding may theoretically provide help to develop potent dual-specific or MDMX inhibitors.

Molecular modeling and molecular dynamics simulation studies on pyrrolopyrimidine-based α-helix mimetic as dual inhibitors of MDM2 and MDMX

Inhibition of the interactions between the tumor suppressor protein p53 and its negative regulators, the MDM2 and MDMX oncogenic proteins, is increasingly gaining interest in cancer therapy and drug design. In this study, we carry out molecular docking, molecular dynamics (MD) simulations, and molecular mechanics Poisson-Boltzmann and generalized Born/surface area (MM-PB/GBSA) binding free energy calculations on an active compound 3a and an inactive compound NC-1, which share a common pyrrolopyrimidine-based scaffold. MD simulations and MM-PB/GBSA calculations show that the compound NC-1 may not bind to MDM2 and MDMX, in agreement with the experimental results. Detailed MM-PB/GBSA calculations on the MDM2-3a and MDMX-3a complexes unravel that the binding free energies are similar for the two complexes. Furthermore, the van der Waals energy is the largest component of the binding free energy for both complexes, which indicates that the interactions between the compound 3a and MDM2 and MDMX are dominated by shape complementarity. In addition, the analysis of individual residue contribution and protein-ligand binding mode show that the three functional groups on R₁, R₂, and R₃ of the compound 3a can mimic the spatial orientation of the side chains of Phe19, Trp23, and Leu26 of p53, respectively. The obtained computational results suggest that the compound 3a can act as a dual inhibitor of MDM2-p53 and MDMX-p53 interactions, consistent with the experimental results.