The papillomavirus E2 proteins: structure, function, and biology

HPV and RNA Binding Proteins: What We Know and What Remains to Be Discovered

Papillomavirus gene regulation is largely post-transcriptional due to overlapping open reading frames and the use of alternative polyadenylation and alternative splicing to produce the full suite of viral mRNAs. These processes are controlled by a wide range of cellular RNA binding proteins (RPBs), including constitutive splicing factors and cleavage and polyadenylation machinery, but also factors that regulate these processes, for example, SR and hnRNP proteins. Like cellular RNAs, papillomavirus RNAs have been shown to bind many such proteins. The life cycle of papillomaviruses is intimately linked to differentiation of the epithelial tissues the virus infects. For example, viral late mRNAs and proteins are expressed only in the most differentiated epithelial layers to avoid recognition by the host immune response. Papillomavirus genome replication is linked to the DNA damage response and viral chromatin conformation, processes which also link to RNA processing. Challenges with respect to elucidating how RBPs regulate the viral life cycle include consideration of the orchestrated spatial aspect of viral gene expression in an infected epithelium and the epigenetic nature of the viral episomal genome. This review discusses RBPs that control viral gene expression, and how the connectivity of various nuclear processes might contribute to viral mRNA production.

Determination and Characterization of Novel Papillomavirus and Parvovirus Associated with Mass Mortality of Chinese Tongue Sole (Cynoglossus semilaevis) in China

A massive mortality event concerning farmed Chinese tongue soles occurred in Tianjin, China, and the causative agent remains unknown. Here, a novel Cynoglossus semilaevis papillomavirus (CsPaV) and parvovirus (CsPV) were simultaneously isolated and identified from diseased fish via electron microscopy, virus isolation, genome sequencing, experimental challenges, and fluorescence in situ hybridization (FISH). Electron microscopy showed large numbers of virus particles present in the tissues of diseased fish. Viruses that were isolated and propagated in flounder gill cells (FG) induced typical cytopathic effects (CPE). The cumulative mortality of fish given intraperitoneal injections reached 100% at 7 dpi. The complete genomes of CsPaV and CsPV comprised 5939 bp and 3663 bp, respectively, and the genomes shared no nucleotide sequence similarities with other viruses. Phylogenetic analysis based on the L1 and NS1 protein sequences revealed that CsPaV and CsPV were novel members of the Papillomaviridae and Parvoviridae families. The FISH results showed positive signals in the spleen tissues of infected fish, and both viruses could co-infect single cells. This study represents the first report where novel papillomavirus and parvovirus are identified in farmed marine cultured fish, and it provides a basis for further studies on the prevention and treatment of emerging viral diseases.

Molecular cloning, biophysical and in silico studies of Human papillomavirus 33 E2 DNA binding domain

Human papillomavirus 33, a high-risk HPV strain, is mainly responsible for HPV infection and cervical cancer in Asian countries. The E2 protein of HPV 33 is a DNA-binding protein that plays a crucial role in viral replication and transcription. We have cloned, overexpressed, and purified the DNA binding domain of the E2 protein. Size exclusion chromatography results suggested that the protein exists in a homodimeric state in the native form. Circular dichroism data showed that the protein has a higher content of β-sheet. The melting temperature obtained from differential scanning calorimetry is 52.59 °C, and the protein is stable at pH 8 and is in a dimeric form at basic pH. The protein is monomeric or unfolded at a very low pH. Chemical denaturation studies suggested that the protein denatured and dissociated simultaneously. The DNA binding activity of the protein was also confirmed and it showed binding affinity in the order of 106 M-1. The protein structure was modeled using homology modeling and other bioinformatic tools. The virtual screening and molecular dynamic simulation studies were performed to find compounds that can act as potent inhibitors against E2 DBD. This study expands the understanding of the conserved structural and binding properties of HPV33 E2 DBD and provides the first report on the characterization of the viral protein.Communicated by Ramaswamy H. Sarma.

Biomolecular Condensation of the Human Papillomavirus E2 Master Regulator with p53: Implications in Viral Replication

p53 exerts its tumour suppressor activity by modulating hundreds of genes and it can also repress viral replication. Such is the case of human papillomavirus (HPV) through targeting the E2 master regulator, but the biochemical mechanism is not known. We show that the C-terminal DNA binding domain of HPV16 E2 protein (E2C) triggers heterotypic condensation with p53 at a precise 2/1 E2C/p53 stoichiometry at the onset for demixing, yielding large regular spherical droplets that increase in size with E2C concentration. Interestingly, transfection experiments show that E2 co-localizes with p53 in the nucleus with a grainy pattern, and recruits p53 to chromatin-associated foci, a function independent of the DNA binding capacity of p53 as judged by a DNA binding impaired mutant. Depending on the length, DNA can either completely dissolve or reshape heterotypic droplets into irregular condensates containing p53, E2C, and DNA, and reminiscent of that observed linked to chromatin. We propose that p53 is a scaffold for condensation in line with its structural and functional features, in particular as a promiscuous hub that binds multiple cellular proteins. E2 appears as both client and modulator, likely based on its homodimeric DNA binding nature. Our results, in line with the known role of condensation in eukaryotic gene enhancement and silencing, point at biomolecular condensation of E2 with p53 as a means to modulate HPV gene function, strictly dependent on host cell replication and transcription machinery.

Protein-DNA Interactions Regulate Human Papillomavirus DNA Replication, Transcription, and Oncogenesis

Human papillomavirus (HPV) is a group of alpha papillomaviruses that cause various illnesses, including cancer. There are more than 160 types of HPV, with many being "high-risk" types that have been clinically linked to cervical and other types of cancer. "Low-risk" types of HPV cause less severe conditions, such as genital warts. Over the past few decades, numerous studies have shed light on how HPV induces carcinogenesis. The HPV genome is a circular double-stranded DNA molecule that is approximately 8 kilobases in size. Replication of this genome is strictly regulated and requires two virus-encoded proteins, E1 and E2. E1 is a DNA helicase that is necessary for replisome assembly and replication of the HPV genome. On the other hand, E2 is responsible for initiating DNA replication and regulating the transcription of HPV-encoded genes, most importantly the E6 and E7 oncogenes. This article explores the genetic characteristics of high-risk HPV types, the roles of HPV-encoded proteins in HPV DNA replication, the regulation of transcription of E6 and E7 oncogenes, and the development of oncogenesis.

