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

Slow motions in A·T rich DNA sequence

In free B-DNA, slow (microsecond-to-millisecond) motions that involve equilibrium between Watson-Crick (WC) and Hoogsteen (HG) base-pairing expand the DNA dynamic repertoire that could mediate DNA-protein assemblies. R_1ρ relaxation dispersion NMR methods are powerful tools to capture such slow conformational exchanges in solution using ¹³C/¹⁵ N labelled DNA. Here, these approaches were applied to a dodecamer containing a TTAAA element that was assumed to facilitate nucleosome formation. NMR data and inferred exchange parameters assign HG base pairs as the minor, transient conformers specifically observed in three successive A·T base pairs forming the TAA·TTA segment. The abundance of these HG A·T base pairs can be up to 1.2% which is high compared to what has previously been observed. Data analyses support a scenario in which the three adenines undergo non-simultaneous motions despite their spatial proximity, thus optimising the probability of having one HG base pair in the TAA·TTA segment. Finally, revisiting previous NMR data on H2 resonance linewidths on the basis of our results promotes the idea of there being a special propensity of A·T base pairs in TAA·TTA tracts to adopt HG pairing. In summary, this study provides an example of a DNA functional element submitted to slow conformational exchange. More generally, it strengthens the importance of the role of the DNA sequence in modulating its dynamics, over a nano- to milli-second time scale.

Innate Molecular and Cellular Signature in the Skin Preceding Long-Lasting T Cell Responses after Electroporated DNA Vaccination

DNA vaccines delivered with electroporation (EP) have shown promising results in preclinical models and are evaluated in clinical trials. In this study, we aim to characterize early mechanisms occurring in the skin after intradermal injection and EP of the auxoGTUmultiSIV DNA vaccine in nonhuman primates. First, we show that EP acts as an adjuvant by enhancing local inflammation, notably via granulocytes, monocytes/macrophages, and CD1a^int-expressing cell recruitment. EP also induced Langerhans cell maturation, illustrated by CD86, CD83, and HLA-DR upregulation and their migration out of the epidermis. Second, we demonstrate the crucial role of the DNA vaccine in soluble factors release, such as MCP-1 or IL-15. Transcriptomic analysis showed that EP played a major role in gene expression changes postvaccination. However, the DNA vaccine is required to strongly upregulate several genes involved in inflammatory responses (e.g., Saa4), cell migration (e.g., Ccl3, Ccl5, or Cxcl10), APC activation (e.g., Cd86), and IFN-inducible genes (e.g., Ifit3, Ifit5, Irf7, Isg15, orMx1), illustrating an antiviral response signature. Also, AIM-2, a cytosolic DNA sensor, appeared to be strongly upregulated only in the presence of the DNA vaccine and trends to positively correlate with several IFN-inducible genes, suggesting the potential role of AIM-2 in vaccine sensing and the subsequent innate response activation leading to strong adaptive T cell responses. Overall, these results demonstrate that a combined stimulation of the immune response, in which EP and the auxoGTUmultiSIV vaccine triggered different components of the innate immunity, led to strong and persistent cellular recall responses.

A dynamic view of DNA structure within the nucleosome: Biological implications

Most of eukaryotic cellular DNA is packed in nucleosome core particles (NCPs), in which the DNA (DNA_NCP) is wrapped around histones. The influence of this organization on the intrinsic local dynamics of DNA is largely unknown, in particular because capturing such information from experiments remains notoriously challenging. Given the importance of dynamical properties in DNA functions, we addressed this issue using CHARMM36 MD simulations of a nucleosome containing the NCP positioning 601 sequence and four related free dodecamers. Comparison between DNA_NCP and free DNA reveals a limited impact of the dense DNA-histone interface on correlated motions of dinucleotide constituents and on fluctuations of inter base pair parameters. A characteristic feature intimately associated with the DNA_NCP super-helical path is a set of structural periodicities that includes a marked alternation of regions enriched in backbone BI and BII conformers. This observation led to uncover a convincing correspondence between the sequence effect on BI/BII propensities in both DNA_NCP and free DNA, strengthening the idea that the histone preference for particular DNA sequences relies on those intrinsic structural properties. These results offer for the first time a detailed view of the DNA dynamical behavior within NCP. They show in particular that the DNA_NCP dynamics is substantial enough to preserve the ability to structurally adjust to external proteins, for instance remodelers. Also, fresh structural arguments highlight the relevance of relationships between DNA sequence and structural properties for NCP formation. Overall, our work offers a more rational framework to approach the functional, biological roles of NCP.

