A single G-to-C change causes human centromere TGGAA repeats to fold back into hairpins

Structure of native four-repeat satellite III sequence with non-canonical base interactions

Tandem-repetitive DNA (where two or more DNA bases are repeated numerous times) can adopt non-canonical secondary structures. Many of these structures are implicated in important biological processes. Human Satellite III (HSat3) is enriched for tandem repeats of the sequence ATGGA and is located in pericentromeric heterochromatin in many human chromosomes. Here, we investigate the secondary structure of the four-repeat HSat3 sequence 5'-ATGGA ATGGA ATGGA ATGGA-3' using X-ray crystallography, NMR, and biophysical methods. Circular dichroism spectroscopy, thermal stability, native PAGE, and analytical ultracentrifugation indicate that this sequence folds into a monomolecular hairpin with non-canonical base pairing and B-DNA characteristics at concentrations below 0.9 mM. NMR studies at 0.05-0.5 mM indicate that the hairpin is likely folded-over into a compact structure with high dynamics. Crystallographic studies at 2.5 mM reveal an antiparallel self-complementary duplex with the same base pairing as in the hairpin, extended into an infinite polymer. The non-canonical base pairing includes a G-G intercalation sandwiched by sheared A-G base pairs, leading to a cross-strand four guanine stack, so called guanine zipper. The guanine zippers are spaced throughout the structure by A-T/T-A base pairs. Our findings lend further insight into recurring structural motifs associated with the HSat3 and their potential biological functions.

STRavinsky STR database and PGTailor PGT tool demonstrate superiority of CHM13-T2T over hg38 and hg19 for STR-based applications

Short-Tandem-Repeats (STRs) have long been studied for possible roles in biological phenomena, and are utilized in multiple applications such as forensics, evolutionary studies and pre-implantation-genetic-testing (PGT). The two reference genomes most used by clinicians and researchers are GRCh37/hg19 and GRCh38/hg38, both constructed using mainly short-read-sequencing (SRS) in which all-STR-containing-reads cannot be assembled to the reference genome. With the introduction of long-read-sequencing (LRS) methods and the generation of the CHM13 reference genome, also known as T2T, many previously unmapped STRs were finally localized within the human genome. We generated STRavinsky, a compact STR database for three reference genomes, including T2T. We proceeded to demonstrate the advantages of T2T over hg19 and hg38, identifying nearly double the number of STRs throughout all chromosomes. Through STRavinsky, providing a resolution down to a specific genomic coordinate, we demonstrated extreme propensity of TGGAA repeats in p arms of acrocentric chromosomes, substantially corroborating early molecular studies suggesting a possible role in formation of Robertsonian translocations. Moreover, we delineated unique propensity of TGGAA repeats specifically in chromosome 16q11.2 and in 9q12. Finally, we harness the superior capabilities of T2T and STRavinsky to generate PGTailor, a novel web application dramatically facilitating design of STR-based PGT tests in mere minutes.

Predicting the aptamer SYL3C-EpCAM complex's structure with the Martini-based simulation protocol

Aptamers (small single strand DNA/RNAs) such as SYL3C are considered as ideal alternatives to antibodies in cancer related research studies. However, 3D structure predictions for aptamers and aptamer-protein complexes are scarce due to the high cost of experimental measurements and unreliable computer-based methods. Thus aptamers' diagnostic and therapeutic applications are severely restricted. To meet the challenge, we proposed a Martini-based aptamer-protein complex prediction protocol. By combining the base-base contact map from simulation and secondary structure prediction from various tools, improved secondary structure predictions can be obtained. This method reduced the risk of providing incorrect or incomplete base pairs in secondary structure prediction. Thus 3D structure modeling based on the secondary structure can be more reliable. We introduced the soft elastic network to the hairpin folded regions of the Martini ssDNAs to preserve their canonical structure. Using our protocol, we predicted the first 3D structure of the aptamer SYL3C and the SYL3C-EpCAM complex. We believe that our work could contribute to the future aptamer-related research studies and medical implications.

Stable inheritance of CENP-A chromatin: Inner strength versus dynamic control

Chromosome segregation during cell division is driven by mitotic spindle attachment to the centromere region on each chromosome. Centromeres form a protein scaffold defined by chromatin featuring CENP-A, a conserved histone H3 variant, in a manner largely independent of local DNA cis elements. CENP-A nucleosomes fulfill two essential criteria to epigenetically identify the centromere. They undergo self-templated duplication to reestablish centromeric chromatin following DNA replication. More importantly, CENP-A incorporated into centromeric chromatin is stably transmitted through consecutive cell division cycles. CENP-A nucleosomes have unique structural properties and binding partners that potentially explain their long lifetime in vivo. However, rather than a static building block, centromeric chromatin is dynamically regulated throughout the cell cycle, indicating that CENP-A stability is also controlled by external factors. We discuss recent insights and identify the outstanding questions on how dynamic control of the long-term stability of CENP-A ensures epigenetic centromere inheritance.

Centromeres under Pressure: Evolutionary Innovation in Conflict with Conserved Function

Centromeres are essential genetic elements that enable spindle microtubule attachment for chromosome segregation during mitosis and meiosis. While this function is preserved across species, centromeres display an array of dynamic features, including: (1) rapidly evolving DNA; (2) wide evolutionary diversity in size, shape and organization; (3) evidence of mutational processes to generate homogenized repetitive arrays that characterize centromeres in several species; (4) tolerance to changes in position, as in the case of neocentromeres; and (5) intrinsic fragility derived by sequence composition and secondary DNA structures. Centromere drive underlies rapid centromere DNA evolution due to the "selfish" pursuit to bias meiotic transmission and promote the propagation of stronger centromeres. Yet, the origins of other dynamic features of centromeres remain unclear. Here, we review our current understanding of centromere evolution and plasticity. We also detail the mutagenic processes proposed to shape the divergent genetic nature of centromeres. Changes to centromeres are not simply evolutionary relics, but ongoing shifts that on one side promote centromere flexibility, but on the other can undermine centromere integrity and function with potential pathological implications such as genome instability.