Role of plasmonics in detection of deadliest viruses: a review

Viruses have threatened animal and human lives since a long time ago all over the world. Some of these tiny particles have caused disastrous pandemics that killed a large number of people with subsequent economic downturns. In addition, the quarantine situation itself encounters the challenges like the deficiency in the online educational system, psychiatric problems and poor international relations. Although viruses have a rather simple protein structure, they have structural heterogeneity with a high tendency to mutation that impedes their study. On top of the breadth of such worldwide worrying issues, there are profound scientific gaps, and several unanswered questions, like lack of vaccines or antivirals to combat these pathogens. Various detection techniques like the nucleic acid test, immunoassay, and microscopy have been developed; however, there is a tradeoff between their advantages and disadvantages like safety in sample collecting, invasiveness, sensitivity, response time, etc. One of the highly resolved techniques that can provide early-stage detection with fast experiment duration is plasmonics. This optical technique has the capability to detect viral proteins and genomes at the early stage via highly sensitive interaction between the biological target and the plasmonic chip. The efficiency of this technique could be proved using commercialized techniques like reverse transcription polymerase chain reaction (RT-PCR) and enzyme-linked immunosorbent assay (ELISA) techniques. In this study, we aim to review the role of plasmonic technique in the detection of 11 deadliest viruses besides 2 common genital viruses for the human being. This is a rapidly moving topic of research, and a review article that encompasses the current findings may be useful for guiding strategies to deal with the pandemics. By investigating the potential aspects of this technique, we hope that this study could open new avenues toward the application of point-of-care techniques for virus detection at early stage that may inhibit the progressively hygienic threats.

Whole Genomic Analysis and Comparison of Two Canine Papillomavirus Type 9 Strains in Malignant and Benign Skin Lesions

Papillomaviruses (PVs) usually cause benign proliferative lesions in the stratified epithelium of various animal species. However, some high-risk types of PVs have been proven to lead to malignant transformations. In dogs, several canine papillomaviruses (CPVs) have been identified in malignant lesions and are suggested as one of the risk factors for the development of squamous cell carcinomas (SCCs). In the present study, the full genomes of two CPV9 strains from recurrent SCCs of Dog 1 and skin viral papilloma (viral plaque) of Dog 2 were sequenced. Alignment of the two CPV9 sequences with the genome of the reference CPV9 strain (accession no. JF800656.1) derived from a solitary pigmented plaque was performed. Compared with the reference strain, a 27 bp in-frame insertion in the E1 gene was identified in both CPV9 strains in this study. In comparison with the CPV9 strains derived from benign lesions, the CPV9 from the SCCs of Dog 1 exhibited a 328 bp deletion at the 3′ end of the E2 and spacer sequence, which encoded a truncated deduced E2 protein and a chimeric E8^E2 protein. However, there was no difference in the mRNA expression levels of viral oncoproteins of E6 and E7 between the two CPV9 cases, suggesting that the oncogenesis of CPV9 for malignant transformation might be different from that of human papillomaviruses. The roles of E2 and E8^E2 deleted CPV9 in the oncogenesis of benign and malignant lesions should be further investigated.

Epigenetic Regulation of the Human Papillomavirus Life Cycle

Persistent infection with certain types of human papillomaviruses (HPVs), termed high risk, presents a public health burden due to their association with multiple human cancers, including cervical cancer and an increasing number of head and neck cancers. Despite the development of prophylactic vaccines, the incidence of HPV-associated cancers remains high. In addition, no vaccine has yet been licensed for therapeutic use against pre-existing HPV infections and HPV-associated diseases. Although persistent HPV infection is the major risk factor for cancer development, additional genetic and epigenetic alterations are required for progression to the malignant phenotype. Unlike genetic mutations, the reversibility of epigenetic modifications makes epigenetic regulators ideal therapeutic targets for cancer therapy. This review article will highlight the recent advances in the understanding of epigenetic modifications associated with HPV infections, with a particular focus on the role of these epigenetic changes during different stages of the HPV life cycle that are closely associated with activation of DNA damage response pathways.

The Role of RNA Splicing Factors in Cancer: Regulation of Viral and Human Gene Expression in Human Papillomavirus-Related Cervical Cancer

The spliceosomal complex components, together with the heterogeneous nuclear ribonucleoproteins (hnRNPs) and serine/arginine-rich (SR) proteins, regulate the process of constitutive and alternative splicing, the latter leading to the production of mRNA isoforms coding multiple proteins from a single pre-mRNA molecule. The expression of splicing factors is frequently deregulated in different cancer types causing the generation of oncogenic proteins involved in cancer hallmarks. Cervical cancer is caused by persistent infection with oncogenic human papillomaviruses (HPVs) and constitutive expression of viral oncogenes. The aberrant activity of hnRNPs and SR proteins in cervical neoplasia has been shown to trigger the production of oncoproteins through the processing of pre-mRNA transcripts either derived from human genes or HPV genomes. Indeed, hnRNP and SR splicing factors have been shown to regulate the production of viral oncoprotein isoforms necessary for the completion of viral life cycle and for cell transformation. Target-therapy strategies against hnRNPs and SR proteins, causing simultaneous reduction of oncogenic factors and inhibition of HPV replication, are under development. In this review, we describe the current knowledge of the functional link between RNA splicing factors and deregulated cellular as well as viral RNA maturation in cervical cancer and the opportunity of new therapeutic strategies.

The HPV E2 Transcriptional Transactivation Protein Stimulates Cellular DNA Polymerase Epsilon

The papillomavirus (PV) protein E2 is one of only two proteins required for viral DNA replication. E2 is the viral transcriptional regulator/activation protein as well as the initiator of viral DNA replication. E2 is known to interact with various cellular DNA replication proteins, including the PV E1 protein, the cellular ssDNA binding complex (RPA), and topoisomerase I. Recently, we observed that cellular DNA polymerase ε (pol ε) interacts with the PV helicase protein, E1. E1 stimulates its activity with a very high degree of specificity, implicating pol ε in PV DNA replication. In this paper, we evaluated whether E2 also shows a functional interaction with pol ε. We found that E2 stimulates the DNA synthesis activity of pol ε, independently of pol ε’ s processivity factors, RFC, PCNA, and RPA, or E1. This appears to be specific for pol ε, as cellular DNA polymerase δ is unaffected by E1. However, unlike other known stimulatory factors of pol ε, E2 does not affect the processivity of pol ε. The domains of E2 were analyzed individually and in combination for their ability to stimulate pol ε. Both the transactivation and hinge domains were found to be important for this stimulation, while the E2 DNA-binding domain was dispensable. These findings support a role for E2 beyond E1 recruitment in viral DNA replication, demonstrate a novel functional interaction in PV DNA replication, and further implicate cellular pol ε in PV DNA replication.