Simulations Meet Experiment to Reveal New Insights into DNA Intrinsic Mechanics

The accurate prediction of the structure and dynamics of DNA remains a major challenge in computational biology due to the dearth of precise experimental information on DNA free in solution and limitations in the DNA force-fields underpinning the simulations. A new generation of force-fields has been developed to better represent the sequence-dependent B-DNA intrinsic mechanics, in particular with respect to the BI ↔ BII backbone equilibrium, which is essential to understand the B-DNA properties. Here, the performance of MD simulations with the newly updated force-fields Parmbsc0εζOLI and CHARMM36 was tested against a large ensemble of recent NMR data collected on four DNA dodecamers involved in nucleosome positioning. We find impressive progress towards a coherent, realistic representation of B-DNA in solution, despite residual shortcomings. This improved representation allows new and deeper interpretation of the experimental observables, including regarding the behavior of facing phosphate groups in complementary dinucleotides, and their modulation by the sequence. It also provides the opportunity to extensively revisit and refine the coupling between backbone states and inter base pair parameters, which emerges as a common theme across all the complementary dinucleotides. In sum, the global agreement between simulations and experiment reveals new aspects of intrinsic DNA mechanics, a key component of DNA-protein recognition.

Automatic workflow for the classification of local DNA conformations

Background

A growing number of crystal and NMR structures reveals a considerable structural polymorphism of DNA architecture going well beyond the usual image of a double helical molecule. DNA is highly variable with dinucleotide steps exhibiting a substantial flexibility in a sequence-dependent manner. An analysis of the conformational space of the DNA backbone and the enhancement of our understanding of the conformational dependencies in DNA are therefore important for full comprehension of DNA structural polymorphism.

Results

A detailed classification of local DNA conformations based on the technique of Fourier averaging was published in our previous work. However, this procedure requires a considerable amount of manual work. To overcome this limitation we developed an automatic classification method consisting of the combination of supervised and unsupervised approaches. A proposed workflow is composed of k-NN method followed by a non-hierarchical single-pass clustering algorithm. We applied this workflow to analyze 816 X-ray and 664 NMR DNA structures released till February 2013. We identified and annotated six new conformers, and we assigned four of these conformers to two structurally important DNA families: guanine quadruplexes and Holliday (four-way) junctions. We also compared populations of the assigned conformers in the dataset of X-ray and NMR structures.

Conclusions

In the present work we developed a machine learning workflow for the automatic classification of dinucleotide conformations. Dinucleotides with unassigned conformations can be either classified into one of already known 24 classes or they can be flagged as unclassifiable. The proposed machine learning workflow permits identification of new classes among so far unclassifiable data, and we identified and annotated six new conformations in the X-ray structures released since our previous analysis. The results illustrate the utility of machine learning approaches in the classification of local DNA conformations.

Persistent immune responses induced by a human immunodeficiency virus DNA vaccine delivered in association with electroporation in the skin of nonhuman primates

