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.     

Major groove recognition by three-stranded beta-sheets: affinity determinants and conserved structural features

We present the results of a rational mutagenesis and binding-affinity study of the three-stranded beta-sheet-DNA interface in the complex formed by the amino-terminal DNA-binding domain of the Tn916 integrase protein and its cognate binding site. The relative importance of interfacial contacts present in its NMR-derived solution structure have been tested through mutagenesis, fluorescence anisotropy, and intrinsic quenching DNA-binding assays. We find that seven protein-DNA hydrogen bonds (two base-specific and five to phosphate groups) significantly contribute to the level of affinity. These interactions span the entire DNA-binding surface on the protein, but primarily originate from residues in only two strands of the sheet and loop L2. Interestingly, we show that highly populated, precisely defined intermolecular hydrogen bonds in the ensemble of conformers are invariably important for DNA-binding, implying that NMR-derived solution structures provide direct insight into the energetics of recognition. Unusual three-stranded beta-sheet-DNA interfaces have recently been discovered in three unrelated protein-DNA complexes. A comparative analysis of these structures reveals similar sheet positioning, the presence of two invariant interfacial contacts to the phosphodiester backbone, and two semi-conserved base-specific hydrogen bonds. Two of these conserved contacts significantly contribute to the affinity of the integrase-DNA complex, suggesting that the three-stranded beta-sheet DNA-binding motif exhibits conserved principles of recognition.     

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

The DNA binding domain of the transposon Tn916 integrase (INT-DBD) binds to its DNA target site by positioning the face of a three-stranded antiparallel beta-sheet within the major groove. Binding of INT-DBD to a 13 base pair duplex DNA target site was studied by isothermal titration calorimetry, differential scanning calorimetry, thermal melting followed by circular dichroism spectroscopy, and fluorescence spectroscopy. The observed heat capacity change accompanying the association reaction (DeltaC(p)) is temperature-dependent, decreasing from -1.4 kJ K(-1) mol(-1) at 4 degrees C to -2.9 kJ K(-1) mol(-1) at 30 degrees C. The reason is that the partial molar heat capacities of the free protein, the free DNA duplex, and the protein-DNA complex are not changing in parallel when the temperature increases and that thermal motions of the protein and the DNA are restricted in the complex. After correction for this effect, DeltaC(p) is -1.8 kJ K(-1) mol(-1) and temperature-independent. However, this value is still higher than DeltaC(p) of -1.2 kJ K(-1) mol(-1) estimated by semiempirical methods from dehydration of surface area buried at the complex interface. We propose that the discrepancy between the measured and the structure-based prediction of binding energetics is caused by incomplete dehydration of polar groups in the complex. In support, we identify cavities at the interface that are large enough to accommodate approximately 10 water molecules. Our results highlight the difficulties of structure-based prediction of DeltaC(p) (and other thermodynamic parameters) and emphasize how important it is to consider changes of thermal motions and soft vibrational modi in protein-DNA association reactions. This requires not only a detailed investigation of the energetics of the complex but also of the folding thermodynamics of the protein and the DNA alone, which are described in the accompanying paper [Milev et al. (2003) Biochemistry 42, 3492-3502].     

Thermodynamics of sequence-specific protein-DNA interactions

The molecular recognition processes in sequence-specific protein-DNA interactions are complex. The only feature common to all sequence-specific protein-DNA structures is a large interaction interface, which displays a high degree of complementarity in terms of shape, polarity and electrostatics. Many molecular mechanisms act in concert to form the specific interface. These include conformational changes in DNA and protein, dehydration of surfaces, reorganization of ion atmospheres, and changes in dynamics. Here we review the current understanding of how different mechanisms contribute to the thermodynamics of the binding equilibrium and the stabilizing effect of the different types of noncovalent interactions found in protein-DNA complexes. The relation to the thermodynamics of small molecule-DNA binding and protein folding is also briefly discussed.     

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.     

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.     

The role of slow and fast protein motions in allosteric interactions

Allostery is fundamentally thermodynamic in nature. Long-range communication in proteins may be mediated not only by changes in the mean conformation with enthalpic contribution but also by changes in dynamic fluctuations with entropic contribution. The important role of protein motions in mediating allosteric interactions has been established by NMR spectroscopy. By using CAP as a model system, we have shown how changes in protein structure and internal dynamics can allosterically regulate protein function and activity. The results indicate that changes in conformational entropy can give rise to binding enhancement, binding inhibition, or have no effect in the expected affinity, depending on the magnitude and sign of enthalpy–entropy compensation. Moreover, allosteric interactions can be regulated by the modulation a low-populated conformation states that serve as on-pathway intermediates for ligand binding. Taken together, the interplay between fast internal motions, which are intimately related to conformational entropy, and slow internal motions, which are related to poorly populated conformational states, can regulate protein activity in a way that cannot be predicted on the basis of the protein’s ground-state structure.     

