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

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

Influence of DNA stiffness in protein-DNA recognition

Protein-DNA recognition plays an essential role in the regulation of gene expression. The protein-DNA binding specificity is based on direct atomic contacts between protein and DNA and/or the conformational properties of DNA. In this work, we have analyzed the influence of DNA stiffness (E) to the specificity of protein-DNA complexes. The average DNA stiffness parameters for several protein-DNA complexes have been computed using the structure based sequence dependent stiffness scale. The relationship between DNA stiffness and experimental protein-DNA binding specificity has been brought out. We have investigated the importance of DNA stiffness with the aid of experimental free energy changes (DeltaDeltaG) due to binding in several protein-DNA complexes, such as, ETS proteins, 434, lambda, Mnt and trp repressors, 434 cro protein, EcoRV endonuclease V and zinc fingers. We found a correlation in the range 0.65-0.97 between DeltaDeltaG and E in these examples. Further, we have qualitatively analyzed the effect of mutations in the target sequence of lambda repressor and we observed that the DNA stiffness could correctly identify 70% of the correct bases among the considered nine positions.     

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.     

Fine structure and activity of discrete RAG-HMG complexes on V(D)J recombination signals

Two lymphoid cell-specific proteins, RAG-1 and RAG-2, initiate V(D)J recombination by introducing DNA breaks at recombination signal sequences (RSSs). Although the RAG proteins themselves bind and cleave DNA substrates containing either a 12-RSS or a 23-RSS, DNA-bending proteins HMG-1 and HMG-2 are known to promote these processes, particularly with 23-RSS substrates. Using in-gel cleavage assays and DNA footprinting techniques, I analyzed the catalytic activity and protein-DNA contacts in discrete 12-RSS and 23-RSS complexes containing the RAG proteins and either HMG-1 or HMG-2. I found that both the cleavage activity and the pattern of protein-DNA contacts in RAG-HMG complexes assembled on 12-RSS substrates closely resembled those obtained from analogous 12-RSS complexes lacking HMG protein. In contrast, 23-RSS complexes containing both RAG proteins and either HMG-1 or HMG-2 exhibited enhanced cleavage activity and displayed an altered distribution of cleavage products compared to 23-RSS complexes containing only RAG-1 and RAG-2. Moreover, HMG-dependent heptamer contacts in 23-RSS complexes were observed. The protein-DNA contacts in RAG-RSS-HMG complexes assembled on 12-RSS or 23-RSS substrates were strikingly similar at comparable positions, suggesting that the RAG proteins mediate HMG-dependent heptamer contacts in 23-RSS complexes. Results of ethylation interference experiments suggest that the HMG protein is positioned 5' of the nonamer in 23-RSS complexes, interacting largely with the side of the duplex opposite the one contacting the RAG proteins. Thus, HMG protein plays the dual role of bringing critical elements of the 23-RSS heptamer into the same phase as the 12-RSS to promote RAG binding and assisting in the catalysis of 23-RSS cleavage.     

Geometric analysis and comparison of protein-DNA interfaces: why is there no simple code for recognition?

Structural studies of protein-DNA complexes have shown that there are many distinct families of DNA-binding proteins, and have shown that there is no simple "code" describing side-chain/base interactions. However, systematic analysis and comparison of protein-DNA complexes has been complicated by the diversity of observed contacts, the sheer number of complexes currently available and the absence of any consistent method of comparison that retains detailed structural information about the protein-DNA interface. To address these problems, we have developed geometric methods for characterizing the local structural environment in which particular side-chain/base interactions are observed. In particular, we develop methods for analyzing and comparing spatial relationships at the protein-DNA interface. Our method involves attaching local coordinate systems to the DNA bases and to the C(alpha) atoms of the peptide backbone (these are relatively rigid structural units). We use these tools to consider how the position and orientation of the polypeptide backbone (with respect to the DNA) helps to determine what contacts are possible at any given position in a protein-DNA complex. Here, we focus on base contacts that are made in the major groove, and we use spatial relationships in analyzing: (i) the observed patterns of side-chain/base interactions; (ii) observed helix docking orientations; (iii) family/subfamily relationships among DNA-binding proteins; and (iv) broader questions about evolution, altered specificity mutants and the limits for the design of new DNA-binding proteins. Our analysis, which highlights differences in spatial relationships in different complexes and at different positions in a complex, helps explain why there is no simple, general code for protein-DNA recognition.     