Strategies to improve vaccine efficacy are still required, especially in the case of chronic infections, including human immunodeficiency virus (HIV). DNA vaccines have potential advantages over conventional vaccines; however, low immunological efficacy has been demonstrated in many experiments involving large animals and in clinical trials. To improve the immunogenicity of DNA vaccines, we have designed a plasmid vector exploiting the binding capacity of the bovine papillomavirus E2 protein and we have used electroporation (EP) to increase DNA uptake after intradermal inoculation. We demonstrated, in nonhuman primates (NHPs), efficient induction of anti-HIV immunity with an improved DNA vaccine vector encoding an artificial fusion protein, consisting of several proteins and selected epitopes from HIV-1. We show that a DNA vaccine delivery method combining intradermal injection and noninvasive EP dramatically increased expression of the vaccine antigen selectively in the epidermis, and our observations strongly suggest the involvement of Langerhans cells in the strength and quality of the anti-HIV immune response. Although the humoral responses to the vaccine were transient, the cellular responses were exceptionally robust and persisted, at high levels, more than 2 years after the last vaccine boost. The immune responses were characterized by the induction of significant proportions of T cells producing both interferon-gamma and interleukin-2 cytokines, in both subpopulations, CD4(+) and CD8(+). This strategy is an attractive approach for vaccination in humans because of its high efficacy and the possible use of newly developed devices for EP.

Investigation of transcription factor Ndt80 affinity differences for wild type and mutant DNA: a molecular dynamics study

Molecular dynamics simulations and free energy calculations have been performed on the transcription factor Ndt80 either in complex with the native DNA sequence or with a mutant DNA with a switched central base pair, C5-G5' to G5-C5'. This mutant has been shown to have a 100-fold decrease in binding affinity of Ndt80, and in this study we explain this both structurally and energetically. The major interactions between the protein and the DNA were maintained in the simulations, apart from around the mutation site. The crystal structure of the Ndt80-DNA complex revealed that R177 makes a base specific bidentate interaction with G5' which is part of a conserved 5'-YpG-3' step. In the simulation with the mutant DNA, the side chain of R177 changes conformation and makes three new stable hydrogen bonds to the DNA backbone. This in turn induces a conformational change in the DNA backbone of the T6'-G5' step from the unusual BII state to the canonical BI state. The affinity difference for the protein-DNA complex with the native DNA compared with the mutant DNA is only about 3 kcal/mol. The free energy calculations of the base pair switch indicated a larger difference than what was found experimentally, about 7.7 kcal/mol, but this is explained in structural terms using the 10 ns simulations of the solvated complexes and the rearrangement of the R177 side chain.

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.

DNA distortion and specificity in a sequence-specific endonuclease

Five new structures of the Q138F HincII enzyme bound to a total of three different DNA sequences and three different metal ions (Ca(2+), Mg(2+), and Mn(2+)) are presented. While previous structures were produced from soaking Ca(2+) into preformed Q138F HincII/DNA crystals, the new structures are derived from cocrystallization with Ca(2+), Mg(2+), or Mn(2+). The Mn(2)(+)-bound structure provides the first view of a product complex of Q138F HincII with cleaved DNA. Binding studies and a crystal structure show how Ca(2+) allows trapping of a Q138F HincII complex with noncognate DNA in a catalytically incompetent conformation. Many Q138F HincII/DNA structures show asymmetry, despite the binding of a symmetric substrate by a symmetric enzyme. The various complexes are fit into a model describing the different conformations of the DNA-bound enzyme and show how DNA conformational energetics determine DNA-cleavage rates by the Q138F HincII enzyme.

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.

Quantification of DNA BI/BII backbone states in solution. Implications for DNA overall structure and recognition

The backbone states of B-DNA influence its helical parameters, groove dimensions, and overall curvature. Therefore, detection and fine characterization of these conformational states are desirable. Using routine NMR experiments on a nonlabeled B-DNA oligomer and analyzing high-resolution X-ray structures, we investigated the relationship between interproton distances and backbone conformational states. The three H2'i-H6/8i+1, H2' 'i-H6/8i+1, and H6/8i-H6/8i+1 sequential distances were found cross-correlated and linearly coupled to epsilon-zeta values in X-ray structures and 31P chemical shifts (deltaP) in NMR that reflect the interconversion between the backbone BI (epsilon-zeta < 0 degrees ) and BII (epsilon-zeta > 0 degrees) states. These relationships provide a detailed check of the NMR data consistency and the possibility to extend the set of restraints for structural refinement through various extrapolations. Furthermore, they allow translation of deltaP in terms of BI/BII ratios. Also, comparison of many published deltaP in solution to crystal data shows that the impact of sequence on the BI/BII propensities is similar in both environments and is therefore an intrinsic and general property of B-DNA. This quantification of the populations of BI and BII is of general interest because these sequence-dependent backbone states act on DNA overall structure, a key feature for DNA-protein-specific recognition.

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

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

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

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