Multiple roles of DNA2 nuclease/helicase in DNA metabolism, genome stability and human diseases

DNA2 nuclease/helicase is a structure-specific nuclease, 5'-to-3' helicase, and DNA-dependent ATPase. It is involved in multiple DNA metabolic pathways, including Okazaki fragment maturation, replication of 'difficult-to-replicate' DNA regions, end resection, stalled replication fork processing, and mitochondrial genome maintenance. The participation of DNA2 in these different pathways is regulated by its interactions with distinct groups of DNA replication and repair proteins and by post-translational modifications. These regulatory mechanisms induce its recruitment to specific DNA replication or repair complexes, such as DNA replication and end resection machinery, and stimulate its efficient cleavage of various structures, for example, to remove RNA primers or to produce 3' overhangs at telomeres or double-strand breaks. Through these versatile activities at replication forks and DNA damage sites, DNA2 functions as both a tumor suppressor and promoter. In normal cells, it suppresses tumorigenesis by maintaining the genomic integrity. Thus, DNA2 mutations or functional deficiency may lead to cancer initiation. However, DNA2 may also function as a tumor promoter, supporting cancer cell survival by counteracting replication stress. Therefore, it may serve as an ideal target to sensitize advanced DNA2-overexpressing cancers to current chemo- and radiotherapy regimens.

TRAF6 mediates human DNA2 polyubiquitination and nuclear localization to maintain nuclear genome integrity

The multifunctional human DNA2 (hDNA2) nuclease/helicase is required to process DNA ends for homology-directed recombination repair (HDR) and to counteract replication stress. To participate in these processes, hDNA2 must localize to the nucleus and be recruited to the replication or repair sites. However, because hDNA2 lacks the nuclear localization signal that is found in its yeast homolog, it is unclear how its migration into the nucleus is regulated during replication or in response to DNA damage. Here, we report that the E3 ligase TRAF6 binds to and mediates the K63-linked polyubiquitination of hDNA2, increasing the stability of hDNA2 and promoting its nuclear localization. Inhibiting TRAF6-mediated polyubiquitination abolishes the nuclear localization of hDNA2, consequently impairing DNA end resection and HDR. Thus, the current study reveals a mechanism for the regulation of hDNA2 localization and establishes that TRAF6-mediated hDNA2 ubiquitination activates DNA repair pathways to maintain nuclear genome integrity.

High Satellite Repeat Turnover in Great Apes Studied with Short- and Long-Read Technologies

Satellite repeats are a structural component of centromeres and telomeres, and in some instances, their divergence is known to drive speciation. Due to their highly repetitive nature, satellite sequences have been understudied and underrepresented in genome assemblies. To investigate their turnover in great apes, we studied satellite repeats of unit sizes up to 50 bp in human, chimpanzee, bonobo, gorilla, and Sumatran and Bornean orangutans, using unassembled short and long sequencing reads. The density of satellite repeats, as identified from accurate short reads (Illumina), varied greatly among great ape genomes. These were dominated by a handful of abundant repeated motifs, frequently shared among species, which formed two groups: 1) the (AATGG)n repeat (critical for heat shock response) and its derivatives; and 2) subtelomeric 32-mers involved in telomeric metabolism. Using the densities of abundant repeats, individuals could be classified into species. However, clustering did not reproduce the accepted species phylogeny, suggesting rapid repeat evolution. Several abundant repeats were enriched in males versus females; using Y chromosome assemblies or Fluorescent In Situ Hybridization, we validated their location on the Y. Finally, applying a novel computational tool, we identified many satellite repeats completely embedded within long Oxford Nanopore and Pacific Biosciences reads. Such repeats were up to 59 kb in length and consisted of perfect repeats interspersed with other similar sequences. Our results based on sequencing reads generated with three different technologies provide the first detailed characterization of great ape satellite repeats, and open new avenues for exploring their functions.

Non-B-Form DNA Is Enriched at Centromeres

Animal and plant centromeres are embedded in repetitive "satellite" DNA, but are thought to be epigenetically specified. To define genetic characteristics of centromeres, we surveyed satellite DNA from diverse eukaryotes and identified variation in <10-bp dyad symmetries predicted to adopt non-B-form conformations. Organisms lacking centromeric dyad symmetries had binding sites for sequence-specific DNA-binding proteins with DNA-bending activity. For example, human and mouse centromeres are depleted for dyad symmetries, but are enriched for non-B-form DNA and are associated with binding sites for the conserved DNA-binding protein CENP-B, which is required for artificial centromere function but is paradoxically nonessential. We also detected dyad symmetries and predicted non-B-form DNA structures at neocentromeres, which form at ectopic loci. We propose that centromeres form at non-B-form DNA because of dyad symmetries or are strengthened by sequence-specific DNA binding proteins. This may resolve the CENP-B paradox and provide a general basis for centromere specification.

Genome-wide Control of Heterochromatin Replication by the Telomere Capping Protein TRF2

Hard-to-replicate regions of chromosomes (e.g., pericentromeres, centromeres, and telomeres) impede replication fork progression, eventually leading, in the event of replication stress, to chromosome fragility, aging, and cancer. Our knowledge of the mechanisms controlling the stability of these regions is essentially limited to telomeres, where fragility is counteracted by the shelterin proteins. Here we show that the shelterin subunit TRF2 ensures progression of the replication fork through pericentromeric heterochromatin, but not centromeric chromatin. In a process involving its N-terminal basic domain, TRF2 binds to pericentromeric Satellite III sequences during S phase, allowing the recruitment of the G-quadruplex-resolving helicase RTEL1 to facilitate fork progression. We also show that TRF2 is required for the stability of other heterochromatic regions localized throughout the genome, paving the way for future research on heterochromatic replication and its relationship with aging and cancer.

