A Molecular Dynamics Study of the Repressor/Operator(OR1,OR3) Complexes from Bacteriophage 434

We have performed three molecular dynamics simulations using the CHARMM molecular modeling program to study the repressor protein from bacteriophage 434 complexed with DNA operators of two different sequences. Two approaches to the modeling of the solvent were used. In the first method, applied to the R1-69/OR1 truncated complex, water molecules were included explicitly in conjunction with a stochastic boundary force to solvate the complex. In the second approach, used for simulations of the R1-69/OR1 and the R1-69/OR3 complexes, the solvent was omitted and implicitly represented by using a distance-dependent dielectric constant and a scaling of the charges on the exposed residues. The simulation with the model which explicitly includes the solvent serves as a validation of the simulations using a simpler solvent representation. In our discussion of the results we focus upon the important interactions between the DNA binding motif of the 434 repressor (motif helix turn helix) and the operators and how the structures of the complexes change with time.Supplementary material to this paper is available in electronic form at http://dx.doi.org/10.1007/s0089460020427     

DNA recognition by the helix-turn-helix motif: investigation by laser Raman spectroscopy of the phage lambda repressor and its interaction with operator sites OL1 and OR3.

The lambda repressor provides a model system for biophysical studies of DNA recognition by the helix-turn-helix motif. We describe laser Raman studies of the lambda operator sites OL1 and OR3 and their interaction with the DNA-binding domain of lambda repressor (residues 1-102). Raman spectra of the two DNA sites exhibit significant differences attributable to interstrand purine-purine steps that differ in the two oligonucleotides. Remarkably, the conformation of each operator is significantly and specifically altered by repressor binding. Protein recognition, which involves hydrogen-bond formation and hydrophobic contacts in the major groove, induces subtle changes in DNA Raman bands of interacting groups. These include (i) site-specific perturbations to backbone phosphodiester geometry at AT-rich domains, (ii) hydrophobic interaction at thymine 5CH3 groups, (iii) hydrogen bonding to guanine 7N and 6C = O acceptors, and (iv) alterations in sugar pucker within the C2'-endo (B-DNA) family. These perturbations differ between aqueous OL1 and OR3 complexes of repressor, indicating that protein binding in solution determines the precise DNA conformation. The overall structure of the lambda domain is not greatly perturbed by binding to either OL1 or OR3, in accord with X-ray studies of other complexes. However, Raman markers indicate a change in hydrogen bonding of the OH group of tyrosine-22, which is a hydrogen-bond acceptor in the absence of DNA but a combined donor and acceptor in the OL1 complex; yet, Y22 hydrogen bonding is not altered in forming the OR3 complex. The present results demonstrate qualitatively different and distinguishable modes of interaction of the lambda repressor DNA-binding domain with operators OL1 and OR3 in solution. This application of laser Raman spectroscopy to a well-characterized system provides a prototype for future Raman studies of other DNA-binding motifs under physiological conditions.     

Monovalent cations regulate DNA sequence recognition by 434 repressor

The bacteriophage 434 repressor distinguishes between its six naturally occurring binding sites using indirect readout. In indirect readout, sequence-dependent differences in the structure and flexibility of non-contacted bases in a protein's DNA-binding site modulate the affinity of DNA for protein. The conformation and flexibility of a DNA sequence can be influenced by the interaction of the DNA bases or backbone with solution components. We examined the effect of changing the cation-type present in solution on the stability and structure of 434 repressor complexes with wild-type and mutant OR1 and OR3, binding sites that differ in their contacted and non-contacted base sequences. We find that the affinity of repressor for OR1, but not for OR3, depends remarkably on the type and concentration of monovalent cation. Moreover, the formation of a stable, specific repressor-OR1 complex requires the presence of monovalent cations; however, repressor-OR3 complex formation has no such requirement. Changing monovalent cation type alters the ability of repressor to protect OR1, but not OR3, from *OH radical cleavage. Altering the relative length of the poly(dA) x poly(dT) tract in the non-contacted regions of the OR1 and OR3 can reverse the cation sensitivity of repressor's affinities for these two sites. Taken together these findings show that cation-dependent alterations in DNA structure underlies indirect readout of DNA sequence by 434 repressor and perhaps other proteins.     

The role of the minor groove substituents in indirect readout of DNA sequence by 434 repressor