Conformations and Dynamics of the Phosphodiester Backbone of a DNA Fragment That Bears a Strong Topoisomerase II Cleavage Site

The dynamics of the DNA phosphodiester backbone conformations have been studied for a strong topoisomerase II cleavage site (site 22) using molecular dynamics simulations in explicit water and in the presence of sodium ions. We investigated the backbone motions and more particularly the BI/BII transitions involving the epsilon and zeta angles. The consensus cleavage site is adjacent to the phosphate which shows the most important phosphodiester backbone flexibility in the sequence. We infer that these latter properties could be responsible for the preferential cleavage at this site possibly through the perturbation of the cleavage/ligation activities of the topoisomerase II. More generally, the steps pur-pur and pyr-pur are those presenting the highest BII contents. Relations are observed between the backbone phosphodiester BI/BII transitions and the flexibility of the deoxyribose sugar and the helical parameters such as roll. The roll is sequence dependent when the related phosphate is in the BI form, whereas this appears not to be true when it is in the BII form. The BI/BII transitions are associated with water migration, and new relations are observed with counterions. Indeed, it is observed that a strong coupling exists between the BII form and the presence of sodium ions near the adjacent sugar deoxyribose. The presence of sodium ions in the O4' surroundings or their binding could assist the BI to BII transition by furnishing energy. The implications of these new findings and, namely, their importance in the context of the sequence-dependent behavior of BI/BII transitions will be investigated in future studies.

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.

Bending in the right direction

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

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

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

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.

Sequence preference for BI/BII conformations in DNA: MD and crystal structure data analysis

Deciphering sequence information from sugar-phosphate backbone is finely tuned through the conformational substates of DNA. BII conformation, one of the conformational substates of B-DNA, is known to play a key role in DNA-protein recognition. BI and BII are identified by the epsilon-zeta difference, which is negative in BI and positive in BII. Our analysis of MD and crystal structures shows that BII conformation is sequence specific and dinucleotides GC, CG, CA, TG, TA show high preference to take up BII conformation, while TT, TC, CT, CC dinucleotides rarely take up this conformation. Significant changes were observed in the dinucleotide parameters viz. twist, roll, and slide for the steps having BII conformation. Interestingly, the magnitude of variation in the dinucleotide parameters is seen to depend mainly on two factors, the magnitude of epsilon-zeta difference and the presence or absence of BII conformation in the second strand, across the WC base-paired dinucleotide step. Based on these two factors, the conformational substate of a dinucleotide step can be further classified as BI.BI (BI conformation in both strands), BI.BII (BI conformation in one strand and BII conformation in the other), and BII.BII (BII conformation in both strands). The occurrence of BII in both strands was found to be quite rare and thus, it can be concluded that BI.BI and BI.BII hybrid steps are more favorable than a BII.BII step. In conformity with the sequence preference seen for dinucleotides in each strand, BII.BII combination of backbone conformation was observed only for GC, CG, CA, and TG containing dinucleotide steps. We further classified BII.BII step as strong BII and weak BII depending on the magnitude of the average epsilon-zeta difference. The dinucleotide steps which belong to the category of strong BII, have large twist, high positive slide and negative roll values, while those in the weak BII group have roll, twist, and slide values similar to that of hybrid BI.BII steps. This conformational property could be contributing to the groove opening/closing and thus can modulate protein-DNA interaction.