The structural basis of DNA target discrimination by papillomavirus E2 proteins

Structural and energetic origins of sequence-specific DNA bending: Monte Carlo simulations of papillomavirus E2-DNA binding sites

DNA bending is an important structural feature for indirect readout in protein-DNA recognition. The binding of papillomavirus E2 transcription factors to their DNA binding sites is associated with DNA bending, providing an attractive model system to study the origins of sequence-specific DNA bending. The consensus E2 target is of the general form ACCGN(4)CGGT with a variable four base pair region. We applied a new all-atom Monte Carlo (MC) algorithm that combines effective sampling with fast conformational equilibration. The resulting MC ensembles resemble the corresponding high-resolution crystal structures very well. Distinct bending is observed for the E2-DNA binding site with a central AATT linker in contrast to an essentially straight DNA with a central ACGT linker. Contributions of specific base pair steps to the overall bending are shown in terms of local structural parameters. The analysis of conformational substates provides new insights into the energetic origins of intrinsic DNA bending.

Molecular dynamics studies on free and bound targets of the bovine papillomavirus type I e2 protein: the protein binding effect on DNA and the recognition mechanism

Molecular dynamics simulations of a total duration of 30 ns in explicit solvent were carried out on the BPV-1-E2 protein complexed to a high-affinity DNA target containing the two hydrogen-bonded ACCG.CGGT half-sites separated by the noncontacted ACGT sequence. The analysis of the trajectories focuses on the DNA structure and on the dynamics. The data are compared to those issued from recent simulations made on three free targets that recognize E2 with different affinities. E2 does not drastically perturb the mechanic properties of the free DNA: the structural relationships between the BI/BII backbone substates and some helical parameters are preserved in the complex despite a severe slowing down of the phosphate group motions. The structures of both free and bound half-sites are very close to each other although the conformational space explored by these regions is narrowed when they are contacted by the protein. The enhanced plasticity found in the best free target spacers, mainly manifested through the backbone motions, allows a clear overlap between several free and bound global DNA features such as the base displacement. Furthermore, this flexibility is preserved in the complex. Our results support the hypothesis that E2 takes advantage of free predistorted structures that may minimize the DNA deformation cost. In addition, we observe that E2 is far from totally stiffening the DNA, suggesting that the entropic penalty inherent in the complex formation could be limited.

Thermodynamic basis for antibody binding to Z-DNA: comparison of a monoclonal antibody and its recombinant derivatives

Antibody engineering represents a promising area in biotechnology. Recombinant antibodies can be easily manipulated generating new ligand and effector activities that can be used as prototype magic bullets. On the other hand, an extensive knowledge of recombinant antibody binding and stability features are essential for an efficient substitution. In this study, we compared the stability and protein binding properties of two recombinant antibody fragments with their parental monoclonal antibody. The recombinant fragments were a monomeric scFv and a dimeric one, harboring human IgG1 CH2-CH3 domains. We have used fluorescence titration quenching to determine the thermodynamics of the interaction between an anti-Z-DNA monoclonal antibody and its recombinant antibody fragments with Z-DNA. All the antibody fragments seemed to bind DNA similarly, in peculiar two-affinity states. Enthalpy-entropy compensation was observed for both affinity states, but a marked entropy difference was observed for the monomeric scFv antibody fragment, mainly for the high affinity binding. In addition, we compared the stability of the dimeric antibody fragment and found differences favoring the monoclonal antibody. These differences seem to derive from the heterologous expression system used.

Human papillomavirus type 18 variant lineages in United States populations characterized by sequence analysis of LCR-E6, E2, and L1 regions

