The structural basis of DNA target discrimination by papillomavirus E2 proteins

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.

Structural predictions of protein-DNA binding: MELD-DNA

Structural, regulatory and enzymatic proteins interact with DNA to maintain a healthy and functional genome. Yet, our structural understanding of how proteins interact with DNA is limited. We present MELD-DNA, a novel computational approach to predict the structures of protein-DNA complexes. The method combines molecular dynamics simulations with general knowledge or experimental information through Bayesian inference. The physical model is sensitive to sequence-dependent properties and conformational changes required for binding, while information accelerates sampling of bound conformations. MELD-DNA can: (i) sample multiple binding modes; (ii) identify the preferred binding mode from the ensembles; and (iii) provide qualitative binding preferences between DNA sequences. We first assess performance on a dataset of 15 protein-DNA complexes and compare it with state-of-the-art methodologies. Furthermore, for three selected complexes, we show sequence dependence effects of binding in MELD predictions. We expect that the results presented herein, together with the freely available software, will impact structural biology (by complementing DNA structural databases) and molecular recognition (by bringing new insights into aspects governing protein-DNA interactions).

Accurate modeling of DNA conformational flexibility by a multivariate Ising model

The sequence-dependent structure and deformability of DNA play a major role for binding of proteins and regulation of gene expression. So far, most efforts to model DNA flexibility are based on unimodal harmonic stiffness models at base-pair resolution. However, multimodal behavior due to distinct conformational substates also contributes significantly to the conformational flexibility of DNA. Moreover, these local substates are correlated to their nearest-neighbor substates. A description for DNA elasticity which includes both multimodality and nearest-neighbor coupling has remained a challenge, which we solve by combining our multivariate harmonic approximation with an Ising model for the substates. In a series of applications to DNA fluctuations and protein-DNA complexes, we demonstrate substantial improvements over the unimodal stiffness model. Furthermore, our multivariate Ising model reveals a mechanical destabilization for adenine (A)-tracts to undergo nucleosome formation. Our approach offers a wide range of applications to determine sequence-dependent deformation energies of DNA and to investigate indirect readout contributions to protein-DNA recognition.

High Levels of Within-Host Variations of Human Papillomavirus 16 E1/E2 Genes in Invasive Cervical Cancer

Human papillomavirus type 16 (HPV16) is the most common HPV genotype found in invasive cervical cancer (ICC). Recent comprehensive genomics studies of HPV16 have revealed that a large number of minor nucleotide variations in the viral genome are present in each infected woman; however, it remains unclear whether such within-host variations of HPV16 are linked to cervical carcinogenesis. Here, by employing next-generation sequencing approaches, we explored the mutational profiles of the HPV16 genome within individual clinical specimens from ICC (n = 31) and normal cervix (n = 21) in greater detail. A total of 367 minor nucleotide variations (167 from ICC and 200 from the normal cervix) were detected throughout the viral genome in both groups, while nucleotide variations at high frequencies (>10% abundance in relative read counts in a single sample) were more prevalent in ICC (10 in ICC versus 1 in normal). Among the high-level variations found in ICC, six were located in the E1/E2 genes, and all of them were non-synonymous substitutions (Q142K, M207I, and L262V for E1; D153Y, R302T, and T357A for E2). In vitro functional analyses of these E1/E2 variants revealed that E1/M207I, E2/D153Y, and E2/R302T had reduced abilities to support viral replication, and that E2/D153Y and E2/R302T failed to suppress the viral early promoter. These results imply that some within-host variations of E1/E2 present at high levels in ICC may be positively selected for and contribute to cervical cancer development through dysfunction or de-stabilization of viral replication/transcription proteins.

Cytotoxicity, Pro-apoptotic Activity and in silico Studies of Dithiocarbamates and their Structure Based Design and SAR Studies

Background:

Cancer is a far-reaching and lethal but curable disease. Researchers have investigated numerous anticancer agents with only a few commercially available effective drugs which are very costly.

Objective:

Herein, we report the synthesis , characterization and anti cancer assays of a series of novel dithiocarbamates derivatives.

Methods:

All compounds were synthesized from different secondary amines and substituted benzyl chlorides in a single step. The structures of newly synthesized dithiocarbamate derivatives were confirmed by spectroscopic techniques (IR, NMR and HR-MS).

