Machine learning shows torsion angle preferences in left-handed and right-handed quadruplex DNAs

Left-handed G quadruplexes (LHG4) have been recently discovered as a new class of G quadruplexes. The biological functions of LHG4s are still unknown, but they share a striking resemblance to Z-DNA in their helicity and jagged phosphate backbone. To further understand structural features of the LHG4s that define their left handedness, we have employed human-interpretable machine-learning methods to classify right- and left-handed G4s purely based on torsional angle analysis. Our results reveal the importance of the α, β, δ, and χ angles in left-handed structuring across both Z-DNAs and LHG4s. Our analysis may serve as the first step to understanding the conditions of formation for LHG4s and their potential biological relevance.     

β‐Bulges: Extensive structural analyses of β‐sheets irregularities

β‐Sheets are quite frequent in protein structures and are stabilized by regular main‐chain hydrogen bond patterns. Irregularities in β‐sheets, named β‐bulges, are distorted regions between two consecutive hydrogen bonds. They disrupt the classical alternation of side chain direction and can alter the directionality of β‐strands. They are implicated in protein‐protein interactions and are introduced to avoid β‐strand aggregation. Five different types of β‐bulges are defined. Previous studies on β‐bulges were performed on a limited number of protein structures or one specific family. These studies evoked a potential conservation during evolution. In this work, we analyze the β‐bulge distribution and conservation in terms of local backbone conformations and amino acid composition. Our dataset consists of 66 times more β‐bulges than the last systematic study (Chan et al. Protein Science 1993, 2:1574–1590). Novel amino acid preferences are underlined and local structure conformations are highlighted by the use of a structural alphabet. We observed that β‐bulges are preferably localized at the N‐ and C‐termini of β‐strands, but contrary to the earlier studies, no significant conservation of β‐bulges was observed among structural homologues. Displacement of β‐bulges along the sequence was also investigated by Molecular Dynamics simulations.     

Mechanical diversity and folding intermediates of parallel-stranded G-quadruplexes with a bulge

A significant number of sequences in the human genome form noncanonical G-quadruplexes (G4s) with bulges or a guanine vacancy. Here, we systematically characterized the mechanical stability of parallel-stranded G4s with a one to seven nucleotides bulge at various positions. Our results show that G4-forming sequences with a bulge form multiple conformations, including fully-folded G4 with high mechanical stability (unfolding forces > 40 pN), partially-folded intermediates (unfolding forces < 40 pN). The folding probability and folded populations strongly depend on the positions and lengths of the bulge. By combining a single-molecule unfolding assay, dimethyl sulfate (DMS) footprinting, and a guanine-peptide conjugate that selectively stabilizes guanine-vacancy-bearing G-quadruplexes (GVBQs), we identified that GVBQs are the major intermediates of G4s with a bulge near the 5′ or 3′ ends. The existence of multiple structures may induce different regulatory functions in many biological processes. This study also demonstrates a new strategy for selectively stabilizing the intermediates of bulged G4s to modulate their functions.     

A left-handed RNA double helix bound by the Z alpha domain of the RNA-editing enzyme ADAR1

The A form RNA double helix can be transformed to a left-handed helix, called Z-RNA. Currently, little is known about the detailed structural features of Z-RNA or its involvement in cellular processes. The discovery that certain interferon-response proteins have domains that can stabilize Z-RNA as well as Z-DNA opens the way for the study of Z-RNA. Here, we present the 2.25 A crystal structure of the Zalpha domain of the RNA-editing enzyme ADAR1 (double-stranded RNA adenosine deaminase) complexed to a dUr(CG)(3) duplex RNA. The Z-RNA helix is associated with a unique solvent pattern that distinguishes it from the otherwise similar conformation of Z-DNA. Based on the structure, we propose a model suggesting how differences in solvation lead to two types of Z-RNA structures. The interaction of Zalpha with Z-RNA demonstrates how the interferon-induced isoform of ADAR1 could be targeted toward selected dsRNAs containing purine-pyrimidine repeats, possibly of viral origin.     

The role of DNA sequence in the formation of Z-DNA versus cruciforms in plasmids