Centromere and Pericentromere Transcription: Roles and Regulation … in Sickness and in Health

The chromosomal loci known as centromeres (CEN) mediate the equal distribution of the duplicated genome between both daughter cells. Specifically, centromeres recruit a protein complex named the kinetochore, that bi-orients the replicated chromosome pairs to the mitotic or meiotic spindle structure. The paired chromosomes are then separated, and the individual chromosomes segregate in opposite direction along the regressing spindle into each daughter cell. Erroneous kinetochore assembly or activity produces aneuploid cells that contain an abnormal number of chromosomes. Aneuploidy may incite cell death, developmental defects (including genetic syndromes), and cancer (>90% of all cancer cells are aneuploid). While kinetochores and their activities have been preserved through evolution, the CEN DNA sequences have not. Hence, to be recognized as sites for kinetochore assembly, CEN display conserved structural themes. In addition, CEN nucleosomes enclose a CEN-exclusive variant of histone H3, named CENP-A, and carry distinct epigenetic labels on CENP-A and the other CEN histone proteins. Through the cell cycle, CEN are transcribed into non-coding RNAs. After subsequent processing, they become key components of the CEN chromatin by marking the CEN locus and by stably anchoring the CEN-binding kinetochore proteins. CEN transcription is tightly regulated, of low intensity, and essential for differentiation and development. Under- or overexpression of CEN transcripts, as documented for myriad cancers, provoke chromosome missegregation and aneuploidy. CEN are genetically stable and fully competent only when they are insulated from the surrounding, pericentromeric chromatin, which must be silenced. We will review CEN transcription and its contribution to faithful kinetochore function. We will further discuss how pericentromeric chromatin is silenced by RNA processing and transcriptionally repressive chromatin marks. We will report on the transcriptional misregulation of (peri)centromeres during stress, natural aging, and disease and reflect on whether their transcripts can serve as future diagnostic tools and anti-cancer targets in the clinic.

Repetitive Fragile Sites: Centromere Satellite DNA As a Source of Genome Instability in Human Diseases

Maintenance of an intact genome is essential for cellular and organismal homeostasis. The centromere is a specialized chromosomal locus required for faithful genome inheritance at each round of cell division. Human centromeres are composed of large tandem arrays of repetitive alpha-satellite DNA, which are often sites of aberrant rearrangements that may lead to chromosome fusions and genetic abnormalities. While the centromere has an essential role in chromosome segregation during mitosis, the long and repetitive nature of the highly identical repeats has greatly hindered in-depth genetic studies, and complete annotation of all human centromeres is still lacking. Here, we review our current understanding of human centromere genetics and epigenetics as well as recent investigations into the role of centromere DNA in disease, with a special focus on cancer, aging, and human immunodeficiency⁻centromeric instability⁻facial anomalies (ICF) syndrome. We also highlight the causes and consequences of genomic instability at these large repetitive arrays and describe the possible sources of centromere fragility. The novel connection between alpha-satellite DNA instability and human pathological conditions emphasizes the importance of obtaining a truly complete human genome assembly and accelerating our understanding of centromere repeats' role in physiology and beyond.

Sequence-Dependent Effects of Monovalent Cations on the Structural Dynamics of Trinucleotide-Repeat DNA Hairpins

Repetitive trinucleotide DNA sequences at specific genetic loci are associated with numerous hereditary, neurodegenerative diseases. The propensity of single-stranded domains containing these sequences to form secondary structure via extensive self-complementarity disrupts normal DNA processing to create genetic instabilities. To investigate these intrastrand structural dynamics, a DNA hairpin system was devised for single-molecule fluorescence study of the folding kinetics and energetics for secondary structure formation between two interacting, repetitive domains with specific numbers of the same trinucleotide motif (CXG), where X = T or A. Single-molecule fluorescence resonance energy transfer (smFRET) data show discrete conformational transitions between unstructured and closed hairpin states. The lifetimes of the closed hairpin states correlate with the number of repeats, with (CTG) _N/(CTG) _N domains maintaining longer-lived, closed states than equivalent-sized (CAG) _N/(CAG) _N domains. NaCl promotes similar degree of stabilization for the closed hairpin states of both repeat sequences. Temperature-based, smFRET experiments reveal that NaCl favors hairpin closing for (CAG) _N/(CAG) _N by preordering single-stranded repeat domains to accelerate the closing transition. In contrast, NaCl slows the opening transition of CTG hairpins; however, it promotes misfolded conformations that require unfolding. Energy diagrams illustrate the distinct folding pathways of (CTG) _N and (CAG) _N repeat domains and identify features that may contribute to their gene-destabilizing effects.

hDNA2 nuclease/helicase promotes centromeric DNA replication and genome stability

DNA2 is a nuclease/helicase that is involved in Okazaki fragment maturation, replication fork processing, and end resection of DNA double‐strand breaks. Similar such helicase activity for resolving secondary structures and structure‐specific nuclease activity are needed during DNA replication to process the chromosome‐specific higher order repeat units present in the centromeres of human chromosomes. Here, we show that DNA2 binds preferentially to centromeric DNA. The nuclease and helicase activities of DNA2 are both essential for resolution of DNA structural obstacles to facilitate DNA replication fork movement. Loss of DNA2‐mediated clean‐up mechanisms impairs centromeric DNA replication and CENP‐A deposition, leading to activation of the ATR DNA damage checkpoints at centromeric DNA regions and late‐S/G2 cell cycle arrest. Cells that escape arrest show impaired metaphase plate formation and abnormal chromosomal segregation. Furthermore, the DNA2 inhibitor C5 mimics DNA2 knockout and synergistically kills cancer cells when combined with an ATR inhibitor. These findings provide mechanistic insights into how DNA2 supports replication of centromeric DNA and give further insights into new therapeutic strategies.

Parity-dependent hairpin configurations of repetitive DNA sequence promote slippage associated with DNA expansion

