Systematic mutational analysis of an ubiquitin ligase (MDM2)-binding peptide: computational studies

To understand the importance of amino acids that comprise the peptide PMI (p53-MDM2/MDMX inhibitor), a p53-mimicking peptide with high affinity for the ubiquitin ligase MDM2, computational alanine scanning has been carried out using various protocols. This approach is very useful for identifying regions of a peptide that can be mutated to yield peptides that bind to their targets with higher affinities. Computational alanine scanning is a very useful technique that involves mutating each amino acid of the peptide in its complex with its target (MDM2 in the current study) to alanine, running short simulations on the mutated complex and computing the difference in interaction energies between the mutant peptides and the target protein (MDM2 in the current study) relative to the interaction energy of the original (wild-type) peptide and the target protein (MDM2 in the current study). We find that running multiple short simulations yield values of computed binding affinities (enthalpies) that are similar to those obtained from a long simulation and are well correlated with the trends in the data available from experiments that used Surface Plasmon Resonance to obtain dissociation constants. The p53-mimicking peptides contain three amino acids (F19, W23 and L26) that are major determinants of the interactions between the peptides and MDM2 and form an essential motif. We find in the current study that the trends amongst the contributions to experimental binding affinities of the hydrophobic residues F19, W23 and L26 are the best reproduced in all the computational protocols examined here. This study suggests that running such short simulations may provide a rapid method to redesign peptides to obtain high-affinity variants against a target protein. We further observe that modelling an extended conformation at the C-terminus of the helical PMI peptides, in accord with the conformation of the p53-peptide complexed to MDM2, reproduces the trends seen amongst the experimental affinities of the peptides that carry the alanine mutations at their C-termini. This suggests that some of the mutant peptides possibly interconvert between helical and extended states and can bind to MDM2 in either conformation. This novel feature, not obvious from the crystallographic data, if factored into modelling protocols, may yield novel high-affinity peptides. Our findings suggest that such protocols may enable rapid investigations of at least certain types of amino acid mutations, notably from large to small amino acids.     

Targeting the conformational transitions of MDM2 and MDMX: insights into dissimilarities and similarities of p53 recognition

MDM2 and MDMX are oncogenic homologue proteins that regulate the activity and stability of p53, a tumor suppressor protein involved in more than 50% of human cancers. While the large body of experiments so far accumulated has validated MDM2 as a therapeutically important target for the development of anticancer drugs, it is only recently that MDMX has also become an attractive target for the treatment of tumor cells expressing wild type p53. The availability of structural information of the N-terminal domain of MDM2 in complex with p53-derived peptides and inhibitors, and the very recent disclosure of the crystal structure of the N-terminal domain of MDMX bound to a p53 peptide, offer an unprecedented opportunity to provide insight into the molecular basis of p53 recognition and the identification of discriminating features affecting the binding of the tumor suppressor protein at MDM2 and MDMX. By using coarse graining simulations, in this study we report the exploration of the conformational transitions featured in the pathway leading from the apo-MDM2 and apo-MDMX states to the p53-bound MDM2 and p53-bound MDMX states, respectively. The results have enabled us to identify a pool of diverse conformational states of the oncogenic proteins that affect the binding of p53 and the presence of conserved and non-conserved interactions along the conformational transition pathway that may be exploited in the design of selective and dual modulators of MDM2 and MDMX activity.     

Use of a retroinverso p53 peptide as an inhibitor of MDM2

An N-terminal helical region of the tumor suppressor p53 binds in a hydrophobic cleft of the oncoprotein MDM2. A retroinverso isomer of the natural N-terminal helical peptide was found to interact with MDM2 using the same hydrophobic residues, Phe, Trp, and Leu. We propose that the retroinverso d-peptide adopts a right-handed helical conformation to achieve functional mimicry of the p53 peptide.     

Quantum mechanical studies of residue-specific hydrophobic interactions in p53-MDM2 binding

Quantum chemistry calculations at the levels of MP2/cc-pVDZ and MP2/cc-PVTZ have been carried out to study residue-specific interactions at the hydrophobic p53-MDM2 binding interface. The result of the calculation, based on structures from nanosecond molecular dynamics simulation, revealed that (19)Phe, (22)Leu, and (23)Trp of p53 have the strongest binding interaction with MDM2 followed by (26)Leu and (27)Pro. The specific residues of MDM2 that have dominant binding interactions with p53 are specifically identified to be (51)Lys, (54)Leu, (62)Met, (67)Tyr, (72)Gln, (94)Lys, (96)His, and (100)Tyr. The p53-MDM2 binding interaction is dominated by van der Waals interaction and to a lesser degree by electrostatic interaction. The MP2 results are in generally good agreement with those from the force field calculation while the DFT/B3LYP calculation failed to give attractive interaction energies for certain residue-residue interactions due to the lack of dispersion energy.     