The capacities of four synthetic sequences containing runs of perfectly alternating purine-pyrimidine base pairs (bp) to adopt left-handed structures were evaluated in a homologous family of recombinant plasmids. All the sequences had the same G+C content (50%) and consisted of simple tetranucleotide repeat units but differed in the relative orientations of these units. For some of the sequences, several alternate secondary structures were theoretically possible; a variety of probes (S1 nuclease, bromoacetaldehyde, OsO4, T7 gene 3 endonuclease, supercoil-induced gel relaxation studies) under a wide range of reaction conditions was used to determine which structures were adopted as a function of superhelical stress. The precise positions at the bp level of reactions with these chemical and enzymatic probes were determined. We conclude that for short (20-24 bp) sequences containing runs of alternating (T-G) and (C-A), the cruciform state is preferred over the similarly allowable left-handed form provided that symmetry constraints allow. However, these sequences can be induced to form a left-handed helix under appropriate conditions. This is the first demonstration of plasmid inserts which will adopt more than one unusual DNA structure in response to negative superhelical stress. The structural properties of a molecule containing a Z-Z junction were studied, and we conclude that the disruption caused by this feature extends over only a few bp although it requires a high energetic penalty.     

Role of amino acid residues at turns in the conformational stability and folding of human lysozyme

To clarify the role of amino acid residues at turns in the conformational stability and folding of a globular protein, six mutant human lysozymes deleted or substituted at turn structures were investigated by calorimetry, GuHCl denaturation experiments, and X-ray crystal analysis. The thermodynamic properties of the mutant and wild-type human lysozymes were compared and discussed on the basis of their three-dimensional structures. For the deletion mutants, Delta47-48 and Delta101, the deleted residues are in turns on the surface and are absent in human alpha-lactalbumin, which is homologous to human lysozyme in amino acid sequence and tertiary structure. The stability of both mutants would be expected to increase due to a decrease in conformational entropy in the denatured state; however, both proteins were destabilized. The destabilizations were mainly caused by the disappearance of intramolecular hydrogen bonds. Each part deleted was recovered by the turn region like the alpha-lactalbumin structure, but there were differences in the main-chain conformation of the turn between each deletion mutant and alpha-lactalbumin even if the loop length was the same. For the point mutants, R50G, Q58G, H78G, and G37Q, the main-chain conformations of these substitution residues located in turns adopt a left-handed helical region in the wild-type structure. It is thought that the left-handed non-Gly residue has unfavorable conformational energy compared to the left-handed Gly residue. Q58G was stabilized, but the others had little effect on the stability. The structural analysis revealed that the turns could rearrange the main-chain conformation to accommodate the left-handed non-Gly residues. The present results indicate that turn structures are able to change their main-chain conformations, depending upon the side-chain features of amino acid residues on the turns. Furthermore, stopped-flow GuHCl denaturation experiments on the six mutants were performed. The effects of mutations on unfolding-refolding kinetics were significantly different among the mutant proteins. The deletion/substitutions in turns located in the alpha-domain of human lysozyme affected the refolding rate, indicating the contribution of turn structures to the folding of a globular protein.     

Pattern recognition in nucleic acid sequences. II. An efficient method for finding locally stable secondary structures.

We present a method for calculating all possible single hairpin loop secondary structures in a nucleic acid sequence by the order of N2 operations where N is the total number of bases. Each structure may contain any number of bulges and internal loops. Most natural sequences are found to be indistinguishable from random sequences in the potential of forming secondary structures, which is defined by the frequency of possible secondary structures calculated by the method. There is a strong correlation between the higher G+C content and the higher structure forming potential. Interestingly, the removal of intervening sequences in mRNAs is almost always accompanied by an increase in the G+C content, which may suggest an involvement of structural stabilization in the mRNA maturation.     

Twist-joints and double twist-joints in RNA structure

Analysis of available RNA crystal structures has allowed us to identify a new family of RNA arrangements that we call double twist-joints, or DTJs. Each DTJ is composed of a double helix that contains two bulges incorporated into different strands and separated from each other by 2 or 3 bp. At each bulge, the double helix is over-twisted, while the unpaired nucleotides of both bulges form a complex network of stacking and hydrogen-bonding with nucleotides of helical regions. In total, we identified 14 DTJ cases, which can be combined in three groups based on common structural characteristics. One DTJ is found in a functional center of the ribosome, another DTJ mediates binding of the pre-tRNA to the RNase P, and two more DTJs form the sensing domains in the glycine riboswitch.     

