Next Article in Journal / Special Issue
Coupled Folding and Specific Binding: Fishing for Amphiphilicity
Previous Article in Journal
Comparison of Free Energy Surfaces Calculations from Ab Initio Molecular Dynamic Simulations at the Example of Two Transition Metal Catalyzed Reactions
Previous Article in Special Issue
Accounting for Large Amplitude Protein Deformation during in Silico Macromolecular Docking
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 12
Issue 2
10.3390/ijms12021410
ijms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Anchoring Intrinsically Disordered Proteins to Multiple Targets: Lessons from N-Terminus of the p53 Protein

by
Yongqi Huang
1,2,3 and
Zhirong Liu
1,2,3,*
1
State Key Laboratory for Structural Chemistry of Unstable and Stable Species, College of Chemistry and Molecular Engineering, Peking University, Beijing 100871, China
2
Center for Theoretical Biology, Peking University, Beijing 100871, China
3
Beijing National Laboratory for Molecular Sciences, Peking University, Beijing 100871, China
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2011, 12(2), 1410-1430; https://doi.org/10.3390/ijms12021410
Submission received: 25 January 2011 / Revised: 10 February 2011 / Accepted: 16 February 2011 / Published: 23 February 2011
(This article belongs to the Special Issue Advances in Molecular Recognition)
Browse Figures
Versions Notes

Abstract

:
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).

Graphical Abstract

1. Introduction

Proteins are the machines of living systems, and their interaction with other molecules is a central step to perform functions. Considerable efforts have been devoted to understanding the principles governing protein–protein interactions, including interface contacts, morphology, residue conservation, and secondary structures [1–5]. In general, the complex interface is not flat, and some residues from one protein deeply insert into the binding groove of the partner, resulting in the greatest changes in the solvent accessible surface area (SASA) among all the interface residues. Such residues are called anchor residues [6]. Anchor residues have been extensively studied, and their critical roles in specific molecular recognition processes have been widely addressed [6–12]. The most remarkable feature of anchor residues is their conformational preference in the unbound state. All-atom molecular dynamics (MD) simulations have shown that anchor residues in the unbound state are in conformations similar to those observed in the bound complexes [6,7]. The conformational preference of the anchor residues has been successfully applied to improve docking efficiency [9,10]. Recently, Csermely et al. [13] figured out an extended view of binding which embraces a repertoire of conformational selection and structural adjustment process, where they highlighted the important role of anchor residues in the binding process. In their mechanism, conformational selection of the anchor residues in the transient encounter process is critical in the stabilization of the encounter complexes due to their large surface area. The encounter complexes then undergo further induced-fit to complete the binding event [14].
Up to now, studies on anchor residues have been restricted to ordered proteins, i.e., proteins that can be described by defined three-dimensional structures. However, not all proteins form unique structures in the free state. There exists another special family of proteins—the intrinsically disordered proteins (IDPs)—which are flexible in the free state and should be described by an ensemble of conformations [15–30]. The sequence composition of IDPs is very different from that of ordered proteins, and therefore IDPs can be reliably predicted through bioinformatics studies. More than 60 computational tools have been developed for disorder prediction and many of them have been reviewed in reference [31]. IDPs are enriched in cellular functions, such as signaling transduction and transcription regulation [26,32,33], and conformational flexibility is extremely important for IDPs to interact with their targets [34–38]. In experiments, NMR techniques, small-angle X-ray scattering, and different spectroscopic and hydrodynamic methods have been widely used to elucidate the structural features of IDPs [39–43], including the molecular sizes, secondary structural elements, coupled folding-binding processes, and aggregation propensities. In particular, combined with experimentally determined restraints, computer simulations have provided the ensemble-level pictures of IDPs [40].
Experimental and simulation studies have provided evidence that specific functional regions of IDPs are spatially exposed and may be the primary contact sites in the binding processes [44]. Through analysis of complex structures and disorder predictions, short binding regions within long disordered sequences were identified and termed as molecular recognition features (MoRFs) [45–48]. MoRFs differ from other disordered regions due to their significant secondary structure propensities, and may possess preformed structures similar to those in the complex state [49–52]. To form the complexes, IDPs use much of their surface to form the interface [47,53,54]; some residues of IDPs insert deeply into the binding partners [47,55]. Consequently, these highly buried residues can also be defined as anchor residues as those defined in ordered proteins. Although the overall structures of MoRFs have been extensively characterized [52,56,57], an atomic-level examination of the MoRFs in terms of anchor residues is still missing. In addition, IDPs are implicated in multiple interactions as their structural plasticity allows them to efficiently interact with different targets. So a residue may be an anchor residue when binding to one target, whereas in a different protein complex this same residue may not be an anchor. Therefore, it would be of significant importance to investigate the behaviors of anchor residues when IDPs bind to different targets.
In this study, we performed all-atom MD simulations on a helical region of the N-terminal transactivation domain (TAD) of the p53 protein (Glu17–Asn29) (abbreviated as p53N hereafter). p53 is a transcription factor and is critical in preventing cancer development. The p53N region is the binding site of multiple targets, such as MDM2, MDMX, CBP/p300, Taz2, and Bcl-XL [58]. Structure analysis has shown that the p53N region is highly disordered with a transient helix structure formed within residues Thr18–Leu26 [59–61]. When binding to its targets, e.g., MDM2, the p53N helix is prolonged and stabilized [62,63]. Therefore, p53N is a MoRF with a preformed helical structure that binds to MDM2. Using p53N as an example of IDPs, we analyzed the side chain conformations of anchor residues in the unbound state, and compared their conformations with those in the p53N-MDM2 and p53N-Taz2 complexes to address the role of anchor residues in the molecular recognition processes.

2. Results and Discussion

2.1. Anchor Residues in the p53N Complexes

p53N is versatile and interacts with multiple targets. The binding profile of p53N in the p53N-MDM2 complex was different from that in the p53N-Taz2 complex (Figure 1). In the p53N-MDM2 complex, Phe19, Trp23, and Leu26 were highly buried in the binding groove (Figure 1a,d). In contrast, in the p53N-Taz2 complex, only Leu22 was highly buried (Figure 1b,e). Structural superposition showed that p53N used different surfaces of a helix structure to bind to these two targets (Figure 1c), i.e., the p53N helix rotated about 90° in these complexes. So the buried residues in the p53N-MDM2 complex (Phe19, Trp23, and Leu26) were exposed to solvent in the p53N-Taz2 complex (Figure 1a,b). In globular protein complexes, anchor residues are identified based on the structure of the complexes. They correspond to solvent exposed residues that are fully buried upon binding to a target, yielding the largest change in SASA [6]. However, for complexes formed by IDPs, the determination of SASA for the unbound state may not be so straightforward, since an ensemble of conformations instead of a unique conformation are needed to describe the unbound state of an IDP. We calculated the SASA for the unbound p53N by two different approaches. The simplest one was removing the Taz2 and MDM2 proteins from the p53N-Taz2 and p53N-MDM2 complexes and then calculating the SASA of the remaining p53N (Figure 1d,e). We further calculated the SASA using the simulated conformations of p53N (Figure S1). In general, different approaches gave similar results. According to the initial definition of anchor residues [6], we suggested that anchor residues of IDPs should expose to solvent in the unbound state and become fully buried (SASA ≤ 15 Å2 was used here) after binding. Therefore, we identified Phe19, Trp23, and Leu26 as the anchor residues in the complex with MDM2, and Leu22 as the anchor residue in the complex with Taz2. It was noted that the side chain conformations of Phe19, Leu22, Trp23, and Leu26 in the MDM2 complex were remarkably different from those in the Taz2 complex (Figure 1c). Consistent with conformation ensemble of the p53 TAD [64], Phe19, Leu22, Trp23, and Leu26 were not buried in our simulated unbound states.