The human papillomavirus replication cycle, and its links to cancer progression: a comprehensive review

HPVs (human papillomaviruses) infect epithelial cells and their replication cycle is intimately linked to epithelial differentiation. There are over 200 different HPV genotypes identified to date and each displays a strict tissue specificity for infection. HPV infection can result in a range of benign lesions, for example verrucas on the feet, common warts on the hands, or genital warts. HPV infects dividing basal epithelial cells where its dsDNA episomal genome enters the nuclei. Upon basal cell division, an infected daughter cell begins the process of keratinocyte differentiation that triggers a tightly orchestrated pattern of viral gene expression to accomplish a productive infection. A subset of mucosal-infective HPVs, the so-called ‘high risk’ (HR) HPVs, cause cervical disease, categorized as low or high grade. Most individuals will experience transient HR-HPV infection during their lifetime but these infections will not progress to clinically significant cervical disease or cancer because the immune system eventually recognizes and clears the virus. Cancer progression is due to persistent infection with an HR-HPV. HR-HPV infection is the cause of >99.7% cervical cancers in women, and a subset of oropharyngeal cancers, predominantly in men. HPV16 (HR-HPV genotype 16) is the most prevalent worldwide and the major cause of HPV-associated cancers. At the molecular level, cancer progression is due to increased expression of the viral oncoproteins E6 and E7, which activate the cell cycle, inhibit apoptosis, and allow accumulation of DNA damage. This review aims to describe the productive life cycle of HPV and discuss the roles of the viral proteins in HPV replication. Routes to viral persistence and cancer progression are also discussed.

Design of peptide inhibitors of human papillomavirus 16 (HPV16) transcriptional regulator E1–E2 formation

Here, we have proposed a new scheme of the computational combinatorial design approach to identify potential inhibitor peptides. It consists of four steps: (i) using “multiple copy simultaneous search” (MCSS) procedure to locate specific functional groups on the protein surface; (ii) the peptide main chain is constructed based on the location of favored N-methylacetamide (NMA) groups; (iii) molecular dynamics simulations of the complex formed between the constructed peptides with the target protein in explicit water molecules are carried to select the peptides with strong binding to the protein and (iv) the sequences of the stable peptides selected from (iii) are aligned and the frequencies of the amino acids at each position of peptide are calculated. Sequence patterns of potential inhibitors are determined based on the frequency of amino acids at each position. It was applied to design peptide inhibitors that bind to the E2 protein of HPV16 so as to disrupt its transcriptional regulator of E1–E2 complex formation. The sequence pattern of these potential inhibitors is in agreement with known inhibitors obtained from phage display, and the MCSS calculations indicate that a hydrophobic pocket on HPV16 E2 plays a significant role in E1–E2 formation and inhibitor-E2 binding.

HPV16-E2 protein modifies self-renewal and differentiation rate in progenitor cells of human immortalized keratinocytes

BACKGROUND

Cervical cancer is the fourth cause of death worldwide by cancer in women and is a disease associated to persistent infection with human papillomavirus (HPV), particularly from two high-risk types HPV16 and 18. The virus initiates its replicative cycle infecting cells located in the basal layer of the epithelium, where a small population of epithelial stem cells is located performing important functions of renewal and maintenance of the tissue. Viral E2 gene is one of the first expressed after infection and plays relevant roles in the replicative cycle of the virus, modifying fundamental processes in the infected cells. Thus, the aim of the present study was to demonstrate the presence of hierarchic subpopulations in HaCaT cell line and evaluate the effect of HPV16-E2 expression, on their biological processes.

METHODS

HaCaT-HPV16-E2 cells were generated by transduction of HaCaT cell line with a lentiviral vector. The α6-integrin-CD71 expression profile was established by immunostaining and flow cytometric analysis. After sorting, cell subpopulations were analyzed in biological assays for self-renewal, clonogenicity and expression of stemness factors (RT-qPCR).

RESULTS

We identified in HaCaT cell line three different subpopulations that correspond to early differentiated cells (α6-integrin^dim), transitory amplifying cells (α6-integrin^bri/CD71^bri) and progenitor cells (α6-integrin^bri/CD71^dim). The last subpopulation showed stem cell characteristics, such as self-renewal ability, clonogenicity and expression of the well-known stem cell factors SOX2, OCT4 and NANOG, suggesting they are stem-like cells. Interestingly, the expression of HPV16-E2 in HaCaT cells changed its α6-integrin-CD71 immunophenotype modifying the relative abundance of the cell subpopulations, reducing significantly the percentage of α6-integrin^bri/CD71^dim cells. Moreover, the expression of the stem cell markers was also modified, increasing the expression of SOX2 and NANOG, but decreasing notably the expression of OCT4.

CONCLUSIONS

Our data demonstrated the presence of a small subpopulation with epithelial "progenitor cells" characteristics in the HaCaT cell line, and that HPV16-E2 expression on these cells induces early differentiation.

Effect of sequence-dependent rigidity on plectoneme localization in dsDNA

We use Monte-Carlo simulations to study the effect of variable rigidity on plectoneme formation and localization in supercoiled double-stranded DNA. We show that the presence of soft sequences increases the number of plectoneme branches and that the edges of the branches tend to be localized at these sequences. We propose an experimental approach to test our results in vitro, and discuss the possible role played by plectoneme localization in the search process of transcription factors for their targets (promoter regions) on the bacterial genome.

Control of human papillomavirus gene expression by alternative splicing

•

Alternative splicing is a key cellular mechanism controlling HPV gene expression.

•

Many cellular SR proteins and hnRNPs have been identified that bind and control production of viral mRNAs.

•

HPV16 E2 protein controls expression of SR proteins and has splicing-related functions.

•

HPV16 infection through its regulatory effects on splicing factors may significantly alter cellular gene expression and cellular metabolism.

