Long-range modulation of chain motions within the intrinsically disordered transactivation domain of tumor suppressor p53

The tumor suppressor p53 is a hub protein with a multitude of binding partners, many of which target its intrinsically disordered N-terminal domain, p53-TAD. Partners, such as the N-terminal domain of MDM2, induce formation of local structure and leave the remainder of the domain apparently disordered. We investigated segmental chain motions in p53-TAD using fluorescence quenching of an extrinsic label by tryptophan in combination with fluorescence correlation spectroscopy (PET-FCS). We studied the loop closure kinetics of four consecutive segments within p53-TAD and their response to protein binding and phosphorylation. The kinetics was multiexponential, showing that the conformational ensemble of the domain deviates from random coil, in agreement with previous findings from NMR spectroscopy. Phosphorylations or binding of MDM2 changed the pattern of intrachain kinetics. Unexpectedly, we found that upon binding and phosphorylation chain motions were altered not only within the targeted segments but also in remote regions. Long-range interactions can be induced in an intrinsically disordered domain by partner proteins that induce apparently only local structure or by post-translational modification.     

Computational strategy for intrinsically disordered protein ligand design leads to the discovery of p53 transactivation domain I binding compounds that activate the p53 pathway

Intrinsically disordered proteins or intrinsically disordered regions (IDPs) have gained much attention in recent years due to their vital roles in biology and prevalence in various human diseases. Although IDPs are perceived as attractive therapeutic targets, rational drug design targeting IDPs remains challenging because of their conformational heterogeneity. Here, we propose a hierarchical computational strategy for IDP drug virtual screening (IDPDVS) and applied it in the discovery of p53 transactivation domain I (TAD1) binding compounds. IDPDVS starts from conformation sampling of the IDP target, then it combines stepwise conformational clustering with druggability evaluation to identify potential ligand binding pockets, followed by multiple docking screening runs and selection of compounds that can bind multi-conformations. p53 is an important tumor suppressor and restoration of its function provides an opportunity to inhibit cancer cell growth. TAD1 locates at the N-terminus of p53 and plays key roles in regulating p53 function. No compounds that directly bind to TAD1 have been reported due to its highly disordered structure. We successfully used IDPDVS to identify two compounds that bind p53 TAD1 and restore wild-type p53 function in cancer cells. Our study demonstrates that IDPDVS is an efficient strategy for IDP drug discovery and p53 TAD1 can be directly targeted by small molecules.

A hierarchical computational strategy for IDP drug virtual screening (IDPDVS) was proposed and successfully applied to identify compounds that bind p53 TAD1 and restore wild-type p53 function in cancer cells.     

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.     

Determination of protein-ligand binding affinity by NMR: observations from serum albumin model systems

The usefulness of bovine serum albumin (BSA) as a model protein for testing NMR methods for the study of protein-ligand interactions is discussed. Isothermal titration calorimetry established the binding affinity and stoichiometry of the specific binding site for L-tryptophan, D-tryptophan, naproxen, ibuprofen, salicylic acid and warfarin. The binding affinities of the same ligands determined by NMR methods are universally weaker (larger KD). This is because the NMR methods are susceptible to interference from additional non-specific binding. The L-tryptophan-BSA and naproxen-BSA systems were the best behaved model systems.     

Multitude of binding modes attainable by intrinsically disordered proteins: a portrait gallery of disorder-based complexes

Proteins are constantly involved in the multitude of various interactions creating sophisticated networks which define and control all (or almost all) the biological processes taking place in any living organism. Intrinsically disordered proteins or regions play a number of crucial roles in mediating protein interactions. The lack of fixed structure protruding to the high level of intrinsic dynamics and almost unrestricted flexibility at various structure levels, being the major characteristics of intrinsically disordered proteins, provides them with unprecedented advantages over the ordered proteins. The binding modes attainable by disordered proteins are highly diverse, creating a multitude of unusual complexes. Although the majority of studied to date intrinsic disorder-based complexes are ordered or static entities originating due to the global or local disorder-to-order transitions, a new development is the discovery of dynamic complexes in which intrinsically disordered proteins continue to sample an ensemble of rapidly interconverting conformations mostly devoid of structure even in their bound state. The goal of this critical review is to illustrate binding plasticity of intrinsically disordered proteins by representing a portrait gallery of the disorder-based complexes (119 references).     

Continuous thermal collapse of the intrinsically disordered protein tau is driven by its entropic flexible domain