2.2. Transient Stable Helix of the p53N

Although the p53N region is rather flexible in solution, experiments indicate that residues Thr18–Leu26 form a transiently stable helix which will be further stabilized in the complex state [59]. To study the conformational preference of unbound p53N, we conducted multiple simulations on the p53N with initial conformations adopted from the p53N-MDM2 and p53N-Taz2 complexes under identical conditions (see the Method Section for details). The helical structures were only transiently stable and unfolded in the simulations under 300 K (Figures 2 and S2). We defined the unfolding time of each trajectory through secondary structure analysis, RMSD relative to the initial helical structure, and inspection of the structures (with detailed results in Table S1). The average unfolding time was 24.5 ± 11.8 ns for the helix from p53N-MDM2 complex and 5.2 ± 3.1 ns for the helix from p53N-Taz2 complex. Clearly, the helix from the p53N-MDM2 complex was much more stable than that from the p53N-Taz2 complex. This may be due to the deformation of the second turn (in the C-terminus) of the p53N helix from the p53N-Taz2 complex, because the unfolding of the p53N helix from either p53N-MDM2 or p53N-Taz2 usually started from the deformation of the second turn during the simulations. We also compared the properties of our ensemble with those from the Daughdrill’s group [64]. Consistent with their results, the distributions of amide nitrogen distances between residues i:i + 5 were bimodal, and Phe19–Trp23 showed greater probabilities of collapsed structures (data not shown).

2.3. Analysis of the Side Chain Conformations of the Anchor Residues

Since the helix of p53N is transiently stable, it is expected that the free energy of the helix state in the free form is higher than the disordered state (Figure S3). The free energy barrier of unfolding is lower than the free energy barrier of folding, resulting in a greater unfolding rate than a folding rate. To obtain conformations of the helix state, we separated the p53N helix from the complex state and carried out simulations. The system quickly relaxed to the free energy basin of the helix state and probably got equilibrium in the basin before it unfolded to the disordered state (Figure 2). To obtain conformations of the disordered state, randomly selected disordered structures were used as initial states from simulations. Then the system got equilibrium in the free energy basin of the disordered state. No disorder-to-helix transition was observed in simulations. It was noted that this strategy did not produce an equilibrium population between the helix state and the disordered state, although it gave the distribution of the side chain conformations within each state (Figure 3).
We analyzed the side chain conformations for the helix state and the disordered state. In our simulations, the helix of p53N from the p53N-Taz2 complex unfolded very quickly (∼5 ns). Within such a short period, conformational sampling of the side chains in the helix state was insufficient. So we analyzed the side chain conformations based on the simulations of p53N from the p53N-MDM2 complex which has a much longer unfolding time. The conformational sampling was found to be more efficient in this case. For example, more than 130 transitions were observed between the two main χ2 conformations of Phe19 and the distribution of the χ2 of Phe19 appeared to be perfectly symmetric (Figure 3a) which is required due to the symmetric nature of Phe.
In this study, we considered two binding targets, so for the four identified anchor residues, each has two bound-like conformations: one is an anchor-type (conformations of Phe19, Trp23, and Leu26 in the p53-MDM2 complex, and that of Leu22 in the p53-Taz2 complex); the other is a non-anchor conformation (conformations of Phe19, Trp23, and Leu26 in the p53-Taz2 complex, and that of Leu22 in the p53-MDM2 complex). Remarkably, the analysis showed that the anchor residues dominantly sampled the anchor-type bound-like conformations rather than the non-anchor bound-like conformations regardless of whether the helix or the disordered states were examined (Figure 3 and Table 1). The anchor residue with the highest population of the non-anchor bound-like conformation was Phe19, but its value was only 12.2% and 6.0% in the helical and disordered states, much smaller than the corresponding value of the anchor-type conformation (59.5% and 19.2%). For the other three anchor residues, the population of non-anchor bound-like conformation was negligible. It was also noted that the formation of a (transient) helical structure enhanced the predominance of the anchor-type bound-like conformations. In the p53N-MDM2 complex, χ1 and χ2 of Phe19 are 177° and 71°. The population of this (anchor-type) rotamer increased by a factor of ∼2 (i.e., 19.2% vs. 59.5%, Table 1), when the p53N transformed from a disordered state to a helix state. Trp23 showed a similar trend. Although the extent of the increase was weaker for Leu22 and Leu26, the formation of a helix still increased the population of the anchor-type conformations.
Although the initial structure for simulations was isolated from the p53N-MDM2 complex, Leu22 sampled conformations similar to that in the p53N-Taz2 complex but not similar to that in the p53N-MDM2 complex (Figure 3b and Table 1). Consequently, the discrimination between the populations of the anchor-type and non-anchor bound-like conformations was not caused by a bias of the initial states. It could not be solely explained in terms of the side-chain rotamer preferences either: rotamer library data [65] indicated that the rotamer of Phe19 preferred the non-anchor bound-like conformation (47.08%) rather than the anchor-type (31.71%); however, the trend was reversed in the helix state (12.2% vs. 59.5%). The preference of Trp23 on the anchor-type conformation was low in the rotamer library (16.21%), which was greatly enhanced in the helix state (62.4%). These observations suggest that the preference of the anchor-type conformations is intrinsic to the transient helical structure of the unbound p53N.
Figure 3 and Table 1 show that the extent of the population shift during the helix-disorder transition was similar for Phe19 and Trp23, suggesting synchronous dynamics between these residues. χ1 of Phe19 sampled three regions (labeled as A, B, and C in Figure 4a. Due to the symmetric nature of Phe, χ2 was not distinguished.). In the χ1–χ2 space, Trp23 sampled six regions (labeled from 1 to 6 in Figure 4b). So there were 18 possible combinations of the conformations for Phe19 and Trp23, where the A1 group was the anchor-type bound-like conformation corresponding to that in the p53N-MDM2 complex. The population analysis showed a remarkable feature that the population of the A1 group (51%) was significantly higher than all other possible combinations (Figure 4c). Furthermore, the population of the A1 group was very close to the individual population of the anchor-type bound-like conformations of Phe19 and Trp23 (59.5% and 62.4%, respectively), showing a strong correlation between these two residues. In the helical conformation of p53N, Phe19 and Trp23 drove (or confined) each other to the anchor-type bound-like conformations.
We also analyzed the side chain conformations of non-anchor residues, e.g., Glu17, Lys24, Leu25, and Glu28 (Figure 5). For Glu17 and Glu28, both bound conformations (in the p53N-MDM2 and p53N-Taz2 complexes) were rarely sampled in the helix and disordered states (data not shown) during the simulations. Therefore, the conformations of these residues in the complex states were induced by interactions with the targets. For Lys24 and Leu25, bound-like conformations were frequently sampled in the simulations and the formation of the helix increased the populations.

