Differences between EcoRI nonspecific and "star" sequence complexes revealed by osmotic stress

The binding of the restriction endonuclease EcoRI to DNA is exceptionally specific. Even a single basepair change ("star" sequence) from the recognition sequence, GAATTC, decreases the binding free energy of EcoRI to values nearly indistinguishable from nonspecific binding. The difference in the number of waters sequestered by the protein-DNA complexes of the "star" sequences TAATTC and CAATTC and by the specific sequence complex determined from the dependence of binding free energy on water activity is also practically indistinguishable at low osmotic pressures from the 110 water molecules sequestered by nonspecific sequence complexes. Novel measurements of the dissociation rates of noncognate sequence complexes and competition equilibrium show that sequestered water can be removed from "star" sequence complexes by high osmotic pressure, but not from a nonspecific complex. By 5 Osm, the TAATTC "star" sequence complex has lost almost 90 of the approximately 110 waters initially present. It is more difficult to remove water from the CAATTC "star" sequence complex. The sequence dependence of water loss correlates with the known sequence dependence of "star" cleavage activity.     

Differences in hydration coupled to specific and nonspecific competitive binding and to specific DNA Binding of the restriction endonuclease BamHI

Using the osmotic stress technique together with a self-cleavage assay we measure directly differences in sequestered water between specific and nonspecific DNA-BamHI complexes as well as the numbers of water molecules released coupled to specific complex formation. The difference between specific and nonspecific binding free energy of the BamHI scales linearly with solute osmolal concentration for seven neutral solutes used to set water activity. The observed osmotic dependence indicates that the nonspecific DNA-BamHI complex sequesters some 120-150 more water molecules than the specific complex. The weak sensitivity of the difference in number of waters to the solute identity suggests that these waters are sterically inaccessible to solutes. This result is in close agreement with differences in the structures determined by x-ray crystallography. We demonstrate additionally that when the same solutes that were used in competition experiments are used to probe changes accompanying the binding of free BamHI to its specific DNA sequence, the measured number of water molecules released in the binding process is strikingly solute-dependent (with up to 10-fold difference between solutes). This result is expected for reactions resulting in a large change in a surface exposed area.     

Solution parameters modulating DNA binding specificity of the restriction endonuclease EcoRV

The DNA binding stringency of restriction endonucleases is crucial for their proper function. The X-ray structures of the specific and non-cognate complexes of the restriction nuclease EcoRV are considerably different suggesting significant differences in the hydration and binding free energies. Nonetheless, the majority of studies performed at pH 7.5, optimal for enzymatic activity, have found a < 10-fold difference between EcoRV binding constants to the specific and nonspecific sequences in the absence of divalent ions. We used a recently developed self-cleavage assay to measure EcoRV-DNA competitive binding and to evaluate the influence of water activity, pH and salt concentration on the binding stringency of the enzyme in the absence of divalent ions. We find the enzyme can readily distinguish specific and nonspecific sequences. The relative specific-nonspecific binding constant increases strongly with increasing neutral solute concentration and with decreasing pH. The difference in number of associated waters between specific and nonspecific DNA-EcoRV complexes is consistent with the differences in the crystal structures. Despite the large pH dependence of the sequence specificity, the osmotic pressure dependence indicates little change in structure with pH. The large osmotic pressure dependence means that measurement of protein-DNA specificity in dilute solution cannot be directly applied to binding in the crowded environment of the cell. In addition to divalent ions, water activity and pH are key parameters that strongly modulate binding specificity of EcoRV.     

Spectroscopic studies on lambda cro protein-DNA interactions

Spectroscopic (circular dichroism and fluorescence) and thermodynamic studies were conducted on lambda Cro-DNA interactions. Some base substitutions were introduced to the operator and the effects on the conformation of the complex and thermodynamic parameters for dissociation of the complex were examined. It was found that, (1) in the specific binding of Cro with DNA which has a (pseudo) consensus sequence, DNA is overwound, while in non-specific binding it is unchanged, or rather unwound; (2) substitution of central base-pairs or the introduction of a mismatched base-pair at the center of the operator reduces the extent of DNA conformational change on Cro binding and lessens the stability of the Cro-DNA complex, even though there is apparently no direct interaction between Cro and DNA at these positions; (3) stability of the complex increases with the degree of DNA conformational change of the same type during binding; (4) in some cases of specific binding, there are three states in the dissociation of the complex as observed by salt titration: two conformational states for the complex depending on salt concentration and, in non-specific binding, dissociation is a two-state transition; (5) the number of ions involved in interactions between Cro and 17 base-pair DNA is about 7.7 for NaCl titrations; (6) dissociation free energy prediction of the Cro-DNA complex by simple addition of the dissociation free energy change of a single base-pair substitution agrees with our experimental results when DNA overwinding occurs during binding, i.e. in specific binding.     

