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

Triple lysine and nucleosome-binding motifs of the viral IE19 protein are required for human cytomegalovirus S-phase infections

Herpesvirus genomes are maintained as extrachromosomal plasmids within the nuclei of infected cells. Some herpesviruses persist within dividing cells, putting the viral genome at risk of being lost to the cytoplasm during mitosis because karyokinesis (nuclear division) requires nuclear envelope breakdown. Oncogenic herpesviruses (and papillomaviruses) avoid genome loss during mitosis by tethering their genomes to cellular chromosomes, thereby ensuring viral genome uptake into newly formed nuclei. These viruses use viral proteins with DNA- and chromatin-binding capabilities to physically link viral and cellular genomes together in a process called tethering. The known viral tethering proteins of human papillomavirus (E2), Epstein-Barr virus (EBNA1), and Kaposi's sarcoma-associated herpesvirus (LANA) each contain two independent domains required for genome tethering, one that binds sequence specifically to the viral genome and another that binds to cellular chromatin. This latter domain is called a chromatin tethering domain (CTD). The human cytomegalovirus UL123 gene encodes a CTD that is required for the virus to productively infect dividing fibroblast cells within the S phase of the cell cycle, presumably by tethering the viral genome to cellular chromosomes during mitosis. The CTD-containing UL123 gene product that supports S-phase infections is the IE19 protein. Here, we define two motifs in IE19 required for S-phase infections: an N-terminal triple lysine motif and a C-terminal nucleosome-binding motif within the CTD.IMPORTANCEThe IE19 protein encoded by human cytomegalovirus (HCMV) is required for S-phase infection of dividing cells, likely because it tethers the viral genome to cellular chromosomes, thereby allowing them to survive mitosis. The mechanism through which IE19 tethers viral genomes to cellular chromosomes is not understood. For human papillomavirus, Epstein-Barr virus, and Kaposi's sarcoma-associated herpesvirus, viral genome tethering is required for persistence (latency) and pathogenesis (oncogenesis). Like these viruses, HCMV also achieves latency, and it modulates the properties of glioblastoma multiforme tumors. Therefore, defining the mechanism through which IE19 tethers viral genomes to cellular chromosomes may help us understand, and ultimately combat or control, HCMV latency and oncomodulation.

A point mutation in the DNA-binding domain of HPV-2 E2 protein increases its DNA-binding capacity and reverses its transcriptional regulatory activity on the viral early promoter

Background

The human papillomavirus (HPV) E2 protein is a multifunctional DNA-binding protein. The transcriptional activity of HPV E2 is mediated by binding to its specific binding sites in the upstream regulatory region of the HPV genomes. Previously we reported a HPV-2 variant from a verrucae vulgaris patient with huge extensive clustered cutaneous, which have five point mutations in its E2 ORF, L118S, S235P, Y287H, S293R and A338V. Under the control of HPV-2 LCR, co-expression of the mutated HPV E2 induced an increased activity on the viral early promoter. In the present study, a series of mammalian expression plasmids encoding E2 proteins with one to five amino acid (aa) substitutions for these mutations were constructed and transfected into HeLa, C33A and SiHa cells.

Results

CAT expression assays indicated that the enhanced promoter activity was due to the co-expressions of the E2 constructs containing A338V mutation within the DNA-binding domain. Western blots analysis demonstrated that the transiently transfected E2 expressing plasmids, regardless of prototype or the A338V mutant, were continuously expressed in the cells. To study the effect of E2 mutations on its DNA-binding activity, a serial of recombinant E2 proteins with various lengths were expressed and purified. Electrophoresis mobility shift assays (EMSA) showed that the binding affinity of E2 protein with A338V mutation to both an artificial probe with two E2 binding sites or HPV-2 and HPV-16 promoter-proximal LCR sequences were significantly stronger than that of the HPV-2 prototype E2. Furthermore, co-expression of the construct containing A338V mutant exhibited increased activities on heterologous HPV-16 early promoter P97 than that of prototype E2.

Conclusions

These results suggest that the mutation from Ala to Val at aa 338 is critical for E2 DNA-binding and its transcriptional regulation.

New human papilloma virus E2 transcription factor mimics: a tripyrrole-peptide conjugate with tight and specific DNA-recognition

Background

Human papillomavirus (HPV) is the main causative agent of cervical cancer, particularly high risk strains such us HPV-16, -18 and -31. The viral encoded E2 protein acts as a transcriptional modulator and exerts a key role in viral DNA replication. Thus, E2 constitutes an attractive target for developing antiviral agents. E2 is a homodimeric protein that interacts with the DNA target through an α-helix of each monomer. However, a peptide corresponding to the DNA recognition helix of HPV-16 E2 binds DNA with lower affinity than its full-length DNA binding domain. Therefore, in an attempt to promote the DNA binding of the isolated peptide, we have designed a conjugate compound of the E2 α-helix peptide and a derivative of the antibiotic distamycin, which involves simultaneous minor- and major-groove interactions.

