The papillomavirus E2 proteins: structure, function, and biology

HPV as a model for the development of prophylactic and therapeutic cancer vaccines

HPV has been linked to many human malignancies and, as such, represents a major public health crisis. The understanding of HPV biology, however, has helped tremendously in developing prophylactic vaccines, which should help in decreasing mortality due to HPV infections. Understanding HPV biology has allowed researchers to use the virus as a model for the development of not only prophylactic vaccines, but also therapeutic ones. The advantages of HPV as a model stem from the limited number of proteins encoded by the HPV genome that can be targeted by vaccines, and also from the restricted expression of certain viral proteins during different stages of infection. In this review, we discuss how HPV can be used as a model for the development of both prophylactic and therapeutic vaccines.

A strained DNA binding helix is conserved for site recognition, folding nucleation, and conformational modulation

Nucleic acid recognition is often mediated by alpha-helices or disordered regions that fold into alpha-helix on binding. A peptide bearing the DNA recognition helix of HPV16 E2 displays type II polyproline (PII) structure as judged by pH, temperature, and solvent effects on the CD spectra. NMR experiments indicate that the canonical alpha-helix is stabilized at the N-terminus, while the PII forms at the C-terminus half of the peptide. Re-examination of the dihedral angles of the DNA binding helix in the crystal structure and analysis of the NMR chemical shift indexes confirm that the N-terminus half is a canonical alpha-helix, while the C-terminal half adopts a 3(10) helix structure. These regions precisely match two locally driven folding nucleii, which partake in the native hydrophobic core and modulate a conformational switch in the DNA binding helix. The peptide shows only weak and unspecific residual DNA binding, 10(4)-fold lower affinity, and 500-fold lower discrimination capacity compared with the domain. Thus, the precise side chain conformation required for modulated and tight physiological binding by HPV E2 is largely determined by the noncanonical strained alpha-helix conformation, "presented" by this unique architecture. (c) 2009 Wiley Periodicals, Inc. Biopolymers 91: 432-443, 2009.

Human papillomavirus type 16 E2 and E6 are RNA-binding proteins and inhibit in vitro splicing of pre-mRNAs with suboptimal splice sites

Human papillomavirus type 16 (HPV16) genome expresses six regulatory proteins (E1, E2, E4, E5, E6, and E7) which regulate viral DNA replication, gene expression, and cell function. We expressed HPV16 E2, E4, E6, and E7 from bacteria as GST fusion proteins and examined their possible functions in RNA splicing. Both HPV16 E2, a viral transactivator protein, and E6, a viral oncoprotein, inhibited splicing of pre-mRNAs containing an intron with suboptimal splice sites, whereas HPV5 E2 did not. The N-terminal half and the hinge region of HPV16 E2 as well as the N-terminal and central portions of HPV16 E6 are responsible for the suppression. HPV16 E2 interacts with pre-mRNAs through its C-terminal DNA-binding domain. HPV16 E6 binds pre-mRNAs via nuclear localization signal (NLS3) in its C-terminal half. Low-risk HPV6 E6, a cytoplasmic protein, does not bind RNA. Notably, both HPV16 E2 and E6 selectively bind to the intron region of pre-mRNAs and interact with a subset of cellular SR proteins. Together, these findings suggest that HPV16 E2 and E6 are RNA binding proteins and might play roles in posttranscriptional regulation during virus infection.

Human papillomavirus type 16 E2 protein transcriptionally activates the promoter of a key cellular splicing factor, SF2/ASF

Human papillomavirus (HPV) gene expression is regulated in concert with the epithelial differentiation program. In particular, expression of the virus capsid proteins L1 and L2 is tightly restricted to differentiated epithelial cells. For HPV16, the capsid proteins are encoded by 13 structurally different mRNAs that are produced by extensive alternative splicing. Previously, we demonstrated that upon epithelial differentiation, HPV16 infection upregulates hnRNP A1 and SF2/ASF, both key factors in alternative splicing regulation. Here we cloned a 1-kb region upstream of and including the transcriptional start site of the SF2ASF gene and used it in in vivo transcription assays to demonstrate that the HPV16 E2 transcription factor transactivates the SF2/ASF promoter. The transactivation domain but not the DNA binding domain of the protein is necessary for this. Active E2 association with the promoter was demonstrated using chromatin immunoprecipitation assays. Electrophoretic mobility shift assays indicated that E2 interacted with a region 482 to 684 bp upstream of the transcription initiation site in vitro. This is the first time that HPV16 E2 has been shown to regulate cellular gene expression and the first report of viral regulation of expression of an RNA processing factor. Such E2-mediated control during differentiation of infected epithelial cells may facilitate late capsid protein expression and completion of the virus life cycle.

Evolutionary and biophysical relationships among the papillomavirus E2 proteins

Infection by human papillomavirus (HPV) may result in clinical conditions ranging from benign warts to invasive cancer. The HPV E2 protein represses oncoprotein transcription and is required for viral replication. HPV E2 binds to palindromic DNA sequences of highly conserved four base pair sequences flanking an identical length variable 'spacer'. E2 proteins directly contact the conserved but not the spacer DNA. Variation in naturally occurring spacer sequences results in differential protein affinity that is dependent on their sensitivity to the spacer DNA's unique conformational and/or dynamic properties. This article explores the biophysical character of this core viral protein with the goal of identifying characteristics that associated with risk of virally caused malignancy. The amino acid sequence, 3d structure and electrostatic features of the E2 protein DNA binding domain are highly conserved; specific interactions with DNA binding sites have also been conserved. In contrast, the E2 protein's transactivation domain does not have extensive surfaces of highly conserved residues. Rather, regions of high conservation are localized to small surface patches. Implications to cancer biology are discussed.

MD simulations of papillomavirus DNA-E2 protein complexes hints at a protein structural code for DNA deformation

