Archaeal histone selection of nucleosome positioning sequences and the procaryotic origin of histone-dependent genome evolution

Archaeal histones and the eucaryal (eucaryotic) nucleosome core histones have almost identical histone folds. Here, we show that DNA molecules selectively incorporated by rHMfB (recombinant archaeal histone B from Methanothermus fervidus) into archaeal nucleosomes from a mixture of approximately 10(14) random sequence molecules contain sequence motifs shown previously to direct eucaryal nucleosome positioning. The dinucleotides GC, AA (=TT) and TA are repeated at approximately 10 bp intervals, with the GC harmonic displaced approximately 5 bp from the AA and TA harmonics [(GCN(3)AA or TA)(n)]. AT and CG were not strongly selected, indicating that TA not equalAT and GC not equalCG in terms of facilitating archaeal nucleosome assembly. The selected molecules have affinities for rHMfB ranging from approximately 9 to 18-fold higher than the level of affinity of the starting population, and direct the positioned assembly of archaeal nucleosomes. Fourier-transform analyses have revealed that AA dinucleotides are much enriched at approximately 10. 1 bp intervals, the helical repeat of DNA wrapped around a nucleosome, in the genomes of Eucarya and the histone-containing Euryarchaeota, but not in the genomes of Bacteria and Crenarchaeota, procaryotes that do not have histones. Facilitating histone packaging of genomic DNA has apparently therefore imposed constraints on genome sequence evolution, and since archaeal histones have no structure in addition to the histone fold, these constraints must result predominantly from histone fold-DNA contacts. Based on the three-domain universal phylogeny, histones and histone-dependent genome sequence evolution most likely evolved after the bacterial-archaeal divergence but before the archaeal-eucaryal divergence, and were subsequently lost in the Crenarchaeota. However, with lateral gene transfer, the first histone fold could alternatively have evolved after the archaeal-eucaryal divergence, early in either the euryarchaeal or eucaryal lineages.     

Histone stoichiometry and DNA circularization in archaeal nucleosomes.

Recombinant (r)HMfB (archaealhistone B fromMethanothermusfervidus) formed complexes with increasing stability with DNA molecules increasing in length from 52 to 100 bp, but not with a 39 bp molecule. By using125I-labeled rHMfB-YY (an rHMfB variant with I31Y and M35Y replacements) and32P-labeled 100 bp DNA, these complexes, designated archaeal nucleosomes, have been shown to contain an archaeal histone tetramer. Consistent with DNA bending and wrapping, addition of DNA ligase to archaeal nucleosomes assembled with 88 and 128 bp DNAs resulted in covalently-closed monomeric circular DNAs which, following histone removal, were positively supercoiled based on their electrophoretic mobilities in the presence of ethidium bromide before and after relaxation by calf thymus topoisomerase I. Ligase addition to mixtures of rHMfB with 53 or 30 bp DNA molecules also resulted in circular DNAs but these were circular dimers and trimers. These short DNA molecules apparently had to be ligated into longer linear multimers for assembly into archaeal nucleosomes and ligation into circles. rHMfB assembled into archaeal nucleosomes at lower histone to DNA ratios with the supercoiled, circular ligation product than with the original 88 bp linear version of this molecule. Archaeal histones are most similar to the globular histone fold region of eukaryal histone H4, and the results reported are consistent with archaeal nucleosomes resembling the structure formed by eukaryal histone (H3+H4)2tetramers.     

Molecular components of the archaeal nucleosome

Here we describe the organization of the archaeal nucleosome, in which four archaeal histones are circumscribed by approximately 80 bp of DNA. Through a combination of sequence comparisons, 3D structural studies, site-directed mutagenesis and assays for DNA binding, we have assigned functions to most of the individual residues in the histone fold of the representative archaeal histone rHMfB. By SELEX selection, the sequences of DNA molecules that are most readily bound and wrapped by rHMfB into archaeal nucleosomes in vitro have been identified, and these define DNA structures that position archaeal nucleosome assembly.     

Archaeal histone tetramerization determines DNA affinity and the direction of DNA supercoiling