The tau protein belongs to the category of Intrinsically Disordered Proteins (IDP), which in their native state lack a folded structure and fluctuate between many conformations. In its physiological state, tau helps nucleating and stabilizing the microtubules' (MTs) surfaces in the axons of the neurons. Tau is mainly composed by two domains: (i) the binding domain that tightly bounds the MT surfaces and (ii) the projection domain that exerts a long-range entropic repulsive force and thus provides the proper spacing between adjacent MTs. Tau is also involved in the genesis and in the development of the Alzheimer disease when it detaches from MT surfaces and aggregates in paired helical filaments. Unfortunately, the molecular mechanisms behind these phenomena are still unclear. Temperature variation, rarely considered in biological studies, is here used to provide structural information on tau correlated to its role as an entropic spacer between adjacent MTs surfaces. In this paper, by means of small-angle X-ray scattering and molecular dynamics simulation, we demonstrate that tau undergoes a counterintuitive collapse phenomenon with increasing temperature. A detailed analysis of our results, performed by the Ensemble Optimization Method, shows that the thermal collapse is coupled to the occurrence of a transient long-range contact between a region encompassing the end of the proline-rich domain P2 and the first part of the repeats domain, and the region of the N-terminal domain entailing residues 80-150. Interestingly these two regions involved in the tau temperature collapse belong to the flexible projection domain that acts as an entropic bristle and regulates the MTs' architecture. Our results show that temperature is an important parameter that influences the dynamics of the tau projection domain, and hence its entropic behavior.     

Structural and binding study of modified siRNAs with the Argonaute 2 PAZ domain by NMR spectroscopy

By using high-resolution NMR spectroscopy, the structures of a natural short interfering RNA (siRNA) and of several altritol nucleic acid (ANA)-modified siRNAs were determined. The interaction of modified siRNAs with the PAZ domain of the Argonaute 2 protein of Drosophila melanogaster was also studied. The structures show that the modified siRNA duplexes (ANA/RNA) adopt a geometry very similar to the naturally occurring A-type siRNA duplex. All ribose residues, except for the 3' overhang, show 3'-endo conformation. The six-membered altritol sugar in ANA occurs in a chair conformation with the nucleobase in an axial position. In all siRNA duplexes, two overhanging nucleotides at the 3' end enhance the stability of the first neighboring base pair by a stacking interaction. The first overhanging nucleotide has a rather fixed position, whereas the second overhanging nucleotide shows larger flexibility. NMR binding studies of the PAZ domain with ANA-modified siRNAs demonstrate that modifications in the double-stranded region of the antisense strand have some small effects on the binding affinity as compared with the unmodified siRNA. Modification of the 3' overhang with thymidine (dTdT) residues shows a sixfold increase in the binding affinity compared with the unmodified siRNA (relative binding affinity of 17% compared with dTdT-modified overhang), whereas modification of the 3' overhang with ANA largely decreases the binding affinity.     

Quantitative lid dynamics of MDM2 reveals differential ligand binding modes of the p53-binding cleft

The oncoprotein MDM2 regulates the activity and stability of the tumor suppressor p53 through protein-protein interaction involving their N-terminal domains. The N-terminal lid of MDM2 has been implicated in p53 regulation; however, due to its flexible nature, limited data are available concerning its role in ligand binding. The quantitative dynamics study using NMR reported here shows, for the first time, that the lid in apo-MDM2 slowly interconverts between a "closed" state that is associated with the p53-binding cleft and an "open" state that is highly flexible. Our results reveal that apo-MDM2 predominantly populates the closed state, whereas the p53-bound MDM2 exclusively populates the open state. Unlike p53 binding, the small molecule MDM2 antagonist nutlin-3 binds to the cleft essentially without perturbing the closed lid state. The lid dynamics thereby represents a signature for the experimental and virtual screening of therapeutic antagonists that target the p53-MDM2 interaction.     

Bicyclic stapled peptides based on p53 as dual inhibitors for the interactions of p53 with MDM2 and MDMX

In recent years, the strategy of inhibiting the interactions of p53 with murine double minute 2(MDM2)and murine double minute X(MDMX) has been proved to be a promising approach for tumor therapy.However, the poor proteolytical stability and low intracellular delivery efficiency of peptide inhibitors limit their clinical application. Here, we designed and synthesized the bicyclic stapled peptides based on p53 by combining all-hydrocarbon stapling and lactam stapling strategies. We demonstrated th...     

Chemical synthesis of transactivation domain (TAD) of tumor suppressor protein p53 by native chemical ligation of three peptide segments