The structural dynamics of the DNA binding domains of the human papillomavirus strain 16 and the bovine papillomavirus strain 1, complexed with their DNA targets, has been investigated by modeling, molecular dynamics simulations, and nuclear magnetic resonance analysis. The simulations underline different dynamical features of the protein scaffolds and a different mechanical interaction of the two proteins with DNA. The two protein structures, although very similar, show differences in the relative mobility of secondary structure elements. Protein structural analyses, principal component analysis, and geometrical and energetic DNA analyses indicate that the two transcription factors utilize a different strategy in DNA recognition and deformation. Results show that the protein indirect DNA readout is not only addressable to the DNA molecule flexibility but it is finely tuned by the mechanical and dynamical properties of the protein scaffold involved in the interaction.

Transition state for protein-DNA recognition

We describe the formation of protein-DNA contacts in the two-state route for DNA sequence recognition by a transcriptional regulator. Surprisingly, direct sequence readout establishes in the transition state and constitutes the bottleneck of complex formation. Although a few nonspecific ionic interactions are formed at this early stage, they mainly play a stabilizing role in the final consolidated complex. The interface is fairly plastic in the transition state, likely because of a high level of hydration. The overall picture of this two-state route largely agrees with a smooth energy landscape for binding that speeds up DNA recognition. This "direct" two-state route differs from the parallel multistep pathway described for this system, which involves nonspecific contacts and at least two intermediate species that must involve substantial conformational rearrangement in either or both macromolecules.

Dissecting direct and indirect readout of cAMP receptor protein DNA binding using an inosine and 2,6-diaminopurine in vitro selection system

The DNA interaction of the Escherichia coli cyclic AMP receptor protein (CRP) represents a typical example of a dual recognition mechanism exhibiting both direct and indirect readout. We have dissected the direct and indirect components of DNA recognition by CRP employing in vitro selection of a random library of DNA-binding sites containing inosine (I) and 2,6-diaminopurine (D) instead of guanine and adenine, respectively. Accordingly, the DNA helix minor groove is structurally altered due to the ‘transfer’ of the 2-amino group of guanine (now I) to adenine (now D), whereas the major groove is functionally intact. The majority of the selected sites contain the natural consensus sequence TGTGAN₆TCACA (i.e. TITIDN₆TCDCD). Thus, direct readout of the consensus sequence is independent of minor groove conformation. Consequently, the indirect readout known to occur in the TG/CA base pair step (primary kink site) in the consensus sequence is not affected by I–D substitutions. In contrast, the flanking regions are selected as I/C rich sequences (mostly I-tracts) instead of A/T rich sequences which are known to strongly increase CRP binding, thereby demonstrating almost exclusive indirect readout of helix structure/flexibility in this region through (anisotropic) flexibility of I-tracts.

Antibody recognition of a flexible epitope at the DNA binding site of the human papillomavirus transcriptional regulator E2

We have obtained a monoclonal antibody (ED15) against the C-terminal DNA-binding domain of the high-risk human papillomavirus strain-16 E2 protein that strongly interferes with its DNA-binding activity. We here characterize the recognition mechanism of this antibody and find that the ED15-E2 interaction has a strong electrostatic component, which correlates with the high proportion of acidic residues found in the antibody combining site. Further circular dichroism experiments in the presence of phosphate show that, in addition to electrostatic screening of key potential interactions, ionic strength affects the conformation of the epitope. In addition, the interaction is strongly modulated by pH, which correlates with the local flexibility of the epitope rather than the presence of pH sensitive residues at the interface. Noticeably, this finding is well correlated with the strong entropic component of the interaction. Site directed mutagenesis indicates that the ED15 epitope involves at least part of the DNA-binding helix of E2, explaining the mAb inhibitory activity. At physiological salt concentrations, the equilibrium dissociation constant of the E2-ED15 interaction is 10(-7) M and the association rate is 10(4) M-1 s-1, at least 1 order of magnitude slower than those generally reported in the most extensively described "nonflexible" antibody-protein interactions, indicating the presence of a slow conformational rearrangement on the antigen as the rate-limiting step. The crucial role of antigen flexibility in antibody-protein recognition is discussed.

Regulation of telomerase and telomeres: human tumor viruses take control

Human tumor viruses are responsible for one-fifth of all cancers worldwide. These viruses have evolved multiple strategies to evade immune defenses and to persist in the host by establishing a latent infection. Proliferation is necessary for pretumor cells to accumulate genetic alterations and to acquire a transformed phenotype. However, each cell division is associated with a progressive shortening of the telomeres, which can suppress tumor development by initiating senescence and irreversible cell cycle arrest. Therefore, the ability of virus-infected cells to circumvent the senescence program is essential for the long-term survival and proliferation of infected cells and the likelihood of transformation. We review the multiple strategies used by human DNA and RNA tumor viruses to subvert telomerase functions during cellular transformation and carcinogenesis. Epstein-Barr virus, Kaposi sarcoma-associated herpesvirus, human papillomavirus, hepatitis B virus, hepatitis C virus, and human T-cell leukemia virus-1 each can increase transcription of the telomerase reverse transcriptase. Several viruses appear to mediate cis-activation or enhance epigenetic activation of telomerase transcription. Epstein-Barr virus and human papillomavirus have each developed posttranscriptional mechanisms to regulate the telomerase protein. Finally, some tumor virus proteins can also negatively regulate telomerase transcription or activity. It is likely that, as future studies further expose the strategies used by viruses to deregulate telomerase activity and control of telomere length, novel mechanisms will emerge and underscore the importance of increased telomerase activity in sustaining virus-infected cells and its potential in therapeutic targeting.

Comprehensive comparison of the interaction of the E2 master regulator with its cognate target DNA sites in 73 human papillomavirus types by sequence statistics