2.4. Transient Formation of Helical Structures Promotes the Binding Process

Because atomic information of the MoRFs is missing and how/why MoRFs and anchor residues initiate the binding process is unclear. To identify the role of anchor residues and the performed structure of MoRFs in the binding processes, we performed binding simulations of p53N to MDM2 and Taz2. Firstly, we randomly selected five conformations from the A1 group (Figure 4) as initial conformations of the p53N. Then, we placed the p53N close to the binding groove of MDM2 with the correct orientation (Figure 6a). This was to mimic the encounter of the two proteins. In these encounter states, the anchor residues did not insert into the binding groove. Based on these in silico encounter states, we conducted MD simulations to track the evolution process. As expected, after the local conformational rearrangements on the binding groove of MDM2, Phe19, Trp23, and Leu26 inserted into the binding groove within ∼1 ns (Figure 6b). This result is consistent with the observations of binding the native p53N to MDM2 [66], because conformations in the A1 group were native-like. Our results confirm the validity of results which are based on the prerequisite that p53N is in a bound conformation during the encounter process [66,67]. Furthermore, we also performed simulations with the p53N in other conformations, e.g., conformations from the B3 group and the disordered states. Within 10 ns of simulations, p53N did not evolve towards the bound conformations but formed non-native interactions with MDM2 (Figure 6c,d). The same conformations of p53N above were also used to simulate binding of p53N to Taz2. For p53N in the helix state (i.e., conformations from the A1 group and B3 group), Leu22 inserted into the binding groove quickly; however, the correct conformations of Phe19 and Trp23 were not observed within 10 ns (Figure 6e,f). This was due to the steric constraints at the binding interface. Therefore, induced formation of the correct conformations of Phe19 and Trp23 may take longer. For p53N in the disordered state, only non-native interactions between p53N and Taz2 were observed.
Although experimental and computational studies have shown that MDM2 undergoes structural rearrangement, in particular in the binding groove, upon the p53N binding [66,68], it does not contradict the concept of a preformed bound-like conformation of the unbound p53N. On the contrary, the bound-like conformation of p53N in the unbound state actually promotes the groove opening of the unbound MDM2 [66]. Furthermore, the important roles of Phe19, Trp23, and Leu26 in the binding process of p53N with MDM2 and Leu22 in the binding process of p53N with Taz2 have been studied thermodynamically [69,70]; however, their roles in the kinetic process is not clear. In this work, through simulations, we found that preformed bound-like conformations of these anchor residues promoted the binding process.

2.5. Discussions: Roles of Anchor Residues in Molecular Recognition

In this study, we tried to extend the concept of anchor residues to IDPs and understand the conformational properties of the anchor residues within a highly flexible context. To this end, we performed atomic MD simulations on an extensively studied system, the p53N region. In the unbound state, p53N is rather flexible; however, once the helix is partially and transiently formed, simulations showed that the anchor residues were restricted to their anchor-type bound-like conformations and were primed for interacting with their targets (Figures 3 and 4, and Table 1).
A comparison of the bound conformations of the anchor residues and non-anchor residues with the rotamer library derived from the Protein Data Bank [65] further supports the concept that anchor residues adopt preformed bound-like conformations (Table 1). In the p53N-MDM2 complex, the conformations of the anchor residues (Phe19, Trp23, and Leu26) represent the major conformations in the rotamer library, whereas the conformations of Glu17, Leu22, and Glu28 exhibit very low values in the rotamer library. Conversely, in the complex p53N-Taz2, the conformation of Leu22 is the major conformation in the library, whereas the conformations of Trp23 and Leu26 exhibit very low populations. The formation of rare rotamers is induced by interactions between p53N and the targets [71].
Significant efforts have been made to identify the structure of MoRFs and understand the molecular recognition processes between IDPs and their binding targets. Experimental and computational studies have shown that preformed structures of MoRFs in the unbound state resemble structures in the bound complexes and therefore facilitate the recognition processes [49,51,52,57,59,72–75]. However, the roles of MoRFs in the molecular recognition processes remain elusive. In this study, we extended the understanding through an investigation of the correlation between the side chain conformations and the overall structure of MoRFs. We suggest that, as in globular proteins, anchor residues also exist in IDPs and are important in the specific molecular recognition processes. Since the side chain conformations are backbone dependent, the preformed structure of MoRFs in the unbound state provides a constraint on the side chain conformations to produce bound-like anchor residues. A recent study showed that the conformational preference of residues in a disordered 20-mer peptide was closely correlated to inhibitory activity [76]. Recently, Kjaergaard et al. determined the core structure of a molten globule by NMR and found that the side chain of some hydrophobic residues had preferred rotamers and made specific interactions [77]. Characterization of the MoRFs of Sendai virus nucleoprotein also showed that the transient formation of a helix optimized the interaction with the negatively charged cleft on the surface of the phosphoprotein PX domain [73]. These results support the concept that particular key residues in IDPs have preferred (function-related) conformations in the unbound state.
Combined with the concepts of MoRFs and anchor residues, a feasible mechanism of the recognition process between a disordered binding region and its target emerges. For example, in the binding process of p53 to MDM2, the binding is initiated by an encounter between the preformed bound-like MoRF of p53 and MDM2, which will produce a transient encounter complex. Insertion of the preformed bound-like anchor residues (Phe19, Trp23, and Leu26), which are located on the MoRF, into the binding groove of MDM2 stabilizes the transient encounter complex (Figure 6a,b). Further folding of the backbone and induced-fit of the side chains take place to finally form the native complex. Therefore, in this mechanism, a binding process between an IDP and its target is initiated by a conformational selection and then proceeds by folding upon binding [51,75]. In the simulated binding process of proline-rich motifs to SH3 domains, electrostatic interactions guide the diffusion to form a nonspecific encounter complex state which is stabilized by subsequent anchoring of an arginine of the peptide into the negatively charged groove of the SH3 domain [78]. This gives some support to our proposed mechanism.
To account for the specific recognition process, various mechanisms have been proposed, including the lock-and-key model [79], the induced-fit mechanism [80], and the conformational selection model [81,82]. Recently, Boehr et al. proposed a general mechanism constituted by a primary conformational selection event followed by an induced-fit of side chains and the backbone to account for the role of dynamics in the biomolecular recognition process [83]. Similarly, Csermely et al. figured out an extended view of binding which embraces a repertoire of selection and adjustment processes [13]. All these mechanisms mainly focus on ordered proteins and the discussions on IDPs are very limited. Our study adds important insights into the molecular recognition mechanisms and extends them to molecular interactions involving IDPs.
IDPs have been proposed to have the ability to bind to multiple targets [38] and p53 is a typical example. The N-terminus of p53 binds MDM2, MDMX, Taz2, and Bcl-XL, while the C-terminus of p53 interacts with S100ββ, Sirtuin, CBP, and Cyclin A2 [58]. SASA analysis shows that p53N uses different anchor residues to bind to different targets (Figure 1). Forming a complex with MDM2, Phe19, Trp23, and Leu26 are the anchor residues, whereas binding with Taz2, Leu22 is the anchor residue. By analyzing the complex structures, Oldfield et al. found that the same residues from the C-terminus of p53 are used to a different extent in binding to different targets and interactions involving the same residue may exclude each other in different complexes [55]. Our simulations showed that, in the helix state of unbound p53N, Phe19, Trp23, and Leu26 significantly sampled conformations similar to those in the complex with MDM2, whereas Leu22 sampled conformations similar to that in the complex with Taz2 (Figures 3 and 4, Table 1); these indicate a new mechanism to account for the binding promiscuity of IDPs. It is possible that all the anchor residues in the MoRFs frequently sample the bound conformations of the corresponding targets and that a particular target selects a particular group of anchor residues in the binding process. This binding mechanism is advantageous in smoothing molecular interactions [84], and reconciles the binding promiscuity and binding kinetics in the binding process and provides a clearer picture for the one-to-many signaling processes where flexibility has been considered as the main source [55,85,86]. A similar mechanism has been proposed to interpret the structural basis of pregnane X receptor binding promiscuity, where the pregnane X receptor has five hot spot regions and, depending on their sizes and shapes, individual PXR ligands extend into two, three, or four hot spot regions [87].
To extend the mechanism obtained from p53N to other IDPs, great efforts are still required. For example, because not all MoRFs adopt helical conformations [46,47], whether the non-helical MoRFs have preformed structures or anchor residues in these non-helical MoRFs also adopt bound-like conformations is unclear. It would also be valuable to experimentally test the side chain conformational preference observed in our simulations.