Human papillomaviruses possess circular double stranded DNA genomes of around 8 kb in size from which multiple mRNAs are synthesized during an infectious life cycle. Although at least three viral promoters are used to initiate transcription, viral mRNAs are largely the product of processing of pre-mRNAs by alternative splicing and polyadenylation. The HPV life cycle and viral gene expression are tightly linked to differentiation of the epithelium the virus infects: there is an orchestrated production of viral mRNAs and proteins. In this review we describe viral mRNA expression and the roles of the SR and hnRNP proteins that respectively positively and negatively regulate splicing. We discuss HPV regulation of splicing factors and detail the evidence that the papillomavirus E2 protein has splicing-related activities. We highlight the possibility that HPV-mediated control of splicing in differentiating epithelial cells may be necessary to accomplish the viral replication cycle.

Human Papillomavirus E2 Protein: Linking Replication, Transcription, and RNA Processing

The human papillomavirus (HPV) life cycle is tightly linked to differentiation of the infected epithelium. This means that viral proteins must exert control over epithelial gene expression in order to optimize viral production. The HPV E2 protein controls replication, transcription, and viral genome partitioning during the viral infectious life cycle. It consists of a nucleic acid-binding domain and a protein-protein interaction domain separated by a flexible serine and arginine-rich hinge region. Over the last few years, mounting evidence has uncovered an important new role for E2 in viral and cellular RNA processing. This Gem discusses the role of E2 in controlling the epithelial cellular environment and how E2 might act to coordinate late events in the viral replication cycle.

Design of peptides as inhibitors of human papillomavirus 16 transcriptional regulator E1-E2

Human papillomavirus 16 (HPV 16) is a DNA virus that is capable of infecting humans and causing cervical cancer. HPV16 E2 plays an important role in viral gene regulation. This work aims to predict the binding conformations and interactions between the dodecapeptides and HPV16 E2 as well as to design novel peptide inhibitors that are capable of binding to HPV16 E2 and disrupt the transcriptional regulator E1-E2 complex formation, using computational protein design techniques. Based on previously reported peptide4 (TWFWPYPYPHLP), novel peptide inhibitors were designed and five peptides that showed lower binding energy to HPV16 E2 than that of peptide4, were selected for in vitro experiments. Enzyme-linked immunosorbent (ELISA) assay showed that Y6R, W4H_Y6R, and W4H peptides bound to HPV16 E2 with higher affinity than peptide4 did. Moreover, Y6R, W4H_Y6R, and W4H peptides more effectively inhibited E1-E2 complex formation than peptide4. This work revealed important interactions between the peptides and E1-E2 complex, suggesting a strategy for development of more potent peptide inhibitors.

Molecular characterization of novel mucosotropic papillomaviruses from a Florida manatee (Trichechus manatus latirostris)

We isolated two new manatee papillomavirus (PV) types, TmPV3 and TmPV4, from a Florida manatee (Trichechus manatus latirostris). Two PV types were previously isolated from this species. TmPV1 is widely dispersed amongst manatees and a close-to-root PV; not much is known about TmPV2. The genomes of TmPV3 and TmPV4 were 7622 and 7771 bp in size, respectively. Both PVs had a genomic organization characteristic of all PVs, with one non-coding region and seven ORFs, including the E7 ORF that is absent in other cetacean PVs. Although these PVs were isolated from separate genital lesions of the same manatee, an enlarged E2/E4 ORF was found only in the TmPV4 genome. The full genome and L1 sequence similarities between TmPV3 and TmPV4 were 63.2 and 70.3 %, respectively. These genomes shared only 49.1 and 50.2 % similarity with TmPV1. The pairwise alignment of L1 nucleotide sequences indicated that the two new PVs nested in a monophyletic group of the genus Rhopapillomavirus, together with the cutaneotropic TmPV1 and TmPV2.

Immunization with Human Papillomavirus 16 L1+E2 Chimeric Capsomers Elicits Cellular Immune Response and Antitumor Activity in a Mouse Model

Development of cervical cancer is associated with persistent infections by high-risk human papillomavirus (HPV). Although current HPV L1-based prophylactic vaccines prevent infection, they do not help to eliminate prevalent infections or lesions. Our aims were (i) to generate a vaccine combining prophylactic and therapeutic properties by producing chimeric capsomers after fusion of the L1 protein to different fragments of E2 from HPV 16, and (ii) to evaluate their capacity to generate an antitumoral cellular response, while conserving L1 neutralizing epitopes. Chimeric proteins were produced in Escherichia coli and purified by glutathione S-transferase (GST)-affinity chromatography. Their structure was characterized using size exclusion chromatography, sucrose gradient centrifugation, electron microscopy, and anti-L1 enzyme-linked immunosorbent assay. All chimeric proteins form capsomers and heterogeneous aggregates. One, containing part of the carboxy-terminal domain of E2 and its hinge region (L1Δ+E2H/NC, aa 206-307), conserved the neutralizing epitope H16.V5. We then evaluated the capacity of this chimeric protein to induce a cytotoxic T-cell response against HPV 16 E2. In (51)Cr release cytotoxicity assays, splenocytes from C57BL/6 immunized mice recognized and lysed TC-1/E2 cells, which express and present endogenously processed E2 peptides. Moreover, this E2-specific cytotoxic response inhibited the growth of tumors of TC-1/E2 cells in mice. Finally, we identified an epitope (aa 292-301) of E2 involved in this cytotoxic response. We conclude that the L1Δ+E2H/NC chimeric protein produced in bacteria can be an effective and economically interesting candidate for a combined prophylactic and therapeutic vaccine that could help eliminating HPV16-positive low-grade cervical lesions and persistent viral infections, thus preventing the development of lesions and, at the same time, the establishment of new infections.

DNA Backbone BI/BII Distribution and Dynamics in E2 Protein-Bound Environment Determined by Molecular Dynamics Simulations

BI and BII conformational substates in the DNA backbone typify canonical B-form DNA. The BI and BII substates are important for structural variation of DNA and have been implicated in protein-nucleic acid recognition mechanisms. Recent refinements have been made to nucleic acid force fields employed in molecular dynamics simulations that demonstrate a better ability to model the BI and BII states, leading to overall improved modeling of DNA structure and dynamics. These force field improvements have yet to be significantly demonstrated in the context of a protein-DNA system extended to long time scales. Our plan was to run molecular dynamics simulations of a well-studied protein-DNA system (E2-DNA) into the microsecond time scale and determine the ability of the force field to populate BII states in the DNA backbone consistent with dinucleotide steps crystallized in the BII conformation. The results showed that the dinucleotide steps in the E2-DNA complex with the highest BII populations from simulation trajectories corresponded to the dinucleotide steps crystallized in the BII state and that decoy BI and BII states converge to the same results within approximately one microsecond.