Methodology/Principal Findings

An E2 α-helix peptide-distamycin conjugate was designed and synthesized. It was characterized by NMR and CD spectroscopy, and its DNA binding properties were investigated by CD, DNA melting and gel shift experiments. The coupling of E2 peptide with distamycin does not affect its structural properties. The conjugate improves significantly the affinity of the peptide for specific DNA. In addition, stoichiometric amounts of specific DNA increase meaningfully the helical population of the peptide. The conjugate enhances the DNA binding constant 50-fold, maintaining its specificity.

Conclusions/Significance

These results demonstrate that peptide-distamycin conjugates are a promising tool to obtain compounds that bind the E2 target DNA-sequences with remarkable affinity and suggest that a bipartite major/minor groove binding scaffold can be a useful approach for therapeutic treatment of HPV infection.

Protein flexibility directs DNA recognition by the papillomavirus E2 proteins

Although DNA flexibility is known to play an important role in DNA–protein interactions, the importance of protein flexibility is less well understood. Here, we show that protein dynamics are important in DNA recognition using the well-characterized human papillomavirus (HPV) type 6 E2 protein as a model system. We have compared the DNA binding properties of the HPV 6 E2 DNA binding domain (DBD) and a mutant lacking two C-terminal leucine residues that form part of the hydrophobic core of the protein. Deletion of these residues results in increased specific and non-specific DNA binding and an overall decrease in DNA binding specificity. Using ¹⁵N NMR relaxation and hydrogen/deuterium exchange, we demonstrate that the mutation results in increased flexibility within the hydrophobic core and loop regions that orient the DNA binding helices. Stopped-flow kinetic studies indicate that increased flexibility alters DNA binding by increasing initial interactions with DNA but has little or no effect on the structural rearrangements that follow this step. Taken together these data demonstrate that subtle changes in protein dynamics have a major influence on protein–DNA interactions.

Mutational analysis of kinetic partitioning in protein folding and protein-DNA binding

Kinetic partitioning between competing routes is present in many biological processes. Here, we propose a methodology to characterize kinetic partitioning through site-directed mutagenesis and apply it to parallel routes for unfolding of the TI I27 protein and for recognition of its target DNA by the human papillomavirus E2 protein. The balance between the two competing reaction routes can be quantified by the partitioning constant K(p). K(p) is easily modulated by point mutations, opening the way for the rational design of kinetic partitioning. Conserved wild-type residues strongly favor one of the two competing reactions, suggesting that in these systems there is an evolutionary pressure to shift partitioning towards a certain route. The mutations with the largest effects on partitioning cluster together in space, defining the protein regions most relevant for the modulation of partitioning. Such regions are neither fully coincident with nor strictly segregated from the regions that are important from each competing reaction. We dissected the mutational effects on partitioning into the contributions from each competing route using a new parameter called pi-value. The results suggest how the design of kinetic partitioning may be approached in each case.

Contributions of the histidine side chain and the N-terminal alpha-amino group to the binding thermodynamics of oligopeptides to nucleic acids as a function of pH

Interactions of histidine with nucleic acid phosphates and histidine pK(a) shifts make important contributions to many protein-nucleic acid binding processes. To characterize these phenomena in simplified systems, we quantified binding of a histidine-containing model peptide HWKK ((+)NH(3)-His-Trp-Lys-Lys-NH(2)) and its lysine analogue KWKK ((+)NH(3)-Lys-Trp-Lys-Lys-NH(2)) to a single-stranded RNA model, polyuridylate (polyU), by changes in tryptophan fluorescence as a function of salt concentration and pH. For both HWKK and KWKK, equilibrium binding constants, K(obs), and magnitudes of log-log salt derivatives, SK(obs) identical with (partial differential logK(obs)/partial differential log[Na(+)]), decreased with increasing pH in the manner expected for a titration curve model in which deprotonation of the histidine and alpha-amino groups weakens binding and reduces its salt-dependence. Fully protonated HWKK and KWKK exhibit the same K(obs) and SK(obs) within uncertainty, and these SK(obs) values are consistent with limiting-law polyelectrolyte theory for +4 cationic oligopeptides binding to single-stranded nucleic acids. The pH-dependence of HWKK binding to polyU provides no evidence for pK(a) shifts nor any requirement for histidine protonation, in stark contrast to the thermodynamics of coupled protonation often seen for these cationic residues in the context of native protein structure where histidine protonation satisfies specific interactions (e.g., salt-bridge formation) within highly complementary binding interfaces. The absence of pK(a) shifts in our studies indicates that additional Coulombic interactions across the nonspecific-binding interface between RNA and protonated histidine or the alpha-amino group are not sufficient to promote proton uptake for these oligopeptides. We present our findings in the context of hydration models for specific vs nonspecific nucleic acid binding.