DNA binding and the topology of DNA have been determined in complexes formed by >20 archaeal histone variants and archaeal histone dimer fusions with residue replacements at sites responsible for histone fold dimer:dimer interactions. Almost all of these variants have decreased affinity for DNA. They have also lost the flexibility of the wild type archaeal histones to wrap DNA into a negative or positive supercoil depending on the salt environment; they wrap DNA into positive supercoils under all salt conditions. The histone folds of the archaeal histones, HMfA and HMfB, from Methanothermus fervidus are almost identical, but (HMfA)(2) and (HMfB)(2) homodimers assemble into tetramers with sequence-dependent differences in DNA affinity. By construction and mutagenesis of HMfA+HMfB and HMfB+HMfA histone dimer fusions, the structure formed at the histone dimer:dimer interface within an archaeal histone tetramer has been shown to determine this difference in DNA affinity. Therefore, by regulating the assembly of different archaeal histone dimers into tetramers that have different sequence affinities, the assembly of archaeal histone-DNA complexes could be localized and used to regulate gene expression.     

Crystal structures of recombinant histones HMfA and HMfB from the hyperthermophilic archaeon Methanothermus fervidus

The hyperthermophilic archaeon Methanothermus fervidus contains two small basic proteins, HMfA (68 amino acid residues) and HMfB (69 residues) that share a common ancestry with the eukaryal nucleosome core histones H2A, H2B, H3, and H4. HMfA and HMfB have sequences that differ at 11 locations, they have different structural stabilities, and the complexes that they form with DNA have different electrophoretic mobilities. Here, crystal structures are documented for recombinant (r) HMfA at a resolution of 1.55 A refined to a crystallographic R-value of 19.8 % (tetragonal form) and at 1.48 A refined to a R-value of 18.8 % (orthorhombic form), and for rHMfB at 1.9 A refined to a R-value of 18.0 %. The rHMfA and rHMfB monomers have structures that are just histone folds in which a long central alpha-helix (alpha2; 29 residues) is separated from shorter N-terminal (alpha1; 11 residues) and C-terminal (alpha3; 10 residues) alpha-helices by two loops (L1 and L2; both 6 residues). Within L1 and L2, three adjacent residues are in extended (beta) conformation. rHMfA and rHMfB assemble into homodimers, with the alpha2 helices anti-parallel aligned and crossing at an angle of close to 35 degrees, and with hydrogen bonds formed between the extended, parallel regions of L1 and L2 resulting in short beta-ladders. Dimerization creates a novel N-terminal structure that contains four proline residues, two from each monomer. As prolines are present at these positions in all archaeal histone sequences, this proline-tetrad structure is likely to be a common feature of all archaeal histone dimers. Almost all residues that participate in monomer-monomer interactions are conserved in HMfA and HMfB, consistent with the ability of these monomers to form both homodimers and (HMfA+HMfB) heterodimers. Differences in side-chain interactions that result from non-conservative residue differences in HMfA and HMfB are identified, and the structure of a (rHMfA)(2)-DNA complex is presented based on the structures documented here and modeled by homology to histone-DNA interactions in the eukaryal nucleosome.     

Conserved eukaryotic histone-fold residues substituted into an archaeal histone increase DNA affinity but reduce complex flexibility

Although the archaeal and eukaryotic nucleosome core histones evolved from a common ancestor, conserved lysine residues are present at DNA-binding locations in all four eukaryotic histones that are not present in the archaeal histones. Introduction of lysine residues at the corresponding locations into an archaeal histone, HMfB, generated a variant with increased affinity for DNA that formed more compact complexes with DNA. However, these complexes no longer facilitated the circularization of short DNA molecules and had lost the flexibility to wrap DNA alternatively in either a negative or positive supercoil.     

Mutational analysis of differences in thermostability between histones from mesophilic and hyperthermophilic archaea

Amino acid residues responsible for the large difference in thermostability between HMfB and HFoB, archaeal histones from the hyperthermophile Methanothermus fervidus and the mesophile Methanobacterium formicicum, respectively, have been identified by site-specific mutagenesis. The thermal denaturation of approximately 70 archaeal histone variants has been monitored by circular dichroism, and the data generated were fit to a two-state unfolding model (dimer-->two random coil monomers) to obtain a standard-state (1M) melting temperature for each variant dimer. The results of single-, double-, and triple-residue substitutions reveal that the much higher stability of rHMfB dimers, relative to rHFoB dimers, is conferred predominantly by improved intermolecular hydrophobic interactions near the center of the histone dimer core and by additional favorable ion pairs on the dimer surface.     

Archaeal chromatin proteins histone HMtB and Alba have lost DNA-binding ability in laboratory strains of <Emphasis Type="Italic">Methanothermobacter thermautotrophicus</Emphasis>

