Skip Navigation

This Article
Right arrow Abstract Freely available
Right arrow Print PDF (394K) Freely available
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in ISI Web of Science
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Add to My Personal Archive
Right arrow Download to citation manager
Right arrow Search for citing articles in:
ISI Web of Science (9)
Right arrowRequest Permissions
Right arrow Commercial Re-use Guidelines
for Open Access NAR Content
Google Scholar
Right arrow Articles by Wright, J. D.
Right arrow Articles by Lim, C.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Wright, J. D.
Right arrow Articles by Lim, C.
Social Bookmarking
Add to CiteULike   Add to Connotea   Add to Del.icio.us  
What's this?

Nucleic Acids Research, 2002, Vol. 30, No. 7 1563-1574
© 2002 Oxford University Press

Factors governing loss and rescue of DNA binding upon single and double mutations in the p53 core domain

Jon D. Wright1, Sergey Yu Noskov1 and Carmay Lim1,2,*

1Institute of Biomedical Sciences, Academia Sinica, 11529 Taipei, Taiwan and 2Department of Chemistry, National Tsing-Hua University, 300 Hsinchu, Taiwan

Received December 14, 2001; Revised and Accepted February 11, 2002.


   ABSTRACT
TOP
 ABSTRACT
INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 REFERENCES
 
The mutation of R273->H in the p53 core domain (p53-CD) is one of the most common mutations found in human cancers. Although the 273H p53-CD retains the wild-type conformation and stability, it lacks sequence-specific DNA binding, a transactivation function and growth suppression. However, mutating T284->R in the 273H p53-CD restores the DNA binding affinity, and transactivation and tumour suppressor functions. Since X-ray/NMR structures of DNA-free or DNA-bound mutant p53-CD molecules are unavailable, the factors governing the loss and rescue of sequence-specific DNA binding in the 273H and 273H+284R p53-CD, respectively, are unclear. Hence, we have carried out molecular dynamics (MD) simulations of the wild-type, single mutant and double mutant p53-CD, free and DNA bound, in the presence of explicit water molecules. Based on the MD structures, the DNA-binding free energy of each p53 molecule has been computed and decomposed into component energies and contributions from the interface residues. The wild-type and mutant p53-CD MD structures were found to be consistent with the antibody-binding, X-ray and NMR data. The predicted DNA binding affinity and specificity of both mutant p53-CDs were also in accord with experimental data. The non-detectable DNA binding of the 273H p53-CD is due mainly to the disruption of a hydrogen-bonding network involving R273, D281 and R280, leading to a loss of major groove binding by R280 and K120. The restoration of DNA binding affinity and specificity of the 273H+284R p53-CD is due mainly to the introduction of another DNA-binding site at position 284, leading to a recovery of major groove binding by R280 and K120. The important role of water molecules and the DNA major groove conformation as well as implications for structure-based linker rescue of the 273H p53-CD DNA-binding affinity are discussed.


   INTRODUCTION
TOP
 ABSTRACT
 INTRODUCTION
MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 REFERENCES
 
Mutations of the tumour supressor protein, p53, are associated with ~50% of human cancers (1) and 95% of lung cancers (2). They are also associated with resistance to therapy (3). The wild-type p53 protein is a key regulator in several cellular processes, including cell cycle control, DNA repair, programmed cell death, differentiation, senescence, angiogenesis and genome stability (4). Following DNA damage, or other genotoxic stress such as hypoxia, oncogene activation and viral infection (5,6), p53 promotes cell growth arrest and/or cell death, thus protecting the genome from accumulating excess mutations. The human p53 protein comprises of 393 amino acids that can be grouped into four functional domains/regions (4). The first 42 amino acids compose a transcriptional transactivation region. Residues 102–292 make up the sequence-specific DNA-binding core/central domain (p53-CD), for which there is a 2.2 Å X-ray structure of a 21 bp DNA duplex and three core domains; two of the core domains bind DNA, while the third does not (7). Residues 323–356 compose the tetramerisation domain for which there are NMR structures (8,9) and an X-ray structure (10). The last 30 residues form a basic regulatory region that can also bind to DNA or RNA, as well as act as a negative regulator of p53 sequence-specific DNA binding (11).

Most of the p53 mutations found in human cancers map to the core domain (12) and ~40% of them occur at only six hot-spots (R175, G245, R248, R249, R273 and R282) (1). The p53-CD·DNA X-ray structure shows the core domain containing a sandwich of a four-stranded and a five-stranded antiparallel ß sheet, and a loop–sheet–helix motif consisting of loop L1 (residues 113–123), a three-stranded ß sheet (124–135 and 271–274) and helix H2 (278–286) (7). Residues from the loop–sheet–helix motif bind to the DNA major groove, and one of the two loops stabilised by Zn2+ (loop L3, residues 237–250) binds to the minor groove. In particular, two of the hot-spot residues, R248 and R273, whose side chains are within 4 Å of phosphate oxygens in the minor and major DNA groove, respectively, make critical contributions to DNA binding. The other four hot-spot residues stabilise the p53 DNA-binding surface but do not contact DNA directly. Based on the p53-CD·DNA X-ray structure and biochemical data, Cho et al. (7) classified the tumour-derived p53 mutants as conformational if they involve residues important for structural integrity but do not directly contact DNA (e.g. R175, G245, R249 and R282), or contact if they involve residues that directly contact DNA (e.g. R248 and R273). Failure of p53 conformational mutants to bind DNA in a sequence-specific manner has been attributed to structural defects, ranging from small structural shifts to local/global unfolding of the core domain. Failure of p53 contact mutants to act as a sequence-specific transactivator has been attributed to loss of critical DNA contacts.

There is considerable interest in the contact mutants because mutations at R248 and R273 account for almost one-fifth of all p53 mutations (1). Here, we have chosen to study the mutation of R273->H in the p53-CD since this is one of the most common mutations found in human cancer (7). The 273H p53 mutant retains the wild-type conformation (13) (i.e. PAb 1620+, PAb 240–, Hsp 70–) and stability (2) and is estimated to be 98% folded at physiological temperature, 37°C (14). However, the mutant 273H p53-CD lacks sequence-specific DNA binding (14) and sequence-specific transactivation function and growth suppression (15). In contrast, the mutation of T284->R significantly enhances not only the DNA binding of the wild-type protein with the last 30 residues deleted (referred to as p53{Delta}364–393), but also the corresponding 273H mutant (16). It can also rescue the DNA-binding affinity and specificity as well as transactivation and tumour suppressor functions of the 273H mutant (16). Since the DNA-free and DNA-bound X-ray/NMR structures of the single mutant and double mutant p53-CD have not been reported, the structural and energetic bases for the loss of DNA binding by the 273H p53-CD and the rescue of sequence-specific DNA binding by mutating T284->R in the single mutant are not clear.

Here, our goals are 2-fold. The first is to identify the factors governing the loss of DNA binding upon mutating R273 to a positively charged His in the p53-CD. The second is to elucidate the factors governing the rescue of sequence-specific DNA binding upon introducing a positively charged Arg at 284 in the mutant 273H p53-CD (see above). First, the three-dimensional (3d) solution models of structures of the wild-type, single mutant and double mutant p53-CD, free and DNA bound, were obtained using molecular dynamics (MD) simulations with explicit water molecules. Based on the MD structures, the free energies of wild-type and mutant p53 molecules binding to DNA were computed and decomposed into component energies and contributions from individual residues. To assess the validity of the theoretical results obtained, the findings were verified by their consistency with available experimental data. The energetic and structural analyses help to elucidate the mechanism of p53-CD reactivation, which in turn would aid the search for suitable ligands that may restore the lost interactions with DNA.


   MATERIALS AND METHODS
TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
RESULTS
 DISCUSSION
 REFERENCES
 