3. Method Section

3.1. Systems Setup

A series of simulations were performed based on the p53 N-terminal domain (Glu17–Asn29, p53N) (Table S1). Initial conformations for MD simulations were taken from the p53N-MDM2 complex (PDB ID 1YCR) [62] (trajectories MS1–MS10), the p53N-Taz2 complex (PDB ID 2K8F) [88] (trajectories TS1–TS5), or the disordered state of unbound p53N (trajectories DS1–DS5). The disordered states of unbound p53N were randomly selected snapshots ranging from 60 ns to 100 ns in the trajectory MS1. The duration for each trajectory varied between 20 and 100 ns (see Table S1). In the binding simulations, we randomly selected five conformations from the A1 group, B3 group (Figure 4), and the disordered state, respectively, as initial conformations of p53N; then we placed p53N close to the binding groove of MDM2 or Taz2 with the correct orientation (Figure 6). All binding simulations lasted for 10 ns.

3.2. Molecular Dynamics Simulations

The MD simulations were performed using the program GROMACS 4.07 [89,90] and the OPLS-AA/L force field [91]. The water molecules were modeled by the SPC/E representation [92]. Each of the starting conformations was placed in the center of a cubic water box with at least 10 Å from the box edge. Periodic boundary conditions were used. Counter ions (Na+ or Cl−) were added to neutralize the net charges. The long-range electrostatic interactions were treated with the particle mesh Ewald method [93]. The cutoff distances were set to 10 Å for short-range coulomb and van der Waals interactions. The bond lengths were fixed by the LINCS algorism [94], and a time step of 2 fs was used. Coordinates were saved every 5 ps.
Each system was first relaxed by 1000 steps of the steepest-descent energy minimization. After the minimization, the system was equilibrated at 300 K by 100 ps under an NVT ensemble and further equilibrated for 200 ps at constant pressure (1 bar). V-rescale [95] and Parrinello-Rahman [96] were used to couple the system to the simulation temperature and pressure with coupling constants of 0.1 ps and 2.0 ps, respectively. Production simulations were performed at constant temperature (300 K) and pressure (1 bar).

3.3. Analysis