Importance of β2-β3 Loop Motion for the Increased Binding and Decreased Selectivity of the ΔLL Mutant of the Human Papillomavirus Type 6 E2 Protein

The binding affinity of the human papillomavirus type 6 E2 protein is strongly mediated by the sequence of the DNA linker region, with high affinity for the AATT linker and low affinity for the CCGG linker. When two terminal leucine residues are removed from the protein, the level of binding to both strands increases, but unequally, resulting in a significant decrease in selectivity for the AATT linker strand. To rationalize this behavior, we performed molecular dynamics simulations of the wild-type and mutant protein in the apo state and bound to DNA with high-affinity AATT and low-affinity CCGG linker strands. While no stable contacts were made between the β2-β3 loop and DNA in the wild type, this loop was repositioned in the mutant complexes and formed electrostatic contacts with the DNA backbone. More contacts were formed when the mutant was bound to the CCGG linker strand than to the AATT linker strand, resulting in a more favorable change in interaction energy for the CCGG strand. In addition, significant differences in correlated motions were found, which further explained the differences in binding. The simulations suggest that β2-β3 loop motions are responsible for the increased affinity and decreased selectivity of the mutant protein.

Characterization of human papillomavirus type 154 and tissue tropism of gammapapillomaviruses

The novel human papillomavirus type 154 (HPV154) was characterized from a wart on the crena ani of a three-year-old boy. It was previously designated as the putative HPV type FADI3 by sequencing of a subgenomic FAP amplicon. We obtained the complete genome by combined methods including rolling circle amplification (RCA), genome walking through an adapted method for detection of integrated papillomavirus sequences by ligation-mediated PCR (DIPS-PCR), long-range PCR, and finally by cloning of four overlapping amplicons. Phylogenetically, the HPV154 genome clustered together with members of the proposed species Gammapapillomavirus 11, and demonstrated the highest identity in L1 to HPV136 (68.6%). The HPV154 was detected in 3% (2/62) of forehead skin swabs from healthy children. In addition, the different detection sites of 62 gammapapillomaviruses were summarized in order to analyze their tissue tropism. Several of these HPV types have been detected from multiple sources such as skin, oral, nasal, and genital sites, suggesting that the gammapapillomaviruses are generalists with a broader tissue tropism than previously appreciated. The study expands current knowledge concerning genetic diversity and tropism among HPV types in the rapidly growing gammapapillomavirus genus.

The papillomavirus E2 proteins

The papillomavirus E2 proteins are pivotal to the viral life cycle and have well characterized functions in transcriptional regulation, initiation of DNA replication and partitioning the viral genome. The E2 proteins also function in vegetative DNA replication, post-transcriptional processes and possibly packaging. This review describes structural and functional aspects of the E2 proteins and their binding sites on the viral genome. It is intended to be a reference guide to this viral protein.

Control of DNA minor groove width and Fis protein binding by the purine 2-amino group

The width of the DNA minor groove varies with sequence and can be a major determinant of DNA shape recognition by proteins. For example, the minor groove within the center of the Fis–DNA complex narrows to about half the mean minor groove width of canonical B-form DNA to fit onto the protein surface. G/C base pairs within this segment, which is not contacted by the Fis protein, reduce binding affinities up to 2000-fold over A/T-rich sequences. We show here through multiple X-ray structures and binding properties of Fis–DNA complexes containing base analogs that the 2-amino group on guanine is the primary molecular determinant controlling minor groove widths. Molecular dynamics simulations of free-DNA targets with canonical and modified bases further demonstrate that sequence-dependent narrowing of minor groove widths is modulated almost entirely by the presence of purine 2-amino groups. We also provide evidence that protein-mediated phosphate neutralization facilitates minor groove compression and is particularly important for binding to non-optimally shaped DNA duplexes.

HPV associated with recurrent respiratory papillomatosis

Papillomaviruses are members of the Papillomaviridae family. Over 150 HPV types have been identified. Recurrent respiratory papillomatosis (RRP) is a chronic condition caused by HPV characterized by recurrent papillomas of the respiratory tract, mainly the larynx. During the early stages, the condition presents with hoarseness, while more advanced disease presents with stridor and respiratory distress. There is no specific cure and treatment consists of repeated surgical procedures to remove the papillomas. Most patients eventually go into remission, but some suffer for many years with this condition, which may be fatal. HPV-6 and HPV-11 are the HPV types most commonly associated with RRP. Although most studies have found RRP due to HPV-11 to be more aggressive than disease due to HPV-6, the variability in disease aggressiveness is probably multifactorial. Information regarding the current epidemiology, molecular diversity and host immune responses is important for strategizing ways to reduce disease. Data on HPV genotypes associated with RRP would provide valuable information for vaccination programs to reduce the incidence of these genotypes in mothers and, in the long term, reduce the incidence of RRP in children. This review focuses on HPV-6 and HPV-11 as the HPV types that cause RRP, and discusses the viral genome and replication, clinical presentation of RRP, current techniques of diagnosis and genotyping, and the molecular diversity of HPV-6 and HPV-11.

Genetic analysis of the E2 transactivation domain dimerization interface from bovine papillomavirus type 1

The bovine papillomavirus type 1 (BPV1) E2 protein binds as a dimer to the viral genome to promote its transcription, replication and maintenance in keratinocytes. Although BPV1 E2 dimerizes primarily through its DNA-binding domain, it was shown previously that its transactivation domain (TAD) can also dimerize in vitro through formation of a disulfide bond between cysteine 57 (C57) of adjacent monomers and of an ion pair between arginine 172 (R172) and aspartic acid 175 (D175). The function of this TAD dimerization interface in vivo remains unknown. Here, we report the effects of substituting C57, R172 and D175 by alanine on the transactivation activity of BPV E2 as well as on its ability to support viral DNA replication using a novel luciferase-based assay. Results for this mutational analysis suggest that the TAD dimerization interface is not essential for either process but may contribute to the DNA replication activity of BPV1 E2.

Treatment failure in patients with HPV 16-induced vulvar intraepithelial neoplasia: understanding different clinical responses to immunotherapy