A computational analysis of the binding model of MDM2 with inhibitors

Abstract  

It is a new and promising strategy for anticancer drug design to block the MDM2-p53 interaction using a non-peptide small-molecule inhibitor. We carry out molecular dynamics simulations to study the binding of a set of six non-peptide small-molecule inhibitors with the MDM2. The relative binding free energies calculated using molecular mechanics Poisson–Boltzmann surface area method produce a good correlation with experimentally determined results. The study shows that the van der Waals energies are the largest component of the binding free energy for each complex, which indicates that the affinities of these inhibitors for MDM2 are dominated by shape complementarity. The A-ligands and the B-ligands are the same except for the conformation of 2,2-dimethylbutane group. The quantum mechanics and the binding free energies calculation also show the B-ligands are the more possible conformation of ligands. Detailed binding free energies between inhibitors and individual protein residues are calculated to provide insights into the inhibitor-protein binding model through interpretation of the structural and energetic results from the simulations. The study shows that G1, G2 and G3 group mimic the Phe19, Trp23 and Leu26 residues in p53 and their interactions with MDM2, but the binding model of G4 group differs from the original design strategy to mimic Leu22 residue in p53.     

Tuning the Binding Affinity and Selectivity of Perfluoroaryl-Stapled Peptides by Cysteine-Editing

A growing number of approaches to “staple” α-helical peptides into a bioactive conformation using cysteine cross-linking are emerging. Here, the replacement of l -cysteine with “cysteine analogues” in combinations of different stereochemistry, side chain length and beta-carbon substitution, is explored to examine the influence that the thiol-containing residue(s) has on target protein binding affinity in a well-explored model system, p53–MDM2/MDMX, which is constituted by the interaction of the tumour suppressor protein p53 and proteins MDM2 and MDMX, which regulate p53 activity. In some cases, replacement of one or more l -cysteine residues afforded significant changes in the measured binding affinity and target selectivity of the peptide. Computationally constructed homology models indicate that some modifications, such as incorporating two d -cysteine residues, favourably alter the positions of key functional amino acid side chains, which is likely to cause changes in binding affinity, in agreement with measured surface plasmon resonance data.     

分子动力学模拟Cu2+对α-突触核蛋白(1-17)肽段构象变化的影响

利用分子动力学模拟研究铜离子(Cu2+)对α-突触核蛋白1-17号氨基酸肽段(α-synuclein(1-17))构象变化的影响,采用GROMOS 43A1力场对Cu2+-α-synuclein(1-17)复合体和α-synuclein(1-17)肽段单体分别进行了6组独立的分子动力学模拟,每组模拟时间为500ns,总模拟时间为3μs.研究结果表明:Cu2+与α-synuclein(1-17)肽段结合使其更易向β折叠片结构折叠,促进了其二级结构的形成,增强了构象的稳定性;Cu2+增大了α-synuclein肽段疏水残基的溶剂可及表面积,增强了其疏水残基的暴露程度.自由能分析指出,Cu2+-α-synuclein(1-17)复合体的自由能比α-synuclein(1-17)肽段低,构象稳定,采样空间紧密,其自由能极小构象为β折叠片结构.构象聚类分析进一步表明,Cu2+使得α-synuclein(1-17)肽段构象趋于稳定.总之,Cu2+诱导固有无序蛋白α-synuclein(1-17)肽段由无序向有序转变,降低了构象的自由能,同时Cu2+增强了α-synuclein(1-17)肽段的疏水性,使得α-synuclein肽段因疏水作用更倾向于形成β折叠片结构,加速其疏水性聚集.     

Molecular Level Characterization of the Structure and Interactions in Peptide‐Functionalized Metal–Organic Frameworks