Some (dis)assembly required: partial unfolding in the Par-6 allosteric switch

Allostery is commonly described as a functional connection between two distant sites in a protein, where a binding event at one site alters affinity at the other. Here, we review the conformational dynamics that encode an allosteric switch in the PDZ domain of Par-6, which is a scaffold protein that organizes other proteins into a complex required to initiate and maintain cell polarity. NMR measurements revealed that the PDZ domain samples an evolutionarily conserved unfolding intermediate allowing rearrangement of two adjacent loop residues that control ligand binding affinity. Cdc42 binding to Par-6 creates a novel interface between the PDZ domain and the adjoining CRIB motif that stabilizes the high-affinity PDZ conformation. Thermodynamic and kinetic studies suggest that partial PDZ unfolding is an integral part of the Par-6 switching mechanism. The Par-6 CRIB-PDZ module illustrates two important structural aspects of protein evolution: the interface between adjacent domains in the same protein can give rise to allosteric regulation, and thermodynamic stability may be sacrificed to increase the sampling frequency of an unfolding intermediate required for conformational switching.     

Molecular dynamics simulation of the RNA complex of a double-stranded RNA-binding domain reveals dynamic features of the intermolecular interface and its hydration

The interaction between double-stranded RNA (dsRNA) and the third double-stranded domain (dsRBD) from Drosophila Staufen protein represents a paradigm to understand how the dsRBD protein family, one of the most common RNA-binding protein units, binds dsRNA. The nuclear magnetic resonance (NMR) structure of this complex and the x-ray structure of another family member revealed the stereochemical basis for recognition, but also raised new questions. Although the crystallographic studies revealed a highly ordered interface containing numerous water-mediated contacts, NMR suggested extensive residual motion at the interface. To address how interfacial motion contributes to molecular recognition in the dsRBD-dsRNA system, we conducted a 2-ns molecular dynamics simulation of the complex derived from Staufen protein and of the separate protein and RNA components. The results support the observation that a high degree of conformational flexibility is retained upon complex formation and that this involves interfacial residues that are critical for dsRBD-dsRNA binding. The structural origin of this residual flexibility is revealed by the analysis of the trajectory of motion. Individual basic side chains switch continuously from one RNA polar group to another with a residence time seldom exceeding 100 ps, while retaining favorable interaction with RNA throughout much of the simulation. Short-lived water molecules mediate some of these interactions for a large fraction of the trajectory studied here. This result indicates that water molecules are not statically associated with the interface, but continuously exchange with the bulk solvent on a 1-10-ps time scale. This work provides new insight into dsRBD-dsRNA recognition and builds upon a growing body of evidence, suggesting that short-lived dynamic interactions play important roles in protein-nucleic acid interactions.     

Entropy Hotspots for the Binding of Intrinsically Disordered Ligands to a Receptor Domain

Proline-rich motifs (PRMs) are widely used for mediating protein-protein interactions with weak binding affinities. Because they are intrinsically disordered when unbound, conformational entropy plays a significant role for the binding. However, residue-level differences of the entropic contribution in the binding of different ligands remain not well understood. We use all-atom molecular dynamics simulation and the maximal information spanning tree formalism to analyze conformational entropy associated with the binding of two PRMs, one from the Abl kinase and the other from the nonstructural protein 1 of the 1918 Spanish influenza A virus, to the N-terminal SH3 (nSH3) domain of the CrkII protein. Side chains of the stably folded nSH3 experience more entropy change upon ligand binding than the backbone, whereas PRMs involve comparable but heterogeneous entropy changes among the backbone and side chains. In nSH3, two conserved nonpolar residues forming contacts with the PRM experience the largest side-chain entropy loss. In contrast, the C-terminal charged residues of PRMs that form polar contacts with nSH3 experience the greatest side-chain entropy loss, although their “fuzzy” nature is attributable to the backbone that remains relatively flexible. Thus, residues that form high-occupancy contacts between nSH3 and PRM do not reciprocally contribute to entropy loss. Furthermore, certain surface residues of nSH3 distal to the interface with PRMs gain entropy, indicating a nonlocal effect of ligand binding. Comparing between the PRMs from cAbl and nonstructural protein 1, the latter involves a larger side-chain entropy loss and forms more contacts with nSH3. Consistent with experiments, this indicates stronger binding of the viral ligand at the expense of losing the flexibility of side chains, whereas the backbone experiences less entropy loss. The entropy “hotspots” as identified in this study will be important for tuning the binding affinity of various ligands to a receptor.     

