Nearest-neighbor non-additivity versus long-range non-additivity in TATA-box structure and its implications for TBP-binding mechanism

TBP recognizes its target sites, TATA boxes, by recognizing their sequence-dependent structure and flexibility. Studying this mode of TATA-box recognition, termed ‘indirect readout’, is important for elucidating the binding mechanism in this system, as well as for developing methods to locate new binding sites in genomic DNA. We determined the binding stability and TBP-induced TATA-box bending for consensus-like TATA boxes. In addition, we calculated the individual information score of all studied sequences. We show that various non-additive effects exist in TATA boxes, dependent on their structural properties. By several criterions, we divide TATA boxes to two main groups. The first group contains sequences with 3–4 consecutive adenines. Sequences in this group have a rigid context-independent cooperative structure, best described by a nearest-neighbor non-additive model. Sequences in the second group have a flexible, context-dependent conformation, which cannot be described by an additive model or by a nearest-neighbor non-additive model. Classifying TATA boxes by these and other structural rules clarifies the different recognition pathways and binding mechanisms used by TBP upon binding to different TATA boxes. We discuss the structural and evolutionary sources of the difficulties in predicting new binding sites by probabilistic weight-matrix methods for proteins in which indirect readout is dominant.     

Molecular dynamics simulations of papilloma virus E2 DNA sequences: dynamical models for oligonucleotide structures in solution

The specificity of papilloma virus E2 protein-DNA binding depends critically upon the sequence of a region of the DNA not in direct contact with the protein, and represents one of the simplest known examples of indirect readout. A detailed characterization of this system in solution is important to the further investigation hypothesis of a structural code for DNA recognition by regulatory proteins. In the crystalline state, the E2 DNA oligonucleotide sequence, d(ACCGAATTCGGT), exhibits three different structural forms. We report herein studies of the structure of E2 DNA in solution based on a series of molecular dynamics (MD) simulations including counterions and water, utilizing both the canonical and various crystallographic structures as initial points of departure. All MDs converged on a single dynamical structure of d(ACCGAATTCGGT) in solution. The predicted structure is in close accord with two of the three crystal structures, and indicates that a significant kink in the double helix at the central ApT step in the other crystal molecule may be a packing effect. The dynamical fine structure was analyzed on the basis of helicoidal parameters. The calculated curvature in the sequence was found to originate primarily from YPR steps in the regions flanking the central AATT tract. In order to study the role of structural adaptation of the DNA in the binding process, a subsequent simulation on the 16-mer cognate sequence d(CAACCGAATTCGGTTG) was initiated from the crystallographic coordinates of the bound DNA in the crystal structure of the protein DNA complex. MD simulations starting with the protein-bound form relaxed rapidly back to the dynamical structure predicted from the previous simulations on the uncomplexed DNA. The MD results show that the bound form E2 DNA is a dynamically unstable structure in the absence of protein, and arises as a consequence of both structural changes intrinsic to the sequence and induced by the interaction with protein.     

Indirect Readout of DNA Sequence by P22 Repressor: Roles of DNA and Protein Functional Groups in Modulating DNA Conformation

The repressor of bacteriophage P22 (P22R) discriminates between its various DNA binding sites by sensing the identity of non-contacted base pairs at the center of its binding site. The “indirect readout” of these non-contacted bases is apparently based on DNA's sequence-dependent conformational preferences. The structures of P22R–DNA complexes indicate that the non-contacted base pairs at the center of the binding site are in the B′ state. This finding suggests that indirect readout and therefore binding site discrimination depend on P22R's ability to either sense and/or impose the B′ state on the non-contacted bases of its binding sites. We show here that the affinity of binding sites for P22R depends on the tendency of the central bases to assume the B′-DNA state. Furthermore, we identify functional groups in the minor groove of the non-contacted bases as the essential modulators of indirect readout by P22R. In P22R–DNA complexes, the negatively charged E44 and E48 residues are provocatively positioned near the negatively charged DNA phosphates of the non-contacted nucleotides. The close proximity of the negatively charged groups on protein and DNA suggests that electrostatics may play a key role in the indirect readout process. Changing either of two negatively charged residues to uncharged residues eliminates the ability of P22R to impose structural changes on DNA and to recognize non-contacted base sequence. These findings suggest that these negatively charged amino acids function to force the P22R-bound DNA into the B′ state and therefore play a key role in indirect readout by P22R.     

