Structural basis of DNA recognition by p53 tetramers

The tumor-suppressor protein p53 is among the most effective of the cell's natural defenses against cancer. In response to cellular stress, p53 binds as a tetramer to diverse DNA targets containing two decameric half-sites, thereby activating the expression of genes involved in cell-cycle arrest or apoptosis. Here we present high-resolution crystal structures of sequence-specific complexes between the core domain of human p53 and different DNA half-sites. In all structures, four p53 molecules self-assemble on two DNA half-sites to form a tetramer that is a dimer of dimers, stabilized by protein-protein and base-stacking interactions. The protein-DNA interface varies as a function of the specific base sequence in correlation with the measured binding affinities of the complexes. The new data establish a structural framework for understanding the mechanisms of specificity, affinity, and cooperativity of DNA binding by p53 and suggest a model for its regulation by regions outside the sequence-specific DNA binding domain.     

The role of flexibility and hydration on the sequence-specific DNA recognition by the Tn916 integrase protein: a molecular dynamics analysis

The N-terminal domain of the Tn916 integrase protein (INT-DBD) is responsible for DNA binding in the process of strand cleavage and joining reactions required for transposition of the Tn916 conjugative transposon. Site-specific association is facilitated by numerous protein-DNA contacts from the face of a three-stranded beta-sheet inserted into the major groove. The protein undergoes a subtle conformational transition and is slightly unfolded in the protein-DNA complex. The conformation of many charged residues is poorly defined by NMR data but mutational studies have indicated that removal of polar side chains decreases binding affinity, while non-polar contacts are malleable. Based on analysis of the binding enthalpy and binding heat capacity, we have reasoned that dehydration of the protein-DNA interface is incomplete. This study presents results from a molecular dynamics investigation of the INT-DBD-DNA complex aimed at a more detailed understanding of the role of conformational dynamics and hydration in site-specific binding. Comparison of simulations (total of 13 ns) of the free protein and of the bound protein conformation (in isolation or DNA-bound) reveals intrinsic flexibility in certain parts of the molecule. Conformational adaptation linked to partial unfolding appears to be induced by protein-DNA contacts. The protein-DNA hydrogen-bonding network is highly dynamic. The simulation identifies protein-DNA interactions that are poorly resolved or only surmised from the NMR ensemble. Single water molecules and water clusters dynamically optimize the complementarity of polar interactions at the 'wet' protein-DNA interface. The simulation results are useful to establish a qualitative link between experimental data on individual residue's contribution to binding affinity and thermodynamic properties of INT-DBD alone and in complex with DNA.     

Tethered-Hopping Model for Protein-DNA Binding and Unbinding Based on Sox2-Oct1-Hoxb1 Ternary Complex Simulations

The sliding and hopping models encapsulate the essential protein-DNA binding process for binary complex formation and dissociation. However, the effects of a cofactor protein on the protein-DNA binding process that leads to the formation of a ternary complex remain largely unknown. Here we investigate the effect of the cofactor Sox2 on the binding and unbinding of Oct1 with the Hoxb1 control element. We simulate the association of Oct1 with Sox2-Hoxb1 using molecular dynamics simulations, and the dissociation of Oct1 from Sox2-Hoxb1 using steered molecular dynamics simulations, in analogy to a hopping event of Oct1. We compare the kinetic and thermodynamic properties of three model complexes (the wild-type and two mutants) in which the Oct1-DNA base-specific interactions or the Sox2-Oct1 protein-protein interactions are largely abolished. We find that Oct1-DNA base-specific interactions contribute significantly to the total interaction energy of the ternary complex, and that nonspecific Oct1-DNA interactions are sufficient for driving the formation of the protein-DNA interface. The Sox2-Oct1 protein-protein binding interface is largely hydrophobic, with remarkable shape complementarity. This interface promotes the formation of the ternary complex and slows the dissociation of Oct1 from its DNA-binding site. We propose a simple two-step reaction model of protein-DNA binding, called the tethered-hopping model, that explains the importance of the cofactor Sox2 and may apply to similar ternary protein-DNA complexes.     

Applying 6-Methylisoxanthopterin-Enhanced Fluorescence To Examine Protein-DNA Interactions in the Picomolar Range