Repetitive DNA sequences are ubiquitous in life, and changes in the number of repeats often have various physiological and pathological implications. DNA repeats are capable of interchanging between different noncanonical and canonical conformations in a dynamic fashion, causing configurational slippage that often leads to repeat expansion associated with neurological diseases. In this report, we used single-molecule spectroscopy together with biophysical analyses to demonstrate the parity-dependent hairpin structural polymorphism of TGGAA repeat DNA. We found that the DNA adopted two configurations depending on the repeat number parity (even or odd). Transitions between these two configurations were also observed for longer repeats. In addition, the ability to modulate this transition was found to be enhanced by divalent ions. Based on the atomic structure, we propose a local seeding model where the kinked GGA motifs in the stem region of TGGAA repeat DNA act as hot spots to facilitate the transition between the two configurations, which may give rise to disease-associated repeat expansion.

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.

Maintenance and replication of the human cytomegalovirus genome during latency

Human cytomegalovirus (HCMV) can establish latent infection in hematopoietic progenitor cells (HPCs) or CD14 (+) monocytes. While circularized viral genomes are observed during latency, how viral genomes persist or which viral factors contribute to genome maintenance and/or replication is unclear. Previously, we identified a HCMV cis-acting viral maintenance element (TR element) and showed that HCMV IE1 exon 4 mRNA is expressed in latently infected HPCs. We now show that a smaller IE1 protein species (IE1x4) is expressed in latently infected HPCs. IE1x4 protein expression is required for viral genome persistence and maintenance and replication of a TR element containing plasmid (pTR). Both IE1x4 and the cellular transcription factor Sp1 interact with the TR, and inhibition of Sp1 binding abrogates pTR amplification. Further, IE1x4 interacts with Topoisomerase IIβ (TOPOIIβ), whose activity is required for pTR amplification. These results identify a HCMV latency-specific factor that promotes viral chromosome maintenance and replication.

The size of the internal loop in DNA hairpins influences their targeting with partially complementary strands

Targeting of noncanonical DNA structures, such as hairpin loops, may have significant diagnostic and therapeutic potential. Oligonucleotides can be used for binding to mRNA, forming a DNA/RNA hybrid duplex that inhibits translation. This kind of modulation of gene expression is called the antisense approach. In order to determine the best strategy to target a common structural motif in mRNA, we have designed a set of stem-loop DNA molecules with sequence: d(GCGCT_nGTAAT₅GTTACT_nGCGC), where n = 1, 3, or 5, “T₅” is an end loop of five thymines. We used a combination of calorimetric and spectroscopy techniques to determine the thermodynamics for the reaction of a set of hairpins containing internal loops with their respective partially complementary strands. Our aim was to determine if internal- and end-loops are promising regions for targeting with their corresponding complementary strands. Indeed, all targeting reactions were accompanied by negative changes in free energy, indicating that reactions proceed spontaneously. Further investigation showed that these negative free energy terms result from a net balance of unfavorable entropy and favorable enthalpy contributions. In particular, unfolding of hairpins and duplexes is accompanied by positive changes in heat capacity, which may be a result of exposure of hydrophobic groups to the solvent. This study provides a new method for the targeting of mRNA in order to control gene expression.

Lentiviral vector integration profiles differ in rodent postmitotic tissues

Lentiviral vectors with self-inactivating (SIN) long terminal repeats (LTRs) are promising for safe and sustained transgene expression in dividing as well as quiescent cells. As genome organization and transcription substantially differs between actively dividing and postmitotic cells in vivo, we hypothesized that genomic vector integration preferences might be distinct between these biological states. We performed integration site (IS) analyses on mouse dividing cells (fibroblasts and hematopoietic progenitor cells (HPCs)) transduced ex vivo and postmitotic cells (eye and brain) transduced in vivo. As expected, integration in dividing cells occurred preferably into gene coding regions. In contrast, postmitotic cells showed a close to random frequency of integration into genes and gene spare long interspersed nuclear elements (LINE). Our studies on the potential mechanisms responsible for the detected differences of lentiviral integration suggest that the lowered expression level of Psip1 reduce the integration frequency in vivo into gene coding regions in postmitotic cells. The motif TGGAA might represent one of the factors for preferred lentiviral integration into mouse and rat Satellite DNA. These observations are highly relevant for the correct assessment of preclinical biosafety studies, indicating that lentiviral vectors are well suited for safe and effective clinical gene transfer into postmitotic tissues.

Nucleic acid folding determined by mesoscale modeling and NMR spectroscopy: solution structure of d(GCGAAAGC)

Determination of DNA solution structure is a difficult task even with the high-sensitivity method used here based on simulated annealing with 35 restraints/residue (Cryoprobe 750 MHz NMR). The conformations of both the phosphodiester linkages and the dinucleotide segment encompassing the sharp turn in single-stranded DNA are often underdetermined. To obtain higher quality structures of a DNA GNRA loop, 5'-d(GCGAAAGC)-3', we have used a mesoscopic molecular modeling approach, called Biopolymer Chain Elasticity (BCE), to provide reference conformations. By construction, these models are the least deformed hairpin loop conformation derived from canonical B-DNA at the nucleotide level. We have further explored this molecular conformation at the torsion angle level with AMBER molecular mechanics using different possible (epsilon,zeta) constraints to interpret the 31P NMR data. This combined approach yields a more accurate molecular conformation, compatible with all the NMR data, than each method taken separately, NMR/DYANA or BCE/AMBER. In agreement with the principle of minimal deformation of the backbone, the hairpin motif is stabilized by maximal base-stacking interactions on both the 5'- and 3'-sides and by a sheared G.A mismatch base pair between the first and last loop nucleotides. The sharp turn is located between the third and fourth loop nucleotides, and only two torsion angles beta6 and gamma6 deviate strongly with respect to canonical B-DNA structure. Two other torsion angle pairs epsilon3,zeta3 and epsilon5,zeta5 exhibit the newly recognized stable conformation BIIzeta+ (-70 degrees, 140 degrees). This combined approach has proven to be useful for the interpretation of an unusual 31P chemical shift in the 5'-d(GCGAAAGC)-3' hairpin.