Experimental snapshots of a protein-DNA binding landscape

Protein recognition of DNA sites is a primary event for gene function. Its ultimate mechanistic understanding requires an integrated structural, dynamic, kinetic, and thermodynamic dissection that is currently limited considering the hundreds of structures of protein-DNA complexes available. We describe a protein-DNA-binding pathway in which an initial, diffuse, transition state ensemble with some nonnative contacts is followed by formation of extensive nonnative interactions that drive the system into a kinetic trap. Finally, nonnative contacts are slowly rearranged into native-like interactions with the DNA backbone. Dissimilar protein-DNA interfaces that populate along the DNA-binding route are explained by a temporary degeneracy of protein-DNA interactions, centered on "dual-role" residues. The nonnative species slow down the reaction allowing for extended functionality.

The basic helix-loop-helix region of human neurogenin 1 is a monomeric natively unfolded protein which forms a "fuzzy" complex upon DNA binding

Neuronal specification is regulated by the activity of transcription factors containing the basic helix-loop-helix motif (bHLH); these regulating proteins include, among others, the neurogenin (Ngn) family, related to the atonal family of genes. Neurogenin 1 (NGN1) is a 237-residue protein that contains a bHLH domain and is involved in neuronal differentiation. In this work, we synthesized the bHLH region of NGN1 (bHLHN) comprising residues 90-150 of the full-length NGN1. The domain is a monomeric natively unfolded protein with a pH-dependent premolten globule conformation, as shown by several spectroscopic techniques (namely, NMR, fluorescence, FTIR, and circular dichroism). The unfolded character of the domain also explains, first, the impossibility of its overexpression in several Escherichia coli strains and, second, its insolubility in aqueous buffers. To the best of our knowledge, this is the first extensive study of the conformational preferences of a bHLH domain under different solution conditions. Upon binding to two DNA E-boxes, the protein forms "fuzzy" complexes (that is, the complexes were not fully folded). The affinities of bHLHN for both DNA boxes were smaller than those of other bHLH domains, which might explain why the protein-DNA complexes were not fully folded.

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.

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.

Dimerization of the human papillomavirus type 16 E2 N terminus results in DNA looping within the upstream regulatory region

Papillomavirus E2 proteins play a central role in regulating viral gene expression and replication. DNA-binding activity is associated with the C-terminal domain of E2, which forms a stable dimer, while the N-terminal domain is responsible for E2's replication and transactivation functions. The crystal structure of the latter domain revealed a second dimerization interface on E2 which may be responsible for DNA loop formation in the regulatory region of the human papillomavirus (HPV) genome. We investigated the biological significance of the N-terminal dimerization by introducing single amino acid substitutions into the dimerization interface. As expected, these substitutions did not influence the C-terminal dimerization and DNA-binding functions of E2. However, the mutations led to reduced transactivation of a synthetic E2-responsive reporter gene, while HPV DNA replication was unaffected. The effect of the mutations on DNA looping was visualized by atomic force microscopy. While wild-type E2 was able to generate DNA loops, all three mutant E2 proteins were defective in this ability. Our results suggest that N-terminal dimerization plays a role in E2-mediated transactivation, probably via DNA looping, a common mechanism for remote regulation of gene transcription.

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.

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.

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.

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.

Papillomavirus proteins and their potential as drug design targets

The papillomaviruses are a family of small, double-stranded DNA viruses that infect the basal cells of cutaneous and mucosal epithelium. While a large percentage of the population is benignly infected with various strains of human papillomavirus (HPV), long-term infection by a subset of HPV strains is associated with malignant transformation. The prospects for prophylaxis against HPV infection have recently received an enormous boost with the approval by the US FDA of a vaccine targeted against the most common cancer-associated HPV strains. However, the large number of people already infected, the high cost of the vaccination regimen (particularly in poorer countries) and the HPV infections that these vaccines do not protect against underscore the need for therapeutic strategies. The elucidation of molecular details underlying fundamental processes in the viral life cycle, such as virus replication, transcription and HPV-induced carcinogenesis, is required to meet this aim. This article provides an overview of high-resolution structures of papillomavirus proteins and their functional complexes, with particular reference to mechanistic and structural features that could be exploited in the rational design of antiviral agents.

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.