Results:

The synthesized compounds showed a significant anti-proliferative effect in cancer cells (HeLa) with the maximum inhibitory activity of compound SHD-2 with an IC50 = 0.31 ± 0.09 μM. However, the same compound exhibited 19.2% inhibition towards Baby Hamster Kidney fibroblasts (BHK-21), normal cell lines. Moreover, quantification of cellular DNA by flow cytometry for the evaluation of pro-apoptotic activity in HeLa cells demonstrates that arrest in cell cycle along with apoptosis advance towards drug cytotoxicity. However, molecular docking studies of the potent compound suggested that it binds to the major groove of the DNA.

Conclusion:

The cytotoxic and pro-apoptotic potential of the potent inhibitor may be further investigated in the animal models to advance their anti-cancer prospective.

E2 protein is the major determinant of specificity at the human papillomavirus origin of replication

The replication of human papillomavirus (HPV) genomes requires E1 and E2 proteins as the viral trans-factors and the replication origin, located in the URR, as a cis-element. The minimal requirements for an HPV replication origin vary among different virus types but always include one or more binding sites for the E2 protein. The requirements for an E1 binding site seem to vary among different HPV genera, with alpha-HPV11 and -18 minimal origins able to replicate without E1 binding site in contrast to beta-HPV8. In the present article, we analysed the sequence requirements for the beta-HPV5 minimal origin of replication. We show that the HPV5 URR is able to replicate in U2OS cells without the sequence proposed as an E1 binding site, albeit at lower levels than wt URR, given that three E2 binding sites are intact and both viral replication proteins are present. The lack of an absolute requirement of the E1 binding site for the origin of replication of HPV5 led us to analyse whether the viral E1 and E2 proteins from other HPV types are competent to support replication from this origin. Surprisingly, the E1 and E2 proteins from beta-HPV types support replication from the origin in contrast to proteins from alpha-HPV types 11, -16, or -18. Furthermore, the replication proteins E1 and E2 of these alpha-HPV types are unable to support the replication of HPV5 URR, even if the E1 binding site is intact. In light of these results, we performed a detailed analysis of the ability of different combinations of E1 and E2 proteins from various alpha- and beta-HPV types to support the replication of URR sequences from the respective HPV types in the U2OS cell line.

Pre- and post-docking sampling of conformational changes using ClustENM and HADDOCK for protein-protein and protein-DNA systems

Incorporating the dynamic nature of biomolecules in the modeling of their complexes is a challenge, especially when the extent and direction of the conformational changes taking place upon binding is unknown. Estimating whether the binding of a biomolecule to its partner(s) occurs in a conformational state accessible to its unbound form (“conformational selection”) and/or the binding process induces conformational changes (“induced‐fit”) is another challenge. We propose here a method combining conformational sampling using ClustENM—an elastic network‐based modeling procedure—with docking using HADDOCK, in a framework that incorporates conformational selection and induced‐fit effects upon binding. The extent of the applied deformation is estimated from its energetical costs, inspired from mechanical tensile testing on materials. We applied our pre‐ and post‐docking sampling of conformational changes to the flexible multidomain protein‐protein docking benchmark and a subset of the protein‐DNA docking benchmark. Our ClustENM‐HADDOCK approach produced acceptable to medium quality models in 7/11 and 5/6 cases for the protein‐protein and protein‐DNA complexes, respectively. The conformational selection (sampling prior to docking) has the highest impact on the quality of the docked models for the protein‐protein complexes. The induced‐fit stage of the pipeline (post‐sampling), however, improved the quality of the final models for the protein‐DNA complexes. Compared to previously described strategies to handle conformational changes, ClustENM‐HADDOCK performs better than two‐body docking in protein‐protein cases but worse than a flexible multidomain docking approach. However, it does show a better or similar performance compared to previous protein‐DNA docking approaches, which makes it a suitable alternative.

Cycloaddition reactions for antiviral compounds