Combinatorial redesign of the DNA binding specificity of a prokaryotic helix-turn-helix repressor

Redesign of the bacteriophage 434 Cro repressor was accomplished by using an in vivo genetic screening system to identify new variants that specifically bound previously unrecognized DNA sequences. Site-directed, combinatorial mutagenesis of the 434 Cro helix-turn-helix (HTH) motif generated libraries of new variants which were screened for binding to new target sequences. Multiple mutations of 434 Cro that functionally converted wild-type (wt) 434 Cro DNA binding-sequence specificity to that of a lambda bacteriophage-specific repressor were identified. The libraries contained variations within the HTH sequence at only three positions. In vivo and in vitro analysis of several of the identified 434 Cro variants showed that the relatively few changes in the recognition helix of the HTH motif of 434 Cro resulted in specific and tight binding of the target DNA sequences. For the best 434 Cro variant identified, an apparent K(d) for lambda O(R)3 of 1 nM was observed. In competition experiments, this Cro variant was observed to be highly selective. We conclude that functional 434 Cro repressor variants with new DNA binding specificities can be generated from wt 434 Cro by mutating just the recognition helix. Important characteristics of the screening system responsible for the successful identifications are discussed. Application of the techniques presented here may allow the identification of DNA binding protein variants that functionally affect DNA regulatory sequences important in disease and industrial and biotechnological processes.     

Functional expression and affinity selection of single-chain cro by phage display: isolation of novel DNA-binding proteins

A robust selection system affording phage display of the DNA-binding helix-turn-helix protein Cro is presented. The aim of the work was to construct an experimental system allowing for the construction and isolation of Cro-derived protein with new DNA-binding properties. A derivative of the phage lambda Cro repressor, scCro8, in which the protein subunits had been covalently connected via a peptide linker was expressed in fusion with the gene 3 protein of Escherichia coli filamentous phage. The phage-displayed single-chain Cro was shown to retain the DNA binding properties of its wild-type Cro counterpart regarding DNA sequence specificity and binding affinity. A kinetic analysis revealed the rate constant of dissociation of the single-chain Cro-phage/DNA complex to be indistinguishable from that of the free single-chain Cro. Affinity selection using a biotinylated DNA with a target consensus operator sequence allowed for a 3000-fold enrichment of phages displaying single-chain Cro over control phages. The selection was based on entrapment of phage/DNA complexes formed in solution on streptavidin-coated paramagnetic beads. The expression system was subsequently used to isolate variant scCro8 proteins, mutated in their DNA-binding residues, that specifically recognized new, unnatural target DNA ligands.     

Break in the heat capacity change at 303 K for complex binding of netropsin to AATT containing hairpin DNA constructs