The sequence of non-contacted bases at the center of the 434 repressor binding site affects the strength of the repressor-DNA complex by influencing the structure and flexibility of DNA (Koudelka, G. B., and Carlson, P. (1992) Nature 355, 89-91). We synthesized 434 repressor binding sites that differ in their central sequence base composition to test the importance of minor groove substituents and/or the number of base pair hydrogen bonds between these base pairs on DNA structure and strength of the repressor-DNA complex. We show here that the number of base pair H-bonds between the central bases apparently has no role in determining the relative affinity of a DNA site for repressor. Instead we find that the affinity of DNA for repressor depends on the absence or presence the N2-NH(2) group on the purine bases at the binding site center. The N2-NH(2) group on bases at the center of the 434 binding site appears to destabilize 434 repressor-DNA complexes by decreasing the intimacy of the specific repressor-DNA contacts, while increasing the reliance on protein contacts to the DNA phosphate backbone. Thus, the presence of an N2-NH(2) group on the purines at the center of a binding site globally alters the precise conformation of the protein-DNA interface.     

Synthesis of nonlinear DNA-binding peptide with binding specificity determinants close to those of 434 Cro-repressor

Design, synthesis and DNA binding activity of a nonlinear 102 residue peptide are reported. The peptide contains four sequence-specific DNA binding domains of 434 Cro protein. These four domains were linked covalently to a symmetrical carboxyterminal crosslinker that contains four arms each ending with an aliphatic aminogroup. From CD studies we have found that in aqueous buffer in the presence of 20% trifluoroethanol the peptide residues assume alpha helical, beta-sheet and random coiled conformations with an alpha helical content of about 16% at room temperature. The alpha helicity is increased up to 40% in the presence of 40% trifluoroethanol. Upon complex formation between the peptide and DNA a change in the peptide conformation takes place which is consistent with an alpha-beta transition in the DNA binding, helix-turn-helix motif of 434 Cro repressor. Evidently residues present in helices alpha(2) and alpha(3) form a beta hairpin which is inserted in the minor DNA groove. The latter inference is supported by our observations that the peptide can displace minor groove binding antibiotic distamycin A from a complex with poly(dA).poly(dT). As revealed from DNase protection studies the peptide exhibits preferences for binding to operator and pseudooperator sites recognized by 434 Cro repressor. It binds strongly to operator sites OR1, OR2 and OR3 and exhibits a greater affinity for pseudooperator site Op1. From analysis of nucleotide sequences in the strong affinity binding sites for the peptide on DNA a conclusion is drawn that it binds to pseudosymmetrical nucleotide sequences 5'-ACAA(W)nCTGT-3', where W is an arbitrary nucleotide. n is equal to six or seven. In the strongest affinity binding site for the peptide on DNA (Op1) motif 5'-ACAA-3' is replaced by sequence 5'-ACCA-3'. A difference in binding specificity shown by the peptide and 434 Cro protein could be attributed to a flexibility of the connecting chains between DNA-binding domains in the peptide molecule as well as to a replacement of Thr - Ala in the alpha 2 helix. Removal of two residues from the N-terminal end of helix alpha 2 in each of the four DNA binding domains of 434 Cro present in the peptide leads to a loss of binding specificity, although the modified peptide binds to DNA unspecifically.     

Expression, purification, and functional characterization of the carboxyl-terminal domain fragment of bacteriophage 434 repressor.

The repressor protein of bacteriophage 434 binds to DNA as a dimer of identical subunits. Its strong dimerization is mediated by the carboxyl-terminal domain. Cooperative interactions between the C-terminal domains of two repressor dimers bound at adjacent sites can stabilize protein-DNA complexes formed with low-affinity binding sites. We have constructed a plasmid, pCT1, which directs the overproduction of the carboxyl-terminal domain of 434 repressor. The protein encoded by this plasmid is called CT-1. Cells transformed with pCT1 are unable to be lysogenized by wild-type 434 phage, whereas control cells are lysogenized at an efficiency of 1 to 5%. The CT-1-mediated interference with lysogen formation presumably results from formation of heteromeric complexes between the phage-encoded repressor and the plasmid-encoded carboxyl-terminal domain fragment. These heteromers are unable to bind DNA and thereby inhibit the repressor's activity in promoting lysogen formation. Two lines of evidence support this conclusion. First, DNase I footprinting experiments show that at a 2:1 ratio of CT-1 to intact 434 repressor, purified CT-1 protein prevents the formation of complexes between 434 repressor and its OR1 binding site. Second, cross-linking experiments reveal that only a specific heterodimeric complex forms between CT-1 and intact 434 repressor. This latter observation indicates that CT-1 interferes with 434 repressor-operator complex formation by preventing dimerization and not by altering the conformation of the DNA-bound repressor dimer. Our other evidence is also consistent with this suggestion. We have used deletion analysis in an attempt to define the region which mediates the 434 repressor-CT-1 interaction. CT-1 proteins which have more than the last 14 amino acids removed are unable to interfere with 434 repressor action in vivo.     