Prominent in the current stage of drug development, antiviral compounds can be efficiently prepared through cycloaddition reactions. The chapter reports the use of classical Diels–Alder and their hetero version for the design and synthesis of compounds that were tested for their antiviral activities against a variety of viruses. Furthermore, 1,3-dipolar cycloaddition reactions of selected 1,3-dipoles, such as azides, nitrones, and nitrile oxides, are reviewed in the light of their application in the preparation of key intermediates for antiviral synthesis. A few examples of [2+2] cycloaddition reactions are also presented. The products obtained from these pericyclic reaction approaches were all tested for their activities in terms of blocking the virus replication, and the relevant biological data are highlighted.

Genetic variations in the DNA replication origins of human papillomavirus family correlate with their oncogenic potential

Human papillomaviruses (HPVs) encompass a large family of viruses that range from benign to highly carcinogenic. The crucial differences between benign and carcinogenic types of HPV remain unknown, except that the two HPV types differ in the frequency of DNA replication. We have systematically analyzed the mechanism of HPV DNA replication initiation in low-risk and high-risk HPVs. Our results demonstrate that HPV-encoded E2 initiator protein and its four binding sites in the replication origin play pivotal roles in determining the destiny of the HPV-infected cell. We have identified strain-specific single nucleotide variations in E2 binding sites found only in the high-risk HPVs. We have demonstrated that these variations result in attenuated formation of the E2-DNA complex. E2 binding to these sites is linked to the activation of the DNA replication origin as well as initiation of DNA replication. Both electrophoretic mobility shift assay and atomic force microscopy studies demonstrated that binding of E2 from either low- or high-risk HPVs with variant binding sequences lacked multimeric E2-DNA complex formation in vitro. These results provided a molecular basis of differential DNA replication in the two types of HPVs and pointed to a correlation with the development of cancer.

Parallels between DNA and collagen - comparing elastic models of the double and triple helix

Multi-stranded helices are widespread in nature. The interplay of polymeric properties with biological function is seldom discussed. This study probes analogies between structural and mechanical properties of collagen and DNA. We modeled collagen with Eulerian rotational and translational parameters of adjacent rungs in the triple-helix ladder and developed statistical potentials by extracting the dispersion of the parameters from a database of atomic-resolution structures. The resulting elastic model provides a common quantitative way to describe collagen deformations upon interacting with integrins or matrix metalloproteinase and DNA deformations upon protein binding. On a larger scale, deformations in Type I collagen vary with a periodicity consistent with the D-periodic banding of higher-order fibers assemblies. This indicates that morphologies of natural higher-order collagen packing might be rooted in the characteristic deformation patterns.

Involvement of Brd4 in different steps of the papillomavirus life cycle

Bromodomain-containing protein 4 (Brd4) is a cellular chromatin-binding factor and transcriptional regulator that recruits sequence-specific transcription factors and chromatin modulators to control target gene transcription. Papillomaviruses (PVs) have evolved to hijack Brd4's activity in order to create a facilitating environment for the viral life cycle. Brd4, in association with the major viral regulatory protein E2, is involved in multiple steps of the PV life cycle including replication initiation, viral gene transcription, and viral genome segregation and maintenance. Phosphorylation of Brd4, regulated by casein kinase II (CK2) and protein phosphatase 2A (PP2A), is critical for viral gene transcription as well as E1- and E2-dependent origin replication. Thus, pharmacological agents regulating Brd4 phosphorylation and inhibitors blocking phospho-Brd4 functions are promising candidates for therapeutic intervention in treating human papillomavirus (HPV) infections as well as associated disease.

BRD4 Phosphorylation Regulates HPV E2-Mediated Viral Transcription, Origin Replication, and Cellular MMP-9 Expression

Post-translational modification can modulate protein conformation and alter binding partner recruitment within gene regulatory regions. Here, we report that bromodomain-containing protein 4 (BRD4), a transcription co-factor and chromatin regulator, uses a phosphorylation-induced switch mechanism to recruit E2 protein encoded by cancer-associated human papillomavirus (HPV) to viral early gene and cellular matrix metalloproteinase-9 (MMP-9) promoters. Enhanced MMP-9 expression, induced upon keratinocyte differentiation, occurs via BRD4-dependent recruitment of active AP-1 and NF-κB to their target sequences. This is triggered by replacement of AP-1 family members JunB and JunD by c-Jun and by re-localization of NF-κB from the cytoplasm to the nucleus. In addition, BRD4 phosphorylation is critical for E2- and origin-dependent HPV DNA replication. A class of phospho-BRD4-targeting compounds, distinct from the BET bromodomain inhibitors, effectively blocks BRD4 phosphorylation-specific functions in transcription and factor recruitment.