Mucosal human papillomaviruses (HPVs) are etiological agents of oral, anal and genital cancer. Properties of high- and low-risk HPV types cannot be reduced to discrete molecular traits. The E2 protein regulates viral replication and transcription through a finely tuned interaction with four sites at the upstream regulatory region of the genome. A computational study of the E2–DNA interaction in all 73 types within the alpha papillomavirus genus, including all known mucosal types, indicates that E2 proteins have similar DNA discrimination properties. Differences in E2–DNA interaction among HPV types lie mostly in the target DNA sequence, as opposed to the amino acid sequence of the conserved DNA-binding alpha helix of E2. Sequence logos of natural and in vitro selected sites show an asymmetric pattern of conservation arising from indirect readout, and reveal evolutionary pressure for a putative methylation site. Based on DNA sequences only, we could predict differences in binding energies with a standard deviation of 0.64 kcal/mol. These energies cluster into six discrete affinity hierarchies and uncovered a fifth E2-binding site in the genome of six HPV types. Finally, certain distances between sites, affinity hierarchies and their eventual changes upon methylation, are statistically associated with high-risk types.

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.

Increased stability and DNA site discrimination of "single chain" variants of the dimeric beta-barrel DNA binding domain of the human papillomavirus E2 transcriptional regulator

Human papillomavirus infects millions of people worldwide and is a causal agent of cervical cancer in women. The HPV E2 protein controls the expression of all viral genes through binding of its dimeric C-terminal domain (E2C) to its target DNA site. We engineered monomeric versions of the HPV16 E2C, in order to probe the link of the dimeric beta-barrel fold to stability, dimerization, and DNA binding. Two single-chain variants, with 6 and 12 residue linkers (scE2C-6 and scE2C-12), were purified and characterized. Spectroscopy and crystallography show that the native structure is unperturbed in scE2C-12. The single chain variants are stabilized with respect to E2C, with effective concentrations of 0.6 to 6 mM. The early folding events of the E2C dimer and scE2C-12 are very similar and include formation of a compact species in the submillisecond time scale and a non-native monomeric intermediate with a half-life of 25 ms. However, monomerization changes the unfolding mechanism of the linked species from two-state to three-state, with a high-energy intermediate. Binding to the specific target site is up to 5-fold tighter in the single chain variants. Nonspecific DNA binding is up to 7-fold weaker in the single chain variants, leading to an overall 10-fold increased site discrimination capacity, the largest described so far for linked DNA binding domains. Titration calorimetric binding analysis, however, shows almost identical behavior for dimer and single-chain species, suggesting very subtle changes behind the increased specificity. Global analysis of the mechanisms probed suggests that the dynamics of the E2C domain, rather than the structure, are responsible for the differential properties. Thus, the plastic and dimeric nature of the domain did not evolve for a maximum affinity, specificity, and stability of the quaternary structure, likely because of regulatory reasons and for roles other than DNA binding played by partly folded dimeric or monomeric conformers.

Sequence-dependent DNA deformability studied using molecular dynamics simulations

Proteins recognize specific DNA sequences not only through direct contact between amino acids and bases, but also indirectly based on the sequence-dependent conformation and deformability of the DNA (indirect readout). We used molecular dynamics simulations to analyze the sequence-dependent DNA conformations of all 136 possible tetrameric sequences sandwiched between CGCG sequences. The deformability of dimeric steps obtained by the simulations is consistent with that by the crystal structures. The simulation results further showed that the conformation and deformability of the tetramers can highly depend on the flanking base pairs. The conformations of xATx tetramers show the most rigidity and are not affected by the flanking base pairs and the xYRx show by contrast the greatest flexibility and change their conformations depending on the base pairs at both ends, suggesting tetramers with the same central dimer can show different deformabilities. These results suggest that analysis of dimeric steps alone may overlook some conformational features of DNA and provide insight into the mechanism of indirect readout during protein-DNA recognition. Moreover, the sequence dependence of DNA conformation and deformability may be used to estimate the contribution of indirect readout to the specificity of protein-DNA recognition as well as nucleosome positioning and large-scale behavior of nucleic acids.

H4 histone tail mediated DNA-DNA interaction and effects on DNA structure, flexibility, and counterion binding. A molecular dynamics study

All-atom molecular dynamics (MD) simulations were performed during 30-45 ns for a system of three identical DNA 22-mers, 14 short fragments of the charged H4 histone tail peptide fragment (amino acids 5-12, KGGKGLGK) with K(+) counterions, and explicit water. The simulation setup mimics the crowded conditions of DNA in eukaryotic chromatin. To assess the influence of tail fragments on DNA structure and dynamics, a "control" 20 ns MD simulation was carried for a system with the same DNA and water content but in the absence of oligopeptides. Results of DNA interaction with the histone tail fragments, K(+), and water is presented. DNA structure and dynamics and its interplay with the histone tail fragments binding are described. The charged side chains of the lysines play a major role in mediating DNA-DNA attraction by forming bridges and coordinating to phosphate groups and electronegative sites in the minor groove. Binding of all species to DNA is dynamic. Some of the tail fragments while being flexible and mobile in each of its functional groups remain associated near certain locations of the DNA oligomer. The present work allows capturing typical features of the histone tail-counterion-DNA structure, interaction, and dynamics.

From rigid base pairs to semiflexible polymers: coarse-graining DNA