Centromere-competent DNA: structure and evolution

Although extant data favour centromere being an epigenetic structure, it is also clear that centromere formation is based on DNA, in particular, tandemly repeated satellite DNA and its transcripts. Presence of conserved structural motifs within satellite DNAs such as periodically distributed AT tracts, protein binding sites, or promoter elements indicate that despite sequence flexibility, there are structural determinants that are prerequisite for centromere function. In addition, existence of functional centromeric DNA transcripts indicates possible importance of structural elements at the level of RNA secondary or tertiary structure. Rapid centromere evolution is explained by homologous recombination followed by extrachromosomal rolling circle replication. This could lead to amplification of different satellite sequences within a genome. However, only those satellites that have inherent centromere-competence in the form of structural requirements necessary for centromere function are after amplification fixed in a population as a new centromere.

Hairpin structures formed by alpha satellite DNA of human centromeres are cleaved by human topoisomerase IIalpha

Although centromere function has been conserved through evolution, apparently no interspecies consensus DNA sequence exists. Instead, centromere DNA may be interconnected through the formation of certain DNA structures creating topological binding sites for centromeric proteins. DNA topoisomerase II is a protein, which is located at centromeres, and enzymatic topoisomerase II activity correlates with centromere activity in human cells. It is therefore possible that topoisomerase II recognizes and interacts with the alpha satellite DNA of human centromeres through an interaction with potential DNA structures formed solely at active centromeres. In the present study, human topoisomerase IIα-mediated cleavage at centromeric DNA sequences was examined in vitro. The investigation has revealed that the enzyme recognizes and cleaves a specific hairpin structure formed by alpha satellite DNA. The topoisomerase introduces a single-stranded break at the hairpin loop in a reaction, where DNA ligation is partly uncoupled from the cleavage reaction. A mutational analysis has revealed, which features of the hairpin are required for topoisomerease IIα-mediated cleavage. Based on this a model is discussed, where topoisomerase II interacts with two hairpins as a mediator of centromere cohesion.

Spectroscopic and structural impact of a stem base-pair change in DNA hairpins: GTTC-ACA-GAAC versus GTAC-ACA-GTAC

Successive investigations over the last decade have revealed and confirmed a stable loop closure in a family of d-[GTAC-5Pur6N7N-GTAC] hairpins, where 5Pur6N7N is a AAA, GAG and AXC loop (X being any nucleotide). The trinucleotide loop is characterized by a well defined 5Pur-7N mispairing mode, and by upfield chemical shifts for three sugar protons of the apical nucleotide 6N. The GTTC-ACA-GAAC DNA hairpin, of interest for its likely involvement in Vibrio cholerae genome mutations, has now been investigated. The GTAC-ACA-GTAC DNA hairpin has also been studied because it is intermediate between the other structures, as it contains the loop of the hairpin under consideration and the stem of the above family. The two hairpins with the ACA loop are stable. They show the same mispairing mode and similar upfield shifts as the previous family, but GTTC-ACA-GAAC seems to be slightly less compact than any other. GTTC-ACA-GAAC is remarkable in that it exhibits a B(II) character for the phosphate-ester conformation at 8Gp9A, together with a swing of the upper hairpin into the major groove that, in particular, brings 6CH1' roughly as close to 7AH2 as to 6CH6. These unexpected structural features are qualitatively deduced from (1)H and (31)P NMR spectra, and confirmed by Raman spectroscopy. This comparative study shows that not only the loop sequence but also the stem sequence may control hairpin structures.

Quadruplex-duplex competition in the nuclease hypersensitive element of human c-myc promoter: C to T mutation in C-rich strand enhances duplex association

The nuclease hypersensitive element NHE III(I) is an important anti-cancer target as the transcription of oncogene c-myc is largely regulated by it. It has been postulated that regulatory control is mediated by G-quadruplex formation in the NHE anti-sense strand through a competition between the duplex and the quadruplex states. A mutation in the NHE has been implicated in cancer. In this study, the reported mutation has been characterized vis-a-vis the kinetics of i-tetraplex formation (in the sense strand) and its effect on duplex formation. We found that i-tetraplex formation was destabilized by approximately 1.4 kcal/mol (DeltaDeltaG at 20 degrees C, pH 5.8). Observed hysteresis allowed us to analyze the kinetics of folding for the mutant (M3). Though we observed higher association (DeltaEon approximately -23.4 kcal/mol) and dissociation (DeltaEoff approximately 22.1 kcal/mol) activation energies (at pH 5.3) for the wild-type (P1) tetraplex folding, the kinetics of folding and unfolding for M3 was somewhat faster at pH 5.3 and 5.8. Interestingly, Surface plasmon resonance (BIAcore) analysis of hybridization at pH 6.6 indicated a higher association constant for M3 (approximately 22.5 x 10(4)M(-1)s(-1)) than P1 (approximately 3.2 x 10(4)M(-1)s(-1)). The equilibrium dissociation constants also indicated favorable duplex association for M3 (approximately 22.2 and approximately 190.6 nM for M3 and P1, respectively). We envisage that the increased affinity for the duplex state due to the mutation could play a functional role in the aberrant regulation of c-myc.

Unusual DNA duplex and hairpin motifs