Incorporation of fluorescent nucleoside analogues into duplex DNA usually leads to a reduction in quantum yield, which significantly limits their potential use and application. We have identified two pentamer DNA sequences containing 6-methylisoxanthopterin (6-MI) (ATFAA and AAFTA, where F is 6-MI) that exhibit significant enhancement of fluorescence upon formation of duplex DNA with quantum yields close to that of monomeric 6-MI. The enhanced fluorescence dramatically increases the utility and sensitivity of the probe and is used to study protein-DNA interactions of nanomolar specificity in this work. The increased sensitivity of 6-MI allows anisotropy binding measurements to be performed at DNA concentrations of 1 nM and fluorescence intensity measurements at 50 pM DNA. The ATFAA sequence was incorporated into DNA constructs to measure the binding affinity of four different protein-DNA interactions that exhibit sequence-specific and non-sequence-specific recognition. In all cases, the K(d) values obtained were consistent with previously reported values measured by other methods. Time-resolved and steady-state fluorescence measurements demonstrate that 6-MI fluorescence is very sensitive to local distortion and reports on different degrees of protein-induced perturbations with single-base resolution, where the largest changes occur at the site of protein binding.     

Rearrangement of side-chains in a Zif268 mutant highlights the complexities of zinc finger-DNA recognition.

Structural and biochemical studies of Cys(2)His(2) zinc finger proteins initially led several groups to propose a "recognition code" involving a simple set of rules relating key amino acid residues in the zinc finger protein to bases in its DNA site. One recent study from our group, involving geometric analysis of protein-DNA interactions, has discussed limitations of this idea and has shown how the spatial relationship between the polypeptide backbone and the DNA helps to determine what contacts are possible at any given position in a protein-DNA complex. Here we report a study of a zinc finger variant that highlights yet another source of complexity inherent in protein-DNA recognition. In particular, we find that mutations can cause key side-chains to rearrange at the protein-DNA interface without fundamental changes in the spatial relationship between the polypeptide backbone and the DNA. This is clear from a simple analysis of the binding site preferences and co-crystal structures for the Asp20-->Ala point mutant of Zif268. This point mutation in finger one changes the specificity of the protein from GCG TGG GCG to GCG TGG GC(G/T), and we have solved crystal structures of the D20A mutant bound to both types of sites. The structure of the D20A mutant bound to the GCG site reveals that contacts from key residues in the recognition helix are coupled in complex ways. The structure of the complex with the GCT site also shows an important new water molecule at the protein-DNA interface. These side-chain/side-chain interactions, and resultant changes in hydration at the interface, affect binding specificity in ways that cannot be predicted either from a simple recognition code or from analysis of spatial relationships at the protein-DNA interface. Accurate computer modeling of protein-DNA interfaces remains a challenging problem and will require systematic strategies for modeling side-chain rearrangements and change in hydration.     

DNA binding of a non-sequence-specific HMG-D protein is entropy driven with a substantial non-electrostatic contribution

The thermal properties of two forms of the Drosophila melanogaster HMG-D protein, with and without its highly basic 26 residue C-terminal tail (D100 and D74) and the thermodynamics of their non-sequence-specific interaction with linear DNA duplexes were studied using scanning and titration microcalorimetry, spectropolarimetry, fluorescence anisotropy and FRET techniques at different temperatures and salt concentrations. It was shown that the C-terminal tail of D100 is unfolded at all temperatures, whilst the state of the globular part depends on temperature in a rather complex way, being completely folded only at temperatures close to 0 degrees C and unfolding with significant heat absorption at temperatures below those of the gross denaturational changes. The association constant and thus Gibbs energy of binding for D100 is much greater than for D74 but the enthalpies of their association are similar and are large and positive, i.e. DNA binding is a completely entropy-driven process. The positive entropy of association is due to release of counterions and dehydration upon forming the protein/DNA complex. Ionic strength variation showed that electrostatic interactions play an important but not exclusive role in the DNA binding of the globular part of this non-sequence-specific protein, whilst binding of the positively charged C-terminal tail of D100 is almost completely electrostatic in origin. This interaction with the negative charges of the DNA phosphate groups significantly enhances the DNA bending. An important feature of the non-sequence-specific association of these HMG boxes with DNA is that the binding enthalpy is significantly more positive than for the sequence-specific association of the HMG box from Sox-5, despite the fact that these proteins bend the DNA duplex to a similar extent. This difference shows that the enthalpy of dehydration of apolar groups at the HMG-D/DNA interface is not fully compensated by the energy of van der Waals interactions between these groups, i.e. the packing density at the interface must be lower than for the sequence-specific Sox-5 HMG box.     

DNA curving and bending in protein–DNA recognition

Most biological events are regulated at the molecular level by site-specific associations between specialized proteins and DNA. These associations may bring distal regions of the genome into functional contact or may lead to the formation of large multisubunit complexes capable of regulating highly site-specific transactional events. It is now believed that sequence-specific protein-DNA recognition and the ability of certain proteins to compete for multiple binding sites is regulated at several levels by the local structure and conformation of the binding partners. These encompass the microstructure of DNA, including its curvature, bending and flexing as well as conformational lability in the DNA-binding domains of the proteins. Possible mechanisms for binding specificity are discussed in the context of specific nucleoprotein systems with particular emphasis given to the roles of DNA conformations in these interactions.     