The elasticity of double-helical DNA on a nm length scale is captured in detail by the rigid base-pair model, whose conformation variables are the relative positions and orientations of adjacent base pairs. Corresponding sequence-dependent elastic potentials have been obtained from all-atom MD simulation and from high-resolution structural data. On the scale of 100 nm, DNA is successfully described by a continuous wormlike chain model with homogeneous elastic properties, characterized by a set of four elastic constants which have been measured in single-molecule experiments. We present here a theory that links these experiments on different scales, by systematically coarse-graining the rigid base-pair model to an effective wormlike chain description. The average helical geometry of the molecule is accounted for exactly, and repetitive as well as random sequences are considered. Structural disorder is shown to produce a small, additive and short-range correction to thermal conformation fluctuations as well as to entropic elasticity. We also discuss the limits of applicability of the homogeneous wormlike chain on short scales, quantifying the anisotropy of bending stiffness, the non-Gaussian bend angle distribution and the variability of stiffness, all of which are noticeable below a helical turn. The coarse-grained elastic parameters show remarkable overall agreement with experimental wormlike chain stiffness. For the best-matching potential, bending persistence lengths of dinucleotide repeats span a range of 37-53 nm, with a random DNA value of 43 nm. While twist stiffness is somewhat underestimated and stretch stiffness is overestimated, the counterintuitive negative sign and the magnitude of the twist-stretch coupling agree with recent experimental findings.

A quasi-spontaneous amyloid route in a DNA binding gene regulatory domain: The papillomavirus HPV16 E2 protein

The DNA binding domain of papillomavirus E2 proteins is at the center of the regulation of gene transcription and replication of the virus. Its unique fold consists of a beta-barrel domain that combines an eight-stranded dimeric beta-barrel core interface with two symmetrical DNA binding alpha-helices and other two helices, packed against the central barrel. Treatment with low amounts of trifluoroethanol readily leads to a mostly beta-sheet oligomeric species, with a loss of near-UV circular dichroism signal and increase in its ANS binding capacity, indicating that buried hydrophobic surfaces become accessible to the solvent. This species subsequently undergoes a slow transition into amyloid aggregates as determined by light scattering and Congo red and thioflavin T binding. Electron microscopy shows short amyloid fibers with a curly aspect as the end product. The amyloid route is completely prevented by addition of stoichiometrical amounts of specific DNA, strongly suggesting that unfolding of the DNA binding alpha-helix is required for the formation of the intermediate. The slow nature of this expanded beta-oligomeric species and the availability of several different conformational probes make it an excellent model for investigating amyloid mechanisms. The mild perturbation required for entering an amyloid route is indicative of a preexisting equilibrium. Oligomerization processes are required for the assembly of transcription initiation and DNA replication machineries, where proteins from different viruses must come together with host cell proteins. The E2 protein is a virus-encoded multifunctional master regulator that may exert one of its multiple functions through its ability to oligomerize.

Recent Advances in the Search for Antiviral Agents against Human Papillomaviruses

Infection by human papillomavirus (HPV) is extremely common and associated with the development of benign warts or malignant lesions of the skin and mucosa. Infection by a high-risk (oncogenic) anogenital HPV type, most often through sexual contacts, is the starting point of virtually all cases of cervical cancers and the majority of anal cancers. The same viral types are also increasingly being linked with a subset of head-and-neck and non-melanoma skin cancers. Although prophylactic vaccines are now available to protect against the four types most commonly found in cervical and anal cancers (HPV16 and HPV18) and anogenital warts (HPV6 and HPV11), these neither protect against all genital HPVs nor are of therapeutic utility for already infected patients. Thus, the need for antiviral agents to treat HPV-associated diseases remains great, but none currently exist. This article reviews the recent progress made towards the development of antiviral agents to treat HPV infections, from target identification and validation to the discovery of lead compounds with therapeutic potential. Emphasis has been placed on novel low-molecular-weight compounds that antagonize HPV proteins or, alternatively, inhibit cellular proteins which have been usurped by papillomaviruses and are mediating their pathogenic effects.

Molecular dynamics of the DNA-binding domain of the papillomavirus E2 transcriptional regulator uncover differential properties for DNA target accommodation

Papillomaviruses are small DNA tumor viruses that infect mammalian hosts, with consequences from benign to cancerous lesions. The Early protein 2 is the master regulator for the virus life cycle, participating in gene transcription, DNA replication, and viral episome migration. All of these functions rely on primary target recognition by its dimeric DNA-binding domain. In this work, we performed molecular dynamics simulations in order to gain insights into the structural dynamics of the DNA-binding domains of two prototypic strains, human papillomavirus strain 16 and the bovine papillomavirus strain 1. The simulations underline different dynamic features in the two proteins. The human papillomavirus strain 16 domain displays a higher flexibility of the beta2-beta3 connecting loop in comparison with the bovine papillomavirus strain 1 domain, with a consequent effect on the DNA-binding helices, and thus on the modulation of DNA recognition. A compact beta-barrel is found in human papillomavirus strain 16, whereas the bovine papillomavirus strain 1 protein is characterized by a loose beta-barrel with a large number of cavities filled by water, which provides great flexibility. The rigidity of the human papillomavirus strain 16 beta-barrel prevents protein deformation, and, as a consequence, deformable spacers are the preferred targets in complex formation. In contrast, in bovine papillomavirus strain 1, a more deformable beta-barrel confers greater adaptability to the protein, allowing the binding of less flexible DNA regions. The flexibility data are confirmed by the experimental NMR S2 values, which are reproduced well by calculation. This feature may provide the protein with an ability to discriminate between spacer sequences. Clearly, the deformability required for the formation of the Early protein 2 C-terminal DNA-binding domain-DNA complexes of various types is based not only on the rigidity of the base sequences in the DNA spacers, but also on the intrinsic deformability properties of each domain.

The DNA binding domain of a papillomavirus E2 protein programs a chimeric nuclease to cleave integrated human papillomavirus DNA in HeLa cervical carcinoma cells