Chemical composition of tumor suppressor protein p53 is altered via multiple post-translational modifications which modulate its cellular lifetime and interactions with other biomolecules. Here we report total chemical synthesis of a 61-residue form of transactivation domain (TAD) of p53 based on native chemical ligation of three peptide segments. The experiments to characterize its binding to nuclear co-activator binding domain (NCBD) of CREB-binding protein confirmed native-like induced folding upon binding to NCBD. Thus, the synthetic approach described herein can be useful for the preparation of various post-translationally modified analogues of TAD-p53 for further functional biochemical and biophysical studies.     

Kinetics of folding and binding of an intrinsically disordered protein: the inhibitor of yeast aspartic proteinase YPrA

The 68 residue peptide IA 3 is an intrinsically unstructured protein that serves as an endogenous inhibitor of the yeast aspartic proteinase A (YPrA). Although unstructured in free solution, IA 3 forms an N-terminal alpha helix as it binds to YPrA, leading to subnanomolar inhibition of the protease. Equilibrium structural and inhibition studies provide little insight into the mechanism and kinetics of the coupled folding and binding interaction. We have used laser temperature jump spectroscopy to study the kinetics of folding of free IA 3 and of the interaction between IA 3 and YPrA. Inducing folding with trifluoroethanol cosolvent allows us to determine the folding rate (kf approximately 0.3 (micros)(-1)) and the unfolding rate (ku approximately 3 (micros)(-1)) for free IA 3 in water at 25 degrees C. A substantially faster relaxation process is observed in the presence of the proteinase; this process appears to be the kinetic signature of an intermediate binding step in the coupled folding and binding interaction of IA 3 and YPrA.     

Virus-ligand interactions: identification and characterization of ligand binding by NMR spectroscopy

We demonstrate the detection and characterization of ligand binding to viruses via NMR. To illustrate the methodology, the interaction of an antiviral compound with human rhinovirus serotype 2 (HRV2) was investigated. Specific interaction of a capsid-binding inhibitor and native HRV2 was monitored utilizing saturation transfer difference (STD) NMR. STD NMR experiments at atomic resolution allowed those regions of the ligand that are involved in the interaction with the virus to be determined. The approach allows for (i) the fast and robust assessment of binding, (ii) the determination of the ligand binding epitope at atomic resolution without the necessity to crystallize virus-ligand complexes, and (iii) the reuse of the virus in subsequent assays. This methodology enables one to easily identify binding of drugs, peptides, and receptor or antibody fragments to the viral capsid.     

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

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

A study of Src SH2 domain protein-phosphopeptide binding interactions by electrospray ionization mass spectrometry

The noncovalent binding of various peptide ligands to pp60^src (Src) SH2 (Src homology 2) domain protein (12.9 ku) has been used as a model system for development of electrospray ionization mass spectrometry (ESI-MS) as a tool to study noncovalently bound complexes. SH2 motifs in proteins are critical in the signal transduction pathways of the tyrosine kinase growth factor receptors and recognize phosphotyrosine-containing proteins and peptides. ESI-MS with a magnetic sector instrument and array detection has been used to detect the protein-peptide complex with low-picomole sensitivity. The relative abundances of the multiply charged ions for the complex formed between Src SH2 protein and several nonphosphorylated and phosphorylated peptides have been compared. The mass spectrometry data correlate well to the measured binding constants derived from solution-based methods, indicating that the mass spectrometry-based method can be used to assess the affinity of such interactions. Solution-phase equilibrium constants may be determined by measuring the amount of bound and unbound species as a function of concentration for construction of a Scatchard graph. ESI-MS of a solution containing Src SH2 with a mixture of phosphopeptides showed the expected protein-phosphopeptide complex as the dominant species in the mass spectrum, demonstrating the method’s potential for screening mixtures from peptide libraries.     

Identifying residual structure in intrinsically disordered systems: a 2D IR spectroscopic study of the GVGXPGVG peptide

The peptide amide-I vibration of a proline turn encodes information on the turn structure. In this study, FTIR, two-dimensional IR spectroscopy and molecular dynamics simulations were employed to characterize the varying turn conformations that exist in the GVGX(L)PGVG family of disordered peptides. This analysis revealed that changing the size of the side chain at the X amino acid site from Gly to Ala to Val substantially alters the conformation of the peptide. To quantify this effect, proline peak shifts and intensity changes were compared to a structure-based spectroscopic model. These simulated spectra were used to assign the population of type-II β turns, bulged turns, and irregular β turns for each peptide. Of particular interest was the Val variant commonly found in the protein elastin, which contained a 25% population of irregular β turns containing two peptide hydrogen bonds to the proline C═O.     

