Modulation of p53 binding to MDM2: computational studies reveal important roles of Tyr100

A high-resolution structural characterization and physicochemical study of how a peptoid binds to an oncoprotein MDM2

Peptoids are a promising drug modality targeting disease-related proteins, but how a peptoid engages in protein binding is poorly understood. This is primarily due to a lack of high-resolution peptoid-protein complex structures and systematic physicochemical studies. Here, we present the first crystal structure of a peptoid bound to a protein, providing high-resolution structural information about how a peptoid binds to a protein. We previously reported a rigid peptoid, oligo(N-substituted alanine) (oligo-NSA), and developed an oligo-NSA-type peptoid that binds to MDM2. X-ray crystallographic analysis of the peptoid bound to MDM2 showed that the peptoid recognizes the MDM2 surface predominantly through the interaction of the N-substituents, while the main chain acts as a scaffold. Additionally, conformational, thermodynamic, and kinetic analysis of the peptoid and its derivatives with a less rigid main chain revealed that rigidification of the peptoid main chain contributes to improving the protein binding affinity. This improvement is thermodynamically attributed to an increased magnitude of the binding enthalpy change, and kinetically to an increased association rate and decreased dissociation rate. This study provides invaluable insights into the design of protein-targeting peptoids.

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.

Hitting on the move: Targeting intrinsically disordered protein states of the MDM2-p53 interaction

Intrinsically disordered proteins are an emerging class of proteins without a folded structure and currently disorder-based drug targeting remains a challenge. p53 is the principal regulator of cell division and growth whereas MDM2 consists its main negative regulator. The MDM2-p53 recognition is a dynamic and multistage process that amongst other, employs the dissociation of a transient α-helical N-terminal ''lid'' segment of MDM2 from the proximity of the p53-complementary interface. Several small molecule inhibitors have been reported to inhibit the formation of the p53-MDM2 complex with the vast majority mimicking the p53 residues Phe19, Trp23 and Leu26. Recently, we have described the transit from the 3-point to 4-point pharmacophore model stabilizing this intrinsically disordered N-terminus by increasing the binding affinity by a factor of 3. Therefore, we performed a thorough SAR analysis, including chiral separation of key compound which was evaluated by FP and 2D NMR. Finally, p53-specific anti-cancer activity towards p53-wild-type cancer cells was observed for several representative compounds.

Role of the N-terminal lid in regulating the interaction of phosphorylated MDMX with p53

Murine double minute 4 protein (MDMX) is crucial for the regulation of the tumor suppressor protein p53. Phosphorylation of the N-terminal domain of MDMX is thought to affect its binding with the transactivation domain of p53, thus playing a role in p53 regulation. In this study, the effects of MDMX phosphorylation on the binding of p53 were investigated using molecular dynamics simulations. It is shown that in addition to the previously proposed mechanism in which phosphorylated Y99 of MDMX inhibits p53 binding through steric clash with P27 of p53, the N-terminal lid of MDMX also appears to play an important role in regulating the phosphorylation-dependent interactions between MDMX and p53. In the proposed mechanism, phosphorylated Y99 aids in pulling the lid into the p53-binding pocket, thus inhibiting the binding between MDMX and p53. Rebinding of p53 appears to be facilitated by the subsequent phosphorylation of Y55, which draws the lid away from the binding pocket by electrostatic attraction of the lid's positively charged N-terminus. The ability to target these mechanisms for the proper regulation of p53 could have important implications for understanding cancer biology and for drug development.

Protein-peptide association kinetics beyond the seconds timescale from atomistic simulations

Understanding and control of structures and rates involved in protein ligand binding are essential for drug design. Unfortunately, atomistic molecular dynamics (MD) simulations cannot directly sample the excessively long residence and rearrangement times of tightly binding complexes. Here we exploit the recently developed multi-ensemble Markov model framework to compute full protein-peptide kinetics of the oncoprotein fragment 25-109Mdm2 and the nano-molar inhibitor peptide PMI. Using this system, we report, for the first time, direct estimates of kinetics beyond the seconds timescale using simulations of an all-atom MD model, with high accuracy and precision. These results only require explicit simulations on the sub-milliseconds timescale and are tested against existing mutagenesis data and our own experimental measurements of the dissociation and association rates. The full kinetic model reveals an overall downhill but rugged binding funnel with multiple pathways. The overall strong binding arises from a variety of conformations with different hydrophobic contact surfaces that interconvert on the milliseconds timescale.