Substituting an alpha-helix switches the sequence-specific DNA interactions of a repressor

It has been suggested that many DNA-binding proteins use an alpha-helix for specific sequence recognition. We have used amino acid sequence homologies to identify the presumptive DNA-recognition helices in two related proteins whose structures are unknown--the repressor and cro protein of bacteriophage 434. The 434 repressor and cro protein each bind to three similar sites in the rightward phage 434 operator, OR, and they make different contacts in each binding site, as revealed by the chemical probe dimethyl sulfate. We substituted the putative recognition alpha-helix of 434 repressor with the putative recognition alpha-helix of 434 cro protein to create a hybrid protein named repressor*. The specific DNA contacts made by repressor* are like those of 434 cro protein.     

Interaction of lambda cro repressor with synthetic operator OR3 studied by competition binding with minor groove binders.

In the present work, we employ a combination of CD spectroscopy and gel retardation technique to characterize thermodynamically the binding of lambda phage cro repressor to a 17 base pair operator OR3. We have found that three minor groove-binding antibiotics, distamycin A, netropsin and sibiromycin, compete effectively with the cro for binding to the operator OR3. Among these antibiotics, sibiromycin binds covalently to DNA in the minor groove at the NH2 of guanine, whereas distamycin A and netropsin interact preferentially with runs of AT base pairs and avoid DNA regions containing guanine bases in the two polynucleotide strands. Only subtle DNA conformation changes are known to take place upon binding of these antibiotics. Both the CD spectral profiles and the results of the gel retardation experiments indicate that distamycin A and netropsin can displace cro repressor from the operator OR3. The binding of cro repressor to the OR3 is accompanied by considerable changes in CD in the far-UV region which appear to be attributed to a DNA-dependent structural transition in the protein. Spectral changes are also induced in the wavelength region of 270-290 nm. The CD spectral profile of the cro-OR3 mixture in the presence of distamycin A can be represented as a sum of the CD spectrum of the repressor-operator complex and spectrum of distamycin-DNA complex at the appropriate molar ratio of the bound antibiotic to the operator DNA (r). When r tends to the saturation level of binding the CD spectrum in the region of 270-360 nm approaches a CD pattern typical of complexes of the antibiotic with the free DNA oligomer. This suggests that simultaneous binding of cro repressor and distamycin A to the same DNA oligomer is not possible and that distamycin A and netropsin can be used to determine the equilibrium affinity constant of cro repressor to the synthetic operator from competition-type experiments. The binding constant of cro repressor to the OR3 is found to be (6 +/- 1).10(6)M-1 at 20 degrees C in 10 mM sodium cacodylate buffer (pH 7.0) in the presence of 0.1 M NH4F.     

Recognition of DNA structure by 434 repressor.

In complexes of bacteriophage 434 binding sites with 434 repressor the central 4 bp of the 14 bp site are not contacted by the protein, although changes in these bases alter binding site affinity for the repressor. Our previous data suggested that the ability of the non-contacted central bases to be overtwisted in repressor-DNA complexes governs affinity of the binding site for 434 repressor. This idea was tested by examining the affinity of two central sequence variant 434 binding sites for 434 repressor as a function of binding site average twist. The 434 repressor preferred the relatively overwound binding site to the two more underwound forms. The greatest affinity enhancement resulting from increasing twist was observed with a binding site that is relatively underwound and more resistant to twisting deformation. Consistent with the idea that 434 repressor overtwists its binding site upon DNA binding, we show that 434 repressor is capable of binding to sites bearing a single base insertion in their center (a 15mer), but binds poorly to binding sites bearing central base deletions (12mer and 13mer). The N-terminal dimer interface plays a large role in determining 434 repressor central base preferences. Mutations in this interface eliminate central base discrimination and/or site size preferences. These mutations also lead to changes in the size of the repressor footprint on the various sized DNA sites that are consistent with their binding characteristics.     

On the importance of van der Waals interaction in the groove binding of DNA with ligands: restrained molecular dynamics study

Comparison of interaction energy between an oligonucleotide and a DNA-binding ligand in the minor and major groove modes was made by use of restrained molecular dynamics. Distortion in DNA was found for the major groove mode whereas less significant changes for both ligand and DNA were detected for the minor groove binding after molecular dynamics simulation. The conformation of the ligand obtained from the major groove mode resembles that computed with the ligand soaked in water. The van der Waals contact energy was found to be as significant as electrostatic energy and more important for difference in binding energy between these two binding modes. The importance of van der Waals force in groove binding was supported by computations on the complex formed by the repressor peptide fragment from the bacteriophage 434 and its operator oligonucleotide.     