Viral DNA binding proteins that direct nucleases or other protein domains to viral DNA in lytically or latently infected cells may provide a novel approach to modulate viral gene expression or replication. Cervical carcinogenesis is initiated by high-risk human papillomavirus (HPV) infection, and viral DNA persists in the cancer cells. To test whether a DNA binding domain of a papillomavirus protein can direct a nuclease domain to cleave HPV DNA in cervical cancer cells, we fused the DNA binding domain of the bovine papillomavirus type 1 (BPV1) E2 protein to the catalytic domain of the FokI restriction endonuclease, generating a BPV1 E2-FokI chimeric nuclease (BEF). BEF introduced DNA double-strand breaks on both sides of an E2 binding site in vitro, whereas DNA binding or catalytic mutants of BEF did not. After expression of BEF in HeLa cervical carcinoma cells, we detected cleavage at E2 binding sites in the integrated HPV18 DNA in these cells and also at an E2 binding site in cellular DNA. BEF-expressing cells underwent senescence, which required the DNA binding activity of BEF, but not its nuclease activity. These results demonstrate that DNA binding domains of viral proteins can target effector molecules to cognate binding sites in virally infected cells.

Shaping the Borrelia burgdorferi genome: crystal structure and binding properties of the DNA-bending protein Hbb

The genome of the Lyme disease-causing spirochete Borrelia burgdorferi encodes only a single polypeptide from the integration host factor (IHF)/HU or 'DNABII' family of nucleoid-associated proteins - Hbb. DNABII proteins induce large bends in DNA and serve as architectural factors in a variety of prokaryotic cellular processes. We have solved the crystal structure of an Hbb-DNA complex in which the DNA is bent by over 180 degrees . We find that like IHF, Hbb relies exclusively on indirect readout to recognize its cognate site. Additional binding studies show that the sequence preferences of Hbb are related to, yet distinct from those of IHF. Defining these binding characteristics may help to uncover additional roles for Hbb in Borrelia DNA metabolism as well as further our understanding of the mechanism of indirect readout.

Stepwise binding and bending of DNA by Escherichia coli integration host factor

Integration host factor (IHF) is a prokaryotic protein required for the integration of lambda phage DNA into its host genome. An x-ray crystal structure of the complex shows that IHF binds to the minor groove of DNA and bends the double helix by 160 degrees [Rice PA, Yang S, Mizuuchi K, Nash HA (1996) Cell 87:1295-1306]. We sought to dissect the complex formation process into its component binding and bending reaction steps, using stopped-flow fluorimetry to observe changes in resonance energy transfer between DNA-bound dyes, which in turn reflect distance changes upon bending. Different DNA substrates that are likely to increase or decrease the DNA bending rate were studied, including one with a nick in a critical kink position, and a substrate with longer DNA ends to increase hydrodynamic friction during bending. Kinetic experiments were carried out under pseudofirst-order conditions, in which the protein concentration is in substantial excess over DNA. At lower concentrations, the reaction rate rises linearly with protein concentration, implying rate limitation by the bimolecular reaction step. At high concentrations the rate reaches a plateau value, which strongly depends on temperature and the nature of the DNA substrate. We ascribe this reaction limit to the DNA bending rate and propose that complex formation is sequential at high concentration: IHF binds rapidly to DNA, followed by slower DNA bending. Our observations on the bending step kinetics are in agreement with results using the temperature-jump kinetic method.

Mechanism of DNA Recognition at a Viral Replication Origin

Recognition of the DNA origin by the Epstein-Barr nuclear antigen 1 (EBNA1) protein is the primary event in latentphase genome replication of the Epstein-Barr virus, a model for replication initiation in eukaryotes. We carried out an extensive thermodynamic and kinetic characterization of the binding mechanism of the DNA binding domain of EBNA1, EBNA1452-641, to a DNA fragment containing a single specific origin site. The interaction displays a binding energy of 12.7 kcal mol-1, with 11.9 kcal mol-1 coming from the enthalpic change with a minimal entropic contribution. Formation of the EBNA1452-641.DNA complex is accompanied by a heat capacity change of -1.22 kcal mol-1 K-1, a very large value considering the surface area buried, which we assign to an unusually apolar protein-DNA interface. Kinetic dissociation experiments, including fluorescence anisotropy and a continuous native electrophoretic mobility shift assay, confirmed that two EBNA1.DNA complex conformers are in slow equilibrium; one dissociates slowly (t1/2 approximately 41 min) through an undissociated intermediate species and the other corresponds to a fast twostep dissociation route (t1/2 approximately 0.8 min). In line with this, at least two parallel association events from two populations of protein conformers are observed, with on-rates of 0.25-1.6x10(8) m-1 s-1, which occur differentially either in excess protein or DNA molecules. Both parallel complexes undergo subsequent firstorder rearrangements of approximately 2.0 s-1 to yield two consolidated complexes. These parallel association and dissociation routes likely allow additional flexible regulatory events for site recognition depending on site availability according to nucleus environmental conditions, which may lock a final recognition event, dissociate and re-bind, or slide along the DNA.

Brd4 links chromatin targeting to HPV transcriptional silencing

The E2 protein encoded by human papillomaviruses (HPVs) inhibits expression of the viral E6 oncoprotein, which, in turn, regulates p53 target gene transcription. To identify cellular proteins involved in E2-mediated transcriptional repression, we isolated an E2 complex from human cells conditionally expressing HPV-11 E2. Surprisingly, the double bromodomain-containing protein Brd4, which is implicated in cell cycle control and viral genome segregation, was found associated with E2 and conferred on E2 the ability to inhibit AP-1-dependent HPV chromatin transcription in an E2-binding site-specific manner as illustrated by in vitro reconstituted chromatin transcription experiments. Knockdown of Brd4 in human cells alleviates E2-mediated repression of HPV transcription. The E2-interacting domain at the extreme C terminus and the chromatin targeting activity of a bromodomain-containing region are both essential for the corepressor activity of Brd4. Interestingly, E2-Brd4 blocks the recruitment of TFIID and RNA polymerase II to the HPV E6 promoter region without inhibiting acetylation of nucleosomal histones H3 and H4, indicating an acetylation-dependent role of Brd4 in the recruitment of E2 for transcriptional silencing of HPV gene activity. Our finding that Brd4 is a component of the virus-assembled transcriptional silencing complex uncovers a novel function of Brd4 as a cellular cofactor modulating viral gene expression.