While HPV 16 variant lineages have been well characterized, the knowledge about HPV 18 variants is limited. In this study, HPV 18 nucleotide variations in the E2 hinge region were characterized by sequence analysis in 47 control and 51 tumor specimens. Fifty of these specimens were randomly selected for sequencing of an LCR-E6 segment and 20 samples representative of LCR-E6 and E2 sequence variants were examined across the L1 region. A total of 2770 nucleotides per HPV 18 variant genome were considered in this study. HPV 18 variant nucleotides were linked among all gene segments analyzed and grouped into three main branches: Asian-American (AA), European (E), and African (Af). These three branches were equally distributed among controls and cases and when stratified by Hispanic and non-Hispanic ethnicities. Among invasive cervical cancer cases, no significant differences in the three HPV variant branches were observed among ethnic groups or when stratified by histopathology (squamous vs. adenocarcinoma). The Af branch showed the greatest nucleotide variability when compared to the HPV 18 reference sequence and was more closely related to HPV 45 than either AA or E branches. Our data also characterize nucleotide and amino acid variations in the L1 capsid gene among HPV 18 variants, which may be relevant to vaccine strategies and subsequent studies of naturally occurring HPV 18 variants. Several novel HPV 18 nucleotide variations were identified in this study.

Free energy contributions to direct readout of a DNA sequence

The energetic contributions of individual DNA-contacting side chains to specific DNA recognition in the human papillomavirus 16 E2C-DNA complex is small (less than 1.0 kcal mol(-1)), independent of the physical and chemical nature of the interaction, and is strictly additive. The sum of the individual contributions differs 1.0 kcal mol(-1) from the binding energy of the wild-type protein. This difference corresponds to the contribution from the deformability of the DNA, known as "indirect readout." Thus, we can dissect the energetic contribution to DNA binding into 90% direct and 10% indirect readout components. The lack of high energy interactions indicates the absence of "hot spots," such as those found in protein-protein interfaces. These results are compatible with a highly dynamic and "wet" protein-DNA interface, yet highly specific and tight, where individual interactions are constantly being formed and broken.

DNA bending by M.EcoKI methyltransferase is coupled to nucleotide flipping

The maintenance methyltransferase M.EcoKI recognizes the bipartite DNA sequence 5′-AACNNNNNNGTGC-3′, where N is any nucleotide. M.EcoKI preferentially methylates a sequence already containing a methylated adenine at or complementary to the underlined bases in the sequence. We find that the introduction of a single-stranded gap in the middle of the non-specific spacer, of up to 4 nt in length, does not reduce the binding affinity of M.EcoKI despite the removal of non-sequence-specific contacts between the protein and the DNA phosphate backbone. Surprisingly, binding affinity is enhanced in a manner predicted by simple polymer models of DNA flexibility. However, the activity of the enzyme declines to zero once the single-stranded region reaches 4 nt in length. This indicates that the recognition of methylation of the DNA is communicated between the two methylation targets not only through the protein structure but also through the DNA structure. Furthermore, methylation recognition requires base flipping in which the bases targeted for methylation are swung out of the DNA helix into the enzyme. By using 2-aminopurine fluorescence as the base flipping probe we find that, although flipping occurs for the intact duplex, no flipping is observed upon introduction of a gap. Our data and polymer model indicate that M.EcoKI bends the non-specific spacer and that the energy stored in a double-stranded bend is utilized to force or flip out the bases. This energy is not stored in gapped duplexes. In this way, M.EcoKI can determine the methylation status of two adenine bases separated by a considerable distance in double-stranded DNA and select the required enzymatic response.

Positive contribution of hydration on DNA binding by E2c protein from papillomavirus

Protein-nucleic acid interactions are responsible for the regulation of key biological events such as genomic transcription and recombination and viral replication. However, the recognition mechanisms involved in these processes are not completely understood. Here, we investigate the dominant forces involved in protein-protein and protein-DNA interactions for the 80-amino-acid C-terminal domain of the E2 protein (E2c) from human papillomavirus (HPV-16). The E2c protein is a homodimer that specifically binds to double-stranded DNA containing the consensus sequence ACCG-N(4)-CGGT, where N is any nucleotide. DNA binding affinity is reduced by lowering water chemical potential, accompanied by an increase in cooperativity. Wyman linkage relations between affinity and water chemical potential indicate that 11 additional water molecules are bound in the formation of the complex between E2c and DNA. Salt dissociation isotherms showed that 10 counterions are released upon association, even at low water activity, indicating that this latter variable does not change the electrostatic component of the interaction. Further analysis demonstrates a strong dependence of cooperativity of binding on the protein concentration. Altogether, these results reveal a novel binding pathway in which the consolidated complex may achieve its final form via a monomer-DNA intermediate, which favors the binding of a second monomer. This molecular mechanism reveals the contributions of multiple conformers in a tight virus genome modulation that seems to be important in the cell infection scenario.