Enthalpic and entropic effects of salt and polyol osmolytes on site-specific protein-DNA association: the integrase Tn916-DNA complex

The effect of low molecular-weight compounds on the equilibrium constant K(A) can be used to explore the energetics and molecular mechanism of protein-DNA interactions. Here we use the complex composed of the integrase Tn916 DNA-binding domain and its target DNA duplex to investigate the effects of salt and the nonionic osmolytes glycerol and sorbitol on sequence-specific protein-DNA association. Increasing Na(+) concentration from 0.12 to 0.32 M weakens the binding affinity by a factor of 20. The decrease of affinity is dominated by a large loss of binding enthalpy but only a small loss of binding entropy. This contrasts the concept that the salt-induced weakening of protein-DNA binding is mainly entropic. The large enthalpy loss is discussed in the light of recent views about the nature of the general salt effect. Addition of up to 2.5 M sorbitol and up to 3.3 M glycerol causes a slight increase of the binding affinity. However, both osmolytes lead to a large enthalpy gain and a similarly large entropy loss. This intriguing enthalpy-entropy compensation can be explained in part by an enthalpic chelate effect: The osmolyte tightens the structure of the protein-DNA complex whereby the formation of enthalpically favorable noncovalent interactions is promoted at the entropic cost of a more rigid complex. The results were obtained by isothermal titration calorimetry. They are supported by kinetic experiments showing that the rate of formation of the complex is reduced by salt, but the rate of complex dissociation is not. Glycerol and sorbitol reduce both rates in line with an only small effect on complex stability. This work clarifies the thermodynamic and kinetic response of a novel protein-DNA complex to increased salt and the presence of two common, nonionic osmolytes.     

Insights into protein-DNA interactions through structure network analysis

Protein-DNA interactions are crucial for many cellular processes. Now with the increased availability of structures of protein-DNA complexes, gaining deeper insights into the nature of protein-DNA interactions has become possible. Earlier, investigations have characterized the interface properties by considering pairwise interactions. However, the information communicated along the interfaces is rarely a pairwise phenomenon, and we feel that a global picture can be obtained by considering a protein-DNA complex as a network of noncovalently interacting systems. Furthermore, most of the earlier investigations have been carried out from the protein point of view (protein-centric), and the present network approach aims to combine both the protein-centric and the DNA-centric points of view. Part of the study involves the development of methodology to investigate protein-DNA graphs/networks with the development of key parameters. A network representation provides a holistic view of the interacting surface and has been reported here for the first time. The second part of the study involves the analyses of these graphs in terms of clusters of interacting residues and the identification of highly connected residues (hubs) along the protein-DNA interface. A predominance of deoxyribose-amino acid clusters in beta-sheet proteins, distinction of the interface clusters in helix-turn-helix, and the zipper-type proteins would not have been possible by conventional pairwise interaction analysis. Additionally, we propose a potential classification scheme for a set of protein-DNA complexes on the basis of the protein-DNA interface clusters. This provides a general idea of how the proteins interact with the different components of DNA in different complexes. Thus, we believe that the present graph-based method provides a deeper insight into the analysis of the protein-DNA recognition mechanisms by throwing more light on the nature and the specificity of these interactions.     

Insights on the acetylated NF-κB transcription factor complex with DNA from molecular dynamics simulations

The nuclear factor-κB (NF-κB) is a DNA sequence-specific regulator of many important biological processes, whose activity is modulated by enzymatic acetylation. In one of the best functionally characterized NF-κB complexes, the p50/p65 heterodimer, acetylation of K221 at p65 causes a decrease of DNA dissociation rate, whilst the acetylation of K122 and K123, also at p65, markedly decreases the binding affinity for DNA. By means of molecular dynamics simulations based on the X-ray structure of the p50/p65 complex with DNA, we provide insights on the structural determinants of the acetylated complexes in aqueous solution. Lysine acetylation involves the loss of favorable electrostatic interactions between DNA and NF-κB, which is partially compensated by the reduction of the desolvation free-energy of the two binding partners. Acetylation at both positions K122 and K123 is associated with a decrease of the electrostatic potential at the p65/DNA interface, which is only partially counterbalanced by an increase of the local Na(+) concentration. It induces the disruption of base-specific and nonspecific interactions between DNA and NF-κB and it is consistent with the observed decrease of binding affinity. In contrast, acetylation at position K221 results in the loss of nonspecific protein-DNA interactions, but the DNA recognition sites are not affected. In addition, the loss of protein-DNA interactions is likely to be counterbalanced by an increase of the configurational entropy of the complex, which provides, at a speculative level, a justification for the observed decrease of NF-κB/DNA dissociation rate.     