Remodeling of the human papillomavirus type 11 replication origin into discrete nucleoprotein particles and looped structures by the E2 protein

The human papillomavirus (HPV) DNA replication origin (ori) shares a common theme with many DNA control elements in having multiple binding sites for one or more proteins spaced over several hundreds of base pairs. The HPV type 11 ori spans 103 bp and contains three palindromic E2 binding sites (E2BS-2, E2BS-3, and E2BS-4) for the dimeric E2 ori binding protein. These sites are separated by 64 and 3 bp. E2BS-1 is located 288 bp upstream of E2BS-2 and is not required for efficient transient or cell-free replication. In this study, electron microscopy was used to visualize complexes of HPV-11 DNA ori bound by purified E2 protein. DNA containing only E2BS-2 showed a single E2 dimer bound. DNA containing E2BS-3 and E2BS-4 showed two side-by-side E2 dimers, while DNA containing E2BS-2, E2BS-3, and E2BS-4 exhibited a large disk/ring-shaped protein particle bound, indicating that the DNA had been remodeled into a discrete complex, likely containing an E2 hexamer. With all four binding sites present, up to 27% of the DNA molecules were arranged into loops by E2, the majority of which spanned E2BS-1 and one of the other three sites. Studies on the dependence of looping on salt, ATP, and DTT using full-length E2 and an E2 protein containing only the carboxyl-terminal DNA binding and protein dimerization domain suggest that looping is dependent on the N-terminal domain and factors that may affect the manner in which E2 scans DNA for binding sites. The role of these structures in the modeling and regulation of the HPV-11 ori is discussed.     

Genetic analysis of high-risk e6 in episomal maintenance of human papillomavirus genomes in primary human keratinocytes

Papillomaviruses possess small DNA genomes that encode five early (E) proteins. Transient DNA replication requires activities of the E1 and E2 proteins and a DNA segment containing their binding sites. The E6 and E7 proteins of cancer-associated human papillomavirus (HPV) transform cells in culture. Recent reports have shown that E6 and E7 are necessary for episomal maintenance of HPV in primary keratinocytes. The functions of E6 necessary for viral replication have not been determined, and to address this question we used a recently developed transfection system based on HPV31. To utilize a series of HPV16 E6 mutations, HPV31 E6 was replaced by its HPV16 counterpart. This chimeric genome was competent for both transient and stable replication in keratinocytes. Four HPV16 E6 mutations that do not stimulate p53 degradation were unable to support stable viral replication, suggesting this activity may be necessary for episomal maintenance. E7 has also been shown to be essential for episomal maintenance of the HPV31 genome. A point mutation in the Rb binding motif of HPV E7 has been reported to render HPV31 unable to stably replicate. Interestingly, HPV31 genomes harboring two of the three p53 degradation-defective E6 mutations combined with this E7 mutation were maintained as replicating episomes. These findings imply that the balance between E6 and E7 functions in infected cells is critical for episomal maintenance of high-risk HPV genomes. This model will be useful to dissect the activities of E6 and E7 necessary for viral DNA replication.     

A C-Terminal Helicase Domain of the Human Papillomavirus E1 Protein Binds E2 and the DNA Polymerase α-Primase p68 Subunit