Intermolecular and intramolecular readout mechanisms in protein-DNA recognition

Protein-DNA recognition plays an essential role in the regulation of gene expression. Regulatory proteins are known to recognize specific DNA sequences directly through atomic contacts (intermolecular readout) and/or indirectly through the conformational properties of the DNA (intramolecular readout). However, little is known about the respective contributions made by these so-called direct and indirect readout mechanisms. We addressed this question by making use of information extracted from a structural database containing many protein-DNA complexes. We quantified the specificity of intermolecular (direct) readout by statistical analysis of base-amino acid interactions within protein-DNA complexes. The specificity of the intramolecular (indirect) readout due to DNA was quantified by statistical analysis of the sequence-dependent DNA conformation. Systematic comparison of these specificities in a large number of protein-DNA complexes revealed that both intermolecular and intramolecular readouts contribute to the specificity of protein-DNA recognition, and that their relative contributions vary depending upon the protein-DNA complexes. We demonstrated that combination of the intermolecular and intramolecular energies derived from the statistical analyses lead to enhanced specificity, and that the combined energy could explain experimental data on binding affinity changes caused by base mutations. These results provided new insight into the relationship between specificity and structure in the process of protein-DNA recognition, which would lead to prediction of specific protein-DNA binding sites.     

Meiotic studies of translocations causing male sterility in the mouse. II. Double heterozygotes for Robertsonian translocations.

Unusual meiotic behavior of the XY chromosome pair was observed in sterile male mice doubly heterozygous for two Robertsonian translocations, Rb(16.17)7Bnr and Rb(8.17)1Iem. Nonrandom association between the X chromosome and the translocation configuration, ascertained from the frequencies of relevant C-band contacts, was found in 9 of 10 sterile males. Besides the nonrandom association, the XY chromosomes showed signs of impaired condensation, as judged by measurement of their lengths at diakinesis/MI of the first meiotic division. In contrast, neither nonrandom contact nor decondensation of the XY chromosomes pair was found in fertile males heterozygous for a single Robertsonian translocation, Rb1Iem or Rb7Bnr. The present observations lend indirect support to the working hypothesis advanced previously, the assumption that interference with X-chromosome inactivation is a possible cause of spermatogenic breakdown in carriers of various male-sterile chromosomal transloations. Alternative explanations of the available data, which cannot be ruled out, are briefly discussed.     

Structure model of a complex between the factor for inversion stimulation (FIS) and DNA: Modeling protein-DNA complexes with dyad symmetry and known protein structures

A method is presented to predict overall conformations of protein-DNA complexes on the basis of the known three-dimensional structures of the proteins. The method is restricted to proteins with a common twofold symmetry axis, which show only minor conformational changes upon binding to DNA. The method uses a numerical finite difference solution of the linearized Poisson-Boltzmann equation and subsequent energy minimization cycles. Structural parameters—the rotation angle of the DNA relative to the protein around the common symmetry axis, the protein-DNA distance, and intermolecular hydrogen-bonding contacts—are presented for two test cases, DNA bound to CAP (catabolite gene activator protein) and to the Cro-repressor of bacteriophage 434. The DNA curvature in the starting model of the docking procedure was chosen as a smoothed approximation of the conformation found in the X-ray structures of these complexes. The method is further used to predict the unknown structure of the complex between the factor for inversion stimulation (FIS) and DNA, which is bent upon binding to FIS. In contrast to the test cases, the unknown curvature of the starting model is derived from a calibration of electrostatic precalculations for different proteins according to crystallographically observed DNA bending. The results of the modeling are in good accordance with the experimentally observed overall structure of protein-DNA complexes for the two test cases; for FIS, they correspond to several of the experimentally proposed protein-DNA contacts. © 1996 Wiley-Liss, Inc.     

Assembly and self-association of oxytricha telomeric nucleoprotein complexes