MD simulations
The core domains of wild-type p53, mutant 273H and 273H+284R p53, free and bound to a 15 bp DNA, were subjected to MD simulation using the CHARMM (17) program at a mean temperature of 300 K. The simulations of p53-CD, 273H p53-CD and 273H+284R p53-CD included a total of 24802, 24358 and 25061 atoms, while the DNA-bound simulations included 58389, 58416 and 58417 atoms.

Starting structures. The initial 3d structure for the wild-type p53-CD·DNA simulation was taken from chain B of PDB entry 1TSR, in which the core domain (residues 96–289) is bound to a consensus DNA binding site, GGGCA AGTCT (7). The DNA strands were truncated to provide a three-base overhang of the binding region, yielding a 15 bp DNA duplex, 5'-3AATTGGGCAAGTCTA17-3'. Initial structures for the mutant p53-CD·DNA simulations were prepared by a side chain rotamer library method (18,19), as described in our previous work (20). Starting structures for the DNA-free simulations were extracted from the averaged MD structure of each complex.

Protonation states. Each simulation was carried out at physiological pH (pH 7) with the following protonation state assignments for the ionisable residues. The three Cys residues that are coordinated to zinc and all Asp, Glu residues and C-terminal groups were deprotonated. The other Cys and all Arg, Lys residues and N-terminal groups were protonated. The His that is coordinated to zinc was deprotonated as well as H214 because of its close proximity to R174. The other His residues were protonated since they are exposed to solvent or in close proximity to carboxylate side chains. In particular, the 273H mutant side chain is close to the negatively charged D281 carboxylate and Thy 11' phosphate. The protonation state of 273H was verified by calculations (data not shown) showing that of all the His residues in p53, the 273H side chain nitrogens have the most negative electrostatic potential, indicating the highest probability of being protonated. With the inclusion of 20 counterions, the net charge of the solvated Zn2+-bound p53-DNA-ion system is –3.

Forcefield and boundary conditions. The simulations employed the all-hydrogen CHARMM forcefield (21) and the TIP3P model of water (22). The vdW interactions were truncated at 15 Å by a potential-shifting function, while the electrostatic interactions were switched from 11 Å to zero at 15 Å by an atom-based force-switching function. The non-bonded interaction list was updated using a heuristic test at a 16-Å cutoff. A spherical boundary potential was applied to the surface water molecules to account for interactions with the implicit bulk solvent (23).

Simulation protocol. The simulation procedure used is described in detail elsewhere (Noskov,S.Y., Wright,J.D. and Lim,C., in preparation.). The fully solvated protein–DNA complex was equilibrated for 60 ps with heavy constraints on only the terminal DNA residues to prevent the DNA strands from unwinding (24). These constraints were reduced and each simulation was continued for 250 ps, during which coordinates were saved every 100 fs.

In the X-ray structure, the root mean square deviation (RMSD) of C{alpha} atoms in the p53-CD·DNA from those in the DNA-free p53-CD is only 0.75 Å, indicating that DNA binding does not alter the protein structure significantly (7). Consequently, the solvated p53 protein was extracted from the average MD structure of each complex and subjected to 250 ps of MD simulation, assuming the same residue protonation state, forcefield and boundary conditions as for the bound-state simulations.

Free energy decomposition
Relative binding free energy. This was based on a thermodynamic cycle (25,26).

{Delta}Ggas

[p53]gas + [DNA]gas->[p53·DNA]gas

{uparrow}–Gsolv(p53){uparrow}–Gsolv(DNA){downarrow}Gsolv(p53·DNA) (Scheme 1)[p53]sln + [DNA]sln-> [p53·DNA]sln

{Delta}G°

The standard free energy change ({Delta}G°) for the binding of the p53-CD to DNA in solution according to Scheme 1 is given by:

{Delta}G° = {Delta}Ggas + Gsolv(p53·DNA) – Gsolv(p53) – Gsolv(DNA) 1

where {Delta}Ggas is the standard free energy change per mole for the non-covalent association of the p53-CD and DNA in the gas phase at 300 K, and –Gsolv corresponds to the work of transferring the molecule in its solution conformation to the same conformation in the gas phase at 300 K.

The focus in this work is not {Delta}G° but {Delta}{Delta}G°, the free energy difference between the wild-type and mutant p53-CD binding to DNA in solution, which is given by:

{Delta}{Delta}G° = ({Delta}Ggasmut – {Delta}Ggaswt) + [Gsolv(p53mut·DNA) – Gsolv(p53wt·DNA)] – [Gsolv(p53mut) – Gsolv(p53wt)] 2

The relative free energy (equation 2) was based on ‘snapshot’ configurations taken every 5 ps from the last 100 ps of the MD trajectories (27,28). All water molecules were removed from the MD trajectories to avoid boundary problems for the different sized systems as the DNA-bound and free proteins had the same number of solute atoms but different numbers of solvent atoms (29). The calculations of the terms in equation 2 are outlined below, and further details can be found in our previous work (26).

Gas-phase binding free energy difference, {Delta}{Delta}Ggas. The free energy difference between wild-type and mutant p53-CD binding to DNA in the gas phase was estimated by:

{Delta}{Delta}Ggas ~ {Delta}{Delta}Egaselec + {Delta}{Delta}EgasvdW = ({Delta}Egaselec,mut – {Delta}Egaselec,wt) + ({Delta}EgasvdW,mut – {Delta}EgasvdW,wt) 3

Equation 3 assumes similar changes in conformation and vibrational entropy between wild-type and mutant p53-CDs binding to DNA. Based on the MD structures of wild-type and mutant p53-CDs (see above), the changes in vdW and electrostatic energies upon binding were calculated using a dielectric constant of one without cut-offs. The energies were computed using the CHARMM program (17) and forcefield (30).

Solvation free energy, Gsolv. The solvation free energy, Gsolv, can be expressed as a sum of three terms:

Gsolv = Gsolvcav + GsolvvdW + Gsolvelec 4

The first two terms in equation 4 were approximated by a linear function of the solvent accessible surface area (SASA) (31):

Gsolvcav + GsolvvdW = ({gamma} cav +{gamma}vdW) x SASA 5

{gamma} vdW was set to –38.8 cal/mol/Å2 (32), while {gamma} cav was set to 46.0 cal/mol/Å2 (33). The absolute SASA of the free or DNA-bound molecule was computed using the GEPOL program (34) and a solvent probe radius of 1.4 Å. The electrostatic contribution to the solvation free energy (Gsolvelec) was calculated using the finite-difference Poisson–Boltzmann method (35), as described in detail elsewhere (Noskov,S.Y., Wright,J.D. and Lim,C., in preparation).


   RESULTS
TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
DISCUSSION
 REFERENCES
 
For each simulation, the RMSD of the protein backbone atoms from the initial coordinates during the simulation (data not shown) shows that it initially rises and then plateaus, fluctuating around means between 1.2 and 1.4 Å during the last 150 ps of the simulation. Hence, the earlier part of each trajectory was included in the equilibration phase, while the final 150 ps were employed in generating an average structure, and in analyses such as computing the average distance between a light/heavy donor and heavy acceptor (see below). The p53 backbone RMSDs are similar to those found in other protein simulations (36–39).

Comparison of MD and X-ray p53 structures
To establish the reliability of the simulation results, hydrogen bonds (defined by a heavy donor–heavy acceptor distance <=3.5 Å, light donor–heavy acceptor distance <=2.5 Å and a deviation of less than ±60° from linearity) made by interface residues (defined by gas-phase interaction energy changes >1 kcal/mol upon binding to DNA) in the p53-CD·DNA simulation were compared to those in the X-ray structure. The hot-spot residues, R175, G245, R249 and R282, which do not contact DNA directly but are believed to play key structural roles in p53, were also used to verify the integrity of the simulation structures of the free and DNA-bound p53-CD. Among the 26 interface residues (see Tables 1 and 2), K120, S241, R273, A276 and R283 are within hydrogen-bonding distance of the major groove phosphate groups, K120, C277 and R280 are within hydrogen-bonding distance of the major groove bases in the X-ray structure, while the R248 side chain has ‘close’ contacts (defined by a heavy donor–heavy acceptor distance >3.5 Å but <=4 Å) to the Thy 12' and Thy 14 phosphate oxygens in the minor groove.