Alignments of the sequences of the all members of the archaeal histone and Alba1 families of chromatin proteins identified isoleucine residues, I19 in HMtB and I39 in MtAlba, in Methanothermobacter thermautotrophicus, at locations predicted to be directly involved in DNA binding. In all other HMfB family members, residue 19 is an arginine (R19), and either arginine or lysine is present in almost all other Alba1 family members at the structural site equivalent to I39 in MtAlba. Electrophoretic mobility shift assays revealed that recombinant HMtB and MtAlba do not bind DNA, but variants constructed with R19 and R39, respectively, bound DNA; and whereas MtAlba(I19) did not bind RNA, MtAlba(R19) bound both single stranded RNA and tRNA. Amplification and sequencing of MT0254 (encodes HMtB) and MT1483 (encodes MtAlba) from several Methanothermobacter thermautotrophicus lineages has revealed that HMtB and MtAlba had arginine residues at positions 19 and 39, respectively, in the original isolate and that spontaneous mutations must have occurred, and been fixed, in some laboratory lineages that now have HMtB(I19) and MtAlba(I39). The retention of these variants suggests some continuing functions and fusion of the HMtB(I19) sequence to HMtA2 resulted in a protein that folds to form a histone fold heterodimer that binds and compacts DNA. The loss of DNA binding by HMtB(I19) does not therefore prevent HMtB from participating in DNA interactions as one partner of an archaeal histone heterodimer.     

An alternative beads‐on‐a‐string chromatin architecture in Thermococcus kodakarensis

We have applied chromatin sequencing technology to the euryarchaeon Thermococcus kodakarensis, which is known to possess histone‐like proteins. We detect positioned chromatin particles of variable sizes associated with lengths of DNA differing as multiples of 30 bp (ranging from 30 bp to >450 bp) consistent with formation from dynamic polymers of the archaeal histone dimer. T. kodakarensis chromatin particles have distinctive underlying DNA sequence suggesting a genomic particle‐positioning code and are excluded from gene‐regulatory DNA suggesting a functional organization. Beads‐on‐a‐string chromatin is therefore conserved between eukaryotes and archaea but can derive from deployment of histone‐fold proteins in a variety of multimeric forms.     

Differences in the binding of H1 variants to DNA. Cooperativity and linker-length related distribution

A study of the complexes formed between short linear DNA and three H1 variants, a typical somatic H1, and the extreme variants H5, from chicken erythrocytes, and spH1 from sea urchin sperm, has revealed differences between H1, H5 and spH1 that have implications for chromatin structure and folding. 1. All three histones bind cooperatively to DNA in 35 mM NaCl forming similar, but not identical, rod-like complexes. With sufficiently long DNA the complexes may be circular, circles forming more easily with H5 and spH1 than with H1. 2. The binding of H5 and spH1 to DNA is cooperative even in 5 mM NaCl, resulting in well-defined thin filaments that appear to contain two DNA molecules bridged by histone molecules. In contrast, H1 binds distributively over all the DNA molecules in 5 mM NaCl, but forms short stretches similar in appearance to the thin filaments formed with H5 and spH1. Rods appear to arise from the intertwining of regular thin filaments containing cooperatively bound histone molecules on raising the NaCl concentration to 35 mM. 3. The compositions of the rods correspond to one histone molecule for about every 47 bp (H1), 81 bp (H5) and 112 bp (spH1), suggesting average spacings of 24 bp (H1), 41 bp (H5) and 56 bp (spH1) in the component thin (double) filaments. Strikingly, these values are proportional to the linker lengths of the chromatins in which the particular H1 variant is the main or sole H1.     

Analysis of DNA compaction profile and intracellular contents of archaeal histones from Pyrococcus kodakaraensis KOD1.

Two histone genes, hpkA and hpkB, from hyperthermophilic archaeon Pyrococcus kodakaraensis KOD1 strain were cloned, sequenced, and expressed in Escherichia coli cells. Both hpkA and hpkB genes encoded a protein of 67 amino acids, however they possessed the different molecular weight (HpkA, 7,378:HpkB, 7,167). Deduced amino acid sequences of HpkA and HpkB were homologous to other archaeal histones and eucaryal core histones (H2A, H4). Gel mobility shift assays by purified proteins demonstrated that HpkB possessed higher affinity to DNA and more extensive ability to compact DNA than HpkA. HpkB prevented double stranded DNA from thermal denaturation in less amount than HpkA in vitro. In order to investigate intracellular contents of HpkA and HpkB in KOD1 cells, immunoblot analysis was performed by using anti-HpkA antisera obtained from immunized BALB/c mice, showing that HpkA was less abundantly expressed than HpkB in KOD1 cells. These results suggest that HpkB plays a major role to protect double stranded DNA from thermal denaturation in vivo.     

