Sequence-dependent Structure Research Articles

Revisiting DNA Sequence-Dependent Deformability in High-Resolution Structures: Effects of Flanking Base Pairs on Dinucleotide Morphology and Global Chain Configuration.

DNA carries more than the list of biochemical ingredients that drive the basic functions of living systems. The sequence of base pairs includes a multitude of structural and energetic signals, which determine the degree to which the long, threadlike molecule moves and how it responds to proteins and other molecules that control its processing and govern its packaging. The chemical composition of base pairs directs the spatial disposition and fluctuations of successive residues. The observed arrangements of these moieties in high-resolution protein–DNA crystal structures provide one of the best available estimates of the natural, sequence-dependent structure and deformability of the double-helical molecule. Here, we update the set of knowledge-based elastic potentials designed to describe the observed equilibrium structures and configurational fluctuations of the ten unique base-pair steps. The large number of currently available structures makes it possible to characterize the configurational preferences of the DNA base-pair steps within the context of their immediate neighbors, i.e., tetrameric context. Use of these knowledge-based potentials shows promise in accounting for known effects of sequence in long chain molecules, e.g., the degree of curvature reported in classic gel mobility studies and the recently reported sequence-dependent responses of supercoiled minicircles to nuclease cleavage.

Sequence Blockiness Controls the Structure of Polyampholyte Necklaces.

A scaling theory of statistical (Markov) polyampholytes is developed to understand how sequence correlations, that is, the blockiness of positive and negative charges, influences conformational behavior. An increase in the charge patchiness leads to stronger correlation attractions between oppositely charged monomers, but simultaneously, it creates a higher charge imbalance in the polyampholyte. A competition between effective short-range attractions and long-range Coulomb repulsions induces globular, pearl-necklace, or fully stretched chain conformations, depending on the average length of the block of like charges. The necklace structure and the underlying distribution of the net charge are also controlled by the sequence. Sufficiently long blocks allow for charge migration from globular beads (pearls) to strings, thereby providing a nonmonotonic change in the number of necklace beads as the blockiness increases. The sequence-dependent structure of polyampholyte necklaces is confirmed by molecular dynamics simulations. The findings presented here provide a framework for understanding the sequence-encoded conformations of synthetic polyampholytes and intrinsically disordered proteins (IDPs).

Insights into Structure and Aggregation Behavior of Elastin-like Polypeptide Coacervates: All-Atom Molecular Dynamics Simulations.

The stimuli-responsive character of elastin-like polypeptides (ELP) has led to their use in a wide range of applications. The temperature-triggered aggregation, or LCST behavior, of ELPs is a complex and multistep phenomenon, which proposed to include the structural transitions, loss of hydrophobic hydration, expulsion of water molecules and physical association of chains. Thus, the origin and detailed mechanism of LCST in ELPs is difficult to elucidate. Here, to gain insights into structure and dynamics of coacervates, we performed all-atom molecular dynamics simulations of 27 90-mer ELPs in explicit water at 350 K. Two sequences, poly(VGPVG)18 and poly(VPGVG)18, were examined due to their experimentally observed differences in thermal hysteresis albeit identical overall composition but different arrangement of amino acids. The simulation results indicate that surface hydrophobicity of poly(VGPVG) aggregate is less than that of the poly(VPGVG) aggregate, and there are marked changes in torsion angles and the propensities of secondary structural motifs during the aggregation process. Moreover, there are significant differences between structure of a single polypeptide in water and structure within the aggregate. Overall, the aggregation process is driven by the formation of peptide-peptide interactions whereas the average hydration of peptides remains almost the same between dissolved and aggregated states. Even though the aggregation is driven by the hydrophobic interactions, ELP coacervate has no hydrophobic core and contains many water molecules. Overall, our findings provide an insight into the sequence-dependent structure of coacervates and molecular behavior of individual peptides during aggregation.

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.

Efficient BiSRU Combined With Feature Dimensionality Reduction for Abnormal Traffic Detection

Abnormal traffic detection is an important network security technology to protect computer systems from malicious attacks. Existing detection methods are usually based on traditional machine learning, such as Support Vector Machine (SVM), Naive Bayes, etc. They rely heavily on manual design of traffic features and usually shallow feature learning, which get a low accuracy for high-dimensional traffic. Although the method based on Long Short-Term Memory (LSTM) has an excellent ability to detect abnormal traffic. The sequence-dependent structure of LSTM cannot realize parallel computation, which leads to slow model training and limits its applicability. To address the above problem, we propose an efficient Bidirectional Simple Recurrent Unit (BiSRU) combined with feature dimensionality reduction for abnormal traffic detection. Specifically, in order to perform feature dimensionality reduction on the original high-dimensional network traffic, we design a stack Sparse Autoencoder (sSAE) to extract the compressed high-level features. For the purpose of realizing efficient parallel computation and accurate feature extraction, a BiSRU is utilized to extract the bidirectional structural features of the traffic. Finally, the experimental results show that our proposed method significantly outperforms existing methods in terms of accuracy and training time. The method we propose can timely and accurately detect various abnormal traffic and achieve effective network security protection.