Studies performed in our laboratory demonstrated the formation of two thermodynamically distinct complexes on binding of netropsin to a number of hairpin-forming DNA sequences containing AATT-binding regions. These two complexes were proposed to differ only by a bridging water molecule between the drug and the DNA in the lower affinity complex. A temperature-dependent isothermal titration calorimetry (ITC)-binding study was performed using one of these constructs (a 20-mer hairpin of sequence 5'-CGAATTCGTCTCCGAATTCG) and netropsin. This study demonstrated a break in the heat capacity change for the formation of the complex containing the bridging water molecule at approximately 303 K. In the plot of the binding enthalpy change versus temperature, the slope (DeltaCp) was -0.67 kcal mol-1 K-1 steeper after the break at 303 K. Because of the relatively low melting temperature of the 20-mer hairpin (341 K (68 degrees C)), the enthalpy change for complex formation might have included some energy of refolding of the partially denatured hairpin, giving the suggestion of a larger DeltaCp. Studies done on the binding of netropsin to similar constructs, a 24-mer and a 28-mer, with added GC basepairs in the hairpin stem to increase thermal stability, exhibit the same nonlinearity in DeltaCp over the temperature range of from 275 to 333 K. The slopes (DeltaCp) were -0.69 and -0.64 kcal mol-1 K-1 steeper after 303 K for the 24-mer and 28-mer, respectively. This observation strengthens the argument regarding the presence of a bridging water molecule in the lower affinity netropsin/DNA complex. The DeltaCp data seem to infer that because the break in the heat capacity change function for the lower affinity binding occurs at the isoequilibrium temperature for water, water may be included or trapped in the complex. The fact that this break does not occur in the heat capacity change function for formation of the higher affinity complex can similarly be taken as evidence that water is not included in the higher affinity complex.     

In vivo non-specific binding of lambda CI and Cro repressors is significant

We propose a thermodynamic model that includes the non-specific binding of the lambda phage regulatory proteins CI and Cro. By fitting the model to experimental in vivo data on activities of the two promoters P(RM) and P(R) versus concentration, we estimate the free energy upon non-specific binding to be -4.1+/-0.9 kcal/mol for CI and -4.2+/-0.8 kcal/mol for Cro. For concentrations >100 nM of CI or Cro, we find that >50% of these proteins are non-specifically bound. In particular, in a lysogen (approximately 250 CI monomeric equivalents per cell) nearly 90% of CI is non-specifically bound.     

Thermodynamic analysis of the structural stability of phage 434 Cro protein

Thermodynamic parameters describing the phage 434 Cro protein have been determined by calorimetry and, independently, by far-UV circular dichroism (CD) measurements of isothermal urea denaturations and thermal denaturations at fixed urea concentrations. These equilibrium unfolding transitions are adequately described by the two-state model. The far-UV CD denaturation data yield average temperature-independent values of 0.99 +/- 0.10 kcal mol(-)(1) M(-)(1) for m and 0.98 +/- 0.05 kcal mol(-)(1) K(-)(1) for DeltaC(p)()(,U), the heat capacity change accompanying unfolding. Calorimetric data yield a temperature-independent DeltaC(p)()(,U) of 0.95 +/- 0.30 kcal mol(-)(1) K(-)(1) or a temperature-dependent value of 1.00 +/- 0.10 kcal mol(-)(1) K(-)(1) at 25 degrees C. DeltaC(p)()(,U) and m determined for 434 Cro are in accord with values predicted using known empirical correlations with structure. The free energy of unfolding is pH-dependent, and the protein is completely unfolded at pH 2.0 and 25 degrees C as judged by calorimetry or CD. The stability of 434 Cro is lower than those observed for the structurally similar N-terminal domain of the repressor of phage 434 (R1-69) or of phage lambda (lambda(6)(-)(85)), but is close to the value reported for the putative monomeric lambda Cro. Since a protein's structural stability is important in determining its intracellular stability and turnover, the stability of Cro relative to the repressor could be a key component of the regulatory circuit controlling the levels and, consequently, the functions of the two proteins in vivo.     

Electrostatic deformation of DNA by a DNA-binding protein

Complementary electrostatic interactions between negatively charged B-DNA and a positively charged array on the lambda Cro repressor protein are shown to substantially contribute to the formation energy of sequence-specific and nonspecific Cro-DNA complexes. The electrostatic interactions favor Cro binding to a bent form of DNA, a geometry which optimizes hydrogen-bonding contacts between Cro and exposed base pair groups in the DNA major groove.     

Methylene blue binding to DNA with alternating AT base sequence: minor groove binding is favored over intercalation

The results presented in this paper on methylene blue (MB) binding to DNA with AT alternating base sequence complement the data obtained in two former modeling studies of MB binding to GC alternating DNA. In the light of the large amount of experimental data for both systems, this theoretical study is focused on a detailed energetic analysis and comparison in order to understand their different behavior. Since experimental high-resolution structures of the complexes are not available, the analysis is based on energy minimized structural models of the complexes in different binding modes. For both sequences, four different intercalation structures and two models for MB binding in the minor and major groove have been proposed. Solvent electrostatic effects were included in the energetic analysis by using electrostatic continuum theory, and the dependence of MB binding on salt concentration was investigated by solving the non-linear Poisson-Boltzmann equation. We find that the relative stability of the different complexes is similar for the two sequences, in agreement with the interpretation of spectroscopic data. Subtle differences, however, are seen in energy decompositions and can be attributed to the change from symmetric 5'-YpR-3' intercalation to minor groove binding with increasing salt concentration, which is experimentally observed for the AT sequence at lower salt concentration than for the GC sequence. According to our results, this difference is due to the significantly lower non-electrostatic energy for the minor groove complex with AT alternating DNA, whereas the slightly lower binding energy to this sequence is caused by a higher deformation energy of DNA. The energetic data are in agreement with the conclusions derived from different spectroscopic studies and can also be structurally interpreted on the basis of the modeled complexes. The simple static modeling technique and the neglect of entropy terms and of non-electrostatic solute-solvent interactions, which are assumed to be nearly constant for the compared complexes of MB with DNA, seem to be justified by the results.     

Structural insight into the kinetics and DeltaCp of interactions between TEM-1 beta-lactamase and beta-lactamase inhibitory protein (BLIP)

In a previous study, we examined thermodynamic parameters for 20 alanine mutants in beta-lactamase inhibitory protein (BLIP) for binding to TEM-1 beta-lactamase. Here we have determined the structures of two thermodynamically distinctive complexes of BLIP mutants with TEM-1 beta-lactamase. The complex BLIP Y51A-TEM-1 is a tight binding complex with the most negative binding heat capacity change (DeltaG = approximately -13 kcal mol(-1) and DeltaCp = approximately -0.8 kcal mol(-1) K(-1)) among all of the mutants, whereas BLIP W150A-TEM-1 is a weak complex with one of the least negative binding heat capacity changes (DeltaG = approximately -8.5 kcal mol(-1) and DeltaCp = approximately -0.27 kcal mol(-1) K(-1)). We previously determined that BLIP Tyr51 is a canonical and Trp150 an anti-canonical TEM-1-contact residue, where canonical refers to the alanine substitution resulting in a matched change in the hydrophobicity of binding free energy. Structure determination indicates a rearrangement of the interactions between Asp49 of the W150A BLIP mutant and the catalytic pocket of TEM-1. The Asp49 of W150A moves more than 4 angstroms to form two new hydrogen bonds while losing four original hydrogen bonds. This explains the anti-canonical nature of the Trp150 to alanine substitution, and also reveals a strong long distance coupling between Trp150 and Asp49 of BLIP, because these two residues are more than 25 angstroms apart. Kinetic measurements indicate that the mutations influence the dissociation rate but not the association rate. Further analysis of the structures indicates that an increased number of interface-trapped water molecules correlate with poor interface packing in a mutant. It appears that the increase of interface-trapped water molecules is inversely correlated with negative binding heat capacity changes.     

1H NMR study of the interaction of bacteriophage lambda Cro protein with the OR3 operator. Evidence for a change of the conformation of the OR3 operator on binding.

The specific complex between the lambda phage OR3 operator and the Cro protein has been studied by proton NMR spectroscopy at 500 MHz. The DNA imino proton resonances of this complex have been assigned to specific base pairs using the known assignments of these resonances for the free operator. Increase of the protein/DNA ratio to complete saturation of the OR3 operator with the Cro protein made it possible to follow the shift changes of the resonances. Ambiguities were resolved by nuclear Overhauser effect measurements on the complex. The shifts of the imino proton resonance positions provide information on the changes induced in the conformation of the operator upon complex formation with a dimer of the Cro protein. The most striking shift occurs for the central (GC 9) base pair, which is known to have no direct contacts with the Cro protein. This shift may be induced by a bend in the OR3 operator DNA at the GC 9 base pair to accommodate the operator for the binding of the Cro protein dimer. The imino proton resonances of two additional base pairs can be observed in the complex, demonstrating an overall stabilization of the DNA structure by the binding of the Cro protein.     

Solvation change and ion release during aminoacylation by aminoacyl-tRNA synthetases

Discrimination between cognate and non-cognate tRNAs by aminoacyl-tRNA synthetases occurs at several steps of the aminoacylation pathway. We have measured changes of solvation and counter-ion distribution at various steps of the aminoacylation pathway of glutamyl- and glutaminyl-tRNA synthetases. The decrease in the association constant with increasing KCl concentration is relatively small for cognate tRNA binding when compared to known DNA–protein interactions. The electro-neutral nature of the tRNA binding domain may be largely responsible for this low ion release stoichiometry, suggesting that a relatively large electrostatic component of the DNA–protein interaction free energy may have evolved for other purposes, such as, target search. Little change in solvation upon tRNA binding is seen. Non-cognate tRNA binding actually increases with increasing KCl concentration indicating that charge repulsion may be a significant component of binding free energy. Thus, electrostatic interactions may have been used to discriminate between cognate and non-cognate tRNAs in the binding step. The catalytic constant of glutaminyl-tRNA synthetase increases with increasing osmotic pressure indicating a water release of 8.4 ± 1.4 mol/mol in the transition state, whereas little change is seen in the case of glutamyl-tRNA synthetase. We propose that the significant amount of water release in the transition state, in the case of glutaminyl-tRNA synthetase, is due to additional contact of the protein with the tRNA in the transition state.     

On the calculation of binding free energies using continuum methods: Application to MHC class I protein-peptide interactions

This paper describes a methodology to calculate the binding free energy (ΔG) of a protein-ligand complex using a continuum model of the solvent. A formal thermodynamic cycle is used to decompose the binding free energy into electrostatic and non-electrostatic contributions. In this cycle, the reactants are discharged in water, associated as purely nonpolar entities, and the final complex is then recharged. The total electrostatic free energies of the protein, the ligand, and the complex in water are calculated with the finite difference Poisson-Boltzmann (FDPB) method. The nonpolar (hydrophobic) binding free energy is calculated using a free energy-surface area relationship, with a single alkane/water surface tension coefficient (γ_aw). The loss in backbone and side-chain configurational entropy upon binding is estimated and added to the electrostatic and the nonpolar components of ΔG. The methodology is applied to the binding of the murine MHC class I protein H-2K^b with three distinct peptides, and to the human MHC class I protein HLA-A2 in complex with five different peptides. Despite significant differences in the amino acid sequences of the different peptides, the experimental binding free energy differences (ΔΔG_exp) are quite small (<0.3 and <2.7 kcal/mol for the H-2K^b and HLA-A2 complexes, respectively). For each protein, the calculations are successful in reproducing a fairly small range of values for ΔΔG_calc (<4.4 and <5.2 kcal/mol, respectively) although the relative peptide binding affinities of H-2K^b and HLA-A2 are not reproduced. For all protein-peptide complexes that were treated, it was found that electrostatic interactions oppose binding whereas nonpolar interactions drive complex formation. The two types of interactions appear to be correlated in that larger nonpolar contributions to binding are generally opposed by increased electrostatic contributions favoring dissociation. The factors that drive the binding of peptides to MHC proteins are discussed in light of our results.     

Repertoire selection of variant single-chain Cro: toward directed DNA-binding specificity of helix-turn-helix proteins

A single-chain derivative of the lambda Cro repressor (scCro) has been randomly mutated in amino acid residues critical for specific DNA recognition to create libraries of protein variants. Utilizing phage display-afforded affinity selection, scCro variants have been isolated for binding to synthetic DNA ligands. Isolated scCro variants were analyzed functionally, both in fusion with phage particles and after expression of the corresponding free proteins. The binding properties with regard to specificity and affinity in binding to different DNA ligands were investigated by inhibition studies and determination of equilibrium dissociation constants for formed complexes. Variant proteins with altered DNA-sequence specificity were identified, which favored binding of targeted synthetic DNA sequences over a consensus operator sequence, bound with high affinity by wild-type Cro. The specificities were relatively modest (2-3-fold, as calculated from K(D) values), which can be attributed to the inherent properties in the design of the selection system; one half-site of the synthetic DNA sequences maintains the consensus operator sequence, and one "subunit" of the variant single-chain Cro dimers was conserved as wild-type sequence. The anticipated interaction between the wild-type subunit and the consensus DNA half-site of target DNA ligands is, hence, expected to contribute to the overlap in sequence discrimination. The binding affinity for the synthetic DNA sequences, however, was improved 10-30-fold in selected variant proteins as compared to "wild-type" scCro.