Different interactions used by Cro repressor in specific and nonspecific DNA binding

The mode of interaction of Cro repressor with specific and nonspecific sites on DNA was explored by chemical modification and protection of lysine and tyrosine residues. Cro has 8 lysines. In the presence of DNA, lysines 32 and 56 are fully protected and lysines 21, 62, and 63 are partially protected from alkylation. However, the terminal amino group and lysines 8, 18, and 39 are not protected. Location of the protected and unprotected lysines on the three-dimensional Cro structure defines a DNA-binding region. The results provide direct experimental support for a mode of interaction between Cro and DNA, in which Cro buries its 2-fold related alpha-helices in consecutive DNA major grooves (Anderson, W. F., Ohlendorf, D. H., Takeda, Y., and Matthews, B. W. (1981) Nature 290, 754-758; Ohlendorf, D. H., Anderson, W. F., Fisher, R. G., Takeda, Y., and Matthews, B. W. (1982) Nature 298, 718-723). In the model, the carboxyl-terminal part of Cro was tentatively presumed to interact with the DNA minor groove. Protection of lysines 62 and 63 confirms the involvement of the carboxyl terminus in DNA binding. Although nonspecific and specific DNA protect the same lysine residues, there are differences in the nature of the interaction of Cro with nonspecific and specific DNA. Cro-nonspecific DNA interaction is salt-sensitive, suggesting that the interaction is predominantly electrostatic. On the other hand, Cro-specific DNA interaction is salt-resistant, suggesting that the interaction may include nonelectrostatic components (hydrogen bonds and hydrophobic interactions) as well. Protection experiments of tyrosine residues (against iodination) suggest that the conformation of Cro repressor changes in two stages: first, when Cro binds at nonspecific sites, and, second, when Cro binds to specific sites on DNA.     

Prediction of 3-D structure of the Cro protein from phage 434

A comparative model building process has been utilized to predict the three-dimensional structure of the bacteriophage 434 Cro protein. Amino acid sequence similarities between the 434 Cro protein and other bacteriophage repressor and Cro proteins have been used, in conjunction with secondary structure prediction and the known structures of other base sequence specific DNA binding proteins, to derive the model. From this model the interactions between the 434 Cro protein and its operator DNA have been deduced. These proposed interactions are consistent with the known properties of the bacteriophage 434 Cro protein.     

Prediction of 3-D Structure of the Cro Protein from Phage 434

Abstract

