Molecular modelling,synthesis and biological evaluation of peptide inhibitors as anti-angiogenic agent targeting neuropilin-1 for anticancer application

Vascular endothelial growth factor (VEGF) and its co-receptor neuropilin-1 (NRP-1) are important targets of many pro-angiogenic factors. In this study, nine peptides were synthesized and evaluated for their molecular interaction with NRP-1 and compared to our previous peptide ATWLPPR. Docking study showed that the investigated peptides shared the same binding region as shown by tuftsin known to bind selectively to NRP-1. Four pentapeptides (DKPPR, DKPRR, TKPPR and TKPRR) and a hexapeptide CDKPRR demonstrated good inhibitory activity against NRP-1. In contrast, peptides having arginine residue at sites other than the C-terminus exhibited low activity towards NRP-1 and this is confirmed by their inability to displace the VEGF₁₆₅ binding to NRP-1. Docking study also revealed that replacement of carboxyl to amide group at the C-terminal arginine of the peptide did not affect significantly the binding interaction to NRP-1. However, the molecular affinity study showed that these peptides have marked reduction in the activity against NRP-1. Pentapeptides having C-terminal arginine showed strong interaction and good inhibitory activity with NRP thus may be a good template for anti-angiogenic targeting agent.     

Mapping the Vif–A3G interaction using peptide arrays: A basis for anti-HIV lead peptides

Human apolipoprotein-B mRNA-editing catalytic polypeptide-like 3G (A3G) is a cytidine deaminase that restricts retroviruses, endogenous retro-elements and DNA viruses. A3G plays a key role in the anti-HIV-1 innate cellular immunity. The HIV-1 Vif protein counteracts A3G mainly by leading A3G towards the proteosomal machinery and by direct inhibition of its enzymatic activity. Both activities involve direct interaction between Vif and A3G. Disrupting the interaction between A3G and Vif may rescue A3G antiviral activity and inhibit HIV-1 propagation. Here, mapping the interaction sites between A3G and Vif by peptide array screening revealed distinct regions in Vif important for A3G binding, including the N-terminal domain (NTD), C-terminal domain (CTD) and residues 83–99. The Vif-binding sites in A3G included 12 different peptides that showed strong binding to either full-length Vif, Vif CTD or both. Sequence similarity was found between Vif-binding peptides from the A3G CTD and NTD. A3G peptides were synthesized and tested for their ability to counteract Vif action. A3G 211–225 inhibited HIV-1 replication in cell culture and impaired Vif dependent A3G degradation. In vivo co-localization of full-length Vif with A3G 211–225 was demonstrated by use of FRET. This peptide has the potential to serve as an anti-HIV-1 lead compound. Our results suggest a complex interaction between Vif and A3G that is mediated by discontinuous binding regions with different affinities.     

Structure-activity relationship study of tetrapeptide inhibitors of the Vascular Endothelial Growth Factor A binding to Neuropilin-1

Neuropilin-1 is considered as one of the key receptors responsible for signaling pathways involved in pathological angiogenesis necessary for tumor progression, therefore targeting of VEGF₁₆₅ binding to NRP-1 could be a relevant strategy for antiangiogenic treatment. It was shown before that the VEGF₁₆₅/NRP-1 interaction can be inhibited by short tetrapeptides with K/RXXR sequence.Here, we present a structure–activity relationship study of the systematic optimization of amino acid residues in positions 1–3 in the above tetrapeptides. All the 13 synthesized analogs possessed C-terminal arginine that is a necessary element for interaction with NRP-1. The obtained results of the inhibitory activity and modeling by molecular dynamics indicate that simultaneous interactions of the basic amino acid residues in position 1 and 4 (Arg) with Neuropilin-1 are crucial and their cooperation strongly affects the inhibitory activity. In addition, the binding strength is modulated by the flexibility of the peptide backbone (in the central part of the peptide), and the nature of the side chain of the amino acids at the second or third position. A dramatic decrease in the activity to the receptor is observed in flexible derivatives that are missing proline residues. The results described in this paper should prove useful for future studies aimed at establishing the best pharmacophore for inhibitors of VEGF₁₆₅ binding to NRP-1.     