The regulatory E2 proteins of human genital papillomaviruses are pro-apoptotic

Cervical carcinomas are frequently associated with infection by human papillomaviruses (HPVs). These viruses encode two oncogenes E6 and E7, which promote cell proliferation and immortalization. The viral E2 protein represses transcription of the E6/E7 oncogenes and activates viral DNA replication together with the viral E1 helicase. The E2 protein is specifically inactivated in HPV18-associated carcinoma, suggesting that it may prevent carcinogenic progression. Indeed, E2 was shown to exhibit a strong anti-proliferative action when ectopically expressed in cervical carcinoma cells, as it induces both G1 cell cycle arrest and cell death by apoptosis. While the cell cycle arrest is due to E2-mediated transcriptional repression of the viral oncogenes, the induction of apoptosis appears to be an autonomous function of E2. The amino-terminal transactivation domain (TAD) of the E2 protein is required for its pro-apoptotic activity, but transcriptional transactivation is not involved. E2 induces apoptosis through the extrinsic pathway, involving the initiator caspase 8. In addition, E2 is cleaved by caspases during apoptosis, providing an example of an apoptotic inducer, which is itself a target for caspase cleavage. The cleaved E2 protein exhibits an enhanced apoptotic activity, suggesting that it may participate in an amplification loop. This article reviews our current knowledge of the pro-apoptotic activity of the oncogenic papillomavirus E2 proteins, and discusses the implications for the viral vegetative cycle.

The role of DNA structure and dynamics in the recognition of bovine papillomavirus E2 protein target sequences

The papillomavirus E2 transcription and replication factors bind to the DNA consensus ACCGN(4)CGGT sequence (E2-BS), through both direct and indirect readout mechanisms. The two symmetric half-sites ACCG.CGGT are highly conserved in the genomes and are hydrogen bound with E2. Although E2 does not contact the N4 spacer, the affinities are modulated by the base composition of this DNA part. Nevertheless, the origin of either the global recognition mechanism or the spacer effect remains unclear, particularly in the case of the bovine papillomavirus type 1 E2 (BPV-1-E2) system, used as model to study the papillomaviruses. We present, herein, studies carried out on oligomers differently recognized by the BPV-1-E2 protein and based on molecular dynamic simulations including counterions and water. The sequences contain the conserved half-sites but three different spacers (CCAT, ACGT and AAAC), resulting in very high, high and low affinity targets for BPV-1-E2. In order to estimate how much the free DNAs resemble the bound conformations, comparisons are made with two DNAs extracted from E2-BS-BPV-1 crystallographic complexes, representative of high and moderate affinity structures. The analysis of 15 ns trajectories reveals that the ACCG/CGGT half-sites, whatever the spacer, have the same behavior and adopt average stable base-pair parameters very close to those of the bound conformations. In contrast, the three different free spacers strongly differ in their BI <--> BII backbone dynamics. The low affinity AAAC spacer exhibits stable BI backbone conformations, the high affinity ACGT spacer is characterized by a dramatic instability of the CpG phosphate groups, and the CpA and GpG backbones in the very high affinity CCAT.ATGG spacer are trapped in BII conformations. All resemble more of the moderate affinity complex DNA than the high affinity one. Nevertheless, the particular behavior of the CCAT and ACGT backbones allows the emergence of BII-rich spacers, a configuration reproducing both local and global helical features of the bound DNA conformation of the high affinity complex and favoring the minor groove curvature required in the complex. In particular, the CCAT-containing site spends almost half of the time in this form that well mimics the bound one. Thus, we propose that the E2 protein could take advantage of the invariant favorable structures of the half-sites to form a pre-complex, but would require a specific spacer intrinsic malleability to lock the interaction. Finally, the backbone conformational states, by their ability to translate information coded in the sequence into structural properties, provide insight into the mechanisms that contribute to fine binding site selection and specific nucleic acid ligand recognition.

Predicting indirect readout effects in protein-DNA interactions