The recognition of local DNA conformation by the human papillomavirus type 6 E2 protein

The E2 proteins are transcription/replication factors from papillomaviruses. Human papillomaviruses (HPVs) can be broadly divided in two groups; low-risk HPV subtypes cause benign warts while high-risk HPVs give rise to cervical cancer. Although a range of crystal structures of E2 DNA-binding domains (DBD) from both high- and low-risk HPV subtypes have been reported previously, structures of E2 DBD:DNA complexes have only been available for high-risk HPV18 and bovine papillomavirus (BPV1). In the present study we report the unliganded and DNA complex structures of the E2 DBD from the low-risk HPV6. As in the previous E2-DNA structures, complex formation results in considerable bending of the DNA, which is facilitated by sequences with A:T-rich spacers that adopt a pre-bent conformation. The low-risk HPV6 E2-DNA complex differs from the earlier structures in that minimal deformation of the protein accompanies complex formation. Stopped-flow kinetic studies confirm that both high- and low-risk E2 proteins adapt their structures on binding to DNA, although this is achieved more readily for HPV6 E2. It therefore appears that the higher selectivity of the HPV6 E2 protein may arise from its limited molecular adaptability, a property that might distinguish the behaviour of E2 proteins from high- and low-risk HPV subtypes.

Structural and thermodynamic basis for the enhanced transcriptional control by the human papillomavirus strain-16 E2 protein

Strain 16 of the human papillomavirus is responsible for the largest number of cases of cervical cancers linked to this virus, and the E2 protein is the transcriptional regulator of all viral genes. We present the first structure for the DNA binding domain of HPV16 E2 bound to DNA, and in particular, to a natural cognate sequence. The NMR structure of the protein backbone reveals that the overall conformation remains virtually unchanged, and chemical shift analysis of the protein bound to a shorter DNA duplex uncovered a contact out of the minimal E2 DNA binding site, made by lysine 349. This contact was confirmed by titration calorimetry and mutagenesis, with a contribution of 1.0 kcal mol(-)(1) to binding energy. HPV16 E2 has the highest DNA binding affinity and exerts a strict transcriptional control, translated into the repression of the E6 and E7 oncogenes. These novel features provide the structural and thermodynamic basis for this tight transcriptional control, the loss of which correlates with carcinogenesis.

Dynamic localization of the human papillomavirus type 11 origin binding protein E2 through mitosis while in association with the spindle apparatus

Papillomaviral DNA replicates as extrachromosomal plasmids in squamous epithelium. Viral DNA must segregate equitably into daughter cells to persist in dividing basal/parabasal cells. We have previously reported that the viral origin binding protein E2 of human papillomavirus types 11 (HPV-11), 16, and 18 colocalized with the mitotic spindles. In this study, we show the localization of the HPV-11 E2 protein to be dynamic. It colocalized with the mitotic spindles during prophase and metaphase. At anaphase, it began to migrate to the central spindle microtubules, where it remained through telophase and cytokinesis. It was additionally observed in the midbody at cytokinesis. A peptide spanning residues 285 to 308 in the carboxyl-terminal domain of HPV-11 E2 (E2C) is necessary and sufficient to confer localization on the mitotic spindles. This region is conserved in HPV-11, -16, and -18 and bovine papillomavirus type 4 (BPV-4) E2 and is also required for the respective E2C to colocalize with the mitotic spindles. The E2 protein of bovine papillomavirus type 1 is tethered to the mitotic chromosomes via the cellular protein Brd4. However, the HPV-11 E2 protein did not associate with Brd4 during mitosis. Lastly, a chimeric BPV-1 E2C containing the spindle localization domain from HPV-11 E2C gained the ability to localize to the mitotic spindles, whereas the reciprocal chimera lost the ability. We conclude that this region of HPV E2C is critical for localization with the mitotic apparatus, enabling the HPV DNA to sustain persistent infections.

Specificity in DNA recognition by a peptide from papillomavirus E2 protein

The E2 proteins of papillomavirus specifically bind to double-stranded DNA containing the consensus sequence ACCG-N4-CGGT, where N is any nucleotide. Here, we show the binding and recognition of dissimilar DNA sequences by an 18 amino-acid peptide (alpha1E2), which corresponds to the DNA-recognition helix, alpha-helix-1. Isothermal DNA binding assays performed with the DNA consensus sequence show saturable curves with alpha1E2 peptide, and the alpha1E2 peptide is converted to an ordered conformation upon complexation. Measurements performed with non-specific DNA sequence fail to saturate, a behavior characteristic of non-specific binding. Binding of the alpha1E2 peptide to these DNA sequences display a different counter-ion dependence, indicating a dissimilar, sequence-dependent mechanism of interaction. Quantitative stoichiometric measurements revealed the specificity in alpha1E2 peptide recognition of the ACCG half-site, demonstrating capacity for discrimination of nucleic acid bases sequences without the need of a whole protein architecture.

Indirect readout of DNA sequence by papillomavirus E2 proteins depends upon net cation uptake