We use density functional theory, newly parameterized molecular dynamics simulations, and last generation ¹⁵N dynamic nuclear polarization surface enhanced solid‐state NMR spectroscopy (DNP SENS) to understand graft–host interactions and effects imposed by the metal–organic framework (MOF) host on peptide conformations in a peptide‐functionalized MOF. Focusing on two grafts typified by MIL‐68‐proline ( ‐Pro ) and MIL‐68‐glycine‐proline ( ‐Gly‐Pro ), we identified the most likely peptide conformations adopted in the functionalized hybrid frameworks. We found that hydrogen bond interactions between the graft and the surface hydroxyl groups of the MOF are essential in determining the peptides conformation(s). DNP SENS methodology shows unprecedented signal enhancements when applied to these peptide‐functionalized MOFs. The calculated chemical shifts of selected MIL‐68‐NH‐ Pro and MIL‐68‐NH‐ Gly‐Pro conformations are in a good agreement with the experimentally obtained ¹⁵N NMR signals. The study shows that the conformations of peptides when grafted in a MOF host are unlikely to be freely distributed, and conformational selection is directed by strong host–guest interactions.     

Azetidine-derived amino acids versus proline derivatives. alternative trends in reverse turn induction

The influence of 2-alkyl-2-carboxyazetidines (Aze) on the 3D structure of model tetrapeptides R2CO-2-R1Aze-l-Ala-NHMe has been analyzed by molecular modeling, 1H NMR, and FT-IR studies. The conformational constraints introduced by the four-membered ring resulted in an effective way to stabilize gamma-turn-like conformations in these short peptides. The conformational preferences of these Aze-containing peptides have been compared to those of the corresponding peptide analogues containing Pro or alpha-MePro in the place of 2-alkyl-Aze residue. In the model studied, both Pro and Aze derivatives are able to induce reverse turns, but the nature of the turn is different as a function of the ring size. While the five-membered ring of Pro tends to induce beta-turns, as previously suggested by different authors, the four-membered ring of Aze residues forces the peptide to preferentially adopt gamma-turn conformations. In both cases, the presence of an alkyl group at the alpha-position of Pro or the azetidine-2-carboxylate ring enhances significantly the turn-inducing ability. These results might open the opportunity of using 2-alkyl-Aze residues as versatile tools in defining the role of gamma-turn structures within the bioactive conformation of selected peptides, and represent an alternative to Pro derivatives as turn inducers.     

Interchain doubly-bridged α-helical peptides for the development of protein binders

This work reported the design and synthesis of interchain doubly-bridged α-helical peptides, involving mutual stabilization of two α -helical peptides crosslinked by two interchain bisthioether crosslinkers.     

Multiple peptide conformations give rise to similar binding affinities: molecular simulations of p53-MDM2

Molecular dynamics simulations, guided by experimental information (Zondlo et al. Biochemistry 2006, 45, 11945-11957) have been used successfully to reproduce experimental trends in binding affinities of variant p53 peptides with MDM2. Simulations reveal how the conformations of the peptides and the receptor modulate each other to optimize interactions. The conformations of the uncomplexed peptides are governed by a combination of helix and intrinsic disorder (in agreement with experiments), while in the complexed state two very different conformations can coexist. This yields very similar binding affinities, driven by either enthalpy or entropy.     

Influence of an Unnatural Amino Acid Side Chain on the Conformational Dynamics of Peptides

In this work, a non‐natural amino acid, H‐propargylglycine‐OH (Pra), is chosen to examine the side‐chain effect on the backbone conformation of small peptides. The conformations of two synthesized Pra‐containing tripeptides, Ac‐Pra‐Pra‐NH₂ (PPTP) and Ac‐Pra‐Ala‐NH₂ (PATP), are examined by infrared (IR) spectroscopy in combination with molecular dynamics (MD) simulations and quantum chemical computations. By analyzing the joint distributions of backbone torsional angles, several significant conformations can be identified for the two tripeptides solvated in D₂O. At room temperature, 44 % of PPTP exists in the α‐α conformation and 33 % of PATP exists in the α‐polyproline‐II conformation. Larger structural inhomogeneity is seen in both cases by MD simulations at elevated temperatures. Thus even a small side chain, such as the propargyl group can significantly alter the peptide backbone conformations. The results suggest that there is no overwhelming conformational propensity of the Pra residue in short peptides. IR spectra simulated in the amide‐I region using two different methods, reasonably reproduce the experimental IR spectra and their temperature dependence.     

An extended planar C5 conformation and a 310-helical structure of peptide foldamer composed of diverse alpha-ethylated alpha,alpha-disubstituted alpha-amino acids

