Associations between nucleosome phasing, sequence asymmetry, and tissue-specific expression in a set of inbred Medaka species.

BACKGROUND: Transcription start sites (TSSs) with pronounced and phased nucleosome arrays downstream and nucleosome-depleted regions upstream of TSSs are observed in various species. RESULTS: We have characterized sequence variation and expression properties of this set of TSSs (which we call "Nucleocyclic TSSs") using germline and somatic cells of three medaka (Oryzias latipes) inbred isolates from different locations. We found nucleocyclic TSSs in medaka to be associated with higher gene expression and characterized by a clear boundary in sequence composition with potentially-nucleosome-destabilizing A/T-enrichment upstream (p < 10(-60)) and nucleosome- accommodating C/G-enrichment downstream (p < 10(-40)) that was highly conserved from an ancestor. A substantial genetic distance between the strains facilitated the in-depth analysis of patterns of fixed mutations, revealing a localization-specific equilibrium between the rates of distinct mutation categories that would serve to maintain the conserved sequence anisotropy around TSSs. Downstream of nucleocyclic TSSs, C to T, T to C, and other mutation rates on the sense strand increased around first nucleosome dyads and decreased around first linkers, which contrasted with genomewide mutational patterns around nucleosomes (p < 5 %). C to T rates are higher than G to A rates around nucleosome associated with germline nucleocyclic TSS sites (p < 5 %), potentially due to the asymmetric effect of transcription-coupled repair. CONCLUSIONS: Our results demonstrate an atypical evolutionary process surrounding nucleocyclic TSSs.

Chromatin-associated periodicity in genetic variation downstream of transcriptional start sites.

Might DNA sequence variation reflect germline genetic activity and underlying chromatin structure? We investigated this question using medaka (Japanese killifish, Oryzias latipes), by comparing the genomic sequences of two strains (Hd-rR and HNI) and by mapping approximately 37.3 million nucleosome cores from Hd-rR blastulae and 11,654 representative transcription start sites from six embryonic stages. We observed a distinctive approximately 200-base pair (bp) periodic pattern of genetic variation downstream of transcription start sites; the rate of insertions and deletions longer than 1 bp peaked at positions of approximately +200, +400, and +600 bp, whereas the point mutation rate showed corresponding valleys. This approximately 200-bp periodicity was correlated with the chromatin structure, with nucleosome occupancy minimized at positions 0, +200, +400, and +600 bp. These data exemplify the potential for genetic activity (transcription) and chromatin structure to contribute to molding the DNA sequence on an evolutionary time scale.

Experimental characterization of the folding kinetics of the notch ankyrin domain.

Proteins constructed from linear arrays of tandem repeats provide a simplified architecture for understanding protein folding. Here, we examine the folding kinetics of the ankyrin repeat domain from the Drosophila Notch receptor, which consists of six folded ankyrin modules and a seventh partly disordered N-terminal ankyrin repeat sequence. Both the refolding and unfolding kinetics are best described as a sum of two exponential phases. The slow, minor refolding phase is limited by prolyl isomerization in the denatured state (D). The minor unfolding phase, which appears as a lag during fluorescence-detected unfolding, is consistent with an on-pathway intermediate (I). This intermediate, although not directly detected during refolding, is shown to be populated by interrupted refolding experiments. When plotted against urea, the rate constants for the major unfolding and refolding phases define a single non-linear v-shaped chevron, as does the minor unfolding phase. These two chevrons, along with unfolding amplitudes, are well-fitted by a sequential three-state model, which yields rate constants for the individual steps in folding and unfolding. Based on these fitted parameters, the D to I step is rate-limiting, and closely matches the major observed refolding phase at low denaturant concentrations. I appears to be midway between N and D in folding free energy and denaturant sensitivity, but has Trp fluorescence properties close to N. Although the Notch ankyrin domain has a simple architecture, folding is slow, with the limiting refolding rate constant as much as seven orders of magnitude smaller than expected from topological predictions.

An experimentally determined protein folding energy landscape.

Energy landscapes have been used to conceptually describe and model protein folding but have been difficult to measure experimentally, in large part because of the myriad of partly folded protein conformations that cannot be isolated and thermodynamically characterized. Here we experimentally determine a detailed energy landscape for protein folding. We generated a series of overlapping constructs containing subsets of the seven ankyrin repeats of the Drosophila Notch receptor, a protein domain whose linear arrangement of modular structural units can be fragmented without disrupting structure. To a good approximation, stabilities of each construct can be described as a sum of energy terms associated with each repeat. The magnitude of each energy term indicates that each repeat is intrinsically unstable but is strongly stabilized by interactions with its nearest neighbors. These linear energy terms define an equilibrium free energy landscape, which shows an early free energy barrier and suggests preferred low-energy routes for folding.

Measuring the stability of partly folded proteins using TMAO.

Standard methods for measuring free energy of protein unfolding by chemical denaturation require complete folding at low concentrations of denaturant so that a native baseline can be observed. Alternatively, proteins that are completely unfolded in the absence of denaturant can be folded by addition of the osmolyte trimethylamine N-oxide (TMAO), and the unfolding free energy can then be calculated through analysis of the refolding transition. However, neither chemical denaturation nor osmolyte-induced refolding alone is sufficient to yield accurate thermodynamic unfolding parameters for partly folded proteins, because neither method produces both native and denatured baselines in a single transition. Here we combine urea denaturation and TMAO stabilization as a means to bring about baseline-resolved structural transitions in partly folded proteins. For Barnase and the Notch ankyrin domain, which both show two-state equilibrium unfolding, we found that DeltaG degrees for unfolding depends linearly on TMAO concentration, and that the sensitivity of DeltaG degrees to urea (the m-value) is TMAO independent. This second observation confirms that urea and TMAO exert independent effects on stability over the range of cosolvent concentrations required to bring about baseline-resolved structural transitions. Thermodynamic parameters calculated using a global fit that assumes additive, linear dependence of DeltaG degrees on each cosolvent are similar to those obtained by standard urea-induced unfolding in the absence of TMAO. Finally, we demonstrate the applicability of this method to measurement of the free energy of unfolding of a partly folded protein, a fragment of the full-length Notch ankyrin domain.