Secondary structure analysis was assigned by the DSSP program [97]. NACCESS ( http://wolf.bms.umist.ac.uk/naccess/) was used to calculate the SASA, using a default solvent probe radius of 1.4 Å. We monitored the side chain dihedral angles (χ1 and χ2) distribution of anchor residues in the helix state and disordered state. χ1 is defined by N-Cα-Cβ-Cγ and χ2 is defined by Cα-Cβ-Cγ-Cδ(1). Protein figures were produced by the VMD program [98].

4. Conclusions

In this study, we identified the anchor residues of the disordered p53N when binding to different targets, i.e., Phe19, Trp23, and Leu26 for the p53N-MDM2 complex, and Leu22 for the p53N-Taz2 complex. From the all-atom simulations of the unbound p53N, we found that, as in ordered proteins, these anchor residues in p53N frequently sampled conformations similar to those in the complex states where anchor residues act as anchors, but seldom sampled those in the alternative complexes. We suggest that the bound-like conformations of anchor residues in the unbound state are an important factor in controlling the specific interaction between IDPs and their targets, in particular, stabilizing the transient encounter complexes. We also propose a new mechanism to account for the binding promiscuity of IDPs.

Supplementary Materials

  • Figure S1.:

  • Figure S2.:

  • Figure S3.:

    • Table
    Table S1.
    Table S1.
    SimulationsInitial Conformation*Simulation Duration (ns)Unfolding Time (ns)
    MS1M10043
    MS2M2011
    MS3M10042
    MS4M10025
    MS5M2017
    MS6M6028
    MS7M2010
    MS8M2015
    MS9M6023
    MS10M6031
    TS1T202
    TS2T204
    TS3T203
    TS4T209
    TS5T208
    DS1D50-
    DS2D50-
    DS3D50-
    DS4D50-
    DS5D50-
    *M denotes the initial p53N structure from p53N-MDM2 complex; T for the initial p53N structure from p53N-Taz2 complex; D for the disordered state.

    Acknowledgments

    This work was supported by the Ministry of Science and Technology of China (No. 2009CB918500) and the National Natural Science Foundation of China (No. 20973016 and No. 10721403).

    References

    1. Larsen, TA; Olson, AJ; Goodsell, DS. Morphology of protein-protein interfaces. Structure 1998, 6, 421–427. [Google Scholar]
    2. Lo Conte, L; Chothia, C; Janin, J. The atomic structure of protein-protein recognition sites. J. Mol. Biol 1999, 285, 2177–2198. [Google Scholar]
    3. Ma, BY; Elkayam, T; Wolfson, H; Nussinov, R. Protein-protein interactions: Structurally conserved residues distinguish between binding sites and exposed protein surfaces. Proc. Natl. Acad. Sci. USA 2003, 100, 5772–5777. [Google Scholar]
    4. Bai, HJ; Ma, WZ; Liu, SY; Lai, LH. Dynamic property is a key determinant for protein-protein interactions. Proteins 2008, 70, 1323–1331. [Google Scholar]
    5. Zhang, QC; Petrey, D; Norel, R; Honig, BH. Protein interface conservation across structure space. Proc. Natl. Acad. Sci. USA 2010, 107, 10896–10901. [Google Scholar]
    6. Rajamani, D; Thiel, S; Vajda, S; Camacho, CJ. Anchor residues in protein-rotein interactions. Proc. Natl. Acad. Sci. USA 2004, 101, 11287–11292. [Google Scholar]
    7. Kimura, SR; Brower, RC; Vajda, S; Camacho, CJ. Dynamical view of the positions of key side chains in protein-protein recognition. Biophys. J 2001, 80, 635–642. [Google Scholar]
    8. Li, X; Keskin, O; Ma, B; Nussinov, R; Liang, J. Protein-protein interactions: Hot spots and structurally conserved residues often locate in complemented pockets that pre-organized in the unbound states: Implications for docking. J. Mol. Biol 2004, 344, 781–795. [Google Scholar]
    9. Camacho, CJ. Modeling side-chains using molecular dynamics improve recognition of binding region in CAPRI targets. Proteins 2005, 60, 245–251. [Google Scholar]
    10. Smith, GR; Sternberg, MJE; Bates, PA. The relationship between the flexibility of proteins and their conformational states on forming protein-protein complexes with an application to protein-protein docking. J. Mol. Biol 2005, 347, 1077–1101. [Google Scholar]
    11. Yogurtcu, ON; Erdemli, SB; Nussinov, R; Turkay, M; Keskin, O. Restricted mobility of conserved residues in protein-protein interfaces in molecular simulations. Biophys. J 2008, 94, 3475–3485. [Google Scholar]
    12. Ben-Shimon, A; Eisenstein, M. Computational mapping of anchoring spots on protein surfaces. J. Mol. Biol 2010, 402, 259–277. [Google Scholar]
    13. Csermely, P; Palotai, R; Nussinov, R. Induced fit, conformational selection and independent dynamic segments: An extended view of binding events. Trends Biochem. Sci 2010, 35, 539–546. [Google Scholar]
    14. Wlodarski, T; Zagrovic, B. Conformational selection and induced fit mechanism underlie specificity in noncovalent interactions with ubiquitin. Proc. Natl. Acad. Sci. USA 2009, 106, 19346–19351. [Google Scholar]
    15. Fink, AL. Compact intermediate states in protein folding. Annu. Rev. Biophys. Biomol. Struct 1995, 24, 495–522. [Google Scholar]
    16. Wright, PE; Dyson, HJ. Intrinsically unstructured proteins: Re-assessing the protein structure-function paradigm. J. Mol. Biol 1999, 293, 321–331. [Google Scholar]
    17. Uversky, VN; Gillespie, JR; Fink, AL. Why are “natively unfolded” proteins unstructured under physiologic conditions? Proteins 2000, 41, 415–427. [Google Scholar]
    18. Dunker, AK; Lawson, JD; Brown, CJ; Williams, RM; Romero, P; Oh, JS; Oldfield, CJ; Campen, AM; Ratliff, CM; Hipps, KW; Ausio, J; Nissen, MS; Reeves, R; Kang, C; Kissinger, CR; Bailey, RW; Griswold, MD; Chiu, W; Garner, EC; Obradovic, Z. Intrinsically disordered protein. J. Mol. Graphics Modell 2001, 19, 26–59. [Google Scholar]
    19. Dunker, AK; Obradovic, Z. The protein trinity–linking function and disorder. Nat. Biotechnol 2001, 19, 805–806. [Google Scholar]
    20. Tompa, P. Intrinsically unstructured proteins. Trends Biochem. Sci 2002, 27, 527–533. [Google Scholar]
    21. Uversky, VN. Natively unfolded proteins: A point where biology waits for physics. Protein Sci 2002, 11, 739–756. [Google Scholar]
    22. Uversky, VN. What does it mean to be natively unfolded? Eur. J. Biochem 2002, 269, 2–12. [Google Scholar]
    23. Dunker, AK; Brown, CJ; Lawson, JD; Iakoucheva, LM; Obradović, Z. Intrinsic disorder and protein function. Biochemistry 2002, 41, 6573–6582. [Google Scholar]
    24. Dunker, AK; Brown, CJ; Obradović, Z. Identification and functions of usefully disordered proteins. Adv. Protein Chem 2002, 62, 25–49. [Google Scholar]
    25. Tompa, P; Csermely, P. The role of structural disorder in the function of RNA and protein chaperones. FASEB J 2004, 18, 1169–1175. [Google Scholar]
    26. Dyson, HJ; Wright, PE. Intrinsically unstructured proteins and their functions. Nat. Rev. Mol. Cell Biol 2005, 6, 197–208. [Google Scholar]
    27. Radivojac, P; Iakoucheva, LM; Oldfield, CJ; Obradovic, Z; Uversky, VN; Dunker, AK. Intrinsic disorder and functional proteomics. Biophys. J 2007, 92, 1439–1456. [Google Scholar]
    28. Huang, YQ; Liu, ZR. Intrinsically disordered proteins: The new sequence-structure-function relations. Acta Phys. Chim. Sin 2010, 26, 2061–2072. [Google Scholar]
    29. Uversky, VN; Dunker, AK. Understanding protein non-folding. Biochim. Biophys. Acta 2010, 1804, 1231–1264. [Google Scholar]
    30. Turoverov, KK; Kuznetsova, IM; Uversky, VN. The protein kingdom extended: Ordered and intrinsically disordered proteins, their folding, supramolecular complex formation, and aggregation. Prog. Biophys. Mol. Biol 2010, 102, 73–84. [Google Scholar]
    31. He, B; Wang, KJ; Liu, YL; Xue, B; Uversky, VN; Dunker, AK. Predicting intrinsic disorder in proteins: An overview. Cell Res 2009, 19, 929–949. [Google Scholar]
    32. Tompa, P; Fuxreiter, M. Fuzzy complexes: Polymorphism and structural disorder in protein-protein interactions. Trends Biochem. Sci 2008, 33, 2–8. [Google Scholar]
    33. Galea, CA; Wang, Y; Sivakolundu, SG; Kriwacki, RW. Regulation of cell division by intrinsically unstructured proteins: Intrinsic flexibility, modularity, and signaling conduits. Biochemistry 2008, 47, 7598–7609. [Google Scholar]
    34. Shoemaker, BA; Portman, JJ; Wolynes, PG. Speeding molecular recognition by using the folding funnel: The fly-casting mechanism. Proc. Natl. Acad. Sci. USA 2000, 97, 8868–8873. [Google Scholar]
    35. Huang, YQ; Liu, ZR. Kinetic advantage of intrinsically disordered proteins in coupled folding-binding process: A critical assessment of the “fly-casting” mechanism. J. Mol. Biol 2009, 393, 1143–1159. [Google Scholar]
    36. Mittag, T; Kay, LE; Forman-Kay, JD. Protein dynamics and conformational disorder in molecular recognition. J. Mol. Recognit 2010, 23, 105–116. [Google Scholar]
    37. Huang, YQ; Liu, ZR. Nonnative interactions in coupled folding and binding processes of intrinsically disordered proteins. PLoS ONE 2010, e15375. [Google Scholar]
    38. Patil, A; Kinoshita, K; Nakamura, H. Hub promiscuity in protein-protein interaction networks. Int. J. Mol. Sci 2010, 11, 1930–1943. [Google Scholar]
    39. Receveur-Bréchot, V; Bourhis, J-M; Uversky, VN; Canard, B; Longhi, S. Assessing protein disorder and induced folding. Proteins 2006, 62, 24–45. [Google Scholar]
    40. Mittag, T; Forman-Kay, JD. Atomic-level characterization of disordered protein ensembles. Curr. Opin. Struct. Biol 2007, 17, 3–14. [Google Scholar]
    41. Chen, HM; Rhoades, E. Fluorescence characterization of denatured proteins. Curr. Opin. Struct. Biol 2008, 18, 516–524. [Google Scholar]
    42. Eliezer, D. Biophysical characterization of intrinsically disordered proteins. Curr. Opin. Struct. Biol 2009, 19, 23–30. [Google Scholar]
    43. Uversky, V; Longhi, S. Instrumental Analysis of Intrinsically Disordered Proteins: Assessing Structure and Conformation; John Wiley & Sons, Inc: Hoboken, NJ, USA, 2010. [Google Scholar]
    44. Csizmók, V; Bokor, M; Bánki, P; Klement, T; Medzihradszky, KF; Friedrich, P; Tompa, KA; Tompa, P. Primary contact sites in intrinsically unstructured proteins: The case of calpastatin and microtubule-associated protein 2. Biochemistry 2005, 44, 3955–3964. [Google Scholar]
    45. Oldfield, CJ; Cheng, Y; Cortese, MS; Brown, CJ; Uversky, VN; Dunker, AK. Comparing and combining predictors of mostly disordered proteins. Biochemistry 2005, 44, 1989–2000. [Google Scholar]
    46. Mohan, A; Oldfield, CJ; Radivojac, P; Vacic, V; Cortese, MS; Dunker, AK; Uversky, VN. Analysis of molecular recognition features (MoRFs). J. Mol. Biol 2006, 362, 1043–1059. [Google Scholar]
    47. Vacic, V; Oldfield, CJ; Mohan, A; Radivojac, P; Cortese, MS; Uversky, VN; Dunker, AK. Characterization of molecular recognition features, MoRFs, and their binding partners. J. Proteome Res 2007, 6, 2351–2366. [Google Scholar]
    48. Cheng, YG; Oldfield, CJ; Meng, JW; Romero, P; Uversky, VN; Dunker, AK. Mining α-helix-forming molecular recognition features with cross species sequence alignments. Biochemistry 2007, 46, 13468–13477. [Google Scholar]
    49. Fuxreiter, M; Simon, I; Friedrich, P; Tompa, P. Preformed structural elements feature in partner recognition by intrinsically unstructured proteins. J. Mol. Biol 2004, 338, 1015–1026. [Google Scholar]
    50. Dancheck, B; Naim, AC; Peti, W. Detailed structural characterization of unbound protein phosphatase 1 inhibitors. Biochemistry 2008, 47, 12346–12356. [Google Scholar]
    51. Yoon, M-K; Venkatachalam, V; Huang, A; Choi, B-S; Stultz, CM; Chou, JJ. Residual structure within the disordered C-terminal segment of p21(Waf1/Cip1/Sdi1) and its implications for molecular recognition. Protein Sci 2009, 18, 337–347. [Google Scholar]
    52. Marsh, JA; Dancheck, B; Ragusa, MJ; Allaire, M; Forman-Kay, JD; Peti, W. Structural diversity in free and bound states of intrinsically disordered protein phosphatase 1 regulators. Structure 2010, 18, 1094–1103. [Google Scholar]
    53. Gunasekaran, K; Tsai, C-J; Nussinov, R. Analysis of ordered and disordered protein complexes reveals structural features discriminating between stable and unstable monomers. J. Mol. Biol 2004, 341, 1327–1341. [Google Scholar]
    54. Mészáros, B; Tompa, P; Simon, I; Dosztányi, Z. Molecular principles of the interactions of disordered proteins. J. Mol. Biol 2007, 372, 549–561. [Google Scholar]
    55. Oldfield, CJ; Meng, JW; Yang, JY; Yang, MQ; Uversky, VN; Dunker, AK. Flexible nets: Disorder and induced fit in the associations of p53 and 14-3-3 with their partners. BMC Genomics 2008, 9, S1. [Google Scholar]
    56. Solt, I; Magyar, C; Simon, I; Tompa, P; Fuxreiter, M. Phosphorylation-induced transient intrinsic structure in the kinase-inducible domain of CREB facilitates its recognition by the KIX domain of CBP. Proteins 2006, 64, 749–757. [Google Scholar]
    57. Ganguly, D; Chen, JH. Atomistic details of the disordered states of KID and pKID. Implications in coupled binding and folding. J. Am. Chem. Soc 2009, 131, 5214–5223. [Google Scholar]
    58. Joerger, AC; Fersht, AR. Structural biology of the tumor suppressor p53. Annu. Rev. Biochem 2008, 77, 557–582. [Google Scholar]
    59. Lee, H; Mok, KH; Muhandiram, R; Park, K-H; Suk, J-E; Kim, D-H; Chang, J; Sung, YC; Choi, KY; Han, K-H. Local structural elements in the mostly unstructured transcriptional activation domain of human p53. J. Biol. Chem 2000, 275, 29426–29432. [Google Scholar]
    60. Espinoza-Fonseca, LM; Trujillo-Ferrara, JG. Transient stability of the helical pattern of region F19-L22 of the N-terminal domain of p53: A molecular dynamics simulation study. Biochem. Biophys. Res. Comm 2006, 343, 110–116. [Google Scholar]
    61. Xu, HB; Ye, H; Osman, NE; Sadler, K; Won, E-Y; Chi, S-W; Yoon, HS. The MDM2-binding region in the transactivation domain of p53 also acts as a Bcl-XL-binding motif. Biochemistry 2009, 48, 12159–12168. [Google Scholar]
    62. Kussie, PH; Gorina, S; Marechal, V; Elenbaas, B; Moreau, J; Levine, AJ; Pavletich, NP. Structure of the MDM2 oncoprotein bound to the p53 tumor suppressor transactivation domain. Science 1996, 274, 948–953. [Google Scholar]
    63. Chen, H-F; Luo, R. Binding induced folding in p53-MDM2 complex. J. Am. Chem. Soc 2007, 129, 2930–2937. [Google Scholar]
    64. Lowry, DF; Stancik, A; Shrestha, RM; Daughdrill, GW. Modeling the accessible conformations of the intrinsically unstructured transactivation domain of p53. Proteins 2008, 71, 587–598. [Google Scholar]
    65. Dunbrack, RL; Cohen, FE. Bayesian statistical analysis of protein side-chain rotamer preferences. Protein Sci 1997, 6, 1661–1681. [Google Scholar]
    66. Dastidar, SG; Lane, DP; Verma, CS. Modulation of p53 binding to MDM2: Computational studies reveal important roles of Tyr100. BMC Bioinformatics 2009, 10, S6. [Google Scholar]
    67. Dastidar, SG; Madhumalar, A; Fuentes, G; Lane, DP; Verma, CS. Forces mediating protein-protein interactions: A computational study of p53 “approaching” MDM2. Theor. Chem. Acc 2010, 125, 621–635. [Google Scholar]
    68. Uhrinova, S; Uhrin, D; Powers, H; Watt, K; Zheleva, D; Fischer, P; McInnes, C; Barlow, PN. Structure of free MDM2 N-terminal domain reveals conformational adjustments that accompany p53-binding. J. Mol. Biol 2005, 350, 587–598. [Google Scholar]
    69. Li, C; Pazgier, M; Li, CQ; Yuan, WR; Liu, M; Wei, G; Lu, W-Y; Lu, W. Systematic mutational analysis of peptide inhibition of the p53-MDM2/MDMX interactions. J. Mol. Biol 2010, 398, 200–213. [Google Scholar]
    70. Zondlo, SC; Lee, AE; Zondlo, NJ. Determinants of specificity of MDM2 for the activation domains of p53 and p65: Proline27 disrupts the MDM2-binding motif of p53. Biochemistry 2006, 45, 11945–11957. [Google Scholar]
    71. Guharoy, M; Janin, J; Robert, CH. Side-chain rotamer transitions at protein-protein interfaces. Proteins 2010, 78, 3219–3225. [Google Scholar]
    72. Sivakolundu, SG; Bashford, D; Kriwacki, RW. Disordered p27Kip1 exhibits intrinsic structure resembling the Cdk2/Cyclin A-bound conformation. J. Mol. Biol 2005, 353, 1118–1128. [Google Scholar]
    73. Jensen, MR; Houben, K; Lescop, E; Blanchard, L; Ruigrok, RWH; Blackledge, M. Quantitative conformational analysis of partially folded proteins from residual dipolar couplings: application to the molecular recognition element of Sendai virus nucleoprotein. J. Am. Chem. Soc 2008, 130, 8055–8061. [Google Scholar]
    74. Ghosh, RP; Nikitina, T; Horowitz-Scherer, RA; Gierasch, LM; Uversky, VN; Hite, K; Hansen, JC; Woodcock, CL. Unique physical properties and interactions of the domains of methylated DNA binding protein 2. Biochemistry 2010, 49, 4395–4410. [Google Scholar]
    75. Gely, S; Lowry, DF; Bernard, C; Jensen, MR; Blackledge, M; Costanzo, S; Bourhis, J-M; Darbon, H; Daughdrill, G; Longhi, S. Solution structure of the C-terminal X domain of the measles virus phosphoprotein and interaction with the intrinsically disordered C-terminal domain of the nucleoprotein. J. Mol. Recognit 2010, 23, 435–447. [Google Scholar]
    76. Jiao, WT; McDonald, DQ; Coxon, JM; Parker, EJ. Molecular modeling studies of peptide inhibitors highlight the importance of conformational prearrangement for inhibition of calpain. Biochemistry 2010, 49, 5533–5539. [Google Scholar]
    77. Kjaergaard, M; Teilum, K; Poulsen, FM. Conformational selection in the molten globule state of the nuclear coactivator binding domain of CBP. Proc. Natl. Acad. Sci. USA 2010, 107, 12535–12540. [Google Scholar]
    78. Ahmad, M; Gu, W; Helms, V. Mechanism of fast peptide recognition by SH3 domains. Angew. Chem. Int. Ed 2008, 47, 7626–7630. [Google Scholar]
    79. Fischer, E. Einfluss der configuration auf die wirkung der enzyme. Ber. Dtsch. Chem. Ges 1894, 27, 2985–2993. [Google Scholar]
    80. Koshland, DE. Application of a theory of enzyme specificity to protein synthesis. Proc. Natl. Acad. Sci. USA 1958, 44, 98–104. [Google Scholar]
    81. Monod, J; Wyman, J; Changeux, J-P. On the nature of allosteric transitions: A plausible model. J. Mol. Biol 1965, 12, 88–118. [Google Scholar]
    82. Tsai, CJ; Kumar, S; Ma, BY; Nussinov, R. Folding funnels, binding funnels, and protein function. Protein Sci 1999, 8, 1181–1190. [Google Scholar]
    83. Boehr, DD; Nussinov, R; Wright, PE. The role of dynamic conformational ensembles in biomolecular recognition. Nat. Chem. Biol 2009, 5, 789–796. [Google Scholar]
    84. Huang, YQ; Liu, ZR. Smoothing molecular interactions: The “kinetic buffer” effect of intrinsically disordered proteins. Proteins 2010, 78, 3251–3259. [Google Scholar]
    85. Dunker, AK; Garner, E; Guilliot, S; Romero, P; Albrecht, K; Hart, J; Obradovic, Z; Kissinger, C; Villafranca, JE. Protein disorder and the evolution of molecular recognition: Theory, predictions and observations. Pac. Symp. Biocomput 1998, 3, 473–484. [Google Scholar]
    86. Dunker, AK; Cortese, MS; Romero, P; Iakoucheva, LM; Uversky, VN. Flexible nets-The roles of intrinsic disorder in protein interaction networks. FEBS J 2005, 272, 5129–5148. [Google Scholar]
    87. Ngan, C-H; Beglov, D; Rudnitskaya, AN; Kozakov, D; Waxman, DJ; Vajda, S. The structural basis of pregnane X receptor binding promiscuity. Biochemistry 2009, 48, 11572–11581. [Google Scholar]
    88. Feng, H; Jenkins, LMM; Durell, SR; Hayashi, R; Mazur, SJ; Cherry, S; Tropea, JE; Miller, M; Wlodawer, A; Appella, E; Bai, Y. Structural basis for p300 Taz2-p53 TAD1 binding and modulation by phosphorylation. Structure 2009, 17, 202–210. [Google Scholar]
    89. Berendsen, HJC; van der Spoel, D; van Drunen, R. GROMACS: A message-passing parallel molecular dynamics implementation. Comput. Phys. Commun 1995, 91, 43–56. [Google Scholar]
    90. Hess, B; Kutzner, C; Spoel, Dvd; Lindahl, E. GROMACS 4: Algorithms for highly efficient, load-balanced, and scalable molecular simulation. J. Chem. Theory Comput 2008, 4, 435–447. [Google Scholar]
    91. Kaminski, GA; Friesner, RA; Tirado-Rives, J; Jorgensen, WL. Evaluation and reparametrization of the OPLS-AA force field for proteins via comparison with accurate quantum chemical calculations on peptides. J. Phys. Chem. B 2001, 105, 6474–6487. [Google Scholar]
    92. Berendsen, HJC; Grigera, JR; Straatsma, TP. The missing term in effective pair potentials. J. Phys. Chem 1987, 91, 6269–6271. [Google Scholar]
    93. Darden, T; York, D; Pedersen, L. Particle mesh Ewald: An N·log(N) method for Ewald sums in large systems. J. Chem. Phys 1993, 98, 10089–10092. [Google Scholar]
    94. Hess, B; Bekker, H; Berendsen, HJC; Fraaije, JGEM. LINCS: A linear constraint solver for molecular simulations. J. Comput. Chem 1997, 18, 1463–1472. [Google Scholar]
    95. Bussi, G; Donadio, D; Parrinello, M. Canonical sampling through velocity rescaling. J. Chem. Phys 2007, 126, 014101. [Google Scholar]
    96. Parrinello, M; Rahman, A. Polymorphic transitions in single crystals: A new molecular dynamics method. J. Appl. Phys 1981, 52, 7182–7190. [Google Scholar]
    97. Kabsch, W; Sander, C. Dictionary of protein secondary structure: Pattern recognition of hydrogen-bonded and geometrical features. Biopolymers 1983, 22, 2577–2637. [Google Scholar]
    98. Humphrey, W; Dalke, A; Schulten, K. VMD: Visual molecular dynamics. J. Mol. Graphics 1996, 14, 33–38. [Google Scholar]
    Ijms 12 01410f1 1024
    Figure 1. The analysis of the structures of p53N when bound to different targets. (a,b) Complex structures of the p53N-MDM2 complex (a) and p53N-Taz2 complex (b). (c) Structure superposition of p53N from p53N-MDM2 and p53N-Taz2 based on the backbone RMSD. In (a–c), the MDM2 and Taz2 are presented as surface, whereas the p53N is presented as ribbon. Phe19, Trp23, and Leu26 are shown as blue balls in (a,b) and sticks in (c), Leu22 is shown in red, whereas other residues of the p53N are colored cyan (when in complex with MDM2) or green (when in complex with Taz2). (d,e) SASA analysis of the p53N in complex with MDM2 (d) and Taz2 (e).
    Figure 1. The analysis of the structures of p53N when bound to different targets. (a,b) Complex structures of the p53N-MDM2 complex (a) and p53N-Taz2 complex (b). (c) Structure superposition of p53N from p53N-MDM2 and p53N-Taz2 based on the backbone RMSD. In (a–c), the MDM2 and Taz2 are presented as surface, whereas the p53N is presented as ribbon. Phe19, Trp23, and Leu26 are shown as blue balls in (a,b) and sticks in (c), Leu22 is shown in red, whereas other residues of the p53N are colored cyan (when in complex with MDM2) or green (when in complex with Taz2). (d,e) SASA analysis of the p53N in complex with MDM2 (d) and Taz2 (e).
    Ijms 12 01410f1
    Ijms 12 01410f2 1024
    Figure 2. RMSD relative to the backbone of the helical region (Phe19–Leu26) for p53N from p53N-MDM2 (a) and from p53-Taz2 (b) during the simulations. Simulation trajectories MS1–MS5 and TS1–TS5 are presented.
    Figure 2. RMSD relative to the backbone of the helical region (Phe19–Leu26) for p53N from p53N-MDM2 (a) and from p53-Taz2 (b) during the simulations. Simulation trajectories MS1–MS5 and TS1–TS5 are presented.
    Ijms 12 01410f2
    Ijms 12 01410f3 1024
    Figure 3. Conformational analysis of the anchor residues in the helix state and disordered state of the unbound p53N. The bound conformation values in the p53N-MDM2 and p53N-Taz2 complexes are denoted by red and green markers, respectively.
    Figure 3. Conformational analysis of the anchor residues in the helix state and disordered state of the unbound p53N. The bound conformation values in the p53N-MDM2 and p53N-Taz2 complexes are denoted by red and green markers, respectively.
    Ijms 12 01410f3
    Ijms 12 01410f4 1024
    Figure 4. Combined analysis of the conformations of Phe19 and Trp23 in the helix state. (a) Side chain dihedral distribution of Phe19; (b) Side chain dihedral distribution of Trp23; (c) The combination of conformations of Phe19 and Trp23.
    Figure 4. Combined analysis of the conformations of Phe19 and Trp23 in the helix state. (a) Side chain dihedral distribution of Phe19; (b) Side chain dihedral distribution of Trp23; (c) The combination of conformations of Phe19 and Trp23.
    Ijms 12 01410f4
    Ijms 12 01410f5 1024
    Figure 5. Conformational analysis of non-anchor residues in the helix state. Bound conformations in the p53N-MDM2 and p53N-Taz2 complexes are denoted by red and green markers, respectively.
    Figure 5. Conformational analysis of non-anchor residues in the helix state. Bound conformations in the p53N-MDM2 and p53N-Taz2 complexes are denoted by red and green markers, respectively.
    Ijms 12 01410f5
    Ijms 12 01410f6 1024
    Figure 6. The roles of anchor residue conformations in the evolution of the transient encounter complex towards the bound state. (a,b) p53N in a bound-like conformation; (c,d) p53N in a distorted conformation; (e,f) p53N in a bound-like conformation. p53N is shown as cyan ribbons and the binding targets are shown as surface. Phe19, Trp23, and Leu26, are shown as blue sticks. Leu22 is shown as red sticks.
    Figure 6. The roles of anchor residue conformations in the evolution of the transient encounter complex towards the bound state. (a,b) p53N in a bound-like conformation; (c,d) p53N in a distorted conformation; (e,f) p53N in a bound-like conformation. p53N is shown as cyan ribbons and the binding targets are shown as surface. Phe19, Trp23, and Leu26, are shown as blue sticks. Leu22 is shown as red sticks.
    Ijms 12 01410f6
    Table
    Table 1. Population of bound-like conformations of anchor residues of the unbound p53N in different states.
    Table 1. Population of bound-like conformations of anchor residues of the unbound p53N in different states.
    Complex comparedState of p53NPopulation of bound-like conformation *
    Phe19Leu22Trp23Leu26
    p53N-MDM2Helix59.5%–62.4%82.7%
    Disordered19.2%–18.7%76.2%
    Rotamer library31.71%3.65%16.21%62.52%
    p53N-Taz2Helix12.2%81.8%––
    Disordered6.0%67.1%––
    Rotamer library47.08%62.52%5.32%<1%
    *“–” indicates the conformation in the bound state does not match a major group in the simulated χ1–χ2 distribution. Data for the rotamer library were adopted from Reference [65].

    Share and Cite

    MDPI and ACS Style

    Huang, Y.; Liu, Z. Anchoring Intrinsically Disordered Proteins to Multiple Targets: Lessons from N-Terminus of the p53 Protein. Int. J. Mol. Sci. 2011, 12, 1410-1430. https://doi.org/10.3390/ijms12021410

    AMA Style

    Huang Y, Liu Z. Anchoring Intrinsically Disordered Proteins to Multiple Targets: Lessons from N-Terminus of the p53 Protein. International Journal of Molecular Sciences. 2011; 12(2):1410-1430. https://doi.org/10.3390/ijms12021410

    Chicago/Turabian Style

    Huang, Yongqi, and Zhirong Liu. 2011. "Anchoring Intrinsically Disordered Proteins to Multiple Targets: Lessons from N-Terminus of the p53 Protein" International Journal of Molecular Sciences 12, no. 2: 1410-1430. https://doi.org/10.3390/ijms12021410

    APA Style

    Huang, Y., & Liu, Z. (2011). Anchoring Intrinsically Disordered Proteins to Multiple Targets: Lessons from N-Terminus of the p53 Protein. International Journal of Molecular Sciences, 12(2), 1410-1430. https://doi.org/10.3390/ijms12021410

    Article Metrics

    Back to TopTop