Recognition of DNA by proteins relies on direct interactions with specific DNA-functional groups, along with indirect effects that reflect variable energetics in the response of DNA sequences to twisting and bending distortions induced by proteins. Predicting indirect readout requires knowledge of the variations in DNA curvature and flexibility in the affected region, which we have determined for a series of DNA-binding sites for the E2 regulatory protein by using the cyclization kinetics method. We examined 16 sites containing different noncontacted spacer sequences, which vary by more than three orders of magnitude in binding affinity. For 15 of these sites, the variation in affinity was predicted within a factor of 3, by using experimental curvature and flexibility values and a statistical mechanical theory. The sole exception was traced to differential magnesium ion binding.

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.

Comparison of the structure and DNA-binding properties of the E2 proteins from an oncogenic and a non-oncogenic human papillomavirus

Human papillomaviruses (HPVS) that infect the genital tract can be divided into two groups: high-risk HPV types, such as HPV 16 and HPV 18, are associated with cancer, low-risk HPV types, such as HPV 6, are associated with benign warts. In both high-risk and low-risk HPV types, the papillomavirus E2 protein binds to four sites within the viral long control region (LCR) and regulates viral gene expression. Here, we present the crystal structure of the minimal DNA-binding domain (DBD) from the HPV 6 E2 protein. We show that the HPV 6 E2 DBD is structurally more similar to the HPV 18 and bovine papillomavirus type 1 (BPV1) E2 proteins than it is to the HPV 16 E2 protein. Using gel retardation assays, we show that the hierarchy of E2 sites within the HPV 16 and HPV 6 LCRs are different. However, despite these differences in structure and site preference, both the HPV 16 and 6 E2 DBDs recognise an extended version of the consensus E2 binding site derived from studies of the BPV1 E2 protein. In both cases, the preferred binding site is 5'AACCGN(4)CGGTT3', where the additional flanking base-pairs are in bold and N(4) represents a four base-pair central spacer. Both of these HPV proteins bind preferentially to E2 sites that contain an A:T-rich central spacer. We show that the preference for an A:T-rich central spacer is due, at least in part, to the need to adopt a DNA conformation that facilitates protein contacts with the flanking base-pairs.

Solution measurement of DNA curvature in papillomavirus E2 binding sites

'Indirect readout' refers to the proposal that proteins can recognize the intrinsic three-dimensional shape or flexibility of a DNA binding sequence apart from direct protein contact with DNA base pairs. The differing affinities of human papillomavirus (HPV) E2 proteins for different E2 binding sites have been proposed to reflect indirect readout. DNA bending has been observed in X-ray structures of E2 protein-DNA complexes. X-ray structures of three different E2 DNA binding sites revealed differences in intrinsic curvature. DNA sites with intrinsic curvature in the direction of protein-induced bending were bound more tightly by E2 proteins, supporting the indirect readout model. We now report solution measurements of intrinsic DNA curvature for three E2 binding sites using a sensitive electrophoretic phasing assay. Measured E2 site curvature agrees well the predictions of a dinucleotide model and supports an indirect readout hypothesis for DNA recognition by HPV E2.

A protein-DNA binding mechanism proceeds through multi-state or two-state parallel pathways

The DNA-binding mechanism of the dimeric C-terminal domain of the papillomavirus E2 protein with its specific DNA target was investigated and shown to proceed through two parallel pathways. A sequential multi-step reaction is initiated by the diffusion-controlled formation of an encounter complex, with no evidence of base sequence discrimination capacity. Following a substantial conformational rearrangement of the protein, a solvent exclusion step leading to the formation of a final protein-DNA complex was identified. This last step involves the largest burial of surface area from the interface and involves the consolidation of the direct readout of the DNA bases. Double-jump stopped-flow experiments allowed us to characterize the sequence of events and demonstrated that a fast-formed consolidated complex can take place through a parallel route. We present the simplest model for the overall mechanism with a description of all the intermediate species in energetic terms.

Inhibition of human papilloma virus E2 DNA binding protein by covalently linked polyamides