Proton-linked contributions to site-specific interactions of lambda cI repressor and OR

The effects of proton activity on the site-specific interactions of cI repressors with operator sites OR were studied by using DNase I footprint titration. Individual-site binding isotherms were obtained for the binding of repressor to each site of wild-type OR and of mutant operators in which binding to some sites is eliminated. The Gibbs energies for binding and for cooperativity (in every operator configuration) were determined at each pH (range 5-8). The proton-linked effects clearly account for a significant fraction of the difference in affinities for the three operator sites. The most dramatic effects on the repressor-operator binding interactions are at acid pH, and therefore do not involve the basic groups in the repressor N-terminal arm known to contact the DNA. Also, the proton-linked effects are different at the three operator sites as indicated by significantly different derivative relationships, partial derivative of ln k versus partial derivative of ln aH = net proton absorption (delta nu bar(H)). These results implicate ionizable repressor groups which may not contact the DNA and conformational differences between the three repressor-operator site complexes as being important components to the mechanism of site specificity. The extensive data base generated by these studies was also used to reevaluate the traditional models used to describe cooperativity in this system. The results confirm the lack of significant cooperative interaction between OR1 and OR3 at all conditions. However, the data for some experimental conditions are clearly inconsistent with the (selection) rule, that cooperative interaction between OR2 and OR3 is eliminated by ligation at OR1.     

Design and synthesis of peptides capable of specific binding to DNA

In the present communication, design, synthesis and DNA binding activities of the following two peptides are reported: Dns-Gly-Ala-Gln-Lys-Leu-Ala-Cly-Lys-Val-Gly-Thr-Lys-Val-Lys-Val-Gl y-Thr-Lys-Thr - Val-OH (I) and [(H-Ala-Lys-Leu-Ala-Thr-Lys-Ala-Gly-Val-Lys-Gln-Gln-Ser-Ile-Gln-Leu-Ile- Thr- Ala-Aca-Lys-Aca)2Lys-Aca]2Lys-Val-OH (II), where Aca = NH(CH2)5CO--; Dns is a residue of 5-dimethylaminonaphtalene-1-sulfonic acid. Peptide I contains a large fraction (ca.30%) of valyl and threonyl residues, which possess a high potential for beta structure formation. Peptide II contains four repeats of the amino acid sequence present in the presumed DNA binding helix-turn-helix unit of 434 Cro repressor. These four domains are linked in such a way that two domains can interact with two halves a 14 base pair long operator site on DNA. From CD studies we have found that peptide I is in a random coil conformation in the aqueous solution in the presence of 20% trifluoroethanol. By contrast, amino acid residues of peptide II assume alpha helical, beta and random coiled conformations under the same conditions. A change in the secondary structure of the two peptides upon binding to DNA is observed. The difference CD spectra obtained by subtracting the spectra of free DNA from the spectra of peptide I--DNA complexes gives rise to a beta-like pattern. The difference CD spectra obtained for complexes of peptide II with various natural and synthetic DNAs suggest that alpha-beta-transition takes place in the presumed helix-turn-helix repeat units of peptide II upon binding to DNA. Peptide I binds more strongly to poly(dG).poly(dC) than to poly(dA).poly(dT) and poly[d(GC)].poly[d(GC)]. The binding takes place in the minor DNA groove because minor groove binding antibiotic sibiromycin can displace peptide I from a complex with poly(dG).poly(dC). Analysis of footprinting diagramms shows that peptide I specifically protects phosphodiester bonds within operator sites OR1 and OR2 of phage lambda from nuclease cleavage. By contrast, peptide II does not react specifically with operators OR1, OR2 and OR3 of phage 434 although it forms very tight complexes with DNA which are stable in the presence of 1M NH4F.     

Dimerization specificity of P22 and 434 repressors is determined by multiple polypeptide segments.