View this table:
[in this window]
[in a new window]
  		Table 1. Changes in DNA-binding free energy components of the mutant 273H p53-CD relative to wild typea

 

View this table:
[in this window]
[in a new window]
  		Table 2. Changes in DNA-binding free energy components of the 273H+284R p53-CD relative to the single 273H mutant

 

Protein–DNA interactions. Among the residues that interact with the DNA in the X-ray structure, R280 makes the most number (six) of contacts (not necessarily hydrogen bonds) within 4 Å of base atoms in the major groove, followed by K120 with five contacts to base or phosphate atoms. R248 has three contacts <=4 Å of a phosphate or sugar oxygen in the minor groove. S241, R273 and A276 each have two contacts <=4 Å of the phosphate oxygen, while R283 has one. The C277 side chain is within hydrogen-bonding distance of the Cyt 9' base N4. All the <=4 Å DNA contacts found in the X-ray structure are preserved in the simulation except the following.

The ‘close’ contact between R283 N{eta}2 and Gua 7 O1P (3.53 Å in the X-ray structure) is extended to 4.68 ± 0.41 Å in the simulation. However, the R283 side chain has a high average X-ray B-factor (69 Å2), indicating a flexible side chain that may not be able to maintain a stable salt bridge with the Gua 7 phosphate oxygen. In the major groove the R273 N{eta}2 is <3.20 Å to two Thy 11' phosphate oxygens in the X-ray structure, but is 3.95 ± 0.40 Å to Thy 11' O5' and 3.91 ± 0.37 Å to Gua 10' O1P in the simulation, while the R280 N{eta}1...Gua 9 O6 and R280 N{eta}2...Gua 10' O6 ‘close’ contacts in the X-ray structure are attenuated by an intervening water molecule (Fig. 1). New hydrogen bonds are formed between a R280 amino proton and the D281 carboxylate oxygen (Fig. 1) as well as between the K120 backbone oxygen and the G279 and R280 backbone amide protons in the simulation (not shown in Fig. 1). In the minor groove R248 N{eta}1 is closer to Thy 12' O3' in the simulation (3.41 ± 0.37 Å) than in the X-ray structure (3.83 Å) but it does not form a hydrogen bond (R248 H{eta}1...Thy 12' O3' = 2.88 ± 0.58 Å), while it is in close contact with Thy 11' O2P (3.87 ± 0.55 Å) instead of Thy 14 O3' and O4' in the X-ray structure (3.83 Å).



View larger version (45K):
[in this window]
[in a new window]
  		Figure 1. Hydrogen-bonding network in the average wild-type p53-CD·DNA structure in the vicinity of the major and minor grooves. Only certain key residues (discussed in text) and their hydrogen bonding partners (including specific DNA base/phosphate groups and water molecules) are shown for the sake of clarity. However, a water molecule hydrogen bonding to the R249 carbonyl oxygen could not be depicted. The loop–sheet–helix motif is in magenta while the L3 loop is in yellow. Oxygen, phosphate, hydrogen and nitrogen atoms are coloured red, orange, green and blue, respectively.

 

Protein–water–DNA interactions. Three water molecules form bridges between the protein and the DNA in the p53-CD·DNA simulation (Fig. 1). The first is a water that bridges the R273 side chain and the Gua 10' phosphate. The second water bridge is between the R280 side chain and the Gua 9 base oxygen (see above). The third water bridge, which is also seen in the crystal structure, links the R248 side chain with the Thy 12' sugar oxygen. Although the oxygen of this water is within hydrogen-bonding distance to the Gua 13 N (3.35 ± 0.34 Å), the water hydrogen to Gua 13 N distance (2.66 ± 0.44 Å) indicates a non-ideal hydrogen bond.

Protein–protein interactions. All the hydrogen bonds made by the four hot-spot ‘structural’ residues, R175, G245, R249 and R282, in the free and DNA-bound p53-CD X-ray structure are preserved in the simulation except the following. The R175 N{eta}1...D184 O{delta}2 salt bridge in the p53-CD X-ray structure is absent from the DNA-bound crystal structure and in the DNA-free simulation. In the p53-CD·DNA X-ray structure, the R175 N{eta}1 is 2.59 Å to S183 O{gamma} and 3.41 Å to P191 O, but in the simulation the hydrogen bond to S183 O{gamma} is maintained (3.10 ± 0.26 Å), but that to P191 O is unstable (4.56 ± 1.13 Å). In the p53-CD·DNA crystal structure the R282 amino nitrogens (N{eta}) hydrogen bond to Y126 O and T118 O{gamma}, while its imino nitrogen (N{varepsilon}) makes bifurcated hydrogen bonds to the E286 and S127 side chain oxygens. However, in the simulation the R282 N{eta} hydrogen bonds to T118 O{gamma} but not to Y126 O, while its N{varepsilon} breaks and forms hydrogen bonds with the E286 carboxylate oxygens (average R282 N{varepsilon}...E286 O{varepsilon} distance = 3.75 ± 0.96 Å).

Dynamical fluctuations. The B-factors (data not shown) averaged over the backbone atoms, N, C{alpha} and C(=O), in the p53-CD·DNA simulation generally reproduce the qualitative features of the X-ray main-chain B-factors with the largest discrepancy found in two loops; the L2 loop (residues 164–194) and the loop connecting ß-strands S7 and S8 (residues 220–229). Apart from these two loops, the computed average B-factors are generally lower than the X-ray values, as found and rationalised in previous studies (36,39–41). The RMSD of the DNA backbone atoms (P, O5', C5', C4', C3' and O3') in the wild-type simulation from the corresponding X-ray structure (data not shown) shows that the major groove region remained relatively stable with an average RMSD of 1.4 Å for base pairs 6–10, while the minor groove was more flexible with an average RMSD of 1.9 Å for base pairs 12–14. The largest deviation was found for base pair 11 and 11' (RMSD = 2.55 Å).

Changes in component free energies upon mutation
In the following, we will first present the free energy decomposition results for the 273H p53-CD single mutant, then the 273H+284R p53-CD double mutant. The contributions to {Delta}{Delta}G° (see Material and Methods) from the protein–DNA interface residues were decomposed according to the electrostatic, vdW and cavity components:

6

The summation in equation 6 is over all interface residues listed in Tables 1 and 2. The net electrostatic and vdW contributions were further broken down into protein–solvent (denoted by subscript ‘solv’) and protein–protein or protein–DNA interactions (collectively known as solute–solute interactions, denoted by subscript ‘gas’); i.e.

7a

7b

The individual terms in equations 6 and 7 for the 273H and 273H+284R p53 mutants are summarised in Tables 1 and 2, respectively.

In the following, a single {Delta} denotes the free energy difference between the DNA-bound and free wild-type or mutant p53-CD. A double {Delta}{Delta} represents the DNA-binding free energy of the single mutant 273H p53 mutant relative to the wild-type protein, whereas {Delta}{Delta}' represents the binding free energy of the double mutant 273H+284R p53 relative to the single mutant. Interface residues with

and changes in the corresponding residue–protein+DNA and residue–solvent interaction energies upon mutation were identified. The energy changes were correlated with changes in the respective residue–protein, residue–DNA and residue–solvent hydrogen-bonding interactions as well as solvent accessibility seen in the simulations. In interpreting the results, we emphasise the trend (sign) of the large free energy changes upon each mutation rather than their absolute values.

R273->H. When R273 is mutated to His, the p53-CD mutant is predicted to bind to DNA far less well than the wild-type protein (total {Delta}{Delta}G° = 21.0 kcal/mol). Not surprisingly, the sum of the binding free energy contributions from all interface residues in the mutant is much less favourable than that in the wild-type protein

Table 1). This is because the net solute–solute electrostatic and vdW contributions of the interface residues with DNA binding in the mutant are unfavourable compared with those in wild-type p53