Structural and biochemical insights into the DNA-binding mode of MjSpt4p:Spt5 complex at the exit tunnel of RNAPII

Spt5 (NusG in bacteria) is the only RNA polymerase-associated factor known to be conserved in all three domains of life. In archaea and eukaryotes, Spt5 associates with Spt4, an elongation factor that is absent in bacteria, to form a functional heterodimeric complex. Previous studies suggest that the Spt4:Spt5 complex interacts directly with DNA at the double-stranded DNA exit tunnel of RNA polymerase to regulate gene transcription. In this study, the DNA-binding ability of Spt4:Spt5 from the archaeon Methanocaldococcus jannaschii was confirmed via nuclear magnetic resonance chemical shift perturbation and fluorescence polarization assays. Crystallographic analysis of the full-length MjSpt4:Spt5 revealed two distinct conformations of the C-terminal KOW domain of Spt5. A similar alkaline region was found on the Spt4:Spt5 surface in both crystal forms, and identified as double-stranded DNA binding patch through mutagenesis-fluorescence polarization assays. Based on these structural and biochemical data, the Spt4:Spt5-DNA binding model was built for the first time.

Human Papillomavirus Type 18 cis-Elements Crucial for Segregation and Latency

Stable maintenance replication is characteristic of the latency phase of HPV infection, during which the viral genomes are actively maintained as extrachromosomal genetic elements in infected proliferating basal keratinocytes. Active replication in the S-phase and segregation of the genome into daughter cells in mitosis are required for stable maintenance replication. Most of our knowledge about papillomavirus genome segregation has come from studies of bovine papillomavirus type 1 (BPV-1), which have demonstrated that the E2 protein cooperates with cellular trans-factors and that E2 binding sites act as cis-regulatory elements in the viral genome that are essential for the segregation process. However, the genomic organization of the regulatory region in HPVs, and the properties of the viral proteins are different from those of their BPV-1 counterparts. We have designed a segregation assay for HPV-18 and used it to demonstrate that the E2 protein performs segregation in combination with at least two E2 binding sites. The cooperative binding of the E2 protein to two E2 binding sites is a major determinant of HPV-18 genome segregation, as demonstrated by the change in spacing between adjacent binding sites #1 and #2 in the HPV-18 Upstream Regulatory Region (URR). Duplication or triplication of the natural 4 bp 5’-CGGG-3’ spacer between the E2 binding sites increased the cooperative binding of the E2 molecules as well as E2-dependent segregation. Removal of any spacing between these sites eliminated cooperative binding of the E2 protein and disabled segregation of the URR and HPV-18 genome. Transfer of these configurations of the E2 binding sites into viral genomes confirmed the role of the E2 protein and binding sites #1 and #2 in the segregation process. Additional analysis demonstrated that these sites also play an important role in the transcriptional regulation of viral gene expression from different HPV-18 promoters.

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.

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.

Replication of Epstein-Barr viral DNA

Epstein-Barr virus (EBV) is a paradigm for human tumor viruses: it is the first virus recognized to cause cancer in people; it causes both lymphomas and carcinomas; yet these tumors arise infrequently given that most people in the world are infected with the virus. EBV is maintained extrachromosomally in infected normal and tumor cells. Eighty-four percent of these viral plasmids replicate each S phase, are licensed, require a single viral protein for their synthesis, and can use two functionally distinct origins of DNA replication, oriP, and Raji ori. Eighty-eight percent of newly synthesized plasmids are segregated faithfully to the daughter cells. Infectious viral particles are not synthesized under these conditions of latent infection. This plasmid replication is consistent with survival of EBV's host cells. Rare cells in an infected population either spontaneously or following exogenous induction support EBV's lytic cycle, which is lethal for the cell. In this case, the viral DNA replicates 100-fold or more, uses a third kind of viral origin of DNA replication, oriLyt, and many viral proteins. Here we shall describe the three modes of EBV's replication as a function of the viral origins used and the viral and cellular proteins that mediate the DNA synthesis from these origins focusing, where practical, on recent advances in our understanding.