The human papillomavirus (HPV) E1 and E2 proteins bind cooperatively to the viral origin of replication (ori), forming an E1-E2-ori complex that is essential for initiation of DNA replication. All other replication proteins, including DNA polymerase α-primase (polα-primase), are derived from the host cell. We have carried out a detailed analysis of the interactions of HPV type 16 (HPV-16) E1 with E2, ori, and the four polα-primase subunits. Deletion analysis showed that a C-terminal region of E1 (amino acids [aa] 432 to 583 or 617) is required for E2 binding. HPV-16 E1 was unable to bind the ori in the absence of E2, but the same C-terminal domain of E1 was sufficient to tether E1 to the ori via E2. Of the polα-primase subunits, only p68 bound E1, and binding was competitive with E2. The E1 region required (aa 397 to 583) was the same as that required for E2 binding but additionally contained 34 N-terminal residues. In confirmation of these differences, we found that a monoclonal antibody, mapping adjacent to the N-terminal junction of the p68-binding region, blocked E1-p68 but not E1-E2 binding. Sequence alignments and secondary-structure prediction for HPV-16 E1 and other superfamily 3 (SF3) viral helicases closely parallel the mapping data in suggesting that aa 439 to 623 constitute a discrete helicase domain. Assuming a common nucleoside triphosphate-binding fold, we have generated a structural model of this domain based on the X-ray structures of the hepatitis C virus and Bacillus stearothermophilus (SF2) helicases. The modelling closely matches the deletion analysis in suggesting that this region of E1 is indeed a structural domain, and our results suggest that it is multifunctional and critical to several stages of HPV DNA replication.     

Signals for TBP/TATA box recognition

The TATA box-binding protein (TBP) recognizes its target sites (TATA boxes) by indirectly reading the DNA sequence through its conformation effects (indirect readout). Here, we explore the molecular mechanisms underlying indirect readout of TATA boxes by TBP by studying the binding of TBP to adenovirus major late promoter (AdMLP) sequence variants, including alterations inside as well as in the sequences flanking the TATA box. We measure here the dissociation kinetics of complexes of TBP with AdMLP targets and, by phase-sensitive assay, the intrinsic bending in the TATA box sequences as well as the bending of the same sequence induced by TBP binding. In these experiments we observe a correlation of the kinetic stability to sequence changes within the TATA recognition elements. Comparison of the kinetic data with structural properties of TATA boxes in known crystalline TBP/TATA box complexes reveals several "signals" for TATA box recognition, which are both on the single base-pair level, as well as larger DNA tracts within the TATA recognition element. The DNA bending induced by TBP on its binding sites is not correlated to the stability of TBP/TATA box complexes. Moreover, we observe a significant influence on the kinetic stability of alteration in the region flanking the TATA box. This effect is limited however to target sites with alternating TA sequences, whereas the AdMLP target, containing an A tract, is not influenced by these changes.     

The E2 binding sites determine the efficiency of replication for the origin of human papillomavirus type 18.

Human papillomaviruses (HPV-s) have been shown to possess transforming and immortalizing activity for many different, mainly keratinocyte cell lines and they have been detected in 90% of anogenital cancer tissues, which suggests a causative role in the induction of anogenital and other tumours. We have exploited a quantitative assay to identify and characterize the origin of replication of the human papillomavirus type 18 (HPV-18), one of the most prevalent types in the high-risk HPV group. Replication of HPV origin fragments was studied transiently by cotransfection with a protein expression vector providing replication proteins E1 and E2. We have localized the HPV-18 origin to nucleotides 7767-119. This region contains three E2 binding sites and an essential A/T rich DNA region (nucleotides 9-35) that is partly homologous to the E1 binding site found in bovine papillomavirus type 1 (BPV-1) genome. At least one of the three E2 binding sites was absolutely required for origin function; addition of other E2 sites had cooperative stimulating effect. This is the first quantitative analysis of the E2 binding sites for papillomavirus replication.     

X-ray crystal structures of half the human papilloma virus E2 binding site: d(GACCGCGGTC).

The X-ray crystal structure of the DNA decamer d(GACCGCGGTC), containing half the human papilloma virus E2 binding site, has been solved from two crystals grown at different ionic conditions (50 mM MgCl2and 50 mM spermine or 1.56 mM MgCl2and 1.56 mM spermine). Despite the variation in salt concentration, the two DNA structures are in a very similar, A-type DNA conformation, with helical axes curving towards the major groove. Although the salt concentrations do not effect the helical parameters or hydration to a large degree, there is a change in the overall helical curvature; 18 degrees and 31 degrees for the low and high salt structures, respectively. This curvature appears to be sequence specific and biologically relevant when compared with similar DNA structures, including the E2 binding site of a protein-DNA complex.     