Failure of the immune system to launch a strong and effective immune response to high-risk HPV is related to viral persistence and the development of anogenital (pre)malignant lesions such as vulvar intraepithelial neoplasia (VIN). Different forms of immunotherapy, aimed at overcoming the inertia of the immune system, have been developed and met with clinical success. Unfortunately these, in principal successful, therapeutic approaches also fail to induce clinical responses in a substantial number of cases. In this review, the authors summarize the traits of the immune response to HPV in healthy individuals and in patients with HPV-induced neoplasia. The potential mechanisms involved in the escape of HPV-induced lesions from the immune system indicate gaps in our knowledge. Finally, the interaction between the immune system and VIN is discussed with a special focus on the different forms of immunotherapy applied to treat VIN and the potential causes of therapy failure. The authors conclude that there are a number of pre-existing conditions that determine the patients' responsiveness to immunotherapy. An immunotherapeutic strategy in which different aspects of immune failure are attacked by complementary approaches, will improve the clinical response rate.

Identification of a novel bat papillomavirus by metagenomics

The discovery of novel viruses in animals expands our knowledge of viral diversity and potentially emerging zoonoses. High-throughput sequencing (HTS) technology gives millions or even billions of sequence reads per run, allowing a comprehensive survey of the genetic content within a sample without prior nucleic acid amplification. In this study, we screened 156 rectal swab samples from apparently healthy bats (n = 96), pigs (n = 9), cattles (n = 9), stray dogs (n = 11), stray cats (n = 11) and monkeys (n = 20) using a HTS metagenomics approach. The complete genome of a novel papillomavirus (PV), Miniopterus schreibersii papillomavirus type 1 (MscPV1), with L1 of 60% nucleotide identity to Canine papillomavirus (CPV6), was identified in a specimen from a Common Bent-wing Bat (M. schreibersii). It is about 7.5kb in length, with a G+C content of 45.8% and a genomic organization similar to that of other PVs. Despite the higher nucleotide identity between the genomes of MscPV1 and CPV6, maximum-likelihood phylogenetic analysis of the L1 gene sequence showed that MscPV1 and Erethizon dorsatum papillomavirus (EdPV1) are most closely related. Estimated divergence time of MscPV1 from the EdPV1/MscPV1 common ancestor was approximately 60.2–91.9 millions of years ago, inferred under strict clocks using the L1 and E1 genes. The estimates were limited by the lack of reliable calibration points from co-divergence because of possible host shifts. As the nucleotide sequence of this virus only showed limited similarity with that of related animal PVs, the conventional approach of PCR using consensus primers would be unlikely to have detected the novel virus in the sample. Unlike the first bat papillomavirus RaPV1, MscPV1 was found in an asymptomatic bat with no apparent mucosal or skin lesions whereas RaPV1 was detected in the basosquamous carcinoma of a fruit bat Rousettus aegyptiacus. We propose MscPV1 as the first member of the novel Dyolambda-papillomavirus genus.

Genital human papillomavirus infections: current and prospective therapies

Infection with human papillomaviruses (HPVs) is very common and associated with benign and malignant epithelial proliferations of skin and internal squamous mucosae. A subset of the mucosal HPVs are oncogenic and associated with 5 % of all cancers in men and women. There are two licensed prophylactic vaccines, both target HPV 16 and 18, the two most pathogenic, oncogenic types and one, additionally, targets HPV 6 and 11 the cause of genital warts. The approach of deliberate immunization with oncogenic HPV E6 and/or E7 proteins and the generation of antigen-specific cytotoxic T-cells as an immunotherapy for HPV-associated cancer and their high-grade pre-cancers has been tested with a wide array of potential vaccine delivery systems in Phase I/II trials with varying success. Understanding local viral and tumour immune evasion strategies is a prerequisite for the rational design of therapeutic vaccines for HPV-associated infection and disease, progress in this is discussed. There are no antiviral drugs for the treatment of HPV infection and disease. Current therapies are not targeted antiviral therapies, but either attempt physical removal of the lesion or induce inflammation and a bystander immune response. There has been recent progress in the identification and characterization of molecular targets for small molecule antagonists of the HPV proteins E1, E2 and E6 or their interactions with their cellular targets. Lead compounds that could disrupt E1-E2 protein-protein interactions have been discovered as have inhibitors of E6-E6-AP-binding interactions. Some of these compounds showed nanomolar affinities and high specificities and demonstrate the feasibility of this approach for HPV infections. These studies are, however, at an early phase and it is unlikely that any specific anti-HPV chemotherapeutic will be in the clinic within the next 10-20 years.

Evolutionary variation of papillomavirus E2 protein and E2 binding sites

Background

In an effort to identify the evolutionary changes relevant to E2 function, within and between papillomavirus genera, we evaluated the E2 binding sites (E2BS)s inside the long-control-region (LCR), and throughout the genomes. We identified E2BSs in the six largest genera of papillomaviruses: Alpha, Beta, Gamma, Delta, Lambda, and Xi-papillomaviruses (128 genomes), by comparing the sequences with a model consensus we created from known functional E2BSs (HPV16, HPV18, BPV1). We analyzed the sequence conservation and nucleotide content of the 4-nucleotide spacer within E2BSs. We determined that there is a statistically significant difference in GC content of the four-nucleotide E2BS spacer, between Alpha and Delta-papillomaviruses, as compared to each of the other groups. Additionally, we performed multiple alignments of E2 protein sequences using members of each genus in order to identify evolutionary changes within the E2 protein.

Results

When a phylogenetic tree was generated from E2 amino acid sequences, it was discovered that the alpha-papillomavirus genera segregates into two distinct subgroups (α1 and α2). When these subgroups were individually analyzed, it was determined that the subgroup α1 consensus E2BS favored a spacer of AAAA, whereas subgroup α2 favored the opposite orientation of the same spacer; TTTT. This observation suggests that these conserved inverted linkers could have functional importance.

Characterization of the complete genomes of Camelus dromedarius papillomavirus types 1 and 2