A structural basis for the assembly and functions of a viral polymer that inactivates multiple tumor suppressors

Evolution of minimal DNA tumor virus' genomes has selected for small viral oncoproteins that hijack critical cellular protein interaction networks. The structural basis for the multiple and dominant functions of adenovirus oncoproteins has remained elusive. E4-ORF3 forms a nuclear polymer and simultaneously inactivates p53, PML, TRIM24, and MRE11/RAD50/NBS1 (MRN) tumor suppressors. We identify oligomerization mutants and solve the crystal structure of E4-ORF3. E4-ORF3 forms a dimer with a central β core, and its structure is unrelated to known polymers or oncogenes. E4-ORF3 dimer units coassemble through reciprocal and nonreciprocal exchanges of their C-terminal tails. This results in linear and branched oligomer chains that further assemble in variable arrangements to form a polymer network that partitions the nuclear volume. E4-ORF3 assembly creates avidity-driven interactions with PML and an emergent MRN binding interface. This reveals an elegant structural solution whereby a small protein forms a multivalent matrix that traps disparate tumor suppressors.

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.

Identification and analysis of papillomavirus E2 protein binding sites in the human genome

Papillomavirus E2 protein is required for the replication and maintenance of viral genomes and transcriptional regulation of viral genes. E2 functions through sequence-specific binding to 12-bp DNA motifs-E2 binding sites (E2BS)-in the virus genome. Papillomaviruses are able to establish persistent infection in their host and have developed a long-term relationship with the host cell in order to guarantee the propagation of the virus. In this study, we have analyzed the occurrence and functionality of E2BSs in the human genome. Our computational analysis indicates that most E2BSs in the human genome are found in repetitive DNA regions and have G/C-rich spacer sequences. Using a chromatin immunoprecipitation approach, we show that human papillomavirus type 11 (HPV11) E2 interacts with a subset of cellular E2BSs located in active chromatin regions. Two E2 activities, sequence-specific DNA binding and interaction with cellular Brd4 protein, are important for E2 binding to consensus sites. E2 binding to cellular E2BSs has a moderate or no effect on cellular transcription. We suggest that the preference of HPV E2 proteins for E2BSs with A/T-rich spacers, which are present in the viral genomes and underrepresented in the human genome, ensures E2 binding to specific binding sites in the virus genome and may help to prevent extensive and possibly detrimental changes in cellular transcription in response to the viral protein.

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.

Design and characterization of an enhanced repressor of human papillomavirus E2 protein

Papillomaviruses are causative agents of cervical and anogenital cancers. The viral E2 protein mediates viral DNA replication and transactivation of viral oncogenes and thus represents a specific target for therapeutic intervention. Short forms of E2, E2R, contain only the C-terminal dimerization domain, and repress the normal function of E2 due to formation of an inactive heterodimer. Using structure-guided design, we replaced conserved residues at the dimer interface to design a heterodimer with increased stability. One E2R mutant in which histidine was replaced by a glutamate residue showed preferential heterodimer formation in vitro, as well as an increase in plasticity at the interface, as a result of histidine-glutamate pair formation, as observed spectroscopically and in the crystal structure, determined to 2.2-Å resolution. In addition, the enhanced E2R showed greater repression of transcription from E2-responsive reporter plasmids in mammalian cell culture. Recent advances in protein delivery into the cell raise the possibility of using exogenously added proteins as therapeutic agents. More generally, this approach may be used to target the subunit interfaces of any multisubunit protein having a similar mechanism of action.

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.

Pushing the limits of what is achievable in protein-DNA docking: benchmarking HADDOCK's performance

The intrinsic flexibility of DNA and the difficulty of identifying its interaction surface have long been challenges that prevented the development of efficient protein-DNA docking methods. We have demonstrated the ability our flexible data-driven docking method HADDOCK to deal with these before, by using custom-built DNA structural models. Here we put our method to the test on a set of 47 complexes from the protein-DNA docking benchmark. We show that HADDOCK is able to predict many of the specific DNA conformational changes required to assemble the interface(s). Our DNA analysis and modelling procedure captures the bend and twist motions occurring upon complex formation and uses these to generate custom-built DNA structural models, more closely resembling the bound form, for use in a second docking round. We achieve throughout the benchmark an overall success rate of 94% of one-star solutions or higher (interface root mean square deviation ≤4 A and fraction of native contacts >10%) according to CAPRI criteria. Our improved protocol successfully predicts even the challenging protein-DNA complexes in the benchmark. Finally, our method is the first to readily dock multiple molecules (N > 2) simultaneously, pushing the limits of what is currently achievable in the field of protein-DNA docking.