Polyamides are a class of heterocyclic small molecules with the potential of controlling gene expression by binding to the minor groove of DNA in a sequence-specific manner. To evaluate the feasibility of this class of compounds as antiviral therapeutics, molecules were designed to essential sequence elements occurring numerous times in the HPV genome. This sequence element is bound by a virus-encoded transcription and replication factor E2, which binds to a 12 bp recognition site as a homodimeric protein. Here, we take advantage of polyamide:DNA and E2:DNA co-crystal structural information and advances in polyamide synthetic chemistry to design tandem hairpin polyamides that are capable of displacing the major groove-binding E2 homodimer from its DNA binding site. The binding of tandem hairpin polyamides and the E2 DNA binding protein to the DNA site is mutually exclusive even though the two ligands occupy opposite faces of the DNA double helix. We show with circular permutation studies that the tandem hairpin polyamide prevents the intrinsic bending of the E2 DNA site important for binding of the protein. Taken together, these results illustrate the feasibility of inhibiting the binding of homodimeric, major groove-binding transcription factors by altering the local DNA geometry using minor groove-binding tandem hairpin polyamides.

Antibody response to a viral transcriptional regulator

The E2 transcriptional activator of the human papillomavirus regulates the expression of most viral transcripts. Its binding to specific target DNA sequences involves large conformational changes in the interacting macromolecules. The high stability of the E2:DNA complex prompted us to analyze the role of macromolecular interactions and adjuvant emulsions in the appearance of conformation-specific antibodies. We demonstrate that immunization with free or DNA-complexed E2 emulsified in an oil-in-water adjuvant elicits a humoral response shifted to the recognition of discontinuous epitopes. Epitope mapping and functional analysis of the generated anti-E2 mAbs reveals that two separate antibodies populations can be obtained: those able to form a stable ternary complex with protein and DNA, and those which recognize the DNA-binding surface of the transcription factor, interfering with E2 binding to DNA.

Transcriptional activity among high and low risk human papillomavirus E2 proteins correlates with E2 DNA binding

The full-length E2 protein, encoded by human papillomaviruses (HPVs), is a sequence-specific transcription factor found in all HPVs, including cancer-causing high risk HPV types 16 and 18 and wart-inducing low risk HPV types 6 and 11. To investigate whether E2 proteins encoded by high risk HPVs may function differentially from E2 proteins encoded by low risk HPVs and animal papillomaviruses, we conducted comparative DNA-binding and transcription studies using electrophoretic mobility shift assays and cell-free transcription systems reconstituted with purified general transcription factors, cofactor, RNA polymerase II, and with E2 proteins encoded by HPV-16, HPV-18, HPV-11, and bovine papillomavirus type 1 (BPV-1). We found that although different types of E2 proteins all exhibited transactivation and repression activities, depending on the sequence context of the E2-binding sites, HPV-16 E2 shows stronger transcription activity and greater DNA-binding affinity than those displayed by the other E2 proteins. Surprisingly, HPV-18 E2 behaves more similarly to BPV-1 E2 than HPV-16 E2 in its functional properties. Our studies thus categorize HPV-18 E2 and BPV-1 E2 in the same protein family, a finding consistent with the available E2 structural data that separate the closely related HPV-16 and HPV-18 E2 proteins but classify together the more divergent BPV-1 and HPV-18 E2 proteins.

The papillomavirus E2 proteins: structure, function, and biology

Nearly twenty years after the first high-resolution crystal structures of specific protein-DNA complexes were determined, the stereo-chemical basis for protein-DNA recognition remains an active area of investigation. One outstanding question is, how are proteins able to detect noncontacted sequences in their binding sites? The papillomavirus E2 proteins represent a particularly suitable group of proteins in which to examine the mechanisms of "indirect readout." Coordinated structural and thermodynamic studies of the E2-DNA interaction conducted over the past five years are summarized in this review. The data support a model in which the electrostatic properties of the individual E2 proteins correlate with their affinities for intrinsically flexible or rigidly prebent DNA targets.

Structure of the retinal determination protein Dachshund reveals a DNA binding motif

The Dachshund proteins are essential components of a regulatory network controlling cell fate determination. They have been implicated in eye, limb, brain, and muscle development. These proteins cannot be assigned to any recognizable structural or functional class based on amino acid sequence analysis. The 1.65 A crystal structure of the most conserved domain of human DACHSHUND is reported here. The protein forms an alpha/beta structure containing a DNA binding motif similar to that found in the winged helix/forkhead subgroup of the helix-turn-helix family. This unexpected finding alters the previously proposed molecular models for the role of Dachshund in the eye determination pathway. Furthermore, it provides a rational framework for future mechanistic analyses of the Dachshund proteins in several developmental contexts.