Discovery of novel inhibitors of vascular endothelial growth factor-A–Neuropilin-1 interaction by structure-based virtual screening

Neuropilin-1 (NRP-1), one of the most important co-receptors of vascular endothelial growth factor-A (VEGF-A), increases its angiogenic action in several chronic diseases including cancer by increasing the activity of associated tyrosine kinase receptors, VEGFR1 and VEGFR2. Binding of VEGF-A to NRP-1 plays a critical role in pathological angiogenesis and tumor progression. Today, targeting this interaction is a validated approach to fight against angiogenesis-dependent diseases. Only anti-NRP-1 antibodies, peptide and peptidomimetic drug-candidates or hits have been developed thus far. In order to identify potent orally active small organic molecules various experimental and in silico approaches can be used. Here we report, novel promising small drug-like molecules disrupting the binding of VEGF-A₁₆₅ to NRP-1. We carried out structure-based virtual screening experiments using the ChemBridge compound collection on the VEGF-A₁₆₅ binding pocket of NRP-1. After docking and two rounds of similarity search computations, we identified 4 compounds that inhibit the biotinylated VEGF-A₁₆₅ binding to recombinant NRP-1 with K_i of about 10 μM. These compounds contain a common chlorobenzyloxy alkyloxy halogenobenzyl amine scaffold that can serve as a base for further development of new NRP-1 inhibitors.     

Identification of single amino acid residues essential for the binding of lipopolysaccharide (LPS) to LPS binding protein (LBP) residues 86-99 by using an Ala-scanning library.

Lipopolysaccharide binding protein (LBP) is a 60 kDa acute phase glycoprotein capable of binding to LPS of Gram-negative bacteria and facilitating its interaction with cellular receptors. This process is thought to be of great importance in systemic inflammatory reactions such as septic shock. A peptide corresponding to residues 86-99 of human LBP (LBP86-99) has been reported to bind specifically with high affinity the lipid A moiety of LPS and to inhibit the interaction of LPS with LBP. We identified essential amino acids in LBP86-99 for binding to LPS by using a peptide library corresponding to the Ala-scanning of human LBP residues 86-99. Amino acids Trp91 and Lys92 were indispensable for peptide-LPS interaction and inhibition of LBP-LPS binding. In addition, several alanine-substituted synthetic LBP-derived peptides inhibited LPS-LBP interaction. Substitution of amino acids Arg94, Lys95 and Phe98 by Ala increased the inhibitory effect. The mutant Lys95 was the most active in blocking LPS binding to LBP. These findings emphasize the importance of single amino acids in the LPS binding capacity of small peptides and may contribute to the development of new drugs for use in the treatment of Gram-negative bacterial sepsis.     

Binding modes of Bcl-2 homology 3 (BH3) peptides with anti-apoptotic protein A1 and redesign of peptide inhibitors: a computational study

The interaction between protein and peptide ligand is a challenging problem in molecular biology and drug design. The binding of the Bcl-2 homology 3 (BH3) peptide to the anti-apoptotic protein A1 was revealed as a critical step in the regulation of apoptosis. These BH3 peptides hold high structural similarity, but are diverse in their regulation abilities. Based on molecular simulations and MM-P(G)BSA methods, this work presented a detailed analysis on binding mechanism of the BH3 peptides derived from PUMA and BMF. Residue-level energy decomposition showed that the core regions of BH3 peptides maintain in stable helical conformations and the four conserved hydrophobic residues together with an invariant aspartic acid contribute the major driving force for binding, whereas their two terminal segments exhibit obvious flexibility and various binding modes. Such kind of behavior was suggested as the reason for binding diversity and selectivity of BH3 peptides. As a further step, several BH3-mimetic peptides have been redesigned by computational mutation. Those new peptides showed not only stronger affinities when binding to protein A1, as well transferable binding patterns at some specific positions. A long-range coupling effect was disclosed for BH3 peptides, side-chain orientation and binding contribution of terminal residues were even affected by mutations at large sequence interval. Overall, this work reports that the binding modes of BH3 peptides are primarily dependent on its two terminal segments. The computational methods applied herein are also demonstrated to be of great assistance in the rational design of peptide inhibitors.     