Drugging challenging targets using fragment-based approaches

Fragment-based approaches have now become firmly established in the drug discovery armoury. After notable early successes against protein kinases, the versatility and power of fragment-based approaches are increasingly being demonstrated on more diverse and difficult protein targets. This review highlights seven examples including targeting protein-protein interactions, a RNA polymerase and a DNA-binding protein. It shows how fragment-based approaches using small libraries have been successful when large HTS screens have failed. It also highlights the range of biophysical approaches being used and the interplay between experimental and in silico screens. The examples all show the iterative way in which potency is built up by synthetic elaboration of the initial fragment hits.

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.

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.

Evolutionary and biophysical relationships among the papillomavirus E2 proteins

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

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

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

Transition state for protein-DNA recognition

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

DBD-Hunter: a knowledge-based method for the prediction of DNA-protein interactions

The structures of DNA–protein complexes have illuminated the diversity of DNA–protein binding mechanisms shown by different protein families. This lack of generality could pose a great challenge for predicting DNA–protein interactions. To address this issue, we have developed a knowledge-based method, DNA-binding Domain Hunter (DBD-Hunter), for identifying DNA-binding proteins and associated binding sites. The method combines structural comparison and the evaluation of a statistical potential, which we derive to describe interactions between DNA base pairs and protein residues. We demonstrate that DBD-Hunter is an accurate method for predicting DNA-binding function of proteins, and that DNA-binding protein residues can be reliably inferred from the corresponding templates if identified. In benchmark tests on ∼4000 proteins, our method achieved an accuracy of 98% and a precision of 84%, which significantly outperforms three previous methods. We further validate the method on DNA-binding protein structures determined in DNA-free (apo) state. We show that the accuracy of our method is only slightly affected on apo-structures compared to the performance on holo-structures cocrystallized with DNA. Finally, we apply the method to ∼1700 structural genomics targets and predict that 37 targets with previously unknown function are likely to be DNA-binding proteins. DBD-Hunter is freely available at http://cssb.biology.gatech.edu/skolnick/webservice/DBD-Hunter/.

Magnitude and direction of DNA bending induced by screw-axis orientation: influence of sequence, mismatches and abasic sites

DNA-bending flexibility is central for its many biological functions. A new bending restraining method for use in molecular mechanics calculations and molecular dynamics simulations was developed. It is based on an average screw rotation axis definition for DNA segments and allows inducing continuous and smooth bending deformations of a DNA oligonucleotide. In addition to controlling the magnitude of induced bending it is also possible to control the bending direction so that the calculation of a complete (2-dimensional) directional DNA-bending map is now possible. The method was applied to several DNA oligonucleotides including A(adenine)-tract containing sequences known to form stable bent structures and to DNA containing mismatches or an abasic site. In case of G:A and C:C mismatches a greater variety of conformations bent in various directions compared to regular B-DNA was found. For comparison, a molecular dynamics implementation of the approach was also applied to calculate the free energy change associated with bending of A-tract containing DNA, including deformations significantly beyond the optimal curvature. Good agreement with available experimental data was obtained offering an atomic level explanation for stable bending of A-tract containing DNA molecules. The DNA-bending persistence length estimated from the explicit solvent simulations is also in good agreement with experiment whereas the adiabatic mapping calculations with a GB solvent model predict a bending rigidity roughly two times larger.

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.

The plasmid replicon of Epstein-Barr virus: mechanistic insights into efficient, licensed, extrachromosomal replication in human cells