The papillomavirus E2 proteins bind with high affinity to palindromic DNA sequences consisting of two highly conserved four base-pair sequences flanking a variable "spacer" of identical length (ACCG NNNN CGGT). While intimate contacts are observed between the bound proteins and conserved DNA in the available co-crystal structures, no contact is seen between the proteins and the spacer DNA. The ability of human papillomavirus strain 16 (HPV-16) E2 and bovine papillomavirus strain 1 (BPV-1) E2 to discriminate among binding sites with different spacer sequences is dependent on their sensitivity to the unique conformational and/or dynamic properties of the spacer DNA in a process termed "indirect readout". Differential sequence-specific K(+) uptake in low ionic strength solutions lacking Mg(2+) is observed upon E2 protein binding to sites containing the AATT, TTAA or ACGT spacer sequences. In contrast, the cation displacement typical of protein-DNA complex formation is observed at high K(+) concentrations or in the presence of Mg(2+). These results are interpreted to reflect the sequence-specific stabilization of bent DNA conformations by cations localized within the narrowed minor grooves of the protein-bound DNA and the intrinsic structure and flexibility of the DNA target. Mg(2+) differentially affects the binding of the HPV-16 E2 DNA binding domain (HPV16-E2/D) and the BPV-1 E2 DNA binding domain (BPV1-E2/D) to sites bearing different spacer sequences. This study suggests that monovalent and divalent cations contribute to the discrimination of DNA structure and flexibility that could in turn contribute to the specificity with which HPV16-E2/D and BPV1-E2/D mediate DNA replication and gene transcription.

Bending in the right direction

The success of a new Monte Carlo algorithm in describing sequence-specific DNA bending is reported. The approach offers the possibility of improved treatments of "indirect readout" effects in the prediction of transcription factor binding sites.

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.

Combined prophylactic and therapeutic cancer vaccine: Enhancing CTL responses to HPV16 E2 using a chimeric VLP in HLA-A2 mice

We identified the strategies to induce a CTL response to human papillomavirus (HPV) 16 E2 in HLA-A2 transgenic mice (AAD). A chimeric HPV16 virus-like particle (VLP) that includes full length HPV16 E7 and E2 (VLP-E7E2) was generated. The combination of E2 and E7 has the advantage that E2 is expressed in early dysplasia and neoplasia lesions, where E7 is expressed in more advance lesions. Since T cell response to E2 is less defined, we first evaluated the strategies to enhancing CD8(+) T cell responses to HPV E7, using different combinations of immune-modulators with VLP-E7E2. Data showed that the CTL response to E7 could be significantly enhanced by coinjection of GM-CSF and anti-CD40 antibodies with chimeric VLP-E7E2 without adjuvant. However, using the same combination, a low level of CD8(+) T cell response to E2 was detected. To enhance the CD8+ T cell response to E2, we analyzed T cell epitopes from E2 sequence. A heterogenous prime-boost with chimeric VLP-E7E2 and E2 peptides was performed. The data showed that the priming with chimeric VLP-E7E2, followed by boosting with E2 peptides, gave a better CTL response than 2 immunizations with E2 peptides. The enhanced immunity is due to the increase of CD11c(+) and CD11c(+) CD40(+) double positive dendritic cells in mice that received immune-modulators, GM-CSF and anti-CD40. Furthermore, the level of anti-L1 antibodies remains similar in mice immunized with chimeric VLP with/without immune-modulators. Thus, the data suggested that the chimeric VLP-E7E2 has a therapeutic potential for the treatment of HPV-associated CINs and cancer without diminishing VLPs potential as a prophylactic vaccine by inducing anti-L1 antibodies against free virus.

Patterns of nucleotide difference in overlapping and non-overlapping reading frames of papillomavirus genomes

On the basis of a phylogenetic analysis of papillomavirus complete genomes, 30 pairs of closely related genomes were chosen for examination of the pattern of synonymous and non-synonymous nucleotide difference in overlapping and non-overlapping reading frames. The proportion of synonymous differences per synonymous site (pS) was generally reduced in overlapping regions in comparison to non-overlapping regions. This reduction in pS appears to result from purifying selection acting on the overlapping reading frame because of functional constraint on the amino acid sequence that it encodes. As a consequence of such purifying selection, a trade-off was observed between synonymous substitution in one reading frame and non-synonymous substitution in the other; this pattern was particularly evident in the overlap between the E2 and E1 reading frames, which encodes a highly conserved portion of the E2 protein. Fifty-three sequences of the hinge region of the HPV16 E2 protein (where the E2 and E4 reading frames overlap) showed an elevated level of non-synonymous nucleotide substitution in the E2 reading frame, particularly in a region believed to be involved in immune recognition. Non-synonymous substitutions in this region occurred disproportionately in such a way as to cause synonymous substitutions in the E4 reading frame. Thus, positive diversifying selection on the E2 hinge region appears to coexist with purifying selection on the overlapping region of E4.

Fine mapping of inherent flexibility variation along DNA molecules: validation by atomic force microscopy (AFM) in buffer

Curvature and flexibility are structural properties of central importance to genome function. However, due to the difficulties in finding suitable experimental conditions, methods for studying one without the interference of the other have proven to be difficult. We propose a new approach that provides a measure of inherent flexibility of DNA by taking advantage of two powerful techniques, X-ray crystallography and nuclear magnetic resonance. Both techniques are able to detect local curvature on DNA fragments but, while the first analyzes DNA in the solid state, the second works on DNA in solution. Comparison of the two data sets allowed us to calculate the relative contribution to flexibility of the three rotations and three translations, which relate successive base pair planes for the ten different dinucleotide steps. These values were then used to compute the variation of flexibility along a given nucleotide sequence. This allowed us to validate the method experimentally through comparisons with maps of local fluctuations in DNA molecule trajectory constructed from atomic force microscopy imaging in solution. We conclude that the six dinucleotide-step parameters defined here provide a powerful tool for the exploration of DNA structure and, consequently will make an important contribution to our understanding of DNA-sequence-dependent biological processes.

Nucleo-cytoplasmic shuttling of high risk human Papillomavirus E2 proteins induces apoptosis