Single-stranded DNA or double-stranded DNA has the potential to adopt a wide variety of unusual duplex and hairpin motifs in the presence (trans) or absence (cis) of ligands. Several principles for the formation of those unusual structures have been established through the observation of a number of recurring structural motifs associated with different sequences. These include: (i) internal loops of consecutive mismatches can occur in a B-DNA duplex when sheared base pairs are adjacent to each other to confer extensive cross- and intra-strand base stacking; (ii) interdigitated (zipper-like) duplex structures form instead when sheared G*A base pairs are separated by one or two pairs of purine*purine mismatches; (iii) stacking is not restricted to base, deoxyribose also exhibits the potential to do so; (iv) canonical G*C or A.T base pairs are flexible enough to exhibit considerable changes from the regular H-bonded conformation. The paired bases become stacked when bracketed by sheared G.A base pairs, or become extruded out and perpendicular to their neighboring bases in the presence of interacting drugs; (v) the purine-rich and pyrimidine-rich loop structures are notably different in nature. The purine-rich loops form compact triloop structures closed by a sheared G*A, A*A, A*C or sheared-like G(anti)*C(syn) base pair that is stacked by a single residue. On the other hand, the pyrimidine-rich loops with a thymidine in the first position exhibit no base pairing but are characterized by the folding of the thymidine residue into the minor groove to form a compact loop structure. Identification of such diverse duplex or hairpin motifs greatly enlarges the repertoire for unusual DNA structural formation.

DNA tri- and tetra-loops and RNA tetra-loops hairpins fold as elastic biopolymer chains in agreement with PDB coordinates

The biopolymer chain elasticity (BCE) approach and the new molecular modelling methodology presented previously are used to predict the tri- dimensional backbones of DNA and RNA hairpin loops. The structures of eight remarkably stable DNA or RNA hairpin molecules closed by a mispair, recently determined in solution by NMR and deposited in the PDB, are shown to verify the predicted trajectories by an analysis automated for large numbers of PDB conformations. They encompass: one DNA tetraloop, -GTTA-; three DNA triloops, -AAA- or -GCA-; and four RNA tetraloops, -UUCG-. Folding generates no distortions and bond lengths and bond angles of main atoms of the sugar-phosphate backbone are well restored upon energy refinement. Three different methods (superpositions, distance of main chain atoms to the elastic line and RMSd) are used to show a very good agreement between the trajectories of sugar-phosphate backbones and between entire molecules of theoretical models and of PDB conformations. The geometry of end conditions imposed by the stem is sufficient to dictate the different characteristic DNA or RNA folding shapes. The reduced angular space, consisting of the new parameter, angle Omega, together with the chi angle offers a simple, coherent and quantitative description of hairpin loops.

Quadruple intercalated G-6 stack: a possible motif in the fold-back structure of the Drosophila centromeric dodeca-satellite?

The purine-rich strand d(GTACGGGACCGA)(n) of the Drosophila centromeric dodeca-satellite sequence is highly conserved and was found to form stable fold-back structures in which the homopurine 5'-GGGA-3' sequence was determined to play a crucial role. Here, we report the stable formation of the d(GGGA)(2) motif in the stem of a DNA hairpin closed by a single-residue d(ACC) loop. Similar to the zipper-like d(GGA)(2) motif observed in the human centromeric (TGGAA)(n) sequence, the central four guanosine bases in the d(GGGA)(2) motif do not pair, but interdigitate to form an elongated zipper-like quadruple-intercalated G-6 stack bracketed by sheared G.A base-pairs. Comparison between the current d(GGGA)(2) structure and the published crystal d(GAAA)(2) structure implies that the alignment of the unpaired purine bases plays an important role in determining the minor groove width of the purine-rich d(GPuPuA)(2) motif. Similarity between the zipper-like motifs possibly present in the Drosophila centromeric dodeca-satellite sequence and in the human centromeric (TGGAA)(n) sequence led us to propose that these special zipper-like motifs may constitute common cores in organizing eukaryotic centromeres.

Solution structure of a DNA double helix incorporating four consecutive non-Watson-Crick base-pairs

A series of DNA 21-mers containing a variety of the 4 x 4 internal loop sequence 5'-CAAG-3'/3'-ACGT-5' were studied using nuclear magnetic resonance (NMR) methodology and distance geometry (DG)/molecular dynamics (MD) approaches. Such oligomers exhibit excellent resolution in the NMR spectra and reveal many unusual NOEs (nuclear Overhauser effect) that allow for the detailed characterization of a DNA hairpin incorporating a track of four different non-Watson-Crick base-pairs in the stem. These include a wobble C.A base-pair, a sheared A.C base-pair, a sheared A.G base-pair, and a wobble G.T base-pair. Significantly different twisting angles were observed between the base-pairs in internal loop that results with excellent intra-strand and inter-strand base stacking within the four consecutive mismatches and the surrounding canonical base-pairs. This explains why it melts at 52 degrees C even though five out of ten base-pairs in the stem adopt non-Watson-Crick pairs. However, the 4 x 4 internal loop still fits into a B-DNA double helix very well without significant change in the backbone torsion angles; only zeta torsion angles between the tandem sheared base-pairs are changed to a great extent from the gauche(-) domain to the trans domain to accommodate the cross-strand base stacking in the internal loop. The observation that several consecutive non-canonical base-pairs can stably co-exist with Watson-Crick base-pairs greatly increases the limited repertoire of irregular DNA folds and reveals the possibility for unusual structural formation in the functionally important genomic regions that have potential to become single-stranded.

Zipper-like Watson-Crick base-pairs

A series of DNA heptadecamers containing the DNA analogues of RNA E-like 5'-d(GXA)/(AYG)-5' motifs (X/Y is complementary T/A, A/T, C/G, or G/C pair) were studied using nuclear magnetic resonance (NMR) methodology and distance geometry (DG)/molecular dynamics (MD) approaches. Such oligomers reveal excellent resolution in NMR spectra and exhibit many unusual nuclear Overhauser effects (NOEs) that allow for good characterization of an unusual zipper-like conformation with zipper-like Watson-Crick base-pairs; the potential canonical X.Y H-bonding is not present, and the central X/Y pairs are transformed instead into inter-strand stacks that are bracketed by sheared G.A base-pairs. Such phenomenal structural change is brought about mainly through two backbone torsional angle adjustments, i.e. delta from C2'-endo to C3'-endo for the sugar puckers of unpaired residues and gamma from gauche(+) to trans for the following 3'-adenosine residues. Such motifs are analogous to the previously studied (GGA)(2) motif presumably present in the human centromeric (TGGAA)(n) tandem repeat sequence. The novel zipper-like motifs are only 4-7 deg. C less stable than the (GGA)(2) motif, suggesting that inter-strand base stacking plays an important role in stabilizing unusual nucleic acid structures. The discovery that canonical Watson-Crick G.C or A.T hydrogen-bonded pairs can be transformed into stacking pairs greatly increases the repertoire for unusual nucleic acid structural motifs.