Camel papillomatosis has been described previously, but the genome of the suspected papillomavirus (PV) has not been identified. An outbreak of papillomatosis occurred in a dromedary farm of 55 animals in Sudan during August 2009. The disease was only present in young animals aged about 3-7 months, of which 44 % (11/25) were affected with lesions, mainly on the lips and lower jaw. This study reports for the first time the complete genomes of Camelus dromedarius papillomavirus types 1 (CdPV1) and 2 (CdPV2), isolated from a cauliflower-like nodule and a round oval raised nodule, respectively. Pairwise comparisons of their L1 nucleotide sequences revealed 69.2 % identity, and phylogenetic analyses suggested that these two PV types are grouped within the genus Deltapapillomavirus. Both viruses were isolated from fibropapillomas, although no putative E5 proteins homologous to that of bovine papillomavirus type 1 were identified. The genetic information will be useful for evolutionary studies of the family Papillomaviridae, as well as for the development of diagnostic methods for surveillance of the disease in dromedaries.

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.

GRIM-19 disrupts E6/E6AP complex to rescue p53 and induce apoptosis in cervical cancers

BACKGROUND

Our previous studies showed a down-regulation of GRIM-19 in primary human cervical cancers, and restoration of GRIM-19 induced tumor regression. The induction of tumor suppressor protein p53 ubiquitination and degradation by E6 oncoportein of high risk-HPV through forming a stable complex with E6AP is considered as a critical mechanism for cervical tumor development. The aims of this study were to determine the potential role of GRIM-19 in rescuing p53 protein and inducing cervical cancer cell apoptosis.

METHODOLOGY/PRINCIPAL FINDINGS

The protein levels of GRIM-19 and p53 were detected in normal cervical tissues from 45 patients who underwent hysterectomy for reasons other than neoplasias of either the cervix or endometrium, and cervical cancer tissues from 60 patients with non-metastatic squamous epithelial carcinomas. Coimmunoprecipitation and GST pull-down assay were performed to examine the interaction of GRIM-19 with 18E6 and E6AP in vivo and in vitro respectively. The competition of 18E6 with E6AP in binding GRIM-19 by performing competition pull-down assays was designed to examine the disruption of E6/E6AP complex by GRIM-19. The augment of E6AP ubiquitination by GRIM-19 was detected in vivo and in vitro ubiquitination assay. The effects of GRIM-19-dependent p53 accumulation on cell proliferation, cell cycle, apoptosis were explored by MTT, flow cytometry and transmission electron microscopy respectively. The tumor suppression was detected by xenograft mouse model.

CONCLUSION/SIGNIFICANCE

The levels of GRIM-19 and p53 were concurrently down regulated in cervical cancers. The restoration of GRIM-19 can induce ubiquitination and degradation of E6AP, and disrupt the E6/E6AP complex through the interaction of N-terminus of GRIM-19 with both E6 and E6AP, which protected p53 from degradation and promoted cell apoptosis. Tumor xenograft studies also revealed the suppression of p53 degradation in presence of GRIM-19. These data suggest that GRIM-19 can block E6/E6AP complex; and synergistically suppress cervical tumor growth with p53.

Small molecule inhibitors of human papillomavirus protein - protein interactions

Human papillomaviruses (HPV) have now been identified as a necessary cause of benign and malignant lesions of the differentiating epithelium, particularly cervical cancer, the second most prevalent cancer in women worldwide. While two prophylactic HPV vaccines and screening programs are available, there is currently no antiviral drug for the treatment of HPV infections and associated diseases. The recent progress toward the identification and characterization of specific molecular targets for small molecule-based approaches provides prospect for the development of effective HPV antiviral compounds. Traditionally, antiviral therapies target viral enzymes. HPV encode for few proteins, however, and rely extensively on the infected cell for completion of their life cycle. This article will review the functions of the viral E1 helicase, which encodes the only enzymatic function of the virus, of the E2 regulatory protein, and of the viral E6 and E7 oncogenes in viral replication and pathogenesis. Particular emphasis will be placed on the recent progress made towards the development of novel small molecule inhibitors that specifically target and inhibit the functions of these viral proteins, as well as their interactions with other viral and/or cellular proteins.

In vivo analysis of the cell cycle dependent association of the bovine papillomavirus E2 protein and ChlR1

It has been shown that the genomes of episomally maintained DNA viruses are tethered to host cell chromosomes during cell division, facilitating maintenance in dividing cells. The papillomavirus E2 protein serves this mechanism of viral genome persistence by simultaneously associating with chromatin and the viral genome during mitosis. Several host cell proteins are reported to be necessary for the association of E2 with chromatin including the cohesion establishment factor ChlR1. Here we use fluorescence resonance energy transfer (FRET) technology to confirm the interaction between BPV-1 E2 and ChlR1. Furthermore, we use synchronised live cells to study the temporal nature of this dynamic protein interaction and show that ChlR1 and E2 interact during specific phases of the cell cycle. These data provide evidence that the association of E2 with ChlR1 contributes to a loading mechanism during DNA replication rather than direct tethering during mitotic division.

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.

Thermodynamics of cooperative DNA recognition at a replication origin and transcription regulatory site

Binding cooperativity guides the formation of protein−nucleic acid complexes, in particular those that are highly regulated such as replication origins and transcription sites. Using the DNA binding domain of the origin binding and transcriptional regulator protein E2 from human papillomavirus type 16 as model, and through isothermal titration calorimetry analysis, we determined a positive, entropy-driven cooperativity upon binding of the protein to its cognate tandem double E2 site. This cooperativity is associated with a change in DNA structure, where the overall B conformation is maintained. Two homologous E2 domains, those of HPV18 and HPV11, showed that the enthalpic−entropic components of the reaction and DNA deformation can diverge. Because the DNA binding helix is almost identical in the three domains, the differences must lie dispersed throughout this unique dimeric β-barrel fold. This is in surprising agreement with previous results for this domain, which revealed a strong coupling between global dynamics and DNA recognition.

Origins of specificity in protein-DNA recognition

Specific interactions between proteins and DNA are fundamental to many biological processes. In this review, we provide a revised view of protein-DNA interactions that emphasizes the importance of the three-dimensional structures of both macromolecules. We divide protein-DNA interactions into two categories: those when the protein recognizes the unique chemical signatures of the DNA bases (base readout) and those when the protein recognizes a sequence-dependent DNA shape (shape readout). We further divide base readout into those interactions that occur in the major groove from those that occur in the minor groove. Analogously, the readout of the DNA shape is subdivided into global shape recognition (for example, when the DNA helix exhibits an overall bend) and local shape recognition (for example, when a base pair step is kinked or a region of the minor groove is narrow). Based on the >1500 structures of protein-DNA complexes now available in the Protein Data Bank, we argue that individual DNA-binding proteins combine multiple readout mechanisms to achieve DNA-binding specificity. Specificity that distinguishes between families frequently involves base readout in the major groove, whereas shape readout is often exploited for higher resolution specificity, to distinguish between members within the same DNA-binding protein family.