Efficient Atomistic Simulation of Pathways and Calculation of Rate Constants for a Protein-Peptide Binding Process: Application to the MDM2 Protein and an Intrinsically Disordered p53 Peptide

The characterization of protein binding processes - with all of the key conformational changes - has been a grand challenge in the field of biophysics. Here, we have used the weighted ensemble path sampling strategy to orchestrate molecular dynamics simulations, yielding atomistic views of protein-peptide binding pathways involving the MDM2 oncoprotein and an intrinsically disordered p53 peptide. A total of 182 independent, continuous binding pathways were generated, yielding a kon that is in good agreement with experiment. These pathways were generated in 15 days using 3500 cores of a supercomputer, substantially faster than would be possible with "brute force" simulations. Many of these pathways involve the anchoring of p53 residue F19 into the MDM2 binding cleft when forming the metastable encounter complex, indicating that F19 may be a kinetically important residue. Our study demonstrates that it is now practical to generate pathways and calculate rate constants for protein binding processes using atomistic simulation on typical computing resources.

Hydrophobic Interactions Are a Key to MDM2 Inhibition by Polyphenols as Revealed by Molecular Dynamics Simulations and MM/PBSA Free Energy Calculations

p53, a tumor suppressor protein, has been proven to regulate the cell cycle, apoptosis, and DNA repair to prevent malignant transformation. MDM2 regulates activity of p53 and inhibits its binding to DNA. In the present study, we elucidated the MDM2 inhibition potential of polyphenols (Apigenin, Fisetin, Galangin and Luteolin) by MD simulation and MM/PBSA free energy calculations. All polyphenols bind to hydrophobic groove of MDM2 and the binding was found to be stable throughout MD simulation. Luteolin showed the highest negative binding free energy value of -173.80 kJ/mol followed by Fisetin with value of -172.25 kJ/mol. It was found by free energy calculations, that hydrophobic interactions (vdW energy) have major contribution in binding free energy.

Recognition Dynamics of p53 and MDM2: Implications for Peptide Design

Peptides that inhibit MDM2 and attenuate MDM2-p53 interactions, thus activating p53, are currently being pursued as anticancer drug leads for tumors harboring wild type p53. The thermodynamic determinants of peptide-MDM2 interactions have been extensively studied. However, a detailed understanding of the dynamics that underlie these interactions is largely missing. In this study, we explore the kinetics of the binding of a set of peptides using Brownian dynamics simulations. We systematically investigate the effect of peptide C-terminal substitutions (Ser, Ala, Asn, Pro) of a Q16ETFSDLWKLLP27 p53-based peptide and a M1PRFMDYWEGLN12 12/1 phage-derived peptide on their interaction dynamics with MDM2. The substitutions modulate peptide residence times around the MDM2 protein. In particular, the highest affinity peptide, Q16ETFSDLWKLLS27, has the longest residence time (t ∼ 25 μs) around MDM2, suggesting its potentially important contribution to binding affinity. The binding of the p53-based peptides appears to be kinetically driven while that of the phage-derived series appears to be thermodynamically driven. The phage-derived peptides were found to adopt distinctly different modes of interaction with the MDM2 protein compared to their p53-based counterparts. The p53-based peptides approach the N-terminal region of the MDM2 protein with the peptide C-terminal end oriented toward the protein, while the M1PRFMDYWEGLN12-based peptides adopt the reverse orientation. To probe the determinants of this switch in orientation, a designed mutant of the phage-derived peptide, R3E (M1PEFMDYWEGLN12), was simulated and found to adopt the orientation adopted by the p53-based peptides and also to result in almost a 5-fold increase in the peptide residence time (∼120 μs) relative to the p53-based peptides. On this basis, we suggest that the R3E mutant phage-derived peptide has a higher affinity for MDM2 than the p53-based peptides and would therefore, competitively inhibit MDM2-p53. The study, therefore, provides a novel computational framework for kinetics-based lead optimization for anticancer drug development strategies.