Two types of specific telomeric protein-DNA complex are reconstituted upon incubation of purified Oxytricha telomere protein with (T4G4)4, an oligodeoxynucleotide of telomeric sequence. The complexes differ in electrophoretic mobility, in protein-DNA contacts, and in the rate of DNA exchange. The patterns of protein-DNA interaction determined by modification interference suggest a model in which the protein can bind either to the two T4G4 repeats at the 3' end or to two internal repeats; in the latter case, it can make a different set of contacts with the terminal repeat to form the more stable complex. Native telomeric chromatin isolated from Oxytricha contains both types of complexes. The reconstituted monomeric complexes associate to give a high molecular weight form that has an altered chemical footprint. Such interactions may mediate the association of chromosomal telomeres in vivo.     

Positioning the 3'-DNA terminus for topoisomerase II-mediated religation

Despite the importance of the topoisomerase II DNA cleavage/rejoining cycle to genomic integrity, the mechanistic details of religation are poorly understood. Topoisomerase II utilizes covalent protein-DNA interactions to align the 5'-termini of cleaved DNA for religation. However, because the enzyme does not form covalent bonds with the 3'-DNA termini, the basis for the alignment of the 3'-ends is less clear. Three major possibilities exist. The 3'-termini may be positioned for religation (i) by base pairing to their complementary DNA strands, (ii) by base stacking to the adjacent residues, or (iii) by noncovalent interactions with topoisomerase II. To distinguish between these possibilities, the ability of human topoisomerase IIalpha to religate a series of oligonucleotides with altered base pairing or base stacking at their 3'-termini was determined. Substrates containing modifications that disrupted terminal base pairing or base stacking with-out affecting the 3'-terminal base were resealed at wild-type rates. In contrast, substrates that lacked the terminal base (or contained an altered base) displayed very low rates of religation. On the basis of these results, we propose that topoisomerase II positions the 3'-DNA termini for religation through noncovalent protein-DNA contacts.     

Determination of the three-dimensional structure of the Mrf2-DNA complex using paramagnetic spin labeling

Understanding the mechanism of protein-DNA interactions at the molecular level is one of the main focuses in structural and molecular biological investigations. At present, NMR spectroscopy is the only approach that can provide atomic details of protein-DNA recognition in solution. However, determining the structures of protein-DNA complexes using NMR spectroscopy has been dependent on the observation of intermolecular nuclear Overhauser effects (NOE) and their assignments, which are difficult to obtain in many cases. In this study, we have shown that intermolecular distance constraints derived from a single spin-label in combination with docking calculations have defined many specific contacts of the complex between the AT-rich interaction domain (ARID) of Mrf2 and its target DNA. Mrf2 contacts DNA mainly using the two flexible loops, L1 and L2. While the L1 loop contacts the phosphate backbone, L2 and several residues in the adjacent helices interact with AT base pairs in the major groove of DNA. Despite the structural diversity in the ARID family of DNA-binding proteins, Mrf2 maintains contacts with DNA similar to those observed in the homologous Dri-DNA complex.     

Interactions between lac repressor protein and site-specific bromodeoxyuridine-substituted operator DNA. Ultraviolet footprinting and protein-DNA cross-link formation

Specific contacts between the lac repressor and operator have been explored using 5-bromodeoxyuridine-substituted DNA. Substitution of BrdU for single thymidine positions in a synthetic 40-base pair operator provides substrate for ultraviolet irradiation; upon irradiation, strand scission occurs at the BrdU residues. When bound, lac repressor protein provides protection against UV-induced breakage depending on the nature of the sites and type of interaction. We have confirmed 13 unique sites of inducer-sensitive protection along the operator sequence using this method compared to complete substitution with BrdU; differences were observed at two positions for singly substituted versus completely substituted DNAs (Ogata, R., and Gilbert, W. (1977) Proc. Natl. Acad. Sci. U.S.A. 74, 4973-4976). The ability of these photosensitive DNAs to form short range cross-links to bound protein has been used to determine the efficiency with which cross-linked protein-DNA complexes are generated at each individual site of BrdU substitution. Five sites of high efficiency cross-linking to the repressor protein have been identified. At one site, cross-linking without protection from strand scission was observed; this result suggests an unusual mechanism of strand scission and/or cross-linking at this site. Comparison of the UV protection results and the cross-linking data show that these processes provide complementary tools for identifying and analyzing individual protein-DNA contacts.