A negative electrostatic determinant mediates the association between the Escherichia coli trp repressor and its operator DNA.

The electrostatic potential surfaces were characterized for trp repressor models that bind to DNA with sequence specificity, without specificity, and not at all. Comparisons among the surfaces were used to isolate protein surface features likely to be important in DNA binding. Models that differ in protein conformation and tryptophan-analogue binding consistently showed positive potential associated with the protein surfaces that interact with the DNA major groove. However, negative potential is associated with the trp repressor surface that contacts the DNA minor groove. This negative potential is significantly neutralized in the protein conformation that is bound to DNA. Positive potential is also associated with the tryptophan binding-site surface, a consequence of the tryptophan- or tryptophan analogue-induced allosteric change. This protein region is complementary to the strongest negative potential associated with the DNA phosphate backbone and is also present in the isolated protein structure from the protein-DNA complex. The effects of charge-change mutation, pH dependence, and salt dependence on the electrostatic potential surfaces were also examined with regard to their effects on protein-DNA binding constants. A consistent model is formed that defines a role for long-range electrostatics early in the protein-DNA association process and complements previous structural, molecular association, and mutagenesis studies.     

Consecutive steps of phosphorylation affect conformation and DNA binding of the chironomus high mobility group A protein

The high mobility group (HMG) proteins of the AT-hook family (HMGA) lie downstream in regulatory networks with protein kinase C, Cdc2 kinase, MAP kinase, and casein kinase 2 (CK2) as final effectors. In the cells of the midge Chironomus, almost all of the HMGA protein (cHMGA) is phosphorylated by CK2 at two adjacent sites. 40% of the protein population is additionally modified by MAP kinase. Using spectroscopic and protein footprinting techniques, we analyzed how individual and consecutive steps of phosphorylation change the conformation of an HMGA protein and affect its contacts with poly(dA-dT).poly(dA-dT) and a fragment of the interferon-beta promoter. We demonstrate that phosphorylation of cHMGA by CK2 alters its conformation and modulates its DNA binding properties such that a subsequent phosphorylation by Cdc2 kinase changes the organization of the protein-DNA complex. In contrast, consecutive phosphorylation by MAP kinase, which results in a dramatic change in cHMGA conformation, has no direct effect on the complex. Because the phosphorylation of the HMGA proteins attenuates binding affinity and reduces the extent of contacts between the DNA and protein, it is likely that this process mirrors the dynamics and diversity of regulatory processes in chromatin.     

Solution structure of ileal lipid binding protein in complex with glycocholate.

Ileal lipid binding protein (ILBP) is a cytosolic lipid-binding protein that binds both bile acids and fatty acids. We have determined the solution structure of porcine ILBP in complex with glycocholate by homonuclear and heteronuclear two-dimensional NMR spectroscopy. The conformation of the protein-ligand complex was determined by restrained energy minimization and simulated annealing calculations after docking the glycocholate ligand into the protein structure. The overall tertiary structure of ILBP is highly analogous to the three-dimensional structures of several other intracellular lipid binding proteins (LBPs). Like the apo-structure, the bile-acid complex of ILBP is composed of 10 anti-parallel beta-strands that form a water-filled clam-shell structure, and two short alpha-helices. Chemical shift data indicated that the bile acid ligand is bound inside the protein cavity. Furthermore, 13C-edited heteronuclear single-quantum correlation-NOESY experiments showed NOE contacts between several aromatic residues located in the proposed bile acid portal region and the 13C-labeled ligand. A single bile acid molecule is bound inside the protein, with the steroid moiety penetrating deep into the water-accessible internal cavity, such that ring A is located right above the plane of the Trp49 indole ring. The carboxylate tail of the ligand is protruding from the proposed bile acid portal into the surrounding aqueous solution. The body of the steroid moiety is oriented with the nonpolar face in contact with the mostly hydrophobic residues of beta-strands C, D and E, while the polar face shows contacts with the side-chains of Tyr97, His99, Glu110 and Arg121 in beta-strands H, I and J. Thus, the conformational arrangement of the ligand complex suggests that the binding affinity of ILBP for bile acid molecules is based mainly on strong hydrophobic interactions inside the protein cavity. Furthermore, this binding mode explains how ILBP can transport unconjugated and conjugated bile acids.