The MECP2-TRD domain interacts with the DNMT3A-ADD domain at the H3-tail binding site

The DNMT3A DNA methyltransferase and MECP2 methylation reader are highly expressed in neurons. Both proteins interact via their DNMT3A-ADD and MECP2-TRD domains, and the MECP2 interaction regulates the activity and subnuclear localization of DNMT3A. Here, we mapped the interface of both domains using peptide SPOT array binding, protein pull-down, equilibrium peptide binding assays, and structural analyses. The region D529-D531 on the surface of the ADD domain was identified as interaction point with the TRD domain. This includes important residues of the histone H3 N-terminal tail binding site to the ADD domain, explaining why TRD and H3 binding to the ADD domain is competitive. On the TRD domain, residues 214–228 containing K219 and K223 were found to be essential for the ADD interaction. This part represents a folded patch within the otherwise largely disordered TRD domain. A crystal structure analysis of ADD revealed that the identified H3/TDR lysine binding pocket is occupied by an arginine residue from a crystallographic neighbor in the ADD apoprotein structure. Finally, we show that mutations in the interface of ADD and TRD domains disrupt the cellular interaction of both proteins in NIH3T3 cells. In summary, our data show that the H3 peptide binding cleft of the ADD domain also mediates the interaction with the MECP2-TRD domain suggesting that this binding site may have a broader role also in the interaction of DNMT3A with other proteins leading to complex regulation options by competitive and PTM specific binding.     

Glycan-dependent binding of galectin-1 to neuropilin-1 promotes axonal regeneration after spinal cord injury

Following spinal cord injury (SCI), semaphorin 3A (Sema3A) prevents axonal regeneration through binding to the neuropilin-1 (NRP-1)/PlexinA4 receptor complex. Here, we show that galectin-1 (Gal-1), an endogenous glycan-binding protein, selectively bound to the NRP-1/PlexinA4 receptor complex in injured neurons through a glycan-dependent mechanism, interrupts the Sema3A pathway and contributes to axonal regeneration and locomotor recovery after SCI. Although both Gal-1 and its monomeric variant contribute to de-activation of microglia, only high concentrations of wild-type Gal-1 (which co-exists in a monomer–dimer equilibrium) bind to the NRP-1/PlexinA4 receptor complex and promote axonal regeneration. Our results show that Gal-1, mainly in its dimeric form, promotes functional recovery of spinal lesions by interfering with inhibitory signals triggered by Sema3A binding to NRP-1/PlexinA4 complex, supporting the use of this lectin for the treatment of SCI patients.     

Binding of a C-end rule peptide to the neuropilin-1 receptor: a molecular modeling approach

Neuropilin-1 (NRP-1) is a receptor that plays an essential role in angiogenesis, vascular permeability, and nervous system development. Previous studies have shown that peptides with an N-terminal Arg, especially peptides with the four-residue consensus sequence R/K/XXR/K, bind to NRP-1 cell surfaces. Peptides containing such consensus sequences promote binding and internalization into cells, while blocking the C-terminal Arg (or Lys) prevents the internalization. In this study, we use molecular dynamics simulations to model the structural properties of the NRP-1 complex with a prototypic CendR peptide, RPAR. Our simulations show that RPAR binds NRP-1 through specific interactions of the RPAR C-terminus: three hydrogen bonds and a salt bridge anchor the ligand in the receptor pocket. The modeling results were used as the starting point for a systematic computational study of new RPAR analogues based on chemical modifications of their natural amino acids. Comparison of the structural properties of the new peptide-receptor complexes with the original organization suggests that some of the analogues can increase the binding affinity while reducing the natural sensitivity of RXXR to endogenous proteases.     