Bending and flexibility of methylated and unmethylated EcoRI DNA

We used cyclization kinetics experiments and Monte Carlo simulations to determine a structural model for a DNA decamer containing the EcoRI restriction site. Our findings agree well with recent crystal and NMR structures of the EcoRI dodecamer, where an overall bend of seven degrees is distributed symmetrically over the molecule. Monte Carlo simulations indicate that the sequence has a higher flexibility, assumed to be isotropic, compared to that of a "generic" DNA sequence. This model was used as a starting point for the investigation of the effect of cytosine methylation on DNA bending and flexibility. While methylation did not affect bend magnitude or direction, it resulted in a reduction in bending flexibility and under-winding of the methylated nucleotides. We demonstrate that our approach can augment the understanding of DNA structure and dynamics by adding information about the global structure and flexibility of the sequence. We also show that cyclization kinetics can be used to study the properties of modified nucleotides.

Structural basis for Ni(2+) transport and assembly of the urease active site by the metallochaperone UreE from Bacillus pasteurii

Bacillus pasteurii UreE (BpUreE) is a putative chaperone assisting the insertion of Ni(2+) ions in the active site of urease. The x-ray structure of the protein has been determined for two crystal forms, at 1.7 and 1.85 A resolution, using SIRAS phases derived from a Hg(2+)-derivative. BpUreE is composed of distinct N- and C-terminal domains, connected by a short flexible linker. The structure reveals the topology of an elongated homodimer, formed by interaction of the two C-terminal domains through hydrophobic interactions. A single Zn(2+) ion bound to four conserved His-100 residues, one from each monomer, connects two dimers resulting in a tetrameric BpUreE known to be formed in concentrated solutions. The Zn(2+) ion can be replaced by Ni(2+) as shown by anomalous difference maps obtained on a crystal of BpUreE soaked in a solution containing NiCl(2). A large hydrophobic patch surrounding the metal ion site is surface-exposed in the biologically relevant dimer. The BpUreE structure represents the first for this class of proteins and suggests a possible role for UreE in the urease nickel-center assembly.

DNA bending by an adenine--thymine tract and its role in gene regulation

To gain insight into the structural basis of DNA bending by adenine-thymine tracts (A-tracts) and their role in DNA recognition by gene-regulatory proteins, we have determined the crystal structure of the high-affinity DNA target of the cancer-associated human papillomavirus E2 protein. The three independent B-DNA molecules of the crystal structure determined at 2.2-A resolution are examples of A-tract-containing helices where the global direction and magnitude of curvature are in accord with solution data, thereby providing insights, at the base pair level, into the mechanism of DNA bending by such sequence motifs. A comparative analysis of E2-DNA conformations with respect to other structural and biochemical studies demonstrates that (i) the A-tract structure of the core region, which is not contacted by the protein, is critical for the formation of the high-affinity sequence-specific protein-DNA complex, and (ii) differential binding affinity is regulated by the intrinsic structure and deformability encoded in the base sequence of the DNA target.

Orientation of a novel DNA binding site affects human papillomavirus-mediated transcription and replication

A consensus binding site for the human papillomavirus (HPV) E2 protein was determined from an unbiased set of degenerate oligonucleotides using cyclic amplification and selection of targets (CASTing). Detectable DNA-protein complexes were formed after six to nine cycles of CASTing. A population of selected binding sites was cloned, and a consensus was determined by statistical analysis of the DNA sequences of individual isolates. Starting from a pool with 20 random bases, a consensus binding site of ACAC-N(5)-GGT was derived. CASTing and electrophoretic mobility shift analyses demonstrate that human but not bovine papillomavirus E2 proteins recognize this sequence. The presence of this sequence in papillomavirus genomes suggests a role for its function. We demonstrate that this site functionally substitutes for the canonical E2 binding site (ACCG-N(4)-CGGT) in both transient-transcription and DNA replication assays. This sequence, in most instances, is interchangeable with the resident E2 binding sites in the context of the HPV type 16 long control region. Where the novel sequence does not support E2-mediated effects on gene expression or DNA replication, we demonstrate that changing the orientation of the novel sequence restores this effect.