Both DNA and histone fold sequences contribute to archaeal nucleosome stability.

The roles and interdependence of DNA sequence and archaeal histone fold structure in determining archaeal nucleosome stability and positioning have been determined and quantitated. The presence of four tandem copies of TTTAAAGCCG in the polylinker region of pLITMUS28 resulted in a DNA molecule with increased affinity (DeltaDeltaG of approximately 700 cal mol(-1)) for the archaeal histone HMfB relative to the polylinker sequence, and the dominant, quantitative contribution of the helical repeats of the dinucleotide TA to this increased affinity has been established. The rotational and translational positioning of archaeal nucleosomes assembled on the (TTTAAAGCCG)(4) sequence and on DNA molecules selectively incorporated into archaeal nucleosomes by HMfB have been determined. Alternating A/T- and G/C-rich regions were located where the minor and major grooves, respectively, sequentially faced the archaeal nucleosome core, and identical positioning results were obtained using HMfA, a closely related archaeal histone also from Methanothermus fervidus. However, HMfA did not have similarly high affinities for the HMfB-selected DNA molecules, and domain-swap experiments have shown that this difference in affinity is determined by residue differences in the C-terminal region of alpha-helix 3 of the histone fold, a region that is not expected to directly interact with DNA. Rather this region is thought to participate in forming the histone dimer:dimer interface at the center of an archaeal nucleosome histone tetramer core. If differences in this interface do result in archaeal histone cores with different sequence preferences, then the assembly of alternative archaeal nucleosome tetramer cores could provide an unanticipated and novel structural mechanism to regulate gene expression.     

Large scale preparation of positively supercoiled DNA using the archaeal histone HMf.

A technique to prepare relatively large quantities (>/=100 microg) of highly positively supercoiled DNA is reported. This uses a recombinant archaeal histone (rHMfB) to introduce toroidal supercoils, and an inexpensive chicken blood extract to relax unrestrained superhelical tension. Preparation of positively supercoiled pUC19 DNA molecules, >50% of which have linking number changes ranging from+8 to+17, is demonstrated. Advantages include the high degree of positive supercoiling that can be achieved, control over the extent of supercoiling, easy production of relatively large quantities of supercoiled DNA, and low cost.     

DNA binding and nuclease protection by the HMf histones from the hyperthermophilic archaeon Methanothermus fervidus

The DNA-binding and nuclease-protection properties of the HMf histones from the hyperthermophilic archaeon Methanothermus fervidus have been shown to be consistent with the formation of nucleosome-like structures (NLS). These proteins bind to DNA molecules as short as 20 bp and form complexes that protect DNA fragments from micrococcal nuclease (MNase) digestion that are 30 bp, ∼ 60 bp and multiples of ∼ 60 bp in length. The sequences of 49 of the ∼ 60-bp DNA fragments protected from MNase digestion by HMfA have been determined and their intrinsic curvatures calculated. A circular permutation gel mobility-shift assay was used to determine directly the curvatures for five of these sequences. HMfA bound to intrinsically curved and noncurved DNAs, but exhibited a slight preference for the model curved DNA in binding competitions with a model noncurved DNA. The results obtained are consistent with the concept that the archaeal NLS is analogous, and possibly homologous, to the central core of the eukaryal nucleosome formed by a histone (H3 + H4)2 tetramer. Received: August 11, 1996 / Accepted: November 12, 1996     

Rice HMGB1 protein recognizes DNA structures and bends DNA efficiently

We analyzed the DNA-binding and DNA-bending properties of recombinant HMGB1 proteins based on a rice HMGB1 cDNA. Electrophoretic mobility shift assay demonstrated that rice HMGB1 can bind synthetic four-way junction (4H) DNA and DNA minicircles efficiently but the binding to 4H can be completed out by HMGA and histone H1. Conformational changes were detected by circular dichroism analysis with 4H DNA bound to various concentrations of HMGB1 or its truncated forms. T4 ligase-mediated circularization assays with short DNA fragments of 123 bp showed that the protein is capable of increasing DNA flexibility. The 123-bp DNA formed closed circular monomers efficiently in its presence, similar to that in an earlier study on maize HMG. Additionally, our results show for the first time that the basic N-terminal domain enhances the affinity of the plant HMGB1 protein for 4H DNA, while the acidic C-terminal domain has the converse effects.