Structure of a C-rich strand fragment of the human centromeric satellite III: a pH-dependent intercalation topology

Repetitive DNA sequences may adopt unusual pairing arrangements. At acid to neutral pH, cytidine-rich DNA oligodeoxynucleotides can form the i-motif structure in which two parallel-stranded duplexes with C.C(+) pairs are intercalated head-to-tail. The i-motif may be formed by multimeric associations or by intra-molecular folding, depending on the number of cytidine tracts, the nucleotide sequences between them, and the experimental conditions. We have found that a natural fragment of the human centromeric satellite III, d(CCATTCCATTCCTTTCC), can form two monomeric i-motif structures that differ in their intercalation topology and that are favored at pH values higher (the eta-form) and lower (the lambda-form) than 4.6. The change in intercalation may be related to adenine protonation in the loops. We studied the uridine derivative methylated on the first cytidine base, d(5mCCATTCCAUTCCUTTCC), whose proton spectrum is better resolved. The intercalation topologies are (C7.C17)/(5mC1.C11)/(C6.C16)/(C2.C12) for form lambda and (5mC1.C11)/(C7.C17)/(C2.C12)/(C6.C16) for form eta. We have solved the structure of the eta-form, and we present a model for the lambda-form. The switch from eta to lambda involves disruption of the i-motif. In both forms, the central AUT linker crosses the wide groove, and the first and the third linkers loop across the minor grooves. The i-motif core is extended in the eta-form by the inter-loop reverse Watson-Crick A3.U13 pair, whose dissociation constant is around 10(-2) at 0 degrees C, and in the lambda-form by the interloop T5.T15 pair. In contrast, d(5mCCATTCCTTACCTTTCC) folds into a pH-independent structure that has the same intercalation topology as the lambda-form. The i-motif core is extended below by the interloop T5.T15 pair and closed on top by the T8.A10 pair.Thus, the C-rich strand of the human satellite III tandem repeats, like the G-rich strand, can fold into various compact structures. The relevance of these features to centromeric function remains unknown.

Conformational analysis of DNA-trinucleotide-hairpin-loop structures using a continuum solvent model

A number of trinucleotide sequences in DNA can form compact and stable hairpin loops that may have significance for DNA replication and transcription. The conformational analysis of these motifs is important for an understanding of the function and design of nucleic acid structures. Extensive conformational searches have been performed on three experimentally known trinucleotide hairpin loops (AGC, AAA, and GCA) closed by a four-base-pair stem. An implicit solvation model based on the generalized Born method has been employed during energy minimization and conformational search. In addition, energy-minimized conformers were evaluated using a finite-difference Poisson-Boltzmann approach. For all three loop sequences, conformations close to experiment were found as lowest-energy structures among several thousand alternative energy minima. The inclusion of reaction-field contributions was found to be important for a realistic conformer ranking. Most generated hairpin loop structures within approximately 5 kcal x mol(-1) of the lowest-energy structure have a similar topology. Structures within approximately 10 kcal x mol(-1) could be classified into about five structural families representing distinct arrangements of loop nucleotides. Although a large number of backbone torsion angle combinations were compatible with each structural class, some specific patterns could be identified. Harmonic mode analysis was used to account for differences in conformational flexibility of low-energy sub-states. Class-specific differences in the pattern of atomic fluctuations along the sequence were observed; however, inclusion of conformational entropy contributions did not change ranking of structural classes. For an additional loop sequence (AAG) with no available experimental structure, the approach suggests a lowest-energy loop topology overall similar to the other three loop sequences but closed by a different non-canonical base-pairing scheme.

Natural abundance heteronuclear NMR studies of the T3 mini-loop hairpin in the terminal repeat of the adenoassociated virus 2

A DNA hairpin containing a T3 loop, as occurs in the terminal repeat of a popular gene therapy vector (Adenoassociated Virus 2, AAV2), has been extensively studied using homo- and heteronuclear NMR experiments. Almost complete assignment of the proton and carbon resonances, including H5'(Pro-S) and H5'(Pro-R) protons, has been accomplished at natural abundance. NOESY spectra in H2O and D2O have revealed many unusual NOEs, which, when combined with the epsilon, beta, gamma, and chi torsion angles determined from heteronuclear 1H-13C, 1H-31P, and 13C-31P coupling constants, have allowed for a more detailed picture of the T3 mini-loop hairpin. The three loop thymidines are all unpaired, yet are highly structured when bracketed by a 5'-GC...GC-3' stem sequence. The structure determined in this manuscript is considerably different from several other structures reported so far. Contrary to an RNA oligomer with a central U3 sequence that has the tendency to form a duplex with three U*U mismatches, the d(GAAGC-TTT-GCTTC) sequence exists mostly as a hairpin under millimolar NMR conditions. Since T3 triloop was found to be an essential element for the site-specific non-homologous integration of the AAV2 virus, and modification of the T3 loop residue abolishes such capability, the structure we report here may be of biological significance.

Stable formation of a pyrimidine-rich loop hairpin in a cruciform promoter