The repressor protein of bacteriophage P22 binds to DNA as a homodimer. This dimerization is absolutely required for DNA binding. Dimerization is mediated by interactions between amino acids in the carboxyl (C)-terminal domain. We have constructed a plasmid, p22CT-1, which directs the overproduction of just the C-terminal domain of the P22 repressor (P22CT-1). Addition of P22CT-1 to DNA-bound P22 repressor causes the dissociation of the complex. Cross-linking experiments show that P22CT-1 forms specific heterodimers with the intact P22 repressor protein, indicating that inhibition of P22 repressor DNA binding by P22CT-1 is mediated by the formation of DNA binding-inactive P22 repressor:P22CT-1 heterodimers. We have taken advantage of the highly conserved amino acid sequences within the C-terminal domains of the P22 and 434 repressors and have created chimeric proteins to help identify amino acid regions required for dimerization specificity. Our results indicate that the dimerization specificity region of these proteins is concentrated in three segments of amino acid sequence that are spread across the C-terminal domain of each of the two phage repressors. We also show that the set of amino acids that forms the cooperativity interface of the P22 repressor may be distinct from those that form its dimer interface. Furthermore, cooperativity studies of the wild-type and chimeric proteins suggest that the location of cooperativity interface in the 434 repressor may also be distinct from that of its dimerization interface. Interestingly, changes in the dimer interface decreases the ability of the 434 repressor to discriminate between its wild-type binding sites, O(R)1, O(R)2, and O(R)3. Since 434 repressor discrimination between these sites depends in large part on the ability of this protein to recognize sequence-specific differences in DNA structure and flexibility, this result indicates that the C-terminal domain is intimately involved in the recognition of sequence-dependent differences in DNA structure and flexibility.     

Site-specific enthalpic regulation of DNA transcription at bacteriophage lambda OR.

Binding of cI repressor to DNA fragments containing the three specific binding sites of the right operator (OR) of bacteriophage lambda was studied in vitro over the temperature range 5-37 degrees C by quantitative footprint titration. The individual-site isotherms, obtained for binding repressor dimers to each site of wild-type OR and to appropriate mutant operator templates, were analyzed for the Gibbs energies of intrinsic binding and pairwise cooperative interactions. It is found that dimer affinity for each of the three sites varies inversely with temperature, i.e., the binding reactions are enthalpy driven, unlike many protein-DNA reactions. By contrast, the magnitude of the pairwise cooperativity terms describing interaction between adjacently site-bound repressor dimers is quite small. This result in combination with the recent finding that repressor monomer-dimer assembly is highly enthalpy driven (with delta H degrees = -16 kcal mol-1) [Koblan, K. S., & Ackers, G. K. (1991) Biochemistry 30, 7817-7821] indicates that the associative contacts between site-bound repressors that mediate cooperativity are unlikely to be the same as those responsible for dimerization. The intrinsic binding enthalpies for all three sites are negative (exothermic) and nearly temperature-invariant, indicating no heat capacity changes on the scale of those inferred in other protein-DNA systems. However, the three operator sites are affected differentially by temperature: the intrinsic binding free energies for sites OR1 and OR3 change in parallel over the entire range, delta H0OR1 = -23.3 +/- 4.0 kcal mol-1 and delta H0OR3 = -22.7 +/- 1.2 kcal mol-1.(ABSTRACT TRUNCATED AT 250 WORDS)     

Protein dynamics and monomer-monomer interactions in AntR activation by electron paramagnetic resonance and double electron-electron resonance

The Anthracis repressor (AntR) is a Mn(II)-activated DNA binding protein that is involved in the regulation of Mn(II) homeostasis in Bacillus anthracis. AntR is structurally and functionally homologous to Mn(II)-activated repressor from Bacillus subtillis (MntR). Our studies on AntR focus on metal-regulated activation of the protein. Line shape analysis of continuous wave electron paramagnetic resonance (EPR) spectra showed that metal binding resulted in a general reduction of backbone dynamics and that there were no further changes in backbone motion upon DNA binding. Double electron-electron resonance (DEER) pulsed EPR spectroscopy was used to measure distances between nitroxide spin labels strategically placed in dimeric AntR. The DEER data were analyzed assuming Gaussian distributions for discrete populations of spins. A structural model for AntR was built from homology to MntR, and the experimentally measured distances were simulated to distinguish between spin label and backbone motions. Together with the computational analysis, the DEER results for apo-AntR indicated relatively narrow conformational distributions for backbone residues at the dimer interface and near the metal binding site. No significant changes were observed on these sites in the presence of metal or DNA. On the other hand, the distribution of the conformers and the distances between the putative DNA binding helices decreased upon metal binding. These results suggest that the DNA binding region of AntR shows large amplitude backbone motions in the absence of metal, which may preclude sequence-specific binding to promoter sites. Metal binding narrows the range of conformations accessible in this region and shortens the mean distance between the DNA binding helices, probably resulting in alignment that optimizes promoter recognition and binding.