Molecular spectroscopy and dynamics of intrinsically fluorescent proteins: coral red (dsRed) and yellow (Citrine)

Gene expression of intrinsically fluorescent proteins in biological systems offers new noninvasive windows into cellular function, but optimization of these probes relies on understanding their molecular spectroscopy, dynamics, and structure. Here, the photophysics of red fluorescent protein (dsRed) from discosoma (coral), providing desired longer emission/absorption wavelengths, and an improved yellow fluorescent protein mutant (Citrine) (S65G/V68L/Q69 M/S72A/T203Y) for significant comparison, are characterized by using fluorescence correlation spectroscopy and time-correlated single-photon counting. dsRed fluorescence decays as a single exponential with a 3.65 +/- 0.07-ns time constant, indicating a single emitting state/species independent of pH 4.4-9.0, in contrast with Citrine. However, laser excitation drives reversible fluorescence flicker at 10(3)-10(4) Hz between dark and bright states with a constant partition fraction f(1) = 0.42 +/- 0.06 and quantum yield of approximately 3 x 10(-3). Unlike Citrine (pKa approximately 5.7), pH-dependent proton binding is negligible (pH 3. 9-11) in dsRed. Time-resolved anisotropy of dsRed reveals rapid depolarization (211 +/- 6 ps) plus slow rotational motion (53 +/- 8 ns), in contrast with a single rotational time (16 +/- 2 ns) for Citrine. The molecular dimensions, calculated from rotational and translational diffusion, indicate that dsRed is hydrodynamically 3.8 +/- 0.4 times larger than predicted for a monomer, which suggests an oligomer (possibly a tetramer) configuration even at approximately 10(-9) M. The fast depolarization is attributed to intraoligomer energy transfer between mobile nonparallel chromophores with the initial anisotropy implying a 24 +/- 3 degrees depolarization angle. Large two-photon excitation cross sections ( approximately 100 GM at 990 nm for dsRed and approximately 50 GM at 970 nm for Citrine), advantageous for two-photon-fluorescence imaging in cells, are measured.

Alignment of (dA).(dT) homopolymer tracts in gene flanking sequences suggests nucleosomal periodicity in D. discoideum DNA

It has been shown that the frequency versus size distribution of A and T overlapping and non-overlapping homopolymer tracts of N > 5 in D. discoideum gene flanking and intron regions are significantly greater than in coding regions(1). In the present report, we demonstrate, that a spatial periodicity exists in long A and T tracts (N > 10) in long flanking sequences by scored alignments of those tracts (N > 10) with the nucleosomal repeat. A tract spacing was found at 185-190 bp that corresponds to a maximum alignment score. This is exactly the average spacing of D. discoideum nucleosomes determined experimentally. A majority of A and T tracts in flanking sequences are often spaced by short DNA stretches and the total length of adjacent A and T tracts plus the interrupting short DNA stretch corresponds closely to the average experimentally measured nucleosomal linker DNA size in D. discoideum-42 bp. These data suggest a model which has A and T runs of N > 10 bp in flanking DNA of D. discoideum organized in a regular phase with nonhomopolymer sequences along the DNA. This model has functional implications for A and T tracts, suggesting that they are found in nucleosomal linker DNA regions of chromatin during some necessary portion(s) of the life of the cell.

Wide variations in neighbor-dependent substitution rates

The pattern of 20,200 point substitutions in the 16 unique neighbor-pair environments has been determined from aligned gene/pseudogene sequences in the current database of human DNA sequences. Substitution rates, representing averages over those for different regions of the genome, are distributed over a 60-fold range with strong biases in particular neighbor-pair environments. The rates for substitutions involving the CG doublet are the most rapid overall, where changes of the C.G pair vary over a tenfold range depending on the type of substitution and the 5' neighbor-pair. In general, the rates are fastest in alternating purine-pyrimidine sequences and slowest in purine.pyrimidine tracts, suggesting that the frequencies of one or both key molecular misadventures that can occur during replication, dNTP misinsertion and transient misalignment, may be associated with structural alternations and flexibility of the backbone. By contrast, purine.pyrimidine tracts are less flexible, less prone to substitution, and therefore their proportions accumulate in sequences over time. Characteristic biases of the content and arrangement of oligonucleotide strings or tuples in all sequence elements, but particularly in non-coding regions, appear to be due to the pattern of different neighbor-dependent substitution rates. Computer simulations of numerous replicative cycles have been carried out with substitutions occurring on the same schedule found in this study for pseudogenes. Statistical analyses of tuple frequencies at periodic intervals during the simulation experiment indicate that sequences slowly change in lexical complexity toward a quasi-equilibrium state that corresponds to that for introns.

Characteristics of the large (dA).(dT) homopolymer tracts in D. discoideum gene flanking and intron sequences

D. discoideum, the slime mold, is one of the most AT rich eukaryotic genomes known. In this paper we examine this organism's database for overlapping N-tuples of high frequency and find A and T tracts possess among the highest frequencies in flanking sequences but not in coding sequences. We examined both overlapping and non-overlapping frequencies of the A, T, G and C homopolymer tracts of 2 < N < 6. Overlapping (dG).(dC) and (dA).(dT) tracts occur at greater frequencies than expected, based on random occurrence. Long (dA).(dT) tracts of N > 10 occur at well above expected frequencies in flanking and intron regions, while (dG).(dC) tracts above N = 5 are rarely found. Some of the implications of these findings for tract origins in slip-strand replication and for chromatin structure are discussed.

The influence of nearest neighbors on the rate and pattern of spontaneous point mutations

The numbers and local sequence environments of the two types of substitution mutation plus additions and deletions have been obtained directly in this study from differences between a large number of extant primate gene and pseudogene sequences. A total of 3786 mutations were scored in regions where similarities between pseudogene and corresponding gene sequences is greater than or equal to 85%, comprising approximately 30% of the pseudogene database of 80,584 bp. The pattern of mutations obtained in this fashion is almost identical to that obtained by Li et al. (1984) using a slightly different, more direct approach and with a smaller database. When mutations were scored, the neighbor pairs on the 5' and 3' sides were also noted, leading to a large 16 x 12 matrix of transitions and transversions. Biases of varying magnitude are found in the rates of substitution of the same base pair in different local sequence environments. The overall order for the effect of the 5' neighbor on the rates of substitution mutation of a pyrimidine is A greater than C much greater than T greater than G, and G greater than A greater than T greater than C for the 3' neighbor; where these results represent the average of substitution rates for the complement purine with complement neighbors of bases ordered above. The order for the 3' neighbor is essentially the same for the two transitions and most of the four transversions as well; however, the order for the 5' neighbor is more variable. The overall rate for the C.G----T.A transition is not unusual, however the presence of a 3' neighboring G.C pair boosts the rate substantially, presumably due to specific cytosine methylation of the CG doublet in primate DNAs. The rate of the T.A----C.G transition is also well above average when the 3' neighbor is an A.T, and to a lesser extent a G.C, pair. The latter bias is typical in that it reflects the association of alternating pyrimidine-purine sequences with increasing mutation rates. The substitution of the pyrimidine in a 5'purine-pyrimidine-purine3' sequence generally occurs much faster than in a pyrimidine tract and points to the local conformation as a major determining factor of the substitution rate. An apparent inverse relationship is found between starting and product doublet frequencies of base pairs undergoing mutations with specific 3' neighbors, indicating that differences in intrinsic substitution rates of base pairs with specific neighbors are a key factor in producing the familiar biases of nearest-neighbor frequencies.