How methyl-sugar interactions determine DNA structure and flexibility.

The sequence dependent structure and flexibility of the DNA double helix is of key importance for gene expression and DNA packing and it can be modulated by DNA modifications. The presence of a C5′-methyl group in thymine or the frequent C5′-methylated-cytosine affects the DNA fine structure, however, the underlying mechanism and steric origins have remained largely unexplained. Employing Molecular Dynamics free energy simulations that allow switching on or off interactions with the methyl groups in several DNA sequences, we systematically identified the physical origin of the coupling between methyl groups and DNA backbone fine structure. Whereas methyl-solvent and methyl–nucleobase interactions were found to be of minor importance, the methyl group interaction with the 5′ neighboring sugar was identified as main cause for influencing the population of backbone substates. The sterical methyl sugar clash prevents the formation of unconventional stabilizing hydrogen bonds between nucleobase and backbone. The technique was also used to study the contribution of methyl groups to DNA flexibility and served to explain why the presence of methyl sugar clashes in thymine and methyl-cytosine can result in an overall local increase of DNA flexibility.

X-ray structure of the MMTV-A nucleosome core

The conformation of DNA bound in nucleosomes depends on the DNA sequence. Questions such as how nucleosomes are positioned and how they potentially bind sequence-dependent nuclear factors require near-atomic resolution structures of the nucleosome core containing different DNA sequences; despite this, only the DNA for two similar α-satellite sequences and a sequence (601) selected in vitro have been visualized bound in the nucleosome core. Here we report the 2.6-Å resolution X-ray structure of a nucleosome core particle containing the DNA sequence of nucleosome A of the 3'-LTR of the mouse mammary tumor virus (147 bp MMTV-A). To our knowledge, this is the first nucleosome core particle structure containing a promoter sequence and crystallized from Mg(2+) ions. It reveals sequence-dependent DNA conformations not seen previously, including kinking into the DNA major groove.

Can nucleosomal DNA be described by an elastic model?: Comment on “Sequence-dependent collective properties of DNAs and their role in biological systems” by Pasquale De Santis and Anita Scipioni

Can nucleosomal DNA be described by an elastic model?: Comment on “Sequence-dependent collective properties of DNAs and their role in biological systems” by Pasquale De Santis and Anita Scipioni

A flexible integrative approach based on random forest improves prediction of transcription factor binding sites

Transcription factor binding sites (TFBSs) are DNA sequences of 6–15 base pairs. Interaction of these TFBSs with transcription factors (TFs) is largely responsible for most spatiotemporal gene expression patterns. Here, we evaluate to what extent sequence-based prediction of TFBSs can be improved by taking into account the positional dependencies of nucleotides (NPDs) and the nucleotide sequence-dependent structure of DNA. We make use of the random forest algorithm to flexibly exploit both types of information. Results in this study show that both the structural method and the NPD method can be valuable for the prediction of TFBSs. Moreover, their predictive values seem to be complementary, even to the widely used position weight matrix (PWM) method. This led us to combine all three methods. Results obtained for five eukaryotic TFs with different DNA-binding domains show that our method improves classification accuracy for all five eukaryotic TFs compared with other approaches. Additionally, we contrast the results of seven smaller prokaryotic sets with high-quality data and show that with the use of high-quality data we can significantly improve prediction performance. Models developed in this study can be of great use for gaining insight into the mechanisms of TF binding.

Click to Join Peptides/Proteins Together

Copper(I) is able to catalyze Huisgen 1,3-dipolar cycloaddition in a "click" fashion. This copper-catalyzed azide-alkyne cycloaddition (CuAAC) reaction presents excellent chemoselectivity and occurs over a wide-range of reaction conditions. It shows tolerance to variation in both pH and solvent polarity, thereby facilitating the ligation of peptides and proteins to produce peptidomimetics and synthetic proteins. In addition, the only product formed is a 1,4-disubstituted-1,2,3-triazole moiety, in many aspects resembling the natural peptide bond, including hydrogen-bonding capability, planarity, distance between the 1 and 4 substituents, and conformational restriction of the peptide backbone; thus the triazole-backbone-modified peptide, in which a triazole replaces the amide bond, may be anticipated to present a secondary structure similar to that of its natural counterpart. This Focus Review describes the scope and applications of copper(I)-catalyzed alkyne–azide cycloaddition in synthetic peptide/protein chemistry.

DNA sequence-dependent variation in nucleosome structure, stability, and dynamics detected by a FRET-based analysisThis paper is one of a selection of papers published in this Special Issue, entitled 29th Annual International Asilomar Chromatin and Chromosomes Conference, and has undergone the Journal’s usual peer review process.

Förster resonance energy transfer (FRET) techniques provide powerful and sensitive methods for the study of conformational features in biomolecules. Here, we review FRET-based studies of nucleosomes, focusing particularly on our work comparing the widely used nucleosome standard, 5S rDNA, and 2 promoter-derived regulatory element-containing nucleosomes, mouse mammary tumor virus (MMTV)-B and GAL10. Using several FRET approaches, we detected significant DNA sequence-dependent structure, stability, and dynamics differences among the three. In particular, 5S nucleosomes and 5S H2A/H2B-depleted nucleosomal particles have enhanced stability and diminished DNA dynamics, compared with MMTV-B and GAL10 nucleosomes and particles. H2A/H2B-depleted nucleosomes are of interest because they are produced by the activities of many transcription-associated complexes. Significant location-dependent (intranucleosomal) stability and dynamics variations were also observed. These also vary among nucleosome types. Nucleosomes restrict regulatory factor access to DNA, thereby impeding genetic processes. Eukaryotic cells possess mechanisms to alter nucleosome structure, to generate DNA access, but alterations often must be targeted to specific nucleosomes on critical regulatory DNA elements. By endowing specific nucleosomes with intrinsically higher DNA accessibility and (or) enhanced facility for conformational transitions, DNA sequence-dependent nucleosome dynamics and stability variations have the potential to facilitate nucleosome recognition and, thus, aid in the crucial targeting process. This and other nucleosome structure and function conclusions from FRET analyses are discussed.

Nearest-neighbor non-additivity versus long-range non-additivity in TATA-box structure and its implications for TBP-binding mechanism

TBP recognizes its target sites, TATA boxes, by recognizing their sequence-dependent structure and flexibility. Studying this mode of TATA-box recognition, termed ‘indirect readout’, is important for elucidating the binding mechanism in this system, as well as for developing methods to locate new binding sites in genomic DNA. We determined the binding stability and TBP-induced TATA-box bending for consensus-like TATA boxes. In addition, we calculated the individual information score of all studied sequences. We show that various non-additive effects exist in TATA boxes, dependent on their structural properties. By several criterions, we divide TATA boxes to two main groups. The first group contains sequences with 3–4 consecutive adenines. Sequences in this group have a rigid context-independent cooperative structure, best described by a nearest-neighbor non-additive model. Sequences in the second group have a flexible, context-dependent conformation, which cannot be described by an additive model or by a nearest-neighbor non-additive model. Classifying TATA boxes by these and other structural rules clarifies the different recognition pathways and binding mechanisms used by TBP upon binding to different TATA boxes. We discuss the structural and evolutionary sources of the difficulties in predicting new binding sites by probabilistic weight-matrix methods for proteins in which indirect readout is dominant.

The stacked-X DNA Holliday junction and protein recognition.

The crystal structure of the four-stranded DNA Holliday junction has now been determined in the presence and absence of junction binding proteins, with the extended open-X form of the junction seen in all protein complexes, but the more compact stacked-X structure observed in free DNA. The structures of the stacked-X junction were crystallized because of an unexpected sequence dependence on the stability of this structure. Inverted repeat sequences that contain the general motif NCC or ANC favor formation of stacked-X junctions, with the junction cross-over occurring between the first two positions of the trinucleotides. This review focuses on the sequence dependent structure of the stacked-X junction and how it may play a role in structural recognition by a class of dimeric junction resolving enzymes that themselves show no direct sequence recognition.

Sequence-dependent structural effects in left-handed DNA

About 25 years ago the DNA sequence d(CGCGCG)2 was crystallized, and the structure solved as a left-handed helix, dubbed Z-DNA. Many crystal structures of variants of this first sequence have since been obtained. In our laboratory we have focused on what happens to the helical structure, inter-helical interactions, and the solvent interactions when a single A:T (adenine–thymine) base pair replaces a G:C (guanine–cytosine) base pair in the hexameric sequence. We present here a review of these studies. The review addresses three topics. In the first part, it discusses the effect of A:T base pairs on the structure of Z-DNA. The crystal structures demonstrate that the degree of structural perturbation that an A:T base pair produces depends on its position in the sequence, as this alters the length of the G:C tract. The presence of four continuous G:C base pairs appears to be necessary and sufficient to nucleate a Z-type structure. In the second part, the review deals with interactions between helices and crystal packing modes in Z-DNA. It examines the consequences to these aspects, of the presence of an A:T base pair in the sequence. Hexameric duplexes may be approximated as cylinders. This leads to degeneracy in crystal packing modes in Z-DNA, particularly at low resolution. The presence of a single A:T base pair breaks the two-fold symmetry in the sequence and gives rise to directionality in the cylinders. This has interesting effects on the inter-helical interactions, and therefore on the crystal packing. Metal ion and solvent interactions are discussed in the third section of the review. Different ions have different modes of binding to the helices. Some of these are specific to the sequence and lead not only to subtle, but well-marked, variations in the helical structure, but also to differences in helical contacts and therefore crystal packing. In some cases they induce disorder in crystal packing. In others, the ions interact with the helices in a non-specific manner. Contents 1. Introduction 338 2. Sequence-dependent structure 339 3. Crystal packing 344 4. Solvent and ion interaction 350 Acknowledgements 353 References 353 Subject index 355

Induced Fit and the Entropy of Structural Adaptation in the Complexation of CAP and λ-Repressor with Cognate DNA Sequences

Molecular dynamics (MD) simulations of 5 ns on protein-DNA complexes of catabolite-activator protein (CAP), λ-repressor, and their corresponding uncomplexed protein and DNA, are reported. These cases represent two extremes of DNA bending, with CAP DNA bent severely and the λ-operator nearly straight when complexed with protein. The calculations were performed using the AMBER suite of programs and the parm94 force field, validated for these studies by good agreement with experimental nuclear magnetic resonance data on DNA. An explicit computational model of structural adaptation and computation of the quasiharmonic entropy of association were obtained from the MD. The results indicate that, with respect to canonical B-form DNA, the extreme bending of the DNA in the complex with CAP is ∼60% protein-induced and 40% intrinsic to the sequence-dependent structure of the free oligomer. The DNA in the complex is an energetically strained form, and the MD results are consistent with a conformational-capture mechanism. The calculated quasiharmonic entropy change accounts for the entropy difference between the two cases. The calculated entropy was decomposed into contributions from protein adaptation, DNA adaptation, and protein-DNA structural correlations. The origin of the entropy difference between CAP and λ-repressor complexation arises more from the additional protein adaptation in the case of λ, than to DNA bending and entropy contribution from DNA bending. The entropy arising from protein DNA cross-correlations, a contribution not previously discussed, is surprisingly large.

Analysis of Chemical and Enzymatic Cleavage Frequencies in Supercoiled DNA

Chemical and enzymatic probing methods are powerful techniques for examining details of sequence-dependent structure in DNA and RNA. Reagents that cleave nucleic acid molecules in a structure-specific, but relatively sequence-non-specific manner, such as hydroxyl radical or DNase I, have been used widely to probe helical geometry in nucleic acid structures, nucleic acid–drug complexes, and in nucleoprotein assemblies. Application of cleavage-based techniques to structures present in superhelical DNA has been hindered by the fact that the cleavage pattern attributable to supercoiling-dependent structures is heavily mixed with non-specific cleavage signals that are inevitable products of multiple cleavage events. We present a rigorous mathematical procedure for extracting the cleavage pattern specific to supercoiled DNA and use this method to investigate the hydroxyl radical cleavage pattern in a cruciform DNA structure formed by a 60 bp inverted repeat sequence embedded in a negatively supercoiled plasmid. Our results support the presence of a stem-loop structure in the expected location and suggest that the helical geometry of the cruciform stem differs from that of the normal duplex form.

Sequence-dependent solution structure and motions of 13 TATA/TBP (TATA-box binding protein) complexes.

The TATA element is a well-known example of a DNA promoter sequence recognized by the TATA box binding protein (TBP) through its intrinsic motion and deformability. Although TBP recognizes the TATA element octamer unusually (through the minor groove, which lacks the distinctive features of the major groove), single base-pair replacements alter transcriptional activity. Recent crystallographic experiments have suggested that TATA/TBP complexes differing by a single base pair retain substantial structural similarity despite their functional differences in activating transcription. To investigate the subtle role of sequence-dependent motion within the TATA element and certain aspects of its effect on assembly of the transcriptional complex, we examine 5-ns dynamics trajectories of 13 variant TATA/TBP complexes differing from each other by a single base pair. They include the wild-type (WT) adenovirus 2 major late promoter (AdMLP) TATA element, TATAAAAG (the octamer specifies positions -31 to -24 with respect to the transcription initiation site), and the variants A31 (i.e., AATAAAAG), T30, A29, C29, G28, T28, T27, G26, T26, C25, T25, and T24. Our simulated TATA/TBP complexes develop sequence-dependent structure and motion trends that may lead to favorable orientations for high-activity variants (with respect to binding TFIIA, TFIIB, and other transcription factors), while conversely, accelerate dissociation of low-activity TATA/TBP complexes. The motions that promote favorable geometries for preinitiation complexes include small rotations between TBP's N- and C-terminal domains, sense strand DNA backbone "slithering," and rotations in TBP's H2 and H2' helices. Low-activity variants tend to translate the H1 and H1' helices and withdraw the intercalating phenylalanines. These cumulative DNA and protein motions lead to a spatial spread of complex orientations up to 4 A; this is associated with an overall bend of the variant TATA/TBP complexes that spans 93 degrees to 110 degrees (107 degrees for the crystal reference). Taken together, our analyses imply larger differences when these local structural and bending changes are extended to longer DNA (upstream and downstream) and suggest that specific local TATA/TBP motions (e.g., shifts in TBP helices and TATA bases and backbone) play a role in modulating the formation and maintenance of the transcription initiation complex.

RADACK, a Stochastic Simulation of Hydroxyl Radical Attack to DNA

RADACK was conceived to simulate the radiation-induced attack to different DNA forms and complexes. It allows to separately calculate the probability of attack to each reactive atom of the sugar and of the base and takes into account the sequence-dependent structure of DNA as known from crystallographic or NMR studies or resulting from molecular modelling. The calculations are aimed to assess sequence-, structure- and ligand-dependent modulation of damages of sugar and bases, leading to single strand breaks (frank strand breaks, FSB) and alkali-labile base modifications (alkali-revealed breaks, ARB), respectively. The modelling procedure and the results of simulations for some representative structures (B, Z and quadruplex forms) are here described and discussed. The calculated relative probabilities of OH. radical attack to all reaction sites are compared to experimental FSB and ARB values. By a fitting procedure, the relative efficiencies of conversion of the C4′ and C5′-centred radicals into FSB, ϵ (C4′): ϵ (C5′), and the relative efficiencies of base radicals—to—ARB conversion, ϵ (T): ϵ(A): ϵ(C): ϵ(G), are then deduced for each DNA form. The ability of the model to account for the distribution of damages in DNA-ligand complexes is proven by its successful application to two DNA-protein systems: the lac repressor-lac operator complex and the nucleosome core.

Sequence-structure relationships in DNA oligomers: a computational approach.

A collective-variable model for DNA structure is used to predict the conformation of a set of 30 octamer, decamer, and dodecamer oligomers for which high-resolution crystal structures are available. The model combines an all-atom base pair representation with an empirical backbone, emphasizing the role of base stacking in fixing sequence-dependent structure. We are able to reproduce trends in roll and twist to within 5 degrees across a large database of both A- and B-DNA oligomers. A genetic algorithm approach is used to search for global minimum structures and this is augmented by a grid search to identify local minimums. We find that the number of local minimums is highly sequence dependent, with certain sequences having a set of minimums that span the entire range between canonical A- and B-DNA conformations. Although the global minimum does not always agree with the crystal structure, for 24 of the 30 oligomers, we find low-energy local minimums that match the experimental step parameters. Discrepancies throw some light on the role of crystal packing in determining the solid-state conformation of double-helical DNA.

The crystal structures of psoralen cross-linked DNAs: drug-dependent formation of Holliday junctions

The single-crystal structures are presented for two DNA sequences with the thymine bases covalently cross-linked across the complementary strands by 4′-hydroxymethyl-4,5′,8-trimethylpsoralen (HMT). The HMT-adduct of d(CCGCTAGCGG) forms a psoralen-induced Holliday junction, showing for the first time the effect of this important class of chemotheraputics on the structure of the recombination intermediate. In contrast, HMT-d(CCGGTACCGG) forms a sequence-dependent junction. In both structures, the DNA duplex is highly distorted at the thymine base linked to the six-member pyrone ring of the drug. The psoralen cross-link defines the intramolecular interactions of the drug-induced junction, while the sequence-dependent structure is nearly identical to the native Holliday junction of d(CCGGTACCGG) alone. The two structures contrast the effects of drug- and sequence-dependent interactions on the structure of a Holliday junction, suggesting a role for psoralen in the mechanism to initiate repair of psoralen-lesions in mammalian DNA.