The shape of the DNA minor groove directs binding by the DNA-bending protein Fis

The bacterial nucleoid-associated protein Fis regulates diverse reactions by bending DNA and through DNA-dependent interactions with other control proteins and enzymes. In addition to dynamic nonspecific binding to DNA, Fis forms stable complexes with DNA segments that share little sequence conservation. Here we report the first crystal structures of Fis bound to high- and low-affinity 27-base-pair DNA sites. These 11 structures reveal that Fis selects targets primarily through indirect recognition mechanisms involving the shape of the minor groove and sequence-dependent induced fits over adjacent major groove interfaces. The DNA shows an overall curvature of approximately 65 degrees , and the unprecedented close spacing between helix-turn-helix motifs present in the apodimer is accommodated by severe compression of the central minor groove. In silico DNA structure models show that only the roll, twist, and slide parameters are sufficient to reproduce the changes in minor groove widths and recreate the curved Fis-bound DNA structure. Models based on naked DNA structures suggest that Fis initially selects DNA targets with intrinsically narrow minor grooves using the separation between helix-turn-helix motifs in the Fis dimer as a ruler. Then Fis further compresses the minor groove and bends the DNA to generate the bound structure.

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.

Anomalous DNA binding by E2 regulatory protein driven by spacer sequence TATA

We have investigated the anomalously weak binding of human papillomavirus (HPV) regulatory protein E2 to a DNA target containing the spacer sequence TATA. Experiments in magnesium (Mg(2+)) and calcium (Ca(2+)) ion buffers revealed a marked reduction in cutting by DNase I at the CpG sequence in the protein-binding site 3' to the TATA spacer sequence, Studies of the cation dependence of DNA-E2 affinities showed that upon E2 binding the TATA sequence releases approximately twice as many Mg(2+) ions as the average of the other spacer sequences. Binding experiments for TATA spacer relative to ATAT showed that in potassium ion (K(+)) the E2 affinity of the two sequences is nearly equal, but the relative dissociation constant (K(d)) for TATA increases in the order K(+ )< Na(+ )< Ca(2+ )< Mg(2+). Except for Mg(2+), K(d) for TATA relative to ATAT is independent of ion concentration, whereas for Mg(2+) the affinity for TATA drops sharply as ion concentration increases. Thus, ions of increasing positive charge density increasingly distort the E2 binding site, weakening the affinity for protein. In the case of Mg(2+), additional ions are bound to TATA that require displacement for protein binding. We suggest that the TATA sequence may bias the DNA structure towards a conformation that binds the protein relatively weakly.

The human papillomavirus E7-E2 interaction mechanism in vitro reveals a finely tuned system for modulating available E7 and E2 proteins

Transcription of the human papillomavirus E7 oncoprotein is negatively controlled by the viral E2 protein, and loss of this repression leads to irreversible transformation and carcinogenesis. Here we show that interaction of the HPV16 E7 protein with the DNA binding domain of the E2 protein (E2C) leads to ionic strength-dependent hetero-oligomerization even at the lowest concentrations measurable. Titration experiments followed by light scattering and native gel electrophoresis show insoluble oligomeric complexes with a >or=2000 nm diameter and intermediate soluble complexes 40 and 115 nm in diameter, respectively, formed in excess of E2C. A discrete oligomeric soluble complex formed in excess of E7 displays a diameter of 12 nm. The N-terminal domain of E7 interacts with E2C with a K(D) of 0.1 muM, where the stretch of residues 25-40 of E7, encompassing both a PEST motif and phosphorylation sites, is sufficient for the interaction. Displacement of the soluble E7-E2C complex by an E2 site DNA duplex and site-directed mutagenesis indicate that the protein-protein interface involves the DNA binding helix of E2. The formation of complexes of different sizes and properties in excess of either of the viral proteins reveals a finely tuned mechanism that could regulate the intracellular levels of both proteins as infection and transformation progress. Sequestering E2 into E7-E2 oligomers provides a possible additional route to uncontrolled E7 expression, in addition and prior to the disruption of the E2 gene during viral integration into the host genome.

Intrinsic flexibility of B-DNA: the experimental TRX scale

B-DNA flexibility, crucial for DNA–protein recognition, is sequence dependent. Free DNA in solution would in principle be the best reference state to uncover the relation between base sequences and their intrinsic flexibility; however, this has long been hampered by a lack of suitable experimental data. We investigated this relationship by compiling and analyzing a large dataset of NMR ³¹P chemical shifts in solution. These measurements reflect the BI ↔ BII equilibrium in DNA, intimately correlated to helicoidal descriptors of the curvature, winding and groove dimensions. Comparing the ten complementary DNA dinucleotide steps indicates that some steps are much more flexible than others. This malleability is primarily controlled at the dinucleotide level, modulated by the tetranucleotide environment. Our analyses provide an experimental scale called TRX that quantifies the intrinsic flexibility of the ten dinucleotide steps in terms of Twist, Roll, and X-disp (base pair displacement). Applying the TRX scale to DNA sequences optimized for nucleosome formation reveals a 10 base-pair periodic alternation of stiff and flexible regions. Thus, DNA flexibility captured by the TRX scale is relevant to nucleosome formation, suggesting that this scale may be of general interest to better understand protein-DNA recognition.

A generalized conformational energy function of DNA derived from molecular dynamics simulations

Proteins recognize DNA sequences by two different mechanisms. The first is direct readout, in which recognition is mediated by direct interactions between the protein and the DNA bases. The second is indirect readout, which is caused by the dependence of conformation and the deformability of the DNA structure on the sequence. Various energy functions have been proposed to evaluate the contribution of indirect readout to the free-energy changes in complex formations. We developed a new generalized energy function to estimate the dependence of the deformability of DNA on the sequence. This function was derived from molecular dynamics simulations previously conducted on B-DNA dodecamers, each of which had one possible tetramer sequence embedded at its center. By taking the logarithm of the probability distribution function (PDF) for the base-step parameters of the central base-pair step of the tetramer, its ability to distinguish the native sequence from random ones was superior to that with the previous method that approximated the energy function in harmonic form. From a comparison of the energy profiles calculated with these two methods, we found that the harmonic approximation caused significant errors in the conformational energies of the tetramers that adopted multiple stable conformations.