Human Papillomavirus (HPV) E2 proteins are the major viral regulators of transcription and replication during the viral life cycle. In addition to these conserved functions, we show that E2 proteins from high risk HPV types 16 and 18, which are associated with cervical cancer, can induce apoptosis. In contrast, E2 proteins from low risk HPV types 6 and 11, which are associated with benign lesions, do not cause cell death. We show that the ability to induce apoptosis is linked to the intracellular localization of the respective E2 proteins rather than to inherent properties of the proteins. Although low risk HPV E2 proteins remain strictly nuclear, high risk HPV E2 proteins are present in both the nucleus and the cytoplasm of expressing cells due to exportin-1 receptor (CRM1)-dependent nucleo-cytoplasmic shuttling. Induction of apoptosis is caused by accumulation of E2 in the cytoplasm and involves caspase 8 activation. We speculate that disruption of the E2 gene during viral genome integration in cervical carcinoma provides a means to avoid E2-induced apoptosis and allow initiation of carcinogenesis.

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.

The Folding Mechanism of a Dimeric β-Barrel Domain

The dimeric beta-barrel domain is an unusual topology, shared only by two viral origin binding proteins, where secondary, tertiary and quaternary structure are coupled, and where the dimerization interface is composed of two four-stranded half-beta-barrels. The folding of the DNA binding domain of the E2 transcriptional regulator from human papillomavirus, strain-16, takes place through a stable and compact monomeric intermediate, with 31% the stability of the folded dimeric domain. Double jump multiple wavelength experiments allowed the reconstruction of the fluorescence spectrum of the monomeric intermediate at 100 milliseconds, indicating that tryptophan residues, otherwise buried in the folded state, are accessible to the solvent. Burial of surface area as well as differential behavior to ionic strength and pH with respect to the native ground state, plus the impossibility of having over 2500 A2 of surface area of the half-barrel exposed to the solvent, indicates that the formation of a non-native compact tertiary structure precedes the assembly of native quaternary structure. The monomeric intermediate can dimerize, albeit with a weaker affinity (approximately 1 microM), to yield a non-native dimeric intermediate, which rearranges to the native dimer through a parallel folding channel, with a unimolecular rate-limiting step. Folding pathways from either acid or urea unfolded states are identical, making the folding model robust. Unfolding takes place through a major phase accounting for apparently all the secondary structure change, with identical rate constant to that of the fluorescence unfolding experiment. In contrast to the folding direction, no unfolding intermediate was found.

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.

Protein-DNA Recognition Patterns and Predictions

Structural data on protein-DNA complexes provide clues for understanding the mechanism of protein-DNA recognition. Although the structures of a large number of protein-DNA complexes are known, the mechanisms underlying their specific binding are still only poorly understood. Analysis of these structures has shown that there is no simple one-to-one correspondence between bases and amino acids within protein-DNA complexes; nevertheless, the observed patterns of interaction carry important information on the mechanisms of protein-DNA recognition. In this review, we show how the patterns of interaction, either observed in known structures or derived from computer simulations, confer recognition specificity, and how they can be used to examine the relationship between structure and specificity and to predict target DNA sequences used by regulatory proteins.

Overlapping reading frames in closely related human papillomaviruses result in modular rates of selection within E2

A core group of four open reading frames (ORFs) is present in all known papillomaviruses (PVs): the E1 and E2 replication/transcription proteins and the L1 and L2 structural proteins. Because they are involved in processes that are essential to PV propagation, the sequences of these proteins are well-conserved. However, sequencing of novel subtypes for human papillomaviruses (HPV) 54 (AE9) and 82 (AE2/IS39), coupled to analysis of four other closely related genital HPV pairs, indicated that E2 has a higher dN/dS ratio than E1, L1 or L2. The elevated ratio is not homogeneous across the length of the ORF, but instead varies with respect to E2's three domains. The E2 hinge region is of particular interest, because its hypervariability (dN/dS>1) differs markedly from the two domains that it joins: the transcription-activation domain and the DNA-binding domain. Deciphering whether the hinge region's high rate of non-synonymous change is the result of positive Darwinian selection or relaxed constraint depends on the evolutionary behaviour of E4, an ORF that overlaps E2. The E2 hinge region is contained within E4 and non-synonymous changes in the hinge are associated with a disproportionate amount of synonymous change in E4, a case of simultaneous positive and purifying selection in overlapping reading frames. Modular rates of selection among E2 domains are a likely consequence of the presence of an embedded E4. E4 appears to be positioned in a part of the HPV genome that can tolerate non-synonymous change and purifying selection of E4 may be indicative of its functional importance.

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 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.

Efficient intracellular assembly of papillomaviral vectors

Although the papillomavirus structural proteins, L1 and L2, can spontaneously coassemble to form virus-like particles, currently available methods for production of L1/L2 particles capable of transducing reporter plasmids into mammalian cells are technically demanding and relatively low-yield. In this report, we describe a simple 293 cell transfection method for efficient intracellular production of papillomaviral-based gene transfer vectors carrying reporter plasmids. Using bovine papillomavirus type 1 (BPV1) and human papillomavirus type 16 as model papillomaviruses, we have developed a system for producing papillomaviral vector stocks with titers of several billion transducing units per milliliter. Production of these vectors requires both L1 and L2, and transduction can be prevented by papillomavirus-neutralizing antibodies. The stocks can be purified by an iodixanol (OptiPrep) gradient centrifugation procedure that is substantially more effective than standard cesium chloride gradient purification. Although earlier data had suggested a potential role for the viral early protein E2, we found that E2 protein expression did not enhance the intracellular production of BPV1 vectors. It was also possible to encapsidate reporter plasmids devoid of BPV1 DNA sequences. BPV1 vector production efficiency was significantly influenced by the size of the target plasmid being packaged. Use of 6-kb target plasmids resulted in BPV1 vector yields that were higher than those with target plasmids closer to the native 7.9-kb size of papillomavirus genomes. The results suggest that the intracellular assembly of papillomavirus structural proteins around heterologous reporter plasmids is surprisingly promiscuous and may be driven primarily by a size discrimination mechanism.