L-nucleotides and 8-methylguanine of d(C1m8G2C3G4C5LG6LC7G8C9G10)2 act cooperatively to promote a left-handed helix under physiological salt conditions

The structure and thermal stability of a hetero chiral decaoligodeoxyribonucleotide duplex d(C1m8 G2C3G4C5LG6LC7G8C9G10)d(C11m8G12C13G14C15LG16LC17G18C19G20) (O1) with two contiguous pairs of enantiomeric 2'-deoxy-L-ribonucleotides (C5LG6L/C15LG16L) at its centre and an 8-methylguanine at position 2/12 was analysed by circular dichroism, NMR and molecular modelling. O1 resolves in a left-handed helical structure already at low salt concentration (0.1 M NaCl). The central L2-sugar portion assumes a B* left-handed conformation (mirror-image of right-handed B-DNA) while its flanking D4-sugar portions adopt the known Z left-handed conformation. The resulting Z4-B2*-Z4 structure (left-handed helix) is the reverse of that of B4-Z2*-B4 (right-handed helix) displayed by the nearly related decaoligodeoxyribonucleotide d(mC1G2mC3G4C5L G6LmC7G8mC9G10)2, at the same low salt concentration (0.1 M NaCl). In the same experimental conditions, d(C1m8G2C3G4C5G6C7G8C9G10)2 (O2), the stereoregular version of O1, resolves into a right-handed B-DNA helix. Thus, both the 8-methylguanine and the enantiomeric step CLpGL at the centre of the molecule are needed to induce left-handed helicity. Remarkably, in the various heterochiral decaoligodeoxyribonucleotides so far analysed by us, when the central CLpGL adopts the B* (respectively Z*) conformation, then the adjacent steps automatically resolves in the Z (respectively B) conformation. This allows a good optimisation of the base-base stackings and base-sugar van der Waals interactions at the ZB*/B*Z (respectively BZ*/Z*B) junctions so that the Z4-B2*-Z4 (respectively B4-Z2*-B4) helix displays a Tm (approximately 65 degrees C) that is only 5 degrees C lower than the one of its homochiral counterpart. Here we anticipate that a large variety of DNA helices can be generated at low salt concentration by manipulating internal factors such as sugar configuration, duplex length, nucleotide composition and base methylation. These helices can constitute powerful tools for structural and biological investigations, especially as they can be used in physiological conditions.     

RNA bulges as architectural and recognition motifs

RNA bulges constitute versatile structural motifs in the assembly of RNA architectures. Three-dimensional structures of RNA molecules and their complexes reveal the role of bulges in RNA architectures and illustrate the molecular mechanisms by which they confer intramolecular interactions and intermolecular recognition.     

The Possible Structural Models for Polyglutamine Aggregation: A Molecular Dynamics Simulations Study

Abstract

Huntington's disease is a neurodegenerative disorder caused by a polyglutamine (polyQ) expansion near the N-terminus of huntingtin. Previous studies have suggested that polyQ aggregation occurs only when the number of glutamine (Q) residues is more than 36-40, the disease threshold. However, the structural characteristics of polyQ nucleation in the very early stage of aggregation still remain elusive. In this study, we designed 18 simulation trials to determine the possible structural models for polyQ nucleation and aggregation with various shapes and sizes of initial β-helical structures, such as left-handed circular, right-handed rectangular, and left- and right-handed triangular. Our results show that the stability of these models significantly increases with increasing the number of rungs, while it is rather insensitive to the number of Qs in each rung. In particular, the 3-rung β-helical models are stable when they adopt the left-handed triangular and right-handed rectangular conformations due to the fact that they preserve high β-turn and β-sheet contents, respectively, during the simulation courses. Thus, we suggested that these two stable β-helical structures with at least 3 rungs might serve as the possible nucleation seeds for polyQ depending on how the structural elements of β-turn and β-sheet are sampled and preserved during the very early stage of aggregation.     

A comprehensive analysis of the Greek key motifs in protein beta-barrels and beta-sandwiches