Optically active peptide foldamers Tfa-[(S)-(alphaEt)Leu]-[(S)-(alphaEt)Nva]-Deg-[(S)-(alphaEt)Nle]-OEt (10) and Tfa-[(S)-(alphaEt)Val]-[(S)-(alphaEt)Leu]-[(S)-(alphaEt)Nva]-Deg-[(S)-(alphaEt)Nle]-OEt (11) composed of diverse alpha-ethylated alpha,alpha-disubstituted alpha-amino acids were synthesized. The dominant conformation of these peptides in solution was an unusual, fully extended planar conformation, and that in the crystal state was both right-handed (P) and left-handed (M) 3(10)-helical structures in 10 and a P 3(10)-helical structure in 11, respectively. The preferred planar C(5) conformation of the peptides prepared from chiral alpha-ethylated alpha,alpha-disubstituted alpha-amino acids was drastically different from the 3(10)-helical structure of the peptides prepared from chiral alpha-methylated alpha,alpha-disubstituted alpha-amino acids.     

Interrogation of MDM2 phosphorylation in p53 activation using native chemical ligation: the functional role of Ser17 phosphorylation in MDM2 reexamined

The E3 ubiquitin ligase MDM2 functions as a crucial negative regulator of the p53 tumor suppressor protein by antagonizing p53 transactivation activity and targeting p53 for degradation. Cellular stress activates p53 by alleviating MDM2-mediated functional inhibition, even though the molecular mechanisms of stress-induced p53 activation still remain poorly understood. Two opposing models have been proposed to describe the functional and structural role in p53 activation of Ser17 phosphorylation in the N-terminal "lid" (residues 1-24) of MDM2. Using the native chemical ligation technique, we synthesized the p53-binding domain (1-109)MDM2 and its Ser17-phosphorylated analogue (1-109)MDM2 pS17 as well as (1-109)MDM2 S17D and (25-109)MDM2, and comparatively characterized their interactions with a panel of p53-derived peptide ligands using surface plasmon resonance, fluorescence polarization, and NMR and CD spectroscopic techniques. We found that the lid is partially structured in apo-MDM2 and occludes p53 peptide binding in a ligand size-dependent manner. Binding of (1-109)MDM2 by the (15-29)p53 peptide fully displaces the lid and renders it completely disordered in the peptide-protein complex. Importantly, neither Ser17 phosphorylation nor the phospho-mimetic mutation S17D has any functional impact on p53 peptide binding to MDM2. Although Ser17 phosphorylation or its mutation to Asp contributes marginally to the stability of the lid conformation in apo-MDM2, neither modification stabilizes apo-MDM2 globally or the displaced lid locally. Our findings demonstrate that Ser17 phosphorylation is functionally neutral with respect to p53 binding, suggesting that MDM2 phosphorylation at a single site is unlikely to play a dominant role in stress-induced p53 activation.     

Interpreting NMR data for beta-peptides using molecular dynamics simulations

NMR is one of the most used techniques to resolve structure of proteins and peptides in solution. However, inconsistencies may occur due to the fact that a polypeptide may adopt more than one conformation. Since the NOE distance bounds and (3)J-values used in such structure determination represent a nonlinear average over the total ensemble of conformers, imposition of NOE or (3)J-value restraints to obtain one unique conformation is not an appropriate procedure in such cases. Here, we show that unrestrained MD simulation of a solute in solution using a high-quality force field yields a conformational ensemble that is largely compatible with the experimental NMR data on the solute. Four 100 ns MD simulations of two forms of a nine-residue beta-peptide in methanol at two temperatures produced conformational ensembles that were used to interpret the NMR data on this molecule and resolve inconsistencies between the experimental NOEs. The protected and unprotected forms of the beta-peptide adopt predominantly a 12/10-helix in agreement with the qualitative interpretation of the NMR data. However, a particular NOE was not compatible with this helix indicating the presence of other conformations. The simulations showed that 3(14)()-helical structures were present in the ensemble of the unprotected form and that their presence correlates with the fulfillment of the particular NOE. Additionally, all inter-hydrogen distances were calculated to compare NOEs predicted by the simulations to the ones observed experimentally. The MD conformational ensembles allowed for a detailed and consistent interpretation of the experimental data and showed the small but specific conformational differences between the protected and unprotected forms of the peptide.     

Conformational profile of a proline-arginine hybrid