A graph theoretical approach to the effect of mutation on the flexibility of the DNA binding domain of p53 protein

Tumor suppressor protein p53 becomes inactive due to mutation on its DNA binding core domain leading to misbehavior of this protein and preventing its interaction with DNA. In the present study, changes of the protein conformation by five hot spot mutations of T-p53C were assessed preventing the mutants wild-type (WT) behavior. While studies of this nature were undertaken both experimentally and theoretically, the focus is fundamentally on the effects of the mutation on the dynamics of the protein. Hence, the basic concept underlying this study is the change in flexibility or rigidity of the protein. It was found that stable variant T-p53C (PDB-ID: 1uol) that is structurally and functionally very close to wild-type p53 is the most rigid structure and each single carcinogenic mutation on it makes the structure more flexible. We hypothesize that these changes of the molecule’s flexibility disrupt the network of hydrogen bonds associated with the interaction of WT not only at interaction but in the internal structures of the mutants as well, which prevents them from interacting in the WT fashion loosing the anti-cancer properties of WT.     

Quantitative analysis of three phases of high-density polyethylene with 2D time-domain proton NMR

With two-dimensional time-domain NMR analysis in the proton spin rotating frame, three phases are identifiable in the solid polyethylene. The major proton magnetization fraction is due to the polymer's crystalline region, where the motion is least isotropic and slowest. A magnetization fraction with intermediate relaxation rate is also intermediate in magnitude. This component is proposed to comprise chain loops on the surfaces of crystallites and effectively entangled chain segments. The most mobile fraction, which is most liquid-like with a T₂ of near 1 ms at 120°C, is also the smallest. It is proposed that it is due to polymer chains in the amorphous phase. In the crystalline phase the chain motion is an unexpectedly effective relaxation mechanism at ∼50 KHz. This process, which involves propagation of a twisted region along the crystallite, as has been suggested before, is either not present in the other two phases of the HDPE or is less efficient because of the lower polymer density. The activation energies for the crystalline, intermediate and amorphous phases are estimated to be 34 ± 3, 13 ± 3, and 8 ± 2 Kcal/mol, respectively. The frequencies of the relaxation process at 120°C are 43 KHz, 350 KHz, and 1.5 MHz, respectively. © 1997 John Wiley & Sons, Inc. J Polym Sci B: Polym Phys 35 : 2551–2558, 1997     

Predicting the coordination geometry for Mg2+ in the p53 DNA-binding domain: insights from computational studies

Zn(2+) in the tumor-suppressor protein p53 DNA-binding domain (DBD) is essential for its structural stability and DNA-binding specificity. Mg(2+) has also been recently reported to bind to the p53DBD and influence its DNA-binding activity. In this contribution, the binding geometry of Mg(2+) in the p53DBD and the mechanism of how Mg(2+) affects its DNA-binding activity were investigated using density functional theory (DFT) calculations and molecular dynamics (MD) simulations. Various possible coordination geometries of Mg(2+) binding to histidines (His), cysteines (Cys), and water molecules were studied at the B3LYP/6-311+g** level of theory. The protonation state of Cys and the environment were taken into account to explore the factors governing the coordination geometry. The free energy of the reaction to form the Mg(2+) complexes was estimated, suggesting that the favorable binding mode changes from a four- to six-coordinated geometry as the number of the protonated Cys increases. Furthermore, MD simulations were employed to explore the binding modes of Mg(2+) in the active site of the p53DBD. The simulation results of the Mg(2+) system and the native Zn(2+) system show that the binding affinity of Mg(2+)to the p53DBD is weaker than that of Zn(2+), in agreement with the DFT calculation results and experiments. In addition, the two metal ions are found to make a significant contribution to maintain a favorable orientation for Arg248 to interact with putative DNA, which is critically important to the sequence-specific DNA-binding activity of the p53DBD. However, the effect of Mg(2+) is less marked. Additionally, analysis of the natural bond orbital (NBO) charge transfer reveals that Mg(2+) has a higher net positive charge than Zn(2+), leading to a stronger electrostatic attractive interaction between Mg(2+) and putative DNA. This may partly explain the higher sequence-independent DNA-binding affinity of p53DBD-Mg(2+) compared to p53DBD-Zn(2+) observed in experiment.