The genome of Epstein-Barr Virus (EBV) and plasmid derivatives of it are among the most efficient extrachromosomal replicons in mammalian cells. The latent origin of plasmid replication (oriP), when supplied with the viral Epstein-Barr Nuclear Antigen 1 (EBNA1) in trans, provides efficient duplication, partitioning and maintenance of plasmids bearing it. In this review, we detail what is known about the viral cis and trans elements required for plasmid replication. In addition, we describe how the cellular factors that EBV usurps are used to complement the functions of the viral constituents. Finally, we propose a model for the sequential assembly of an EBNA1-dependent origin of DNA synthesis into a pre-Replicative Complex (pre-RC), which functions by making use only of cellular enzymatic activities to carry out the replication of the viral plasmid.

Characterization of the nuclear localization signal of high risk HPV16 E2 protein

The E2 protein of high risk human papillomavirus type 16 (HPV16) contains an amino-terminal (N) domain, a hinge (H) region and a carboxyl-terminal (C) DNA-binding domain. Using enhanced green fluorescent protein (EGFP) fusions with full length E2 and E2 domains in transfection assays in HeLa cells, we found that the C domain is responsible for the nuclear localization of E2 in vivo, whereas the N and H domains do not contain additional nuclear localization signals (NLSs). Deletion analysis of EGFP-E2 and EGFP-cE2 determined that the C domain contains an alpha helix cNLS that overlaps with the DNA-binding region. Mutational analysis revealed that the arginine and lysine residues in this cNLS are essential for nuclear localization of HPV16 E2. Interestingly, these basic amino acid residues are well conserved among the E2 proteins of BPV-1 and some high risk HPV types but not in the low risk HPV types, suggesting that there are differences between the NLSs and corresponding nuclear import pathways between these E2 proteins.

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.

Insights into protein-protein interfaces using a Bayesian network prediction method

Identifying the interface between two interacting proteins provides important clues to the function of a protein, and is becoming increasing relevant to drug discovery. Here, surface patch analysis was combined with a Bayesian network to predict protein-protein binding sites with a success rate of 82% on a benchmark dataset of 180 proteins, improving by 6% on previous work and well above the 36% that would be achieved by a random method. A comparable success rate was achieved even when evolutionary information was missing, a further improvement on our previous method which was unable to handle incomplete data automatically. In a case study of the Mog1p family, we showed that our Bayesian network method can aid the prediction of previously uncharacterised binding sites and provide important clues to protein function. On Mog1p itself a putative binding site involved in the SLN1-SKN7 signal transduction pathway was detected, as was a Ran binding site, previously characterized solely by conservation studies, even though our automated method operated without using homologous proteins. On the remaining members of the family (two structural genomics targets, and a protein involved in the photosystem II complex in higher plants) we identified novel binding sites with little correspondence to those on Mog1p. These results suggest that members of the Mog1p family bind to different proteins and probably have different functions despite sharing the same overall fold. We also demonstrated the applicability of our method to drug discovery efforts by successfully locating a number of binding sites involved in the protein-protein interaction network of papilloma virus infection. In a separate study, we attempted to distinguish between the two types of binding site, obligate and non-obligate, within our dataset using a second Bayesian network. This proved difficult although some separation was achieved on the basis of patch size, electrostatic potential and conservation. Such was the similarity between the two interacting patch types, we were able to use obligate binding site properties to predict the location of non-obligate binding sites and vice versa.

Human papillomaviruses: basic mechanisms of pathogenesis and oncogenicity

Human papillomaviruses (HPVs) are small double-stranded DNA viruses that infect the cutaneous and mucosal epithelium. Infection by specific HPV types has been linked to the development of cervical carcinoma. HPV infects epithelial cells that undergo terminal differentiation and so encode multiple mechanisms to override the normal regulation of differentiation to produce progeny virions. Two viral proteins, E6 and E7, alter cell cycle control and are the main arbitrators of HPV-induced oncogenesis. Recent data suggest that E6 and E7 also play a major role in the inhibition of the host cell innate immune response to HPV. The E1 and E2 proteins, in combination with various cellular factors, mediate viral replication. In addition, E2 has been implicated in both viral and cellular transcriptional control. Despite decades of research, the function of other viral proteins still remains unclear. While prophylactic vaccines to block genital HPV infection will soon be available, the widespread nature of HPV infection requires greater understanding of both the HPV life cycle as well as the mechanisms underlying HPV-induced carcinogenesis.

Specificity in DNA recognition by a peptide from papillomavirus E2 protein

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

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

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

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

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

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

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