A comparative model building process has been utilized to predict (he three-dimensional structure of the bacteriophage 434 Cro protein, Amino acid sequence similarities between the 434 Cro protein and other bacteriophage repressor and Cro proteins have been used, in conjunction with secondary structure prediction and the known structures of other base sequence specific DNA binding proteins, to derive the model. From this model the interactions between the 434 Cro protein and its operator DNA have been deduced. These proposed interactions are consistent with the known properties of the bacteriophage 434 Cro protein.     

Helix-capping interaction in λ cro protein: A free energy simulation analysis

The stability mutant Tyr-26 → Asp was studied in the Cro protein from bacteriophage λ using free energy molecular dynamics simulations. The mutant was calculated to be more stable than the wild type by 3.0 ± 1.7 kcal/mol/monomer, in reasonable agreement with experiment (1.4 kcal/mol/monomer). Moreover, the aspartic acid in the mutant was found to form a capping interation with the amino terminus of the third α-helix of Cro. The simulations were analyzed to understand better the source of the stability of this helix-capping interaction and to examine the results in light of previous explanations of stabilizing helix caps-namely, a model of local unsatisfied hydrogen bonds at the helix termini and the helix macro dipole model. Analysis of the simulations shows that the stabilizing effect of this charged helical cap is due both to favorable hydrogen bonds with backbone NH groups at the helix terminus and to favorable electrostatic interactions (but not hydrogen bonds) with their carbonyls (effectively the next row of local dipoles in the helix). However, electrostatic interactions are weak or negligible with backbone dipolar groups in the helix further away from the terminus. Moreover, the importance of other local electrostatic interactions with polar side chains near the helix terminus, which are neglected in most treatments of this effect, are shown to be important. Thus, the results support a model that is intermediate between the two previous explanations: both unsatisfied hydrogen bonds at the helix terminus and other, local preoriented dipolar groups stabilize the helix cap. These findings suggest that similar interactions with preoriented dipolar groups may be important for cooperativity in other charge–dipole interactions and may be employed to advantage for molecular design. © 1994 Wiley-Liss, Inc.     

Electrostatic background of chromatin fiber stretching

We have carried out an investigation of the electrostatic forces involved in gradual removal of the DNA from the histone proteins in chromatin. Two simple models of DNA-histone core dissociation were considered. Calculations of the electrostatic free energy within the Poisson-Boltzmann theory gave similar results for the both models, which turned out to be in a qualitative agreement with recent optical tweezers stretching experiments measuring the force necessary to unwrap DNA from the histone core. Our analysis shows that the electrostatic interactions between the highly negatively charged polymeric DNA and the positively charged histones play a determining role in stabilizing the nucleosomes at physiological conditions.     

Electrostatic Exclusion of Neutral Solutes from Condensed DNA and Other Charged Phases

Motivated by experiments on condensed DNA phases in binary mixtures of water and a low-dielectric solute, we develop a theory for the electrostatic contribution to solute exclusion from a highly charged phase, within the continuum approximation of the medium. Because the electric field is maximum at the surface of each ion, the electrostatic energy is dominated by the Born energy; interactions between charges are of secondary importance. Neglecting interactions and considering only the competition between the Born energy and the free energy of mixing, we predict that low dielectric solutes are excluded from condensed DNA phases in water-cosolvent mixtures. This suggests that the traditional continuum electrostatic approach of modeling binary mixtures with a uniform dielectric constant needs to be modified. The linking of solute exclusion to solute dielectric properties also suggests a mechanism for predicting the electrostatic contribution to preferential hydration of polar and charged surfaces.     

High Resolution Structural Studies of Cro Repressor Protein and Implications for DNA Recognition

Abstract

Cro repressor is a small dimeric protein that binds to specific sites on the DNA of bacteriophage λ. The structure of Cro has been determined and suggests that the protein binds to its sequence-specific sites with a pair of two-fold related α-helices of the protein located within successive major grooves of the DNA.

From the known three-dimensional structure of the repressor, model building and energy refinement have been used to develop a detailed model for the presumed complex between Cro and DNA. Recognition of specific DNA binding sites appears to occur via multiple hydrogen bonds between amino acid side chains of the protein and base pair atoms exposed within the major groove of DNA. The Cro:DNA model is consistent with the calculated electrostatic potential energy surface of the protein.

From a series of amino acid sequence and gene sequence comparisons, it appears that a number of other DNA-binding proteins have an α-helical DNA-binding region similar to that seen in Cro. The apparent sequence homology includes not only DNA-binding proteins from different bacteriophages, but also gene-regulatory proteins from bacteria and yeast. It has also been found that the conformations of part of the presumed DNA-binding regions of Cro repressor, λ repressor and CAP gene activator proteins are strikingly similar. Taken together, these results strongly suggest that a two-helical structural unit occurs in the DNA-binding region of many proteins that regulate gene expression. However, the results to date do not suggest that there is a simple one-to-one recognition code between amino acids and bases.

Crystals have been obtained of complexes of Cro with six-base-pair and nine-basepair DNA oligomers, and X-ray analysis of these co-crystals is in progress.     

Molecular Interactions and Residues Involved in Force Generation in the T4 Viral DNA Packaging Motor

Many viruses utilize molecular motors to package their genomes into preformed capsids. A striking feature of these motors is their ability to generate large forces to drive DNA translocation against entropic, electrostatic, and bending forces resisting DNA confinement. A model based on recently resolved structures of the bacteriophage T4 motor protein gp17 suggests that this motor generates large forces by undergoing a conformational change from an extended to a compact state. This transition is proposed to be driven by electrostatic interactions between complementarily charged residues across the interface between the N- and C-terminal domains of gp17. Here we use atomistic molecular dynamics simulations to investigate in detail the molecular interactions and residues involved in such a compaction transition of gp17. We find that although electrostatic interactions between charged residues contribute significantly to the overall free energy change of compaction, interactions mediated by the uncharged residues are equally if not more important. We identify five charged residues and six uncharged residues at the interface that play a dominant role in the compaction transition and also reveal salt bridging, van der Waals, and solvent hydrogen-bonding interactions mediated by these residues in stabilizing the compact form of gp17. The formation of a salt bridge between Glu309 and Arg494 is found to be particularly crucial, consistent with experiments showing complete abrogation in packaging upon Glu309Lys mutation. The computed contributions of several other residues are also found to correlate well with single-molecule measurements of impairments in DNA translocation activity caused by site-directed mutations.     

Crystal structure of an engineered Cro monomer bound nonspecifically to DNA: possible implications for nonspecific binding by the wild-type protein.

The structure has been determined at 3.0 A resolution of a complex of engineered monomeric Cro repressor with a seven-base pair DNA fragment. Although the sequence of the DNA corresponds to the consensus half-operator that is recognized by each subunit of the wild-type Cro dimer, the complex that is formed in the crystals by the isolated monomer appears to correspond to a sequence-independent mode of association. The overall orientation of the protein relative to the DNA is markedly different from that observed for Cro dimer bound to a consensus operator. The recognition helix is rotated 48 degrees further out of the major groove, while the turn region of the helix-turn-helix remains in contact with the DNA backbone. All of the direct base-specific interactions seen in the wild-type Cro-operator complex are lost. Virtually all of the ionic interactions with the DNA backbone, however, are maintained, as is the subset of contacts between the DNA backbone and a channel on the protein surface. Overall, 25% less surface area is buried at the protein DNA interface than for half of the wild-type Cro-operator complex, and the contacts are more ionic in character due to a reduction of hydrogen bonding and van der Waals interactions. Based on this crystal structure, model building was used to develop a possible model for the sequence-nonspecific interaction of the wild-type Cro dimer with DNA. In the sequence-specific complex, the DNA is bent, the protein dimer undergoes a large hinge-bending motion relative to the uncomplexed form, and the complex is twofold symmetric. In contrast, in the proposed nonspecific complex the DNA is straight, the protein retains a conformation similar to the apo form, and the complex lacks twofold symmetry. The model is consistent with thermodynamic, chemical, and mutagenic studies, and suggests that hinge bending of the Cro dimer may be critical in permitting the transition from the binding of protein at generic sites on the DNA to binding at high affinity operator sites.     

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.     

Role of the Cro repressor carboxy-terminal domain and flexible dimer linkage in operator and nonspecific DNA binding

A series of mutations comprising single and multiple substitutions, deletions, and extensions within the carboxy-terminal domain of the bacteriophage lambda Cro repressor have been constructed. These mutations generally affect the affinity of repressor for specific and nonspecific DNA. Additionally, substitution of the carboxy-terminal alanine with several amino acids capable of hydrogen-bonding interactions leads to improved specific binding affinities. A mutation is also described whereby cysteine links the two Cro monomers by a disulfide bond. As a consequence, a significant improvement in nonspecific binding and a concomitant reduction in specific binding are observed with this mutant. These results provide evidence that the carboxy terminus of Cro repressor is an important DNA binding domain and that a flexible connection between the two repressor monomers is a critical factor in modulating the affinity of wild-type repressor for DNA.     

Peptide nucleic acid (PNA) amphiphiles: synthesis, self-assembly, and duplex stability

Peptide amphiphiles comprising a class of conjugates of peptide nucleic acid (PNA), natural amino acids, and n-alkanes were synthesized and studied. These PNA amphiphiles (PNAA) self-assemble at concentrations between 10 and 50 muM and exhibit water solubilities above 500 muM. The highly specific, stable DNA binding properties of PNAs are preserved by these modifications, with no significant differences between the thermodynamics of DNA binding of the PNA peptide and the PNA amphiphile. Proper solubilization of the PNAA required the attachment of (Lys)(2) and (Glu)(4) peptides to PNAs, which affected the PNAA-DNA duplex stability by electrostatic interactions between these charged amino acids and the negatively charged DNA backbone. These electrostatic effects did not affect the specificity of DNA binding, however. Electrostatic effects are screened with added salt, in a manner consistent with previous studies of PNA-DNA duplex stability and predictions from a charged-cylinder model for the duplex.     

A comparative study of dynamic structures between phage 434 Cro and repressor proteins by normal mode analysis

Two DNA binding proteins, Cro and the amino-terminal domain of the repressor of bacteriophage 434 (434 Cro and 434 repressor) that regulate gene expression and contain a helix-turn-helix (HTH) motif responsible for their site-specific DNA recognition adopt very similar three-dimensional structures when compared to each other. To reveal structural differences between these two similar proteins, their dynamic structures, as examined by normal mode analysis, are compared in this paper. Two kinds of structural data, one for the monomer and the other for a complex with DNA, for each protein, are used in the analyses. From a comparison between the monomers it is found that the interactions of Ala-24 in 434 Cro or Val-24 in 434 repressor, both located in the HTH motif, with residues 44, 47, 48, and 51 located in the domain facing the motif, and the interactions between residues 17, 18, 28, and 32, located in the HTH motif, cause significant differences in the correlative motions of these residues. From the comparison between the monomer and the complex with DNA for each protein, it was found that the first helix in the HTH motif is distorted in the complex form. While the residues in the HTH motif in 434 Cro have relatively larger positive correlation coefficients of motions with other residues within the HTH motif, such correlations are not large in the HTH motif of 434 repressor. It is suggestive to their specificity because the 434 repressor is less specific than 434 Cro. Although a structural comparison of proteins has been performed mainly from a static or geometrical point of view, this study demonstrates that the comparison from a dynamic point of view, using the normal mode analysis, is useful and convenient to explore a difference that is difficult to find only from a geometrical point of view, especially for proteins very similar in structure. © 1996 Wiley-Liss, Inc.     

Single molecule force spectroscopy reveals that electrostatic interactions affect the mechanical stability of proteins

It is well known that electrostatic interactions play important roles in determining the thermodynamic stability of proteins. However, the investigation into the role of electrostatic interactions in mechanical unfolding of proteins has just begun. Here we used single molecule atomic force microscopy techniques to directly evaluate the effect of electrostatic interactions on the mechanical stability of a small protein GB1. We engineered a bi-histidine motif into the force-bearing region of GB1. By varying the pH, histidine residues can switch between protonated and deprotonated states, leading to the change of the electrostatic interactions between the two histidine residues. We found that the mechanical unfolding force of the engineered protein decreased by ∼34% (from 115 pN to 76 pN) on changing the pH from 8.5 to 3, due to the increased electrostatic repulsion between the two positively charged histidines at acidic pH. Our results demonstrated that electrostatic interactions can significantly affect the mechanical stability of elastomeric proteins, and modulating the electrostatic interactions of key charged residues can become a promising method for regulating the mechanical stability of elastomeric proteins.     

Salt Concentration Effects on Equilibrium Melting Curves from DNA Microarrays

DNA microarrays find applications in an increasing number of domains where more quantitative results are required. DNA being a charged polymer, the repulsive interactions between the surface of the microarray and the targets in solution are increasing upon hybridization. Such electrostatic penalty is generally reduced by increasing the salt concentration. In this article, we present equilibrium-melting curves obtained from dedicated physicochemical experiments on DNA microarrays in order to get a better understanding of the electrostatic penalty incurred during the hybridization reaction at the surface. Various salt concentrations have been considered and deviations from the commonly used Langmuir adsorption model are experimentally quantified for the first time in agreement with theoretical predictions.     

Solvent density and long-range dipole field around a DNA-binding protein studied by molecular dynamics

The distribution and orientation of solvent around a DNA-binding protein, 434 Cro, were investigated by molecular dynamics simulations with a periodic-boundary condition. The protein was treated in two states: charged and neutral. The computed high-density sites of the solvent around the protein correlated well with the experimentally determined crystal-water sites, in both the charged and neutral states. A local density map, introduced to investigate the solvent density around the highly mobile regions of the protein, showed a hydration shell around hydrophobic sidechains and hydrogen-bondable sites around hydrophilic sidechains, and also showed that the solvent density is sensitive to the slight concaves of the sidechain surface. The long-range solvent-dipole field was observed around the protein, where the pattern of the dipole ordering was considerably different between the charged and neutral states. A local solvent-dipole field was introduced, and the pattern of the dipole ordering was different between the hydrophobic and hydrophilic sidechains. The dipole field from the charged state provided a higher correlation to the electrostatic field obtained from the Poisson-Boltzmann's equation than that from the neutral state, although the correlation become weak quickly for the both states with increasing the protein-solvent distance.     

Electrostatic Potential of B-DNA: Effect of Interionic Correlations

Modified Poisson-Boltzmann (MPB) equations have been numerically solved to study ionic distributions and mean electrostatic potentials around a macromolecule of arbitrarily complex shape and charge distribution. Results for DNA are compared with those obtained by classical Poisson-Boltzmann (PB) calculations. The comparisons were made for 1:1 and 2:1 electrolytes at ionic strengths up to 1 M. It is found that ion-image charge interactions and interionic correlations, which are neglected by the PB equation, have relatively weak effects on the electrostatic potential at charged groups of the DNA. The PB equation predicts errors in the long-range electrostatic part of the free energy that are only ∼1.5 kJ/mol per nucleotide even in the case of an asymmetrical electrolyte. In contrast, the spatial correlations between ions drastically affect the electrostatic potential at significant separations from the macromolecule leading to a clearly predicted effect of charge overneutralization.     

Electrostatic contributions to the binding free energy of the lambdacI repressor to DNA.

A model based on the nonlinear Poisson-Boltzmann (NLPB) equation is used to study the electrostatic contribution to the binding free energy of the lambdacI repressor to its operator DNA. In particular, we use the Poisson-Boltzmann model to calculate the pKa shift of individual ionizable amino acids upon binding. We find that three residues on each monomer, Glu34, Glu83, and the amino terminus, have significant changes in their pKa and titrate between pH 4 and 9. This information is then used to calculate the pH dependence of the binding free energy. We find that the calculated pH dependence of binding accurately reproduces the available experimental data over a range of physiological pH values. The NLPB equation is then used to develop an overall picture of the electrostatics of the lambdacI repressor-operator interaction. We find that long-range Coulombic forces associated with the highly charged nucleic acid provide a strong driving force for the interaction of the protein with the DNA. These favorable electrostatic interactions are opposed, however, by unfavorable changes in the solvation of both the protein and the DNA upon binding. Specifically, the formation of a protein-DNA complex removes both charged and polar groups at the binding interface from solvent while it displaces salt from around the nucleic acid. As a result, the electrostatic contribution to the lambdacI repressor-operator interaction opposes binding by approximately 73 kcal/mol at physiological salt concentrations and neutral pH. A variety of entropic terms also oppose binding. The major force driving the binding process appears to be release of interfacial water from the protein and DNA surfaces upon complexation and, possibly, enhanced packing interactions between the protein and DNA in the interface. When the various nonelectrostatic terms are described with simple models that have been applied previously to other binding processes, a general picture of protein/DNA association emerges in which binding is driven by the nonpolar interactions, whereas specificity results from electrostatic interactions that weaken binding but are necessary components of any protein/DNA complex.