The Greek key motifs are the topological signature of many beta-barrels and a majority of beta-sandwich structures. An updated survey of these structures integrates many early observations and newly emerging patterns and provides a better understanding of the unique role of Greek keys in protein structures. A stereotypical Greek key beta-barrel accommodates five or six strands and can have 12 possible topologies. All except one six-stranded topologies have been observed, and only one five-stranded topologies have been seen in actual structures. Of the representative beta-barrel structures analyzed here, half have left-handed Greek keys. This result challenges the empirical claim of the handedness regularity of Greek keys in beta-barrels. One of the five-stranded topologies that has not been observed in beta-barrels comprises two overlapping Greek keys. The two three-dimensional forms of this topology constitute a structural unit that is present in a vast majority of known beta-sandwich structures. Using this unit as the root, we have built a new taxonomy tree for the beta-sandwich folds and deduced a set of rules that appear to constrain how other beta-strands adjoin the unit to form a larger double-layered structure. These rules, though derived from a larger data set, are essentially the same as those drawn from earlier studies, suggesting that they may reflect the true topological constraints in the design of beta-sandwich structures. Finally, a novel variant of the Greek key motif (defined here as the twisted Greek key) has emerged which introduces loop crossings into the folded structures. Proteins 2000;40:409-419.     

Sheared-type G(anti).C(syn) base-pair: a unique d(GXC) loop closure motif

Stable DNA loop structures closed by a novel G.C base-pair have been determined for the single-residue d(GXC) loops (X=A, T, G or C) in low-salt solution by high-resolution nuclear magnetic resonance (NMR) techniques. The closing G.C base-pair in these loops is not of the canonical Watson-Crick type, but adopts instead a unique sheared-type (trans Watson-Crick/sugar-edge) pairing, like those occurring in the sheared mismatched G.A or A.C base-pair, to draw the two opposite strands together. The cytidine residue in the closing base-pair is transformed into the rare syn domain to form two H-bonds with the guanine base and to prevent the steric clash between the G 2NH(2) and the C H-5 protons. Besides, the sugar pucker of the syn cytidine is still located in the regular C2'-endo domain, unlike the C3'-endo domain adopted for the pyrimidines of the out-of-alternation left-handed Z-DNA structure. The facile formation of the compact d(GXC) loops closed by a unique sheared-type G(anti).C(syn) base-pair demonstrates the great potential of the single-stranded d(GXC) triplet repeats to fold into stable hairpins.     

Selection on the structural stability of a ribosomal RNA expansion segment in Daphnia obtusa

The high rate of sequence divergence in nuclear ribosomal RNA (rRNA) expansion segments offers a unique opportunity to study the importance of natural selection in their evolution. To this end, we polymerase chain reaction amplified and cloned a 589-nt fragment of the 18S rRNA gene containing expansion segments 43/e1 and 43/e4 from six individual Daphnia obtusa from four populations. We screened 2,588 clones using single-stranded conformation polymorphism analysis and identified 103 unique haplotype sequences. We detected two pairs of indel sites in segment 43/e4 that complement each other when the secondary structure of the linear sequence is formed. Seven of the 12 observed combinations of length variants at these four sites (haplotypes) are shared between individuals from different populations, which may suggest that some of the length variation was present in their common ancestor. Haplotypes with uncompensated indels were only observed at low frequencies, while compensated indel haplotypes were found at a wide range of frequencies, supporting the hypothesis that the energetic stability of expansion segments is a trait under natural selection. In addition, there was strong linkage disequilibrium between the four complementary indel sites, particularly those that pair with one another in the secondary structure. Despite selection against unpaired bulges at these four indel sites, some nucleotides that form unpaired bulges are highly conserved in segment 43/e4, indicating that they are under a different selective constraint, possibly due to their role in higher level structural interactions.     

Chirality and Handedness of Protein Structures

In proteins, the polypeptide chain forms a number of right-and left-handed helices and superhelices, right-and left-turned hairpins, and some other structures that are nonsuperimposable, although they are not mirror images of each other as the Lamino acids are not converted to the Damino acids. This property of protein structures will be referred to here as pseudo-chirality–or handedness. It has been shown that there are two kinds of handedness in proteins–helical handedness and handedness of arrangement. Some protein structures exhibit both the kinds of handedness. Handedness is observed at all levels of protein structural organization–from α-helices, β-strands, hairpins, βαβ-units up to complex structural motifs, superhelices, and supramolecular structures in fibrous and polymer proteins. There are several structures that have unique handedness in proteins, for example, α-helices, αα-corners, βαβ-units, abcd-units, and so on. This property of the polypeptide chain is of particular value in protein folding and protein modeling, because it drastically reduces the number of possible folds.     