A spatiotemporal characterization of the effect of p53 phosphorylation on its interaction with MDM2

The interaction of p53 and MDM2 is modulated by the phosphorylation of p53. This mechanism is key to activating p53, yet its molecular determinants are not fully understood. To study the spatiotemporal characteristics of this molecular process we carried out Brownian dynamics simulations of the interactions of the MDM2 protein with a p53 peptide in its wild type state and when phosphorylated at Thr18 (pThr18) and Ser20 (pSer20). We found that p53 phosphorylation results in concerted changes in the topology of the interaction landscape in the diffusively bound encounter complex domain. These changes hinder phosphorylated p53 peptides from binding to MDM2 well before reaching the binding site. The underlying mechanism appears to involve shift of the peptide away from the vicinity of the MDM2 protein, peptide reorientation, and reduction in peptide residence time relative to wild-type p53 peptide. pThr18 and pSr20 p53 peptides experience reduction in residence times by factors of 13.6 and 37.5 respectively relative to the wild-type p53 peptide, indicating a greater role for Ser20 phosphorylation in abrogating p53 MDM2 interactions. These detailed insights into the effect of phosphorylation on molecular interactions are not available from conventional experimental and theoretical approaches and open up new avenues that incorporate molecular interaction dynamics, for stabilizing p53 against MDM2, which is a major focus of anticancer drug lead development.

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.

Structure of a stapled peptide antagonist bound to nutlin-resistant Mdm2

As key negative regulator of the p53 tumour suppressor, Mdm2 is an attractive therapeutic target. Small molecules such as Nutlin have been developed to antagonise Mdm2, resulting in p53-dependent death of tumour cells. We have recently described a mutation in Mdm2 (M62A), which precludes binding of Nutlin, but not p53. This Nutlin-resistant variant is not, however, refractory to binding and inhibition by stapled peptide antagonists targeting the same region of Mdm2. A detailed understanding of how stapled peptides are recalcitrant to Mdm2 mutations conferring Nutlin-resistance will aid in the further development of potent Mdm2 antagonists. Here, we report the 2.00 Å crystal structure of a stapled peptide antagonist bound to Nutlin resistant Mdm2. The stapled peptide relies on an extended network of interactions along the hydrophobic binding cleft of Mdm2 for high affinity binding. Additionally, as seen in other stapled peptide structures, the hydrocarbon staple itself contributes to binding through favourable interactions with Mdm2. The structure highlights the intrinsic plasticity present in both Mdm2 and the hydrocarbon staple moiety, and can be used to guide future iterations of both small molecules and stapled peptides for improved antagonists of Mdm2.

Identification of a second Nutlin-3 responsive interaction site in the N-terminal domain of MDM2 using hydrogen/deuterium exchange mass spectrometry