We have determined the solution structure of a TCC-loop hairpin in the cruciform promoter for the bacteriophage N4 virion RNA polymerase (N4 vRNAP). This hairpin and its complementary GGA-loop hairpin are extruded at physiological superhelical density and are required for vRNAP recognition. Contrary to its complementary GGA-loop, the three pyrimidines in the TCC-loop are all unpaired. However, with the help of two juxtaposed stem Watson-Crick G.C base-pairs, each nucleotide in the loop employs a special method to stabilize the hairpin structure. The resulting structures display extensive loop base-stacking rearrangement yet minor backbone distortion, which is largely accomplished through some loop zeta and alpha torsional angle changes. Consistent with the structural studies, UV melting of the GAAGCTCCGCTTC hairpin revealed a higher melting temperature (66 degrees C) than that of the GAACGTCCCGTTC hairpin (58 degrees C) with reversed stem G.C base-pairs, indicating significant contribution from the extra three loop-stem H-bonds. Thermodynamic parameters DeltaG degrees 25of the GAAGCTCCGCTTC hairpin and its complementary GAAGCGGAGCTTC hairpin are -4.1 and -4. 3 kcal/mol respectively, indicating approximately equal contribution of each hairpin to the cruciform formation of the N4 virion RNA polymerase promoter. No significant loop dynamics in the microsecond to millisecond NMR time-scale was observed, and the abundant well-defined exchangeable and non-exchangeable proton NOEs allowed us to efficiently determine a well-converged family for the final structures of the TCC-loop hairpin.

Stable sheared A.C pair in DNA hairpins

Single-residue d(Pu1NPu2) (Pu1.Pu2=G.A, G.G or A.A) hairpin loops can be stably closed by sheared purine.purine pairs. These special motifs have been found in several important biological systems. We now extend these loop-closing base-pairs to a sheared purine. pyrimidine (A.C) pair at a neutral pH condition. High-resolution NMR spectroscopy, distance geometry, and molecular dynamics methods were used to study d(GTACANCGTAC) oligomers. Numerous idiosyncratic nuclear Overhauser enhancements, especially those across the A.C base-pair between C4NH2left and right arrow AH1', C4NH2left and right arrow AH2, and CH5left and right arrow AH2 proton pairs, clearly define the novel sheared nature of the closing A.C base-pair. This novel base-pair is possibly present in several biological systems and in two single-stranded DNA aptamers selected from oligonucleotide libraries.

The telomeres of Streptomyces chromosomes contain conserved palindromic sequences with potential to form complex secondary structures

The chromosomes of the gram-positive soil bacteria Streptomyces are linear DNA molecules, usually of about 8Mb, containing a centrally located origin of replication and covalently bound terminal proteins (which are presumably involved in the completion of replication of the telomeres). The ends of the chromosomes contain inverted repeats of variable lengths. The terminal segments of five Streptomyces chromosomes and plasmids were cloned and sequenced. The sequences showed a high degree of conservation in the first 166-168bp. Beyond the terminal homology, the sequences diverged and did not generally cross-hybridize. The homologous regions contained seven palindromes with a few nucleotide differences. Many of these differences occur in complementary pairs, such that the palindromicity is preserved. Energy-optimized modelling predicted that the 3' strand of the terminal palindromes can form extensive hairpin structures that are similar to the 3' ends of autonomous parvovirus genomes. Most of the putative hairpins have a GCGCAGC sequence at the loop, with the potential to form a stable single C-residue loop closed by a sheared G:A pairing. The similarity between the terminal structures of the Streptomyces replicons and the autonomous parvoviral genomes suggests that they may share some structural and/or replication features.

Centromeric pyrimidine strands fold into an intercalated motif by forming a double hairpin with a novel T:G:G:T tetrad: solution structure of the d(TCCCGTTTCCA) dimer

The solution structures of the oligodeoxynucleotides d(CCCGTTTCC) and d(TCCCGTTTCCA) have been determined by two-dimensional NMR spectroscopy. These oligomers are part of a DNA box in human centromeric alpha satellite targeted by the centromere protein B (CENP-B). Both CENP-B and its recognition box in alphoid DNA are conserved in mammals, suggesting an important biological role. At acidic pH, d(CCCGTTTCC), d(TCCCGTTTCCA) and the full d(TCCCGTTTCCAACGAAG) CENP-B box strand all fold and dimerize in solution forming a stable bimolecular structure containing two GTTT hairpin loops that interact through a novel T : G : G : T tetrad. The stem region of the dimer is a four-stranded intercalated motif in which the hairpin monomers are parallel and held together by C : C+ hydrogen-bonding and intercalation. The loops are at the same end of the dimer and lie across the narrow grooves of the tetraplex. They are remarkably structured and stabilized by base-base cross-stacking, sugar-base stacking, and parallel G:G and antiparallel G:T pairing. In the d(TCCCGTTTCCA)2 structure, the intercalated motif is continued at the other end of the dimer with unpaired but stacked adenine and thymine bases. The possible biological implications of these structures are discussed.

Unusual DNA conformations

DNA is on the move across conformational space. Duplexes diversity and, joined by triplexes, quadruplexes, loops, bulges and multiarmed junctions, open the route to a bewildering array of increasingly complex conformations. In addition to this structural growth, DNA has come under increasing scrutiny thanks to the development of chemical and physical techniques for deforming its conformation and probing its properties. These investigations help us to learn more about the mechanics and the activity of this remarkably versatile macromolecule.

Sheared purine x purine pairing in biology

The Watson-Crick G x C and A x T base-paired DNA duplex has been the single most important milestone in modem molecular biology. However, it is possible that other types of stable DNA structures besides the double helix might exist, since only about 5% of the human chromosome is transcribed and expressed. Stable, four-stranded G-tetraplex DNA structures occur in the extensive tandem repeated sequences at the telomeres of chromosome. Formation of stable triplexes of the Py x Pu x Py or Pu x Pu x Py type have been implicated at the control regions of certain human genes. We review and discuss the various types of DNA duplex structures containing stable sheared base-pairs and compare their structural characteristics with that of B-DNA. Pu x Pu structural motifs are found in the highly conserved sequences at the replication origins of several single-stranded DNA viruses and in the peri-centromeric regions of human chromosomes, and may be involved in important biological functions, such as viral DNA replication and centromere formation.