and

Table 1). The latter term is roughly cancelled by favourable solute–solvent electrostatic contributions

Table 1). The remaining solute–solvent vdW and solvent–solvent cavity contributions of the interface residues to DNA binding in the mutant and wild-type protein are similar (Table 1).

The key residues contributing to the net loss of DNA binding across the protein–DNA interface are K120, R280 and D281 in the major groove region and R249 in the minor groove region ({Delta}{Delta}Gi >= 5.5 kcal/mol for i = K120, R249, R280 and D281; Table 1). These four residues as well as the mutant 273H residue have less favourable electrostatic interactions with the protein and/or DNA in the p53 mutant compared with wild-type p53 (positive {Delta}{Delta}Gi,gaselec for i = K120, R249, 273H, R280 and D281; Table 1). The positive {Delta}{Delta}Ggaselec for the mutant 273H, R280 and D281 is due to the disruption of a hydrogen-bonding network involving the R273, D281 and R280 side chains, the Thy 11' phosphate and the Gua 10' base upon mutating R273 to His in the wild-type p53-CD (compare Figs 2 and 1). This has three manifestations. First, the observed loss of protein–protein and protein–DNA interactions in the major groove of the mutant complex is partially compensated by a gain in protein–solvent interactions: the D281 side chain gains three water hydrogen bonds, while the D281 and K120 backbone amide and the T284 side chain oxygen each gain one solvent hydrogen bond (Fig. 2). Second, the loss of protein–protein/DNA hydrogen bonds but gain in protein–solvent interactions in the major groove correlates with an increase in the backbone RMSDs of the major groove base pairs 6–10 and B-factors in the mutant complex simulation relative to the wild-type one (Fig. 3, dotted curves). Third, the behaviour of the RMSD as a function of base pair number in the mutant complex simulation differs from that in the wild-type X-ray structure; e.g., base pair 7 has the highest average RMSD in the mutant simulation but the lowest one in the wild-type X-ray structure (Fig. 3A). This finding, in conjunction with the high major groove backbone RMSDs in the mutant simulation, indicates a change in the major groove conformation (but not its width) in the 273H p53 complex.



View larger version (44K):
[in this window]
[in a new window]
  		Figure 2. Hydrogen-bonding network in the average 273H p53-CD·DNA structure in the vicinity of the major and minor grooves. See also the legend to Figure 1. The T284 HN...R280 O and D281 HN...O (water) hydrogen bonds are not shown for the sake of clarity.

 


View larger version (21K):
[in this window]
[in a new window]
  		Figure 3. (A) The RMSD of the DNA backbone atoms (averaged over base pair) from the wild-type simulation structure as a function of base pair number for wild-type crystal structure (solid curve), 273H p53 (dotted curve) and 273H+284R p53 (dashed curve). (B) The B-factors averaged over the protein backbone atoms as a function of residue number in the simulations of DNA bound to wild-type p53 (solid curve), 273H p53 (dotted curve) and 273H+284R p53 (dashed curve).

 
The distortion of the major groove conformation (in particular, base pairs 7–9) in the mutant complex is likely to cause K120 to lose three hydrogen bonds to the DNA; one from its mainchain nitrogen to the phosphate oxygen of Gua 7, and two from its side chain nitrogens to the base atoms of Gua 8 and Gua 9 (compare Figs 2 and 1), hence {Delta}{Delta}GK120,gaselec is positive. The positive {Delta}{Delta}Ggaselec for R249 is due to the loss of salt bridges between the R249 and E171 side chains upon mutating R273 to His (compare Figs 2 and 1).

Although the mutation of R273 to His results in the loss of direct charge–charge interactions with the DNA phosphate group and D281 carboxylate side chain (Fig. 2), the mutant 273H contributes more favourably than the wild-type Arg to DNA binding (negative {Delta}{Delta}G273H; Table 1). This is due to the large, negative {Delta}{Delta}G273H,solvelec (–18.2 kcal/mol), which more than compensates for the positive {Delta}{Delta}G273H,gaselec (3.6 kcal/mol; Table 1). The negative {Delta}{Delta}G273H,solvelec is because the mutant His has a smaller desolvation cost upon DNA binding compared with the wild-type Arg (151.4 versus 169.6 kcal/mol; Table 3).


View this table:
[in this window]
[in a new window]
  		Table 3. Electrostatic contributions to the DNA-binding free energy of wild-type and mutant p53-CD

 
In addition to the mutant residue, the net binding free energy contributions of V122, R248 and R283 in the p53 mutant are more favourable than those in wild-type p53 (negative {Delta}{Delta}Gi for i = V122, R248 and R283; Table 1). Unlike 273H, these three residues have more favourable electrostatic interactions with the protein and/or DNA in the mutant compared with the wild-type protein (negative {Delta}{Delta}Gi,gaselec for i = V122, R248 and R283; Table 1). The negative {Delta}{Delta}Ggaselec for V122, R248 and R283 is mainly due to hydrogen bonds formed in the mutant DNA simulation that are absent from the X-ray and MD structures of the wild-type complex. The interactions exclusive to the 273H p53 complex correspond to V122 N...K120 O (3.16 ± 0.20 Å), V122 O...Q136 N{varepsilon}2 (3.03 ± 0.32 Å) and R283 O...N288 N (3.03 ± 0.18 Å) hydrogen bonds, as well as salt bridges between the R248 guanidinium side chain and the Thy 12' base oxygen (Fig. 2). The negative {Delta}{Delta}Ggaselec for R248 and R283 are partially offset by positive {Delta}{Delta}Gsolvelec contributions (Table 1) probably due to the hydrogen bonds formed by these two residues, which reduce their solvent accessibility (Table 4), and hence their solvation free energies in the DNA-bound state become more unfavourable (Table 3).


View this table:
[in this window]
[in a new window]
  		Table 4. %SASA in wild-type and mutant p53-CD, free and DNA bounda

 

R273->H+T284->R. When T284 is mutated to Arg in the 273H p53-CD, the DNA binding ability of the double mutant becomes comparable to the wild-type p53-CD (total {Delta}{Delta}G° = –2.8 kcal/mol), and the sum of the binding free energy contributions from all interface residues in the double mutant is similar to that in the wild-type protein

. The improved DNA binding affinity of the 273H+284R p53-CD relative to the 273H p53-CD

is mainly because the net solute–solute electrostatic contribution of the interface residues to DNA binding in the double mutant is more favourable than that in the single mutant

although it is partially cancelled by an unfavourable solute–solvent electrostatic free energy

(Table 2). The net non-electrostatic solvation free energy contribution of the interface residues to DNA binding in the double mutant is similar to that in the single mutant as the positive solute–solvent vdW free energy

is compensated by a negative solvent–solvent cavity term

(Table 2).