The M.EcoRV DNA-(adenine N6)-methyltransferase uses DNA bending for recognition of an expanded EcoDam recognition site

The M.EcoRV DNA methyltransferase recognizes GATATC sites. It is related to EcoDam, which methylates GATC sites. The DNA binding domain of M.EcoRV is similar to that of EcoDam suggesting a similar mechanism of DNA recognition. We show that amino acid residue Lys11 of M.EcoRV is involved in recognition of Gua1 and Arg128 contacts the Gua in base pair 6. These residues correspond to Lys9 and Arg124 in EcoDam, which recognize the Gua residues in both strands of the Dam recognition sequence, indicating that M.EcoRV and EcoDam make similar contacts to outermost base pairs of their recognition sequences and M.EcoRV recognizes its target site as an expanded GATC site. In contrast to EcoDam, M.EcoRV considerably bends the DNA (59+/-4 degrees) suggesting indirect readout of the AT-rich inner sequence. Recognition of an expanded target site by DNA bending is a new principle for changing DNA recognition specificity of proteins during molecular evolution. R128A is inefficient in DNA bending and binding, whereas K11A bends DNA with relaxed sequence specificity. These results suggest a temporal order of the formation of protein-DNA contacts in which the Gua6-Arg128 contact forms early followed by DNA bending and, finally, the formation of the Lys11-Gua1 contact.     

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.     

Indirect DNA Readout on the Protein Side: Coupling between Histidine Protonation, Global Structural Cooperativity, Dynamics, and DNA Binding of the Human Papillomavirus Type 16 E2C Domain

DNA sequence recognition by the homodimeric C-terminal domain of the human papillomavirus type 16 E2 protein (E2C) is known to involve both direct readout and DNA-dependent indirect readout mechanisms, while protein-dependent indirect readout has been deduced but not directly observed. We have investigated coupling between specific DNA binding and the dynamics of the unusual E2C fold, using pH as an external variable. Nuclear magnetic resonance and isothermal titration calorimetry show that pH titration of His318 in the complex interface and His288 in the core of the domain is coupled to both binding and the dynamics of the β-barrel core of E2C, with a tradeoff between dimer stability and function. Specific DNA binding is, in turn, coupled to the slow dynamics and amide hydrogen exchange in the entire β-barrel, reaching residues far apart from the DNA recognition elements but not affecting the two helices of each monomer. The changes are largest in the dimerization interface, suggesting that the E2C β-barrel acts as a hinge that regulates the relative position of the DNA recognition helices. In conclusion, the cooperative dynamics of the human papillomavirus type 16 E2C β-barrel is coupled to sequence recognition in a protein-dependent indirect readout mechanism. The patterns of residue substitution in genital papillomaviruses support the importance of the protonation states of His288 and His318 and suggest that protein-dependent indirect readout and histidine pH titration may regulate DNA binding in the cell.     

Human Papillomavirus DNA Replication Compartments in a Transient DNA Replication System

Many DNA viruses replicate their genomes at nuclear foci in infected cells. Using indirect immunofluorescence in combination with fluorescence in situ hybridization, we colocalized the human papillomavirus (HPV) replicating proteins E1 and E2 and the replicating origin-containing plasmid to nuclear foci in transiently transfected cells. The host replication protein A (RP-A) was also colocalized to these foci. These nuclear structures were identified as active sites of viral DNA synthesis by bromodeoxyuridine (BrdU) pulse-labeling. Unexpectedly, the great majority of RP-A and BrdU incorporation was found in these HPV replication domains. Furthermore, E1, E2, and RP-A were also colocalized to nuclear foci in the absence of an origin-containing plasmid. These observations suggest a spatial reorganization of the host DNA replication machinery upon HPV DNA replication or E1 and E2 expression. Alternatively, viral DNA replication might be targeted to host nuclear domains that are active during the late S phase, when such domains are limited in number. In a fraction of cells expressing E1 and E2, the promyelocytic leukemia protein, a component of nuclear domain 10 (ND10), was either partially or completely colocalized with E1 and E2. Since ND10 structures were recently hypothesized to be sites of bovine papillomavirus virion assembly, our observation suggests that HPV DNA amplification might be partially coupled to virion assembly.