Energetics of sequence-specific protein-DNA association: computational analysis of integrase Tn916 binding to its target DNA

The N-terminal domain of the bacterial integrase Tn916 specifically recognizes the 11 bp DNA target site by positioning the face of a three-stranded beta-sheet into the major groove. Binding is linked to structural adaptation. We have characterized INT-DBD binding to DNA in detail by calorimetry [Milev, S., Gorfe, A., Karshikoff, A., Clubb, R. T., Bosshard, H. R., and Jelesarov, I. (2003) Biochemistry 42, 3481-3491]. Our thermodynamic analysis has indicated that the major driving force of association is the hydrophobic effect while polar interactions contribute less. To gain more comprehensive information about the binding process, we performed a computational analysis of the binding free energy and report here the results. A hybrid molecular mechanics/continuum approach was followed. The total binding free energy is predicted with reasonable accuracy. The calculations confirm that nonpolar effects stabilize the protein-DNA complex while electrostatics opposes binding. Structural changes optimizing surface complementarity are costly in terms of energy. The energetic consequences from the replacement of nine DNA-contacting residues by alanine were investigated. The calculations correctly predict the binding affinity decrease of eight mutations and the destabilizing effect of one wild-type residue. Bulky side chains stabilize the wild-type complex through packing interactions and favorable nonpolar dehydration, but the net nonpolar energy changes do not correlate with the relative affinity loss upon mutation. Discrete protein-DNA electrostatic interactions may be net stabilizing or net destabilizing depending on the local environment. In contrast to nonpolar energy changes, the magnitude of the electrostatic free energy ranks the mutations according to the experimentally measured DeltaDeltaG. Free energy decomposition analysis from a structural perspective leads to detailed information about the thermodynamic strategy used by INT-DBD for sequence-specific DNA binding.     

Molecular architecture of protein-RNA recognition sites

The molecular architecture of protein-RNA interfaces are analyzed using a non-redundant dataset of 152 protein-RNA complexes. We find that an average protein-RNA interface is smaller than an average protein-DNA interface but larger than an average protein–protein interface. Among the different classes of protein-RNA complexes, interfaces with tRNA are the largest, while the interfaces with the single-stranded RNA are the smallest. Significantly, RNA contributes more to the interface area than its partner protein. Moreover, unlike protein–protein interfaces where the side chain contributes less to the interface area compared to the main chain, the main chain and side chain contributions flipped in protein-RNA interfaces. We find that the protein surface in contact with the RNA in protein-RNA complexes is better packed than that in contact with the DNA in protein-DNA complexes, but loosely packed than that in contact with the protein in protein–protein complexes. Shape complementarity and electrostatic potential are the two major factors that determine the specificity of the protein-RNA interaction. We find that the H-bond density at the protein-RNA interfaces is similar with that of protein-DNA interfaces but higher than the protein–protein interfaces. Unlike protein-DNA interfaces where the deoxyribose has little role in intermolecular H-bonds, due to the presence of an oxygen atom at the 2′ position, the ribose in RNA plays significant role in protein-RNA H-bonds. We find that besides H-bonds, salt bridges and stacking interactions also play significant role in stabilizing protein-nucleic acids interfaces; however, their contribution at the protein–protein interfaces is insignificant.     

Quantitative parameters for amino acid-base interaction: implications for prediction of protein-DNA binding sites.

Inspection of the amino acid-base interactions in protein-DNA complexes is essential to the understanding of specific recognition of DNA target sites by regulatory proteins. The accumulation of information on protein-DNA co-crystals challenges the derivation of quantitative parameters for amino acid-base interaction based on these data. Here we use the coordinates of 53 solved protein-DNA complexes to extract all non-homologous pairs of amino acid-base that are in close contact, including hydrogen bonds and hydrophobic interactions. By comparing the frequency distribution of the different pairs to a theoretical distribution and calculating the log odds, a quantitative measure that expresses the likelihood of interaction for each pair of amino acid-base could be extracted. A score that reflects the compatibility between a protein and its DNA target can be calculated by summing up the individual measures of the pairs of amino acid-base involved in the complex, assuming additivity in their contributions to binding. This score enables ranking of different DNA binding sites given a protein binding site and vice versa and can be used in molecular design protocols. We demonstrate its validity by comparing the predictions using this score with experimental binding results of sequence variants of zif268 zinc fingers and their DNA binding sites.