The net free energy contribution of the 284R mutant residue to the improved DNA binding affinity of the double mutant is not significant ({Delta}{Delta}'G284R, = –1.5 kcal/mol; Table 2). Although the 284R side chain forms a salt bridge with the Thy 12' phosphate oxygen (Fig. 4) resulting in a large negative {Delta}{Delta}'G284R,gaselec (–192.8 kcal/mol; Table 2), the latter is offset by an even larger but positive {Delta}{Delta}'G284R,solvelec (195.7 kcal/mol; Table 2) due to the greater desolvation penalty of the positively charged 284R upon binding DNA compared with neutral T284 in the single mutant (207.6 versus 12.0 kcal/mol; Table 3). The net vdW contribution of the 284R mutant residue to the improved DNA binding is negligible ({Delta}{Delta}'G284RvdW = –0.1 kcal/mol; Table 2). In contrast, the solvent–solvent cavity {Delta}{Delta}'G284R,solvcav term is favourable (–4.2 kcal/mol; Table 2) due to the large SASA change at position 284 upon DNA binding ({Delta}SASA = –70 Å2 for 284R compared with –1.8 Å2 for T284 in the single mutant using the average MD structures).



View larger version (42K):
[in this window]
[in a new window]
  		Figure 4. Hydrogen-bonding network in the average 273H+284R p53-CD·DNA structure in the vicinity of the major and minor grooves. See also the legend to Figure 1. As for the wild-type and 273H p53-CD·DNA simulations, a 284R HN...R280 O hydrogen bond is not shown for the sake of clarity.

 
Instead of the mutant 284R, the residues that contribute to the loss of DNA binding upon mutating R273 to His in the wild-type p53-CD [namely K120, R249, C277 and R280 (see above and Table 1)] are now responsible for the gain in DNA binding upon mutating T284 to Arg in the 273H p53-CD ({Delta}{Delta}'Gi negative and {Delta}{Delta}'Gi >6 kcal/mol for i = K120, R249, C277 and R280; Table 2). The negative {Delta}{Delta}'Gi values for these four residues stem mainly from the negative {Delta}{Delta}'Gi,gaselec (Table 2), which is due to the restoration of hydrogen bonds formed by these four side chains in the wild-type protein (compare Figs 4 and 1). In the 273H+284R p53-CD·DNA simulation, the K120 and R280 side chains hydrogen bond to Gua 8, Gua 9 and Gua 10' base atoms in the major groove (Fig. 4), the R249 and E171 side chains form salt bridges (average distance = 3.08 ± 0.39 Å) in the minor groove, and C277 S{gamma} is close (4.00 ± 0.46 Å) to Cyt 9' N4 (not shown in Fig. 4). These interactions are lost in the single mutant complex simulation, but are present in the wild-type complex simulation and/or X-ray structure. However, the negative {Delta}{Delta}'Gi,gaselec values for K120, R249 and R280 are partially offset by positive {Delta}{Delta}'Gsolvelec (Table 2) because these three positively charged residues are better solvated in the DNA-free double mutant than in the unbound single mutant, hence their desolvation penalties in the double mutant are greater than those in the single mutant (Table 3).

The creation of a new hydrogen bond to the DNA phosphate by mutating T284->R in 273H p53 and the restoration of native hydrogen bonds to the major groove bases by K120, C277 and R280 has three other manifestations. First, the observed gain in protein–DNA and protein–protein interactions in the major groove of the 273H+284R p53 complex is at the expense of protein–solvent interactions. Relative to the single mutant complex, 273H, C277 and D281, respectively, lose 2, 1 and 1 water hydrogen bonds (see also %SASA in Table 4). Second, the recovery of protein–DNA/protein–protein hydrogen bonds and loss of protein–solvent interactions in the major groove correlates with a decrease in the backbone RMSDs of the major groove base pairs 6–10 and B-factors in the double mutant complex simulation compared with the single mutant one (compare dashed and dotted curves in Fig. 3). Third, the behaviour of the RMSD versus base pair number in the double mutant complex simulation is similar to that in the wild-type X-ray structure for the major groove base pairs (compare dashed and solid curves in Fig. 3A), suggesting similar major groove conformation in the double mutant and in the wild-type complex.

The residues that contribute favourably to DNA binding in the 273H p53-CD relative to the wild-type protein (namely R248, 273H and R283) now contribute unfavourably to DNA binding in the double mutant compared with the single mutant. The positive {Delta}{Delta}'G R248 and {Delta}{Delta}'G273H stems mainly from the positive {Delta}{Delta}'Gi,gaselec (Table 2). In the case of R248, this is probably because upon mutating T284 to Arg in the 273H p53-CD a direct hydrogen bond from the R248 side chain to the Thy 12' base oxygen in the single mutant complex (Fig. 2) is screened by a water in the double mutant complex (Fig. 4), which shields the electrostatic interactions. Although the 273H side chain forms a salt bridge with the D281 side chain in the double mutant (Fig. 4) but not in the single mutant (Fig. 2), its interaction with the DNA phosphate backbone is attenuated (273H N{varepsilon}...Thy 11' OP = 5.15 ± 0.86 Å compared with 4.14 ± 0.56 Å in the single mutant). Unlike R248 and 273H, the positive {Delta}{Delta}'GR283 stems from a large, positive {Delta}{Delta}'G283,solvelec (16.2 kcal/mol; Table 2), which more than offsets the negative {Delta}{Delta}'GR283,gaselec term (–13.1 kcal/mol; Table 2).

D281 contributes unfavourably to DNA binding upon mutating R273 to His in the wild-type protein ({Delta}{Delta}GD281 = 6.9 kcal/mol; Table 1) as well as upon mutating T284 to Arg in the single mutant ({Delta}{Delta}'GD281 = 11.1 kcal/mol; Table 2). The positive {Delta}{Delta}'GD281 term is dominated by a positive {Delta}{Delta}'GD281,gaselec (8.7 kcal/mol; Table 2) due probably to the greater solvent accessibility of the major groove in the single mutant compared with the double mutant. D281 hydrogen bonds to four and two water molecules in the single mutant and double p53 mutant complexes, respectively. The extra water molecules in the single mutant helps to shield the negatively charged D281 side chain from the negatively charged phosphate backbone, thus reducing the repulsive electrostatic interactions in the single mutant relative to the double mutant complex.


   DISCUSSION
TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
REFERENCES
 
Comparison with available experimental data
Wild-type structure. A comparison of the hydrogen bonds in the simulation and X-ray structures of the free and DNA-bound p53-CD shows that most of the hydrogen bonds observed in  the crystal structure are preserved in the simulation (see  Results). Although minor differences between the p53-CD·DNA X-ray and simulation structures are observed for the side chain interactions of R248, R273, R280, R282 and R283, they do not seem to be inconsistent with the relatively high average side chain X-ray B-factors (~30 Å2 for R273 and R280, 46 Å2 for R248 and R282, and 69 Å2 for R283). Thus, the free and DNA-bound p53-CD simulations have preserved the key protein–protein (and protein–DNA) contacts as well as the qualitative features of the X-ray B-factors.

Mutant structures. Although X-ray/NMR structures of the single mutant and double mutant p53-CD, free or DNA bound, are not available for comparison with the respective simulation structures, we can employ the fact that 273H p53 is recognised by the monoclonal antibody PAb 1620, but not by the monoclonal antibody PAb 240 to verify the mutant structure. PAb 1620 is specific for the wild-type conformation and particularly recognises R209 and N210 (42), while PAb 240 is specific for the ‘denatured’ conformation and recognises residues 212–217 on the S7 ß-strand (43). In the average 273H p53 structure, the average %SASA for R209 and N210 (46%; Table 4) indicates that they are solvent exposed, and can thus be recognised by PAb 1620, while the average %SASA for residues 212–217 (22%; Table 4) indicates that they are relatively buried, and are thus inaccessible to PAb 240. In analogy to the 273H p53 mutant, the average %SASA values for residues 209–210 (58%) and 212–217 (13%) in the average 273H+284R structure are similar to those in the respective wild-type structure (58 and 16%), indicating a wild-type conformation for the double mutant (Table 4).

A NMR study of the free 273H p53-CD has shown that structural changes in the mutant protein are found mainly in the loop–sheet–helix motif and the L3 loop, in particular residues D281, R282, E285 and E286 (44). This is consistent with the change in the hydrogen-bonding network involving these four residues in the simulation of the 273H p53-CD relative to the wild-type protein. Upon mutating R273 to His in the DNA-free p53-CD, the D281 carboxylate oxygens lose salt bridges to the 273H and R280 side chains, whereas the R282 and E286 side chains form salt bridges while the E285 carboxylate group forms new hydrogen bonds with K132. Thus, the simulation of the free 273H p53-CD seems to be consistent with both the antibody-binding and NMR data.

DNA binding affinity. Sequence-specific DNA binding to activate transcription of genes containing p53-binding DNA sites appears to be essential for normal p53 function (45,46). Such p53-binding DNA sites contain a repeat of two consensus pentamer sequences: 5'-Pu-Pu-Pu-C-(A/T) (T/A)-G-Py-Py-Py-3', where Pu represents a purine (Gua or Ade) and Py denotes a pyrimidine (Cyt or Thy) (47). Clearly, the binding affinity ({Delta}G°) of wild-type p53 for DNA depends not only on the DNA sequence, but also on the p53 molecule itself (i.e. if it is full length, minus the last 30 amino acids, or only the core domain) as well as on the temperature (46). For the consensus site studied here, 7GGGCA11 12AGTCT16, DNA binding specificity stems from hydrogen bonding interactions between the side chains of K120, C277 and R280 and the base oxygen or nitrogen of Gua 8, Gua 9, Cyt 9' and Gua 10' in the major groove, while non-specific binding affinity is due to salt bridges between K120, S241, R273, A276 and R283 and the DNA phosphate. Note that K120 and C277 could also hydrogen bond to proton-acceptor base atoms of Ade 8 and Thy 9', respectively, consistent with the variability of the base sequence at these two positions of the consensus site (7).

The calculations here predict that the DNA binding affinity of the mutant 273H p53-CD is much weaker than that of its wild-type counterpart (positive {Delta}{Delta}G°) and that it may not be detectable experimentally, as indicated by the relatively large magnitude of {Delta}{Delta}G° (21 kcal/mol). Furthermore, they predict that most of the lost interactions are between the K120, C277 and R280 side chains and the DNA major groove bases, in particular the invariant Gua 10', suggesting a loss of DNA binding specificity in the single mutant. These results are consistent with experiments showing that at 25°C the proteolytically excised core domain of wild-type, but not 273H p53, can bind specifically to a synthetic 46 bp DNA containing a site, 5'-TGGCA AGCCT-3', similar to the consensus DNA site used here (except at the underlined positions) (15). They are also consistent with experiments showing that the p53{Delta}364–393 and the p53-CD (at 20°C) bind to the gadd45 oligonucleotide containing the consensus site, 5'-GAACA TGTCT-3', but the respective proteins with R273 mutated to His do not (14,16).

The calculations also predict that the DNA binding affinity and specificity of the 273H+284R p53-CD is similar to that of the wild-type and is enhanced relative to the single mutant (see Results). This is consistent with experiments showing that the 273H+284R p53{Delta}364–393 binds to the gadd45 oligonucleotide, exhibited wild-type transcriptional activity and suppressed proliferation of tumour cells.

Factors affecting loss of DNA binding upon mutating R273->H in p53
In previous work (7), the 273H mutant was believed not to bind significantly to DNA due to the loss of the hydrogen bond between the wild-type R273 side chain and the Thy 11' phosphate upon mutation. In contrast, the calculations show that although the mutant H273 does not hydrogen bond to the Thy 11' phosphate as anticipated, it actually contributes to DNA binding more favourably than R273 in wild-type p53. This is because the mutant His has a smaller desolvation cost upon DNA binding compared with the wild-type Arg ({Delta}G273H,solvelec < {Delta}GR273,solvelec; Table 3). The calculations suggest that the loss of the salt bridge from the wild-type R273 to the D281 side chain (rather than to the Thy 11' phosphate) probably triggers the loss of DNA binding by the mutant 273H p53-CD (see below).

The lower DNA binding affinity of the mutant 273H p53-CD stems from the disruption of a hydrogen-bonding network involving R273, D281 and R280 in wild-type p53, leading to a loss of major groove binding by R280 and K120, which is partly compensated by a gain in minor groove binding by R248. When R273 is mutated to His, the shorter imidazolium side chain loses a close contact to the Thy 11' phosphate, and also a salt bridge to the D281 side chain, which in turn loses hydrogen bonds to the R280 guanidinium group (compare Figs 2 and 1). Consequently, R280 is not optimally oriented to bind into the major groove and loses its native hydrogen bonds with the invariant Gua 10' base. The loss of these protein–protein and protein–DNA hydrogen bonds in the 273H mutant complex is accompanied by (i) an overall increase in the flexibility of the mutant p53 (dotted curves in Fig. 3), (ii) an increase in the solvent accessibility of the major groove (Table 4) and (iii) a change in the major groove conformation (compare dotted and solid curves in Fig. 3A), resulting in a loss of DNA contacts by K120 and a gain in DNA contacts by R248 at the expense of its neighbour R249, which loses salt bridges with the E171 side chain (Fig. 2). The net effect is a loss of favourable protein–protein and protein–DNA electrostatic

positive) and vdW

positive) interactions as the protein–DNA interface becomes less well packed in the 273H p53 complex (Table 1).

Thus, the calculations suggest a central role for D281 in helping to orient R273 and R280 for hydrogen bonding to the DNA, consistent with the observation that D281 is strictly conserved across 27 species that p53 has been sequenced for (48). Mutation of D281 to Asn or the charge-preserving Glu in p53 has been found in tumour cells, implying that the mutant cannot bind to DNA specifically and function as a tumour suppressor gene. Furthermore, mutations of D281 to Tyr or Gly show a complete loss of DNA binding (49,50).

Factors governing rescue of DNA binding by mutating T284->R in 273H p53
In previous work (16), the mutant 273H+284R p53-CD was believed to restore DNA binding owing to the introduction of a positively charged residue in the vicinity of the DNA-binding interface that could establish novel protein–DNA contacts. As anticipated, the calculations show that the mutant 284R side chain forms a new hydrogen bond to the Thy 12' phosphate in the double mutant complex simulation. However, the favourable electrostatic interactions between the 284R side chain and the Thy 12' phosphate are offset by the higher cost of desolvating the positively charged 284R upon binding DNA compared with the neutral wild-type T284 ({Delta}G284R,solvelec >> {Delta}GT284,solvelec; Table 3). Consequently, the mutant 284R residue does not contribute as favourably as other residues (see below) to the improved DNA binding affinity of the double mutant relative to the single mutant (Table 2).

The introduction of another DNA-binding site at position 284 leads to a recovery of major groove binding by R280 and K120, and hence to a ‘wild-type’ DNA binding affinity and specificity of the double mutant. When T284 is mutated to Arg in the 273H p53-CD, the positively charged guanidinium side chain establishes a novel DNA contact that has the apparent effect of squeezing some water molecules out of the DNA binding surface (compare Figs 4 and 2; see also Table 4) and restoring the wild-type major groove conformation (compare dashed and solid curves in Fig. 3A). This enables R280 and K120 to regain native major groove contacts (compare Figs 4 and 1), which are essential for sequence-specific binding affinity. The net effect is a gain in favourable electrostatic and cavity free energies

and

negative; Table 2) in the double mutant compared with the single mutant.

Thus, the calculations suggest that the 284R mutant does not simply play a role by acting in lieu of R280 to bind to the DNA major groove in the 273H p53-CD. Rather, it plays a key role in helping R280 and K120 to regain their major groove contacts by helping to restore the ‘native’ conformation of the protein–DNA interface. Furthermore, since R280 and K120 contributes significantly to the improved DNA binding affinity of the double mutant relative to the single mutant (Table 2), a triple R273H+T284R+R280X p53-CD mutant would likely bind DNA with reduced affinity relative to the double mutant or wild-type p53-CD. This is consistent with the observation that the charge-preserving R280K p53 mutant is known not to bind to DNA (51), and that mutations of R280 to Lys, Ser, Thr, Ile, Gly and Pro in p53 are present in tumour cells (1).

Implications for rescue using small linker molecules
The mutant structures generated in this study can be used with molecular docking methods [such as those described in our previous work (52)] to screen databases of small molecule compounds to identify potential ligands that may restore the binding affinity of mutant p53-CD towards DNA. The energetic and structural analyses suggest that it may be possible to rescue the DNA binding affinity of the 273H mutant using an appropriate small ligand in the vicinity of the 284 site instead of the mutation site at position 273. The ligand should provide a short bridge to DNA in order to bring the DNA closer to the protein and force the removal of solvent water molecules from the DNA-binding interface. This may allow K120 and R280 to re-establish their native major groove contacts. The suitability of such a ligand could be verified using a combination of MD simulations and free energy decomposition calculations as presented in this work.


   ACKNOWLEDGEMENTS
 
We thank Drs C. S. Babu, T. W. Chang, Y. -S. Lin and H. Yuan for helpful discussions. We are grateful to Prof. M. Karplus for the CHARMM program. S.N. is supported by a postdoctoral fellowship from Academia Sinica. This work is supported by grant number NSC90-2311-b-001-155 from the National Science Council and the Institute of Biomedical Sciences, Academia Sinica, Taiwan, R.O.C.


   FOOTNOTES
 
* To whom correspondence should be addressed at: Institute of Biomedical Sciences, Academia Sinica, 11529 Taipei, Taiwan. Tel: +886 2 2652 3031; Fax: +886 2 2788 7641; Email: carmay{at}gate.sinica.edu.tw Permanent address: Sergey Yu Noskov, Institute of Solution Chemistry, Russian Academy of Sciences, Akademicheskaya str., 1, 153045 Ivanovo, Russia The authors wish it to be known that, in their opinion, the first two authors should be regarded as joint First Authors Back


   REFERENCES
TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 REFERENCES
 

  1. Hollstein,M., Rice,K., Soussi,T., Fuchs,R., Sorlie,T., Hovig,E., Smith-Sorensen,B., Montesano,R. and Harris,C.C. (1994) Database of p53 gene somatic mutations in human tumors and cell lines. Nucleic Acids Res., 22, 3551–3555.

  2. Bullock,A.N., Henckel,J., DeDecker,B.S., Johnson,C.M., Nikolova,P.V., Proctor,M.R., Lane,D.P. and Fersht,A.R. (1997) Thermodynamic stability of wild-type and mutant p53 core domain. Proc. Natl Acad. Sci. USA, 94, 14338–14342.[Abstract/Free Full Text]

  3. Fisher,D.E. (1994) Apoptosis in cancer therapy: crossing the threshold. Cell, 78, 539–542.[Web of Science][Medline]

  4. May,P. and May,E. (1999) Twenty years of p53 research: structural and functional aspects of the p53 protein. Oncogene, 18, 7621–7636.[Web of Science][Medline]

  5. Ko,L.J. and Prives,C. (1996) p53: puzzle and paradigm. Genes Dev., 10, 1054–1072.[Free Full Text]

  6. Levine,A. (1997) p53, the cellular gatekeeper for growth and division. Cell, 88, 323–331.[Web of Science][Medline]

  7. Cho,Y., Gorina,S., Jeffrey,P.D. and Pavletich,N.P. (1994) Crystal structure of a p53 tumor suppressor-DNA complex: understanding tumorigenic mutations. Science, 265, 346–355.[Abstract/Free Full Text]

  8. Lee,W., Harvey,T.S., Yin,Y., Yau,P., Litchfield,D. and Arrowsmith,C.H. (1994) Solution structure of the tetrameric minimum transforming domain of p53. Nature Struct. Biol., 1, 877–890.[Web of Science][Medline]

  9. Clore,G.M., Omichinski,J.G., Sakaguchi,K., Zambrano,N., Sakamoto,H., Apella,E. and Gronenborn,A.M. (1995) Interhelical angles in the solution structure of the oligomerization domain of p53. Science, 267, 1515–1516.[Free Full Text]

  10. Jeffrey,P.D., Gorina,S. and Pavletich,N.P. (1995) Crystal structure of the tetramerization domain of the p53 tumor suppressor at 1.7 Angstroms. Science, 267, 1498–1502.[Abstract/Free Full Text]

  11. Wang,Y. and Prives,C. (1995) Increased and altered DNA binding of human p53 by S and G2/M but not G1 cyclin-dependent kinases. Nature, 376, 88–91.[Medline]

  12. Hollstein,M., Sidransky,D., Vogelstein,B. and Harris,C.C. (1991) p53 mutations in human cancers. Science, 253, 49–53.[Abstract/Free Full Text]

  13. Ory,K., Legros,Y., Auguin,C. and Soussi,T. (1994) Analysis of the most representative tumour-derived p53 mutants reveals that changes in the protein conformation are not correlated with loss of transactivation or inhibition of cell prolifereation. EMBO J., 13, 3496–3504.[Web of Science][Medline]

  14. Bullock,A.N., Henckel,J. and Fersht,A.R. (2000) Quantitative analysis of residual folding and DNA binding in mutant p53 core domain: definition of mutant states for rescue in cancer therapy. Oncogene, 19, 1245–1256.[Web of Science][Medline]

  15. Bargonetti,J., Manfredi,J., Chen,X., Marshak,D.R. and Prives,C. (1993) A proteolytic fragment from the central region of p53 has marked sequence-specific DNA-binding activity when generated from wild-type but not from oncogenic mutant p53 protein. Genes Dev., 7, 2565–2574.[Abstract/Free Full Text]

  16. Wieczorek,A.M., Waterman,J.L., Waterman,M.J. and Halazonetis,T.D. (1996) Structure-based rescue of common tumor-derived p53 mutants. Nature Med., 2, 1143–1146.[Web of Science][Medline]

  17. Brooks,B.R., Bruccoleri,R.E., Olafson,B.D., States,D.J., Swaminathan,S. and Karplus,M. (1983) CHARMM: A program for macromolecular energy, minimization and dynamics calculations. J. Comp. Chem., 4, 187–217.

  18. Dunbrack,R.L.J. and Karplus,M. (1994) Conformational analyses of the backbone-dependent rotamer preferences of protein sidechains. Nature Struct. Biol., 1, 334–340.[Web of Science][Medline]

  19. Bower,M.J., Cohen,F.E. and Dunbrack,R.L. (1997) Prediction of protein side-chain rotamers from a backbone-dependent rotamer library: a new homology modeling tool. J. Mol. Biol., 267, 1268–1282.[Web of Science][Medline]

  20. Wright,J.D. and Lim,C. (2001) A fast method for predicting amino acid mutations that lead to unfolding. Protein Eng., 14, 479–486.[Abstract/Free Full Text]

  21. Mackerell,A.D.J. and Banavali,N. (2000) All-atom empirical force field for nucleic acids: II. application to molecular dynamics simulations of DNA and RNA in solution. J. Comp. Chem., 21, 105–120.

  22. Jorgensen,W.L., Chandrasekhar,J., Madura,J.D., Impey,R.W. and Klein,M.L. (1983) Comparison of simple potentials for simulating liquid water. J. Chem. Phys., 79, 923–926.

  23. Berkowitz,M. and McCammon,J.A. (1982) Molecular dynamics with stochastic boundary conditions. Chem. Phys. Lett., 90, 215–217.

  24. Eriksson,M.A.L., Hard,T. and Nilsson,L. (1995) Molecular dynamics simulations of the glucocorticoid receptor DNA-binding domain in complex with DNA and free in solution. Biophys. J., 68, 402–426.[Web of Science][Medline]

  25. Sharp,K. (1996) Electrostatic interactions in hirudin-thrombin binding. Biophys. Chem., 61, 37–49.[Web of Science][Medline]

  26. Noskov,S. and Lim,C. (2001) Free energy decomposition of protein-protein interactions. Biophys. J., 81, 737–750.[Web of Science][Medline]

  27. Nifosi,R., Reyes,C.M. and Kollman,P.A. (2000) Molecular dynamics studies of the HIV-1 TAR and its complex with argininamide. Nucleic Acids Res., 28, 4944–4955.[Abstract/Free Full Text]

  28. Wang,W., Lim,W.A., Jakalian,A., Wang,J., Wang,J., Luo,R., Bayly,C.I. and Kollman,P.A. (2001) An analysis of the interactions between the Sem-5 SH3 domain and its ligands using molecular dynamics, free energy calculations and sequence analysis. J. Am. Chem. Soc., 123, 3986–3994.[Web of Science][Medline]

  29. Archontis,G., Simonson,T. and Karplus,M. (2001) Binding free energies and free energy components from molecular dynamics and Poisson-Boltzmann calculations. Application to amino acid recognition by aspartyl-tRNA synthetase. J. Mol. Biol., 306, 307–327.[Web of Science][Medline]

  30. MacKerell,J.A.D., Bashford,D., Bellott,M., Dunbrack,R., Evanseck,J.D., Field,M.J., Fischer,S., Gao,J., Guo,H., Ha,S. et al. (1998) All-hydrogen empirical potential for molecular modeling and dynamics studies of proteins using the CHARMM22 force field. J. Phys. Chem. B., 102, 3586–3616.

  31. Lee,B. and Richards,F.M. (1971) The interpretation of protein structures: estimation of static accessibility. J. Mol. Biol., 55, 379–400.[Web of Science][Medline]

  32. Jayaram,B., McConnell,K.J., Surjit,B.D. and Beveridge,D.L. (1999) Free energy analysis of protein-DNA binding: the EcoRI endonuclease-DNA complex. J. Comp. Phys., 151, 333–357.

  33. Sharp,K.A., Nicholls,A., Fine,R.F. and Honig,B. (1991) Reconciling the magnitude of the macroscopic and microscopic hydrophobic effects. Science, 252, 106–109.[Abstract/Free Full Text]

  34. Pascual-Ahuir,J.L. and Silla,E. (1990) GEPOL: An improved description of molecular surfaces. I. Building the spherical surface set. J. Comp. Chem., 11, 1047–1060.

  35. Gilson,M.K. and Honig,B.H. (1988) Calculation of the electrostatic potential in solution: method and error assessment. J. Comp. Chem., 9, 327–335.

  36. Chandrasekhar,I., Clore,G.M., Szabo,A., Gronenborn,A.M. and Brooks,B.R. (1992) A 500 ps molecular dynamics simulation study of interleukin-1ß in water. J. Mol. Biol., 226, 239–250.[Web of Science][Medline]

  37. Daggett,V. and Levitt,M. (1993) Protein unfolding pathways explored through molecular dynamics simulations. J. Mol. Biol., 232, 600–619.[Web of Science][Medline]

  38. Smith,P.E., van Schaik,R.C., Szyperski,T., Wuthrich,K. and van Gunsteren,W.F. (1995) Internal mobility of the basic trypsin inhibitor in solution: a comparison of NMR spin relaxation measurements and molecular dynamics simulations. J. Mol. Biol., 246, 356–365.[Web of Science][Medline]

  39. Philippopoulos,M. and Lim,C. (1995) Molecular dynamics simulation of E. coli ribonuclease H1 in solution: correlation with NMR and X-ray data and insights into biological function. J. Mol. Biol., 254, 771–792.[Web of Science][Medline]

  40. Petsko,G. and Ringe,D. (1984) Fluctuations in protein structure from X-ray diffraction. Annu. Rev. Biophys. Bioeng., 13, 331–371.[Web of Science][Medline]

  41. Post,C.B., Dobson,C.M. and Karplus,M. (1989) A molecular dynamics analysis of protein structural elements. Proteins, 5, 337–354.[Web of Science][Medline]

  42. Wang,P.L., Sait,F. and Winter,G. (2001) The ‘wildtype’ conformation of p53: epitope mapping using hybrid proteins. Oncogene, 20, 2318–2324.[Web of Science][Medline]

  43. Gannon,J.V., Greaves,R. and Lane,D.P. (1990) Activating mutations in p53 produce a common conformational effect. A monoclonal antibody specific for the mutant form. EMBO J., 9, 1595–1602.[Web of Science][Medline]

  44. Wong,K.-B., DeDecker,B.S., Freund,M.V., Proctor,M.R., Bycroft,M. and Fersht,A.R. (1999) Hot-spot mutants of p53 core domain evince characteristic local structural chnages. Proc. Natl Acad. Sci. USA, 96, 8438–8442.[Abstract/Free Full Text]

  45. Pietenpol,J.A., Tokino,T., Thiagalingam,S., El-Deiry,W.S., Kinzler,K.W. and Vogelstein,B. (1994) Sequence-specific transcriptional activation is essential for growth suppression by p53. Proc. Natl Acad. Sci. USA, 91, 1998–2002.[Abstract/Free Full Text]

  46. Friedlander,P., Legros,Y., Soussi,T. and Prives,C. (1996) Regulation of mutant p53 temperature-sensitive DNA binding. J. Biol. Chem., 271, 25468–25478.[Abstract/Free Full Text]

  47. El-Deiry,A.A., Kern,S.E., Pietenpol,J.A., Kinzler,K.W. and Vogelstein,B. (1992) Definition of a consensus binding site for p53. Nature Genet., 1, 45–49.[Web of Science][Medline]

  48. Trulson,J.A. and Millhauser,G.L. (1999) The effects of mutations on peptide models of the DNA binding helix of p53: evidence for a correlation between structure and tumorigenesis. Biopolymers, 49, 215–224.[Web of Science][Medline]

  49. Rolley,N., Butcher,S. and Milner,J. (1995) Specific DNA binding by different classes of human p53 mutants. Oncogene, 11, 763–770.[Web of Science][Medline]

  50. Zhang,W., Funk,W.D., Wright,W.E., Shay,J.W. and Deisseroth,A.B. (1993) Novel DNA binding of p53 mutants and their role in transcription activation. Oncogene, 8, 2555–2559.[Web of Science][Medline]

  51. Park,J.D., Nakamura,H., Chumakov,A.M., Said,J.W., Miller,C.W., Chen,D.L. and Koeffler,H.P. (1994) Transactivation and DNA binding abilities of endogenous p53 mutant cell lines. Oncogene, 9, 1899–1906.[Web of Science][Medline]

  52. Wright,J.D. and Lim,C. (1998) Prediction of an anti-IgE binding site on IgE. Protein Eng., 11, 421–427.[Abstract/Free Full Text]

  53. Koradi,R., Billeter,M. and Wuthrich,K. (1996) MOLMOL: a program for display and analysis of macromoleular structures. J. Mol. Graphics, 14, 51–55.[Web of Science][Medline]

  54. Gromiha,M.M., Oobatake,M., Kono,H., Uedaira,H. and Sarai,A. (1999) Role of the structural and sequence information in the prediction of protein stability changes: comparison between buried and partially buried mutations. Protein Eng., 12, 549–555.[Abstract/Free Full Text]


Add to CiteULike CiteULike   Add to Connotea Connotea   Add to Del.icio.us Del.icio.us    What's this?


This article has been cited by other articles:


Home page
 Nucleic Acids ResHome page
B. Ma and A. J. Levine
Probing potential binding modes of the p53 tetramer to DNA based on the symmetries encoded in p53 response elements
Nucleic Acids Res., December 3, 2007; 35(22): 7733 - 7747.
[Abstract] [Full Text] [PDF]

Home page
 Proc. Natl. Acad. Sci. USAHome page
B. Ma, Y. Pan, K. Gunasekaran, R. B. Venkataraghavan, A. J. Levine, and R. Nussinov
Comparison of the protein-protein interfaces in the p53-DNA crystal structures: Towards elucidation of the biological interface
PNAS, March 15, 2005; 102(11): 3988 - 3993.
[Abstract] [Full Text] [PDF]

This Article
Right arrow Abstract Freely available
Right arrow Print PDF (394K) Freely available
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in ISI Web of Science
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Add to My Personal Archive
Right arrow Download to citation manager
Right arrow Search for citing articles in:
ISI Web of Science (9)
Right arrowRequest Permissions
Right arrow Commercial Re-use Guidelines
for Open Access NAR Content
Google Scholar
Right arrow Articles by Wright, J. D.
Right arrow Articles by Lim, C.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Wright, J. D.
Right arrow Articles by Lim, C.
Social Bookmarking
Add to CiteULike   Add to Connotea   Add to Del.icio.us  
What's this?