Peptides derived from a secretory yeast library restore factor VIII activity in the presence of an inhibitory antibody

The development of autoantibodies against factor VIII represents one of the major complications in the treatment of hemophilia A patients. We have employed a novel library system to obtain peptides that specifically neutralize the interaction between factor VIII and these inhibitors. The random peptides are presented as carboxy-terminal extensions of the eukaryotic initiation factor 5a, an intracellular protein with a molecular mass of 18 kDa. These random peptides formed an unique binding site, as demonstrated by molecular simulations using the computer programs InsightII and GROMACS. The library was screened to identify peptides binding to the murine monoclonal anti-factor VIII antibody ESH8 and to inhibitors derived from patients with factor VIII antibodies. Ten peptides binding to ESH8 were identified. Their specificity was confirmed by displacement assays. Two peptides with the sequences STKTLGRPLHGPAGPVEGGALAGVAEDADLVTAVSGR and YHCKREDLTDRDATCALRQPPQAVRGLGPRVTAVSGR showed the ability to restore the factor VIII activity from 33% up to approximately 90% in functional tests performed in vitro. Three candidates for binding to factor VIII antibodies derived from four different patient's sera were achieved. Three fusion proteins with the peptide sequences PQLGSRRSTTPSLTFQNASWFPAGGPCARSNRG, SGSRQVCKLARSLQPF and WERGRRVGAQVRHARHLVARVLDGAGHQARLTAVNGP bound to inhibitors derived from different patients. Furthermore, two of the obtained fusion proteins with the peptide sequences RHWTALGPAPTHTCADLNYPLLS and WERGRRVGAQVRHARHLVARVLDGAGHQARLTAVNGP did also bind to the monoclonal antibody ESH8. This study demonstrates the potential of this system to identify peptides that inhibit the activity of potent inhibitory antibodies and also shows potential as a method for screening of bioactive peptides.     

Sub‐angstrom modeling of complexes between flexible peptides and globular proteins

A wide range of regulatory processes in the cell are mediated by flexible peptides that fold upon binding to globular proteins. Computational efforts to model these interactions are hindered by the large number of rotatable bonds in flexible peptides relative to typical ligand molecules, and the fact that different peptides assume different backbone conformations within the same binding site. In this study, we present Rosetta FlexPepDock, a novel tool for refining coarse peptide–protein models that allows significant changes in both peptide backbone and side chains. We obtain high resolution models, often of sub‐angstrom backbone quality, over an extensive and general benchmark that is based on a large nonredundant dataset of 89 peptide–protein interactions. Importantly, side chains of known binding motifs are modeled particularly well, typically with atomic accuracy. In addition, our protocol has improved modeling quality for the important application of cross docking to PDZ domains. We anticipate that the ability to create high resolution models for a wide range of peptide–protein complexes will have significant impact on structure‐based functional characterization, controlled manipulation of peptide interactions, and on peptide‐based drug design. Proteins 2010. © 2010 Wiley‐Liss, Inc.     

Molecular inhibition of telomerase recruitment using designer peptides: an in silico approach

Telomere holds special mechanism for solving end repair problems and maintaining genomic stability. Protection of telomeres 1 (POT1) which belongs to shelterin family is identified as a key protein that recruits telomerase by interacting with telomere repeat binding factors (TRB1-3). Since, deciphering the mechanism through which POT assembles telomerase is of great interest, computational approaches have been undertaken to understand the mechanism in a well- developed model system – Arabidopsis thaliana. As a first step, an untraditional approach was mediated to locate the active site residues on modeled AtPOT1b protein by interaction studies using molecular docking. To keep in trend with the recent developments, peptide construction and validation was promoted as the next step via molecular dynamics simulation studies. Finally, the validated peptides based on propensity score was evaluated for its efficacy as a potent inhibitor for POT and TRB1-3 interactions. The best peptide, namely, (1-2-d) out of 30 designed peptides, was proved to be vital inhibitor by weakening the interacting complexes.     

Binding site on human C-reactive Protein (CRP) recognized by the Leukocyte CRP-receptor

C-reactive protein (CRP), the prototypical inflammatory acute phase reactant in humans, interacts with monocytes and neutrophils via a specific receptor. To map the site on CRP recognized by the CRP receptor (CRP-R), synthetic peptides corresponding to the surface region on each of the five identical subunits were tested as competitors vs. [¹²⁵l]-CRP for cell binding. A peptide of residues 27–38 (TKPLKAFTVCLH) efficiently inhibited CRP binding when compared to other nonoverlapping peptides. This peptide was termed the cell-binding peptide (CB-Pep). The F(ab′)₂ of an IgG Ab to the CB-Pep specifically inhibited CRP binding upon reacting with the ligand. Competitive binding studies with synthetic peptides truncated from either the NH₂- or COOH-terminus of the CB-Pep revealed that the minimum length recognized by the CRP-R consisted of residues 31–36: KAFTVC. Conservative substitutions of residues within the CB-Pep indicated that the four residues AFTV were critical for CRP-R binding. The CB-Pep also inhibited induced superoxide generation by HL-60 granulocytes. The minimum length required for the inhibition was also KAFTVC; however, only Phe-33 and Leu-37 were critical residues in this assay. Anti-CB-Pep IgG Ab reacted more extensively with heat-modified CRP, suggesting that an altered conformation of CRP is preferentially recognized by the CRP-R. The results suggest that this contiguous sequence on a β-strand on one face of each of five subunits of the CRP pentamer serves as a unique recognition motif for inflammatory leukocytes. J. Cell. Biochem. 64:140–151. © 1997 Wiley-Liss, Inc.     

Structural insights into the role of the WW2 domain on tandem WW–PPxY motif interactions of oxidoreductase WWOX

Class I WW domains are present in many proteins of various functions and mediate protein interactions by binding to short linear PPxY motifs. Tandem WW domains often bind peptides with multiple PPxY motifs, but the interplay of WW–peptide interactions is not always intuitive. The WW domain–containing oxidoreductase (WWOX) harbors two WW domains: an unstable WW1 capable of PPxY binding and stable WW2 that cannot bind PPxY. The WW2 domain has been suggested to act as a WW1 domain chaperone, but the underlying mechanism of its chaperone activity remains to be revealed. Here, we combined NMR, isothermal calorimetry, and structural modeling to elucidate the roles of both WW domains in WWOX binding to its PPxY-containing substrate ErbB4. Using NMR, we identified an interaction surface between these two domains that supports a WWOX conformation compatible with peptide substrate binding. Isothermal calorimetry and NMR measurements also indicated that while binding affinity to a single PPxY motif is marginally increased in the presence of WW2, affinity to a dual-motif peptide increases 10-fold. Furthermore, we found WW2 can directly bind double-motif peptides using its canonical binding site. Finally, differential binding of peptides in mutagenesis experiments was consistent with a parallel N- to C-terminal PPxY tandem motif orientation in binding to the WW1–WW2 tandem domain, validating structural models of the interaction. Taken together, our results reveal the complex nature of tandem WW-domain organization and substrate binding, highlighting the contribution of WWOX WW2 to both protein stability and target binding.     

The C-terminal regulatory domain of p53 contains a functional docking site for cyclin A

Radiation injury to cells enhances C-terminal phosphorylation of p53 at both Ser315 and Ser392 in vivo, suggesting the existence of two cooperating DNA damage-responsive pathways that play a role in stimulating p53-dependent gene expression. Our previous data has shown that cyclin A-cdk2 is the major enzyme responsible for modifying p53 at Ser315 in vivo after irradiation damage and in this report we dissect the mechanism of cyclinA-cdk2 binding to and phosphorylation of p53. Although cyclin B(1)-dependent protein kinases can phosphorylate small peptides containing the Ser315 site, cyclin A-cdk2 does not phosphorylate such small peptides suggesting that additional determinants are required for cyclin A-cdk2 interaction with p53. Peptide competition studies have localized a cyclin A interaction site to a Lys381Lys382Leu383Met384Phe385 sequence within C-terminal negative regulatory domain of human p53. An alanine mutation at any one of four key positions abrogates the efficacy of a synthetic peptide containing this motif as an inhibitor of cyclin A-cdk2 phosphorylation of p53 protein. Single amino acid mutations of full-length p53 protein at Lys382, Leu383, or Phe385 decreases cyclin A-cdk2 dependent phosphorylation at Ser315. Cyclin B(1)-cdk2 complexes are not inhibited by KKLMF motif-containing peptides nor is p53 phosphorylation by cyclin B-cdk2 reduced by mutation of the cyclin A interaction site. These data identifying a KKLMF cyclin A docking site on p53 protein highlight a common cyclin A interaction motif that is shared between the tumour suppressor proteins pRb and p53.     

Mapping discontinuous protein‐binding sites via structure‐based peptide libraries: combining in silico and in vitro approaches

To perform their various functions, protein surfaces often have to interact with each other in a specific way. Usually, only parts of a protein are accessible and can act as binding sites. Because proteins consist of polypeptide chains that fold into complex three‐dimensional shapes, binding sites can be divided into two different types: linear sites that follow the primary amino acid sequence and discontinuous binding sites, which are made up of short peptide fragments that are adjacent in spatial proximity. Such discontinuous binding sites dominate protein–protein interactions, but are difficult to identify. To meet this challenge, we combined a computational, structure‐based approach and an experimental, high‐throughput method. SUPERFICIAL is a program that uses protein structures as input and generates peptide libraries to represent the protein's surface. A large number of the predicted peptides can be simultaneously synthesised applying the SPOT technology. The results of a binding assay subsequently help to elucidate protein–protein interactions; the approach is applicable to any kind of protein. The crystal structure of the complex of hen egg lysozyme with the well‐characterised murine IgG1 antibody HyHEL‐5 is available, and the complex is known to have a discontinuous binding site. Using SUPERFICIAL, the entire surface of lysozyme was translated into a peptide library that was synthesised on a cellulose membrane using the SPOT technology and tested against the HyHEL‐5 antibody. In this way, it was possible to identify two peptides (longest common sequence and peptide 19) that represented the discontinuous epitope of lysozyme. Copyright © 2012 John Wiley & Sons, Ltd.     

Synthetic peptides mimicking the binding site of human acetylcholinesterase for its inhibitor fasciculin 2

Molecules capable of mimicking protein binding and/or functional sites present useful tools for a range of biomedical applications, including the inhibition of protein–ligand interactions. Such mimics of protein binding sites can currently be generated through structure‐based design and chemical synthesis. Computational protein design could be further used to optimize protein binding site mimetics through rationally designed mutations that improve intermolecular interactions or peptide stability. Here, as a model for the study, we chose an interaction between human acetylcholinesterase (hAChE) and its inhibitor fasciculin‐2 (Fas) because the structure and function of this complex is well understood. Structure‐based design of mimics of the hAChE binding site for Fas yielded a peptide that binds to Fas at micromolar concentrations. Replacement of hAChE residues known to be essential for its interaction with Fas with alanine, in this peptide, resulted in almost complete loss of binding to Fas. Computational optimization of the hAChE mimetic peptide yielded a variant with slightly improved affinity to Fas, indicating that more rounds of computational optimization will be required to obtain peptide variants with greatly improved affinity for Fas. CD spectra in the absence and presence of Fas point to conformational changes in the peptide upon binding to Fas. Furthermore, binding of the optimized hAChE mimetic peptide to Fas could be inhibited by hAChE, providing evidence for a hAChE‐specific peptide–Fas interaction. Copyright © 2015 European Peptide Society and John Wiley & Sons, Ltd.     

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.     

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.