MDM2 is a multidomain protein that functions as an E3 ubiquitin ligase, transcription repressor, mRNA-binding protein, translation factor, and molecular chaperone. The small molecule Nutlin-3 has been engineered to bind to the N-terminal hydrophobic pocket domain of MDM2. This binding of Nutlin-3 has two consequences: (i) antagonistic effects through competitive disruption of the MDM2-p53 complex and (ii) agonist effects that allosterically stabilize MDM2 protein-protein interactions that increase p53 ubiquitination as well as nucleophosmin deoligomerization. We present a methodology using a hydrogen/deuterium (H/D) exchange platform that measures Nutlin-3 binding to the N-terminal domain of MDM2 (MDM2(1-126)) in order to begin to develop dynamic assays that evaluate MDM2 allostery. In order to localize the regions in MDM2 being suppressed by Nutlin-3, MDM2 was incubated with the ligand and H/D amide exchange was measured after pepsin digestion. One dynamic segment containing amino acids 55-60 exhibited slower deuterium exchange after Nutlin-3 binding, reflecting ligand binding within the hydrophobic pocket. However, another dominant suppression of H/D exchange was observed in a motif from amino acids 103-107 that reflects surface hydrophobic residues surrounding the hydrophobic pocket of MDM2. In order to explore the consequences of this latter Nutlin-3 interaction site on MDM2, the Y104G and L107G mutant series was constructed. The MDM2(Y104G) and MDM2(L107G) mutants were fully active in p53 binding. However, the authentic p53-derived peptide:MDM2(Y104G) complex exhibited partial resistance to Nutlin-3 inhibition, while the p53-mimetic 12.1 peptide:MDM2(Y104G) complex retained normal Nutlin-3 responsiveness. These data reveal the existence of a second functional Nutlin-3-binding site in a surface hydrophobic patch of MDM2, flanking the hydrophobic pocket. This reveals two modes of peptide binding by MDM2 and highlights the utility of H/D exchange as an assay for measuring allosteric effects in MDM2.

Adaptability in protein structures: structural dynamics and implications in ligand design

The basic framework of understanding the mechanisms of protein functions is achieved from the knowledge of their structures which can model the molecular recognition. Recent advancement in the structural biology has revealed that in spite of the availability of the structural data, it is nontrivial to predict the mechanism of the molecular recognition which progresses via situation-dependent structural adaptation. The mutual selectivity of protein-protein and protein-ligand interactions often depends on the modulations of conformations empowered by their inherent flexibility, which in turn regulates the function. The mechanism of a protein's function, which used to be explained by the ideas of 'lock and key' has evolved today as the concept of 'induced fit' as well as the 'population shift' models. It is felt that the 'dynamics' is an essential feature to take into account for understanding the mechanism of protein's function. The design principles of therapeutic molecules suffer from the problems of plasticity of the receptors whose binding conformations are accurately not predictable from the prior knowledge of a template structure. On the other hand, flexibility of the receptors provides the opportunity to improve the binding affinity of a ligand by suitable substitution that will maximize the binding by modulating the receptors surface. In this paper, we discuss with example how the protein's flexibility is correlated with its functions in various systems, revealing the importance of its understanding and for making applications. We also highlight the methodological challenges to investigate it computationally and to account for the flexible nature of the molecules in drug design.

On the origin of the stereoselective affinity of Nutlin-3 geometrical isomers for the MDM2 protein

The stereoselective affinity of small-molecule binding to proteins is typically broadly explained in terms of the thermodynamics of the final bound complex. Using Brownian dynamics simulations, we show that the preferential binding of the MDM2 protein to the geometrical isomers of Nutlin-3, an effective anticancer lead that works by inhibiting the interaction between the proteins p53 and MDM2, can be explained by kinetic arguments related to the formation of the MDM2:Nutlin-3 encounter complex. This is a diffusively bound state that forms prior to the final bound complex. We find that the MDM2 protein stereoselectivity for the Nutlin-3a enantiomer stems largely from the destabilization of the encounter complex of its mirror image enantiomer Nutlin-3b, by the K70 residue that is located away from the binding site. On the other hand, the trans-Nutlin-3a diastereoisomer exhibits a shorter residence time in the vicinity of MDM2 compared with Nutlin-3a due to destabilization of its encounter complex by the collective interaction of pairs of charged residues on either side of the binding site: Glu25 and Lys51 on one side, and Lys94 and Arg97 on the other side. This destabilization is largely due to the electrostatic potential of the trans-Nutlin-3a isomer being largely positive over extended continuous regions around its structure, which are otherwise well-identified into positive and negative regions in the case of the Nutlin-3a isomer. Such rich insight into the binding processes underlying biological selectivity complements the static view derived from the traditional thermodynamic analysis of the final bound complex. This approach, based on an explicit consideration of the dynamics of molecular association, suggests new avenues for kinetics-based anticancer drug development and discovery.

Transient protein states in designing inhibitors of the MDM2-p53 interaction

Reactivation of p53 by release of the functional protein from its inhibition by MDM2 provides an efficient, nongenotoxic approach to a wide variety of cancers. We present the cocrystal structures of two complexes of MDM2 with inhibitors based on 6-chloroindole scaffolds. Both molecules bound to a distinct conformational state of MDM2 with nM-μM affinities. In contrast to other structurally characterized antagonists, which mimic three amino acids of p53 (Phe19, Trp23, and Leu26), the compounds induced an additional hydrophobic pocket on the MDM2 surface and unveiled a four-point binding mode. The enlarged interaction interface of the inhibitors resulted in extension of small molecules binding toward the "lid" segment of MDM2 (residues 19-23)--a nascent element that interferes with p53 binding. As supported by protein engineering and molecular dynamics studies, employing these unstable elements of MDM2 provides an efficient and yet unexplored alternative in development of MDM2-p53 association inhibitors.

On the interaction mechanisms of a p53 peptide and nutlin with the MDM2 and MDMX proteins: a Brownian dynamics study

The interaction of p53 with its regulators MDM2 and MDMX plays a major role in regulating the cell cycle. Inhibition of this interaction has become an important therapeutic strategy in oncology. Although MDM2 and MDMX share a very high degree of sequence/structural similarity, the small-molecule inhibitor nutlin appears to be an efficient inhibitor only of the p53-MDM2 interaction. Here, we investigate the mechanism of interaction of nutlin with these two proteins and contrast it with that of p53 using Brownian dynamics simulations. In contrast to earlier attempts to examine the bound states of the partners, here we locate initial reaction events in these interactions by identifying the regions of space around MDM2/MDMX, where p53/nutlin experience associative encounters with prolonged residence times relative to that in bulk solution. We find that the initial interaction of p53 with MDM2 is long-lived relative to nutlin, but, unlike nutlin, it takes place at the N- and C termini of the MDM2 protein, away from the binding site, suggestive of an allosteric mechanism of action. In contrast, nutlin initially interacts with MDM2 directly at the clefts of the binding site. The interaction of nutlin with MDMX, however, is very short-lived compared with MDM2 and does not show such direct initial interactions with the binding site. Comparison of the topology of the electrostatic potentials of MDM2 and MDMX and the locations of the initial encounters with p53/nutlin in tandem with structure-based sequence alignment revealed that the origin of the diminished activity of nutlin toward MDMX relative to MDM2 may stem partly from the differing topologies of the electrostatic potentials of the two proteins. Glu25 and Lys51 residues underpin these topological differences and appear to collectively play a key role in channelling nutlin directly toward the binding site on the MDM2 surface and are absent in MDMX. The results, therefore, provide new insight into the mechanism of p53/nutlin interactions with MDM2 and MDMX and could potentially have a broader impact on anticancer drug optimization strategies.

An iTRAQ proteomics screen reveals the effects of the MDM2 binding ligand Nutlin-3 on cellular proteostasis

Mouse double minute 2 (MDM2) participates in protein synthesis, folding, and ubiquitin-mediated degradation and is therefore a proteostasic hub protein. The MDM2 interactome contains over 100 proteins, yet stratification of dominant MDM2-interacting proteins has not been achieved. 8-plex iTRAQ (nanoLC-MS/MS) of MCF7 cells treated with the MDM2-binding ligand Nutlin-3 identified the most abundant cellular protein changes over early time points; 1,323 unique proteins were identified including 35 with altered steady-state levels within 2 h of Nutlin-3 treatment, identifying a core group of MDM2 related proteins. Six of these proteins were previously identified MDM2 interactors, and the effects of Nutlin-3 on the MDM2-nucleophosmin interaction (NPM) was further validated. This revealed that Nutlin-3 mediates the in vivo conversion of NPM from an oligomer to a monomer as an MDM2-dependent phenomenon, with Nutlin-3 stimulating MDM2 binding to a peptide motif derived from the oligomerization interface of NPM. These data form the first proteomic screen of Nutlin-3 in cells whereby we (i) identify the most abundant MDM2-interacting proteins whose steady-state levels change early after Nutlin-3 treatment; (ii) identify the first protein apart from p53, nucleophosmin (NPM), whose interaction with MDM2 can be stimulated allosterically by Nutlin-3; and (iii) raise the possibility that Nutlin-3 might act as a general agonist of other MDM2 protein-protein interactions.

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.

Why is F19Ap53 unable to bind MDM2? Simulations suggest crack propagation modulates binding

Why doesn't the F19A mutant of p53 bind to MDM2? Binding thermodynamics have suggested that the loss of packing interactions upon mutating Phe into Ala sidechain results in destabilizing the binding free energy between p53 and MDM2. Does this mutation also modulate the initial recognition between p53 and MDM2? We look at atomistic computer simulations of the process of the initial encounter between wild type p53 peptide and its F19A mutant with the N-terminal domain of MDM2. These simulations show that binding is characterized by a complex multistep process. It starts with the capture of F19 of wild type p53 by certain residues in the MDM2 binding pocket. This initial step anchors the peptide onto the surface of MDM2, and with the consequent reduction in the search space of the peptide, the peptide docks into the partially occluded surface of MDM2. This is similar to a crack forming in an otherwise occluded hydrophobic cavity in MDM2, and the peptide, docked through F19, modulates the propagation of this crack, which subsequently results in the stepwise docking of the rest of the peptide through insertions of W23 and L26. The lack of the bulky sidechain of F in the F19A mutant results in the absence of the initial "grasp" complex, and hence the mutant peptide diffuses randomly on the surface of MDM2 without binding. This is the first such demonstration of the possibility that a "kinetic" effect may partly underlie the destabilized thermodynamics of binding of F19A and is a feature that appears to be conserved in evolution. The observations by Wallace et al. (Mol Cell 2006; 23:251-63) that despite the inability of F19A to bind at the N-terminal domain of MDM2, it gets ubiquitinated, can now be partly understood based on a mechanism whereby the occupation of the binding pocket by ligands/peptides induces, via crack propagation and the dynamics of gatekeeper Y100, the ubiquitination signal for interactions between the acidic domain of MDM2 and the DNA binding domain of p53.

Characterization of the Ligand Receptor Encounter Complex and Its Potential for in Silico Kinetics-Based Drug Development

The study of drug-receptor interactions has largely been framed in terms of the equilibrium thermodynamic binding affinity, an in vitro measure of the stability of the drug-receptor complex that is commonly used as a proxy measure of in vivo biological activity. In response to the growing realization of the importance of binding kinetics to in vivo drug activity we present a computational methodology for the kinetic characterization of drug-receptor interactions in terms of the encounter complex. Using trajectory data from multiple Brownian dynamics simulations of ligand diffusion, we derive the spatial density of the ligand around the receptor and show how it can be quantitatively partitioned into different basins of attraction. Numerical integration of the ligand densities within the basins can be used to estimate the residence time of the ligand within these diffusive binding sites. Simulations of two structurally similar inhibitors of Hsp90 exhibit diffusive binding sites with similar spatial structure but with different ligand residence times. In contrast, a pair of structurally dissimilar inhibitors of MDM2, a peptide and a small molecule, exhibit spatially distinct basins of attraction around the receptor, which in turn reveal differences in ligand orientational order. Thus, our kinetic approach provides microscopic details of drug-receptor dynamics that provide novel insight into the observed differences in the thermodynamic binding affinities for the two inhibitors, such as the differences in the entropic contributions to binding. The characterization of the encounter complex, in terms of the structure, topology, and dynamics of diffusive binding sites, offers a new perspective on ligand-receptor interactions and the potential for greater insight into drug action. The method, which requires no prior knowledge of the bound state, is a first step toward the incorporation of ligand kinetics into in silico drug development protocols.

C-terminal substitution of MDM2 interacting peptides modulates binding affinity by distinctive mechanisms

The complex between the proteins MDM2 and p53 is a promising drug target for cancer therapy. The residues 19–26 of p53 have been biochemically and structurally demonstrated to be a most critical region to maintain the association of MDM2 and p53. Variation of the amino acid sequence in this range obviously alters the binding affinity. Surprisingly, suitable substitutions contiguous to this region of the p53 peptides can yield tightly binding peptides. The peptide variants may differ by a single residue that vary little in their structural conformations and yet are characterized by large differences in their binding affinities. In this study a systematic analysis into the role of single C-terminal mutations of a 12 residue fragment of the p53 transactivation domain (TD) and an equivalent phage optimized peptide (12/1) were undertaken to elucidate their mechanistic and thermodynamic differences in interacting with the N-terminal of MDM2. The experimental results together with atomistically detailed dynamics simulations provide insight into the principles that govern peptide design protocols with regard to protein-protein interactions and peptidomimetic design.

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.

Anchoring intrinsically disordered proteins to multiple targets: lessons from N-terminus of the p53 protein

Anchor residues, which are deeply buried upon binding, play an important role in protein–protein interactions by providing recognition specificity and facilitating the binding kinetics. Up to now, studies on anchor residues have been focused mainly on ordered proteins. In this study, we investigated anchor residues in intrinsically disordered proteins (IDPs) which are flexible in the free state. We identified the anchor residues of the N-terminus of the p53 protein (Glu17–Asn29, abbreviated as p53N) which are involved in binding with two different targets (MDM2 and Taz2), and analyzed their side chain conformations in the unbound states. The anchor residues in the unbound p53N were found to frequently sample conformations similar to those observed in the bound complexes (i.e., Phe19, Trp23, and Leu26 in the p53N-MDM2 complex, and Leu22 in the p53N-Taz2 complex). We argue that the bound-like conformations of the anchor residues in the unbound state are important for controlling the specific interactions between IDPs and their targets. Further, we propose a mechanism to account for the binding promiscuity of IDPs in terms of anchor residues and molecular recognition features (MoRFs).

Modeling the relationship between the p53 C-terminal domain and its binding partners using molecular dynamics

Fifty percent of all cancer cases result from mutations of the TP53 gene, which encodes the tumor suppressor p53, and it is hypothesized that the p53-mediated checkpoint pathway is compromised in most of the remaining cases. The p53 C-terminal domain (CTD) is an important site of p53 regulation but by nature is difficult to study, as it is intrinsically disordered. In this study, we performed molecular dynamics simulations on the p53 CTD and five known regulatory binding partners. We identified distinct trends in fluctuation within and around the p53 CTD binding site on each partner demonstrating a behavior that facilitates association. Further, we present evidence that the size of the hydrophobic pocket in each p53 CTD binding site governs the secondary structure of the p53 CTD when in the bound state. This information will be useful for predicting new binding partners for the p53 CTD, identifying interacting regions within other known partners, and discovering inhibitors that provide additional points of control over p53 activity.

Towards a career in bioinformatics

The 2009 annual conference of the Asia Pacific Bioinformatics Network (APBioNet), Asia's oldest bioinformatics organisation from 1998, was organized as the 8^th International Conference on Bioinformatics (InCoB), Sept. 9-11, 2009 at Biopolis, Singapore. InCoB has actively engaged researchers from the area of life sciences, systems biology and clinicians, to facilitate greater synergy between these groups. To encourage bioinformatics students and new researchers, tutorials and student symposium, the Singapore Symposium on Computational Biology (SYMBIO) were organized, along with the Workshop on Education in Bioinformatics and Computational Biology (WEBCB) and the Clinical Bioinformatics (CBAS) Symposium. However, to many students and young researchers, pursuing a career in a multi-disciplinary area such as bioinformatics poses a Himalayan challenge. A collection to tips is presented here to provide signposts on the road to a career in bioinformatics. An overview of the application of bioinformatics to traditional and emerging areas, published in this supplement, is also presented to provide possible future avenues of bioinformatics investigation. A case study on the application of e-learning tools in undergraduate bioinformatics curriculum provides information on how to go impart targeted education, to sustain bioinformatics in the Asia-Pacific region. The next InCoB is scheduled to be held in Tokyo, Japan, Sept. 26-28, 2010.