The intrinsic conformational preferences of a new nonproteinogenic amino acid have been explored by computational methods. This tailored molecule, named ((β)Pro)Arg, is conceived as a replacement for arginine in bioactive peptides when the stabilization of folded turn-like conformations is required. The new residue features a proline skeleton that bears the guanidilated side chain of arginine at the C(β) position of the five-membered pyrrolidine ring, in either a cis or a trans orientation with respect to the carboxylic acid. The conformational profiles of the N-acetyl-N'-methylamide derivatives of the cis and trans isomers of ((β)Pro)Arg have been examined in the gas phase and in solution by B3LYP/6-31+G(d,p) calculations and molecular dynamics simulations. The main conformational features of both isomers represent a balance between geometric restrictions imposed by the five-membered pyrrolidine ring and the ability of the guanidilated side chain to interact with the backbone through hydrogen bonds. Thus, both cis- and trans-((β)Pro)Arg exhibit a preference for the α(L) conformation as a consequence of the interactions established between the guanidinium moiety and the main-chain amide groups.     

Establishing effective simulation protocols for beta- and alpha/beta-peptides. II. Molecular mechanical (MM) model for a cyclic beta-residue

All-atom molecular mechanical (MM) force field parameters are developed for a cyclic beta-amino acid, amino-cyclo-pentane-carboxylic acid (ACPC), using a multi-objective evolutionary algorithm. The MM model is benchmarked using several short, ACPC-containing alpha/beta-peptides in water and methanol with SCC-DFTB (self consistent charge-density functional tight binding)/MM simulations as the reference. Satisfactory agreements are found between the MM and SCC-DFTB/MM results regarding the distribution of key dihedral angles for the tetra-alpha/beta-peptide in water. For the octa-alpha/beta-peptide in methanol, the MM and SCC-DFTB/MM simulations predict the 11- and 14/15-helical form as the more stable conformation, respectively; however, the two helical forms are very close in energy (2-4 kcal/mol) at both theoretical levels, which is also the conclusion from recent NMR experiments. As the first application, the MM model is applied to an alpha/beta-pentadeca-peptide in water with both explicit and implicit solvent models. The stability of the peptide is sensitive to the starting configuration in the explicit solvent simulations due to their limited length ( approximately 10-40 ns). Multiple ( approximately 20 x 20 ns) implicit solvent simulations consistently show that the 14/15-helix is the predominant conformation of this peptide, although substantially different conformations are also accessible. The calculated nuclear Overhauser effect (NOE) values averaged over different trajectories are consistent with experimental data, which emphasizes the importance of considering conformational heterogeneity in such comparisons for highly dynamical peptides.     

Complex ion effects on polypeptide conformational stability: chloride and sulfate salts of guanidinium and tetrapropylammonium

The effects of chloride and sulfate salts of tetrapropylammonium (TPA(+)) and guanidinium (Gdm(+)) on the conformational stabilities of tryptophan zipper (trpzip) and α-helical (alahel) peptides were measured by circular dichroism spectroscopy. Like Gdm(+), TPA(+) interacts with the planar tryptophan indole group, perturbing the conformational stability of trpzip peptides. TPA(+) effects are largely unaffected by sulfate, indicating an absence of the heteroion pairing that is observed in concentrated Gdm(2)SO(4) solutions. TPA(+) stabilizes helical conformations in alahel peptides, indicating exclusion from the peptide bond. The observations are broadly consistent with predictions of molecular dynamics simulations [Mason, P. E.; et al. J. Phys. Chem. B2009, 113, 3227-3234], indicating that the effects of complex ions on proteins are increasingly predictable in terms of ion hydration, complementary interactions with specific protein groups, and ion-pairing contributions.     

Structure-based design of potent non-peptide MDM2 inhibitors

A successful structure-based design of a class of non-peptide small-molecule MDM2 inhibitors targeting the p53-MDM2 protein-protein interaction is reported. The most potent compound 1d binds to MDM2 protein with a Ki value of 86 nM and is 18 times more potent than a natural p53 peptide (residues 16-27). Compound 1d is potent in inhibition of cell growth in LNCaP prostate cancer cells with wild-type p53 and shows only a weak activity in PC-3 prostate cancer cells with a deleted p53. Importantly, 1d has a minimal toxicity to normal prostate epithelial cells. Our studies provide a convincing example that structure-based strategy can be employed to design highly potent, non-peptide, cell-permeable, small-molecule inhibitors to target protein-protein interaction, which remains a very challenging area in chemical biology and drug design.     

Design and synthesis of functionalized trisaccharides as p53-peptide mimics

Oligosaccharides represent potentially useful scaffolds for the development of peptidomimetics. We report here the design and synthesis of functionalized trisaccharides modeled after an α-helical 15-mer peptide region of p53 which binds to its cellular regulator MDM2. The trisaccharide scaffold was obtained efficiently by applying the sulfoxide glycosylation reaction as a key methodology.