Effect of single-residue bulges on RNA double-helical structures: crystallographic database analysis and molecular dynamics simulation studies

Asymmetric bulge loop motifs are widely dispersed in all types of functional RNAs. They are frequently occurring structural motifs in folded RNA structures and appear commonly in pre-microRNA and ribosomes, where they are involved in specific RNA–RNA and RNA–protein interactions. It is therefore necessary to understand such motifs from a structural point of view. We analyzed all available RNA structures and identified quite a few fragments of double helices that contain bulges. We found that these discontinuities often introduce kinks into the double helices, which also affects the stacking overlap between the base pairs across the irregularity. In order to understand the influence of these bulges on stability and flexibility, we carried out molecular dynamics simulations of three different single-residue bulge-containing RNA helices using the CHARMM36 force field. The structural variability at the junctions of RNA bulges is expected to differ from that in continuous double-helical stretches. The structural features of the junction region were observed to vary noticeably depending on the orientation of the bulge residue. When the base of the bulge residue is looped out, the RNA stretch behaves like a standard long A-form RNA double helix, whereas the entire RNA behaves differently when the base of the bulge residue is intercalated between base pairs inside the RNA stem. Such single-base intercalation was found to introduce a permanent kink into the composite double helix, which could be a recognition element for Dicer during the maturation of miRNA.     

A novel structural motif and structural trees for proteins containing it

In the present study, a novel structural motif that can be represented as a combination of the known βαβ-unit and ψ-motif is described and analyzed. In theory, there are four possible combinations of the motifs since each of them can exist in two forms, left-handed and right-handed. For this study, we have selected 140 nonhomologous proteins in which 158 combinations of such types have been found. The combination of the right-handed ψ-motif and the right-handed βαβ-unit has been shown to occur most often (87 cases out of 158) and the combination of the left-handed βαβ-unit and the left-handed ψ-motif does not occur at all. Three novel structural trees in which the commonly occurring combinations are taken as the root structures have been constructed.     

NMR solution structures of bistranded abasic site lesions in DNA

Ionizing radiation produces clustered lesions in DNA. Since the orientation of bistranded lesions affects their recognition by DNA repair enzymes, clustered damages are more difficult to process and thus more toxic than single oxidative lesions. In order to understand the structural determinants that lead to differential recognition, we used NMR spectroscopy and restrained molecular dynamics to solve the structure of two DNA duplexes, each containing two stable abasic site analogues positioned on opposite strands of the duplex and staggered in the 3' (-1 duplex, (AP) 2-1 duplex) or 5' (+1 duplex, (AP) 2+1 duplex) direction. Cross-peak connectivities observed in the nonexchangeable NOESY spectra indicate compression of the helix at the lesion site of the duplexes, resulting in the formation of two abasic bulges. The exchangeable proton spectra show the AP site partner nucleotides forming interstrand hydrogen bonds that are characteristic of a Watson-Crick G.C base pairs, confirming the extra helical nature of the AP residues. Restrained molecular dynamics simulations generate a set of converging structures in full agreement with the spectroscopic data. In the (AP) 2-1 duplex, the extra helical abasic site residues reside in the minor groove of the helix, while they appear in the major groove in the (AP) 2+1 duplex. These structural differences are consistent with the differential recognition of bistranded abasic site lesions by human AP endonuclease.     

AA.AG@helix.ends: A:A and A:G base-pairs at the ends of 16 S and 23 S rRNA helices.

This study reveals that AA and AG oppositions occur frequently at the ends of helices in RNA crystal and NMR structures in the PDB database and in the 16 S and 23 S rRNA comparative structure models, with the G usually 3' to the helix for the AG oppositions. In addition, these oppositions are frequently base-paired and usually in the sheared conformation, although other conformations are present in NMR and crystal structures. These A:A and A:G base-pairs are present in a variety of structural environments, including GNRA tetraloops, E and E-like loops, interfaced between two helices that are coaxially stacked, tandem G:A base-pairs, U-turns, and adenosine platforms. Finally, given structural studies that reveal conformational rearrangements occurring in regions of the RNA with AA and AG oppositions at the ends of helices, we suggest that these conformationally unique helix extensions might be associated with functionally important structural rearrangements.