Protein Sliding along DNA: Dynamics and Structural Characterization

Efficient search of DNA by proteins is fundamental to the control of cellular regulatory processes. It is currently believed that protein sliding, hopping, and transfer between adjacent DNA segments, during which the protein nonspecifically interacts with DNA, are central to the speed of their specific recognition. In this study, we focused on the structural and dynamic features of proteins when they scan the DNA. Using a simple computational model that represents protein-DNA interactions by electrostatic forces, we identified that the protein makes use of identical binding interfaces for both nonspecific and specific DNA interactions. Accordingly, in its one-dimensional diffusion along the DNA, the protein is bound at the major groove and performs a helical motion, which is stochastic and driven by thermal diffusion. A microscopic structural insight into sliding from our model, which is governed by electrostatic forces, corroborates previous experimental studies suggesting that the active site of some regulatory proteins continually faces the interior of the DNA groove while sliding along sugar-phosphate rails. The diffusion coefficient of spiral motion along the major groove of the DNA is not affected by salt concentration, but the efficiency of the search can be significantly enhanced by increasing salt concentration due to a larger number of hopping events. We found that the most efficient search comprises ∼ 20% sliding along the DNA and ∼ 80% hopping and three-dimensional diffusion. The presented model that captures various experimental features of facilitated diffusion has the potency to address other questions regarding the nature of DNA search, such as the sliding characteristics of oligomeric and multidomain DNA-binding proteins that are ubiquitous in the cell.     

Diffusion of Human Replication Protein A along Single-Stranded DNA

Replication protein A (RPA) is a eukaryotic single-stranded DNA (ssDNA) binding protein that plays critical roles in most aspects of genome maintenance, including replication, recombination and repair. RPA binds ssDNA with high affinity, destabilizes DNA secondary structure and facilitates binding of other proteins to ssDNA. However, RPA must be removed from or redistributed along ssDNA during these processes. To probe the dynamics of RPA–DNA interactions, we combined ensemble and single-molecule fluorescence approaches to examine human RPA (hRPA) diffusion along ssDNA and find that an hRPA heterotrimer can diffuse rapidly along ssDNA. Diffusion of hRPA is functional in that it provides the mechanism by which hRPA can transiently disrupt DNA hairpins by diffusing in from ssDNA regions adjacent to the DNA hairpin. hRPA diffusion was also monitored by the fluctuations in fluorescence intensity of a Cy3 fluorophore attached to the end of ssDNA. Using a novel method to calibrate the Cy3 fluorescence intensity as a function of hRPA position on the ssDNA, we estimate a one-dimensional diffusion coefficient of hRPA on ssDNA of D₁ ~ 5000 nt² s^− 1 at 37 °C. Diffusion of hRPA while bound to ssDNA enables it to be readily repositioned to allow other proteins access to ssDNA.     

Substitution Rate Heterogeneity and the Male Mutation Bias

Germline mutation rates have been found to be higher in males than in females in many organisms, a likely consequence of cell division being more frequent in spermatogenesis than in oogenesis. If the majority of mutations are due to DNA replication error, the male-to-female mutation rate ratio (α_m) is expected to be similar to the ratio of the number of germ line cell divisions in males and females (c), an assumption that can be tested with proper estimates of α_m and c. α_m is usually estimated by comparing substitution rates in putatively neutral sequences on the sex chromosomes. However, substantial regional variation in substitution rates across chromosomes may bias estimates of α_m based on the substitution rates of short sequences. To investigate regional substitution rate variation, we estimated sequence divergence in 16 gametologous introns located on the Z and W chromosomes of five bird species of the order Galliformes. Intron ends and potentially conserved blocks were excluded to reduce the effect of using sequences subject to negative selection. We found significant substitution rate variation within Z chromosome (G₁₅ = 37.6, p = 0.0010) as well as within W chromosome introns (G₁₅ = 44.0, p = 0.0001). This heterogeneity also affected the estimates of α_m, which varied significantly, from 1.53 to 3.51, among the introns (ANOVA: F_13,14 =2.68, p = 0.04). Our results suggest the importance of using extensive data sets from several genomic regions to avoid the effects of regional mutation rate variation and to ensure accurate estimates of α_m. Electronic Supplementary Material Electronic Supplementary material is available for this article at and accessible for authorised users. [Reviewing Editor: Mr. Martin Kreitman] Nick G.C. Smith Deceased     

Two Dimensional Yau-Hausdorff Distance with Applications on Comparison of DNA and Protein Sequences

Comparing DNA or protein sequences plays an important role in the functional analysis of genomes. Despite many methods available for sequences comparison, few methods retain the information content of sequences. We propose a new approach, the Yau-Hausdorff method, which considers all translations and rotations when seeking the best match of graphical curves of DNA or protein sequences. The complexity of this method is lower than that of any other two dimensional minimum Hausdorff algorithm. The Yau-Hausdorff method can be used for measuring the similarity of DNA sequences based on two important tools: the Yau-Hausdorff distance and graphical representation of DNA sequences. The graphical representations of DNA sequences conserve all sequence information and the Yau-Hausdorff distance is mathematically proved as a true metric. Therefore, the proposed distance can preciously measure the similarity of DNA sequences. The phylogenetic analyses of DNA sequences by the Yau-Hausdorff distance show the accuracy and stability of our approach in similarity comparison of DNA or protein sequences. This study demonstrates that Yau-Hausdorff distance is a natural metric for DNA and protein sequences with high level of stability. The approach can be also applied to similarity analysis of protein sequences by graphic representations, as well as general two dimensional shape matching.     

Fine Mapping of Satellite DNA Sequences along the Y Chromosome of Drosophila Melanogaster: Relationships between Satellite Sequences and Fertility Factors

The entirely heterochromatic Y chromosome of Drosophila melanogaster contains a series of simple sequence satellite DNAs which together account for about 80% of its length. Molecular cloning of the three simple sequence satellite DNAs of D. melanogaster (1.672, 1.686 and 1.705 g/ml) revealed that each satellite comprises several distinct repeat sequences. Together 11 related sequences were identified and 9 of them were shown to be located on the Y chromosome. In the present study we have finely mapped 8 of these sequences along the Y by in situ hybridization on mitotic chromosome preparations. The hybridization experiments were performed on a series of cytologically determined rearrangements involving the Y chromosome. The breakpoints of these rearrangements provided an array of landmarks along the Y which have been used to localize each sequence on the various heterochromatic blocks defined by Hoechst and N-banding techniques. The results of this analysis indicate a good correlation between the N-banded regions and 1.705 repeats and between the Hoechst-bright regions and the 1.672 repeats. However, the molecular basis for banding does not appear to depend exclusively on DNA content, since heterochromatic blocks showing identical banding patterns often contain different combinations of satellite repeats. The distribution of satellite repeats has also been analyzed with respect to the male fertility factors of the Y chromosome. Both loop-forming (kl-5, kl-3 and ks-1) and non-loop-forming (kl-2 and ks-2) fertility genes contain substantial amounts of satellite DNAs. Moreover, each fertility region is characterized by a specific combination of satellite sequences rather than by an homogeneous array of a single type of repeat.(ABSTRACT TRUNCATED AT 250 WORDS)     

Rates of Nucleotide Substitution in Angiosperm Mitochondrial DNA Sequences and Dates of Divergence Between Brassica and Other Angiosperm Lineages

We obtained 16 nucleotide sequences (∼1400 bp each) of the first intron of the mitochondrial (mt) gene for NADH subunit 4 (nad4) from 10 species of Brassicaceae. Using these new sequences and five published sequences from GenBank, we constructed a phylogenetic tree of the Brassicaceae species under study and showed that the rate of nucleotide substitution in the first intron of nad4 is very low, about 0.16–0.23 × 10⁻⁹ substitution per site per year, which is about half of the silent rate in exons of nad4. The ratios of substitution rates in this intron, ITS, and IGS are approximately 1:23:73, where ITS is the nuclear intergenic spacer between 18S and 25S rRNA genes and IGS is the intergenic spacer of 5S rRNA genes. A segment (335 bp) in the first intron of nad4 in Brassicaceae species that is absent in wheat was considered as a nonfunctional sequence and used to estimate the neutral rate (the rate of mutation) in mtDNA to be 0.5–0.7 × 10⁻⁹ substitution per site per year, which is about three times higher than the substitution rate in the rest of the first intron of nad4. We estimated that the dates of divergence are 170–235 million years (Myr) for the monocot–dicot split, 112–156 Myr for the Brassicaceae–Lettuce split, 14.5–20.4 Myr for the Brassica–Arabidopsis split, and 14.5–20.4 Myr for the Arabidopsis–Arabideae split. Received: 14 July 1998 / Accepted: 1 October 1998     

Detecting Site-Specific Biochemical Constraints Through Substitution Mapping

The neutral theory of molecular evolution states that most mutations are deleterious or neutral. It results that the evolutionary rate of a given position in an alignment is a function of the level of constraint acting on this position. Inferring evolutionary rates from a set of aligned sequences is hence a powerful method to detect functionally and/or structurally important positions in a protein. Some positions, however, may be constrained while having a high substitution rate, providing these substitutions do not affect the biochemical property under constraint. Here, I introduce a new evolutionary rate measure accounting for the evolution of specific biochemical properties (e.g., volume, polarity, and charge). I then present a new statistical method based on the comparison of two rate measures: a site is said to be constrained for property X if it shows an unexpectedly high conservation of X knowing its total evolutionary rate. Compared to single-rate methods, the two-rate method offers several advantages: it (i) allows assessment of the significance of the constraint, (ii) provides information on the type of constraint acting on each position, and (iii) detects positions that are not proposed by previous methods. I apply this method to a 200-sequence data set of triosephosphate isomerase and report significant cases of positions constrained for polarity, volume, or charge. The three-dimensional localization of these positions shows that they are of potential interest to the molecular evolutionist and to the biochemist.     

Latent Periodicity of Protein Sequences

This article is in the area of protein sequence investigation. It studies protein sequence periodicity. The notion of latent periodicity is introduced. A mathematical method for searching for latent periodicity in protein sequences is developed. Implementation of the method developed for known cases of perfect and imperfect periodicity is demonstrated. Latent periodicity of many protein sequences from the SWISS-PROT data bank is revealed by the method and examples of latent periodicity of amino acid sequences are demonstrated for: the translation initiation factor EIF-2B (epsilon subunit) of Saccharomyces cerevisiae from the E2BE_YEAST sequence; the E.coli ferrienterochelin receptor from the FEPA_ECOLI sequence; the lysozyme of Bacteriophage SF6 from the LY_BPSF6 sequence; lipoamide dehydrogenase of Azotobacter vinelandii from the DLDH_AZOVI sequence. These protein sequences have latent periods equal to six, two, seven and 19 amino acids, respectively. We propose that a possible purpose of the amino acid sequence latent periodicity is to determine certain protein structures.     

Protein displacement by an assembly of helicase molecules aligned along single-stranded DNA

Helicases are molecular motors that unwind double-stranded DNA or RNA. In addition to unwinding nucleic acids, an important function of these enzymes seems to be the disruption of protein-nucleic acid interactions. Bacteriophage T4 Dda helicase can displace proteins bound to DNA, including streptavidin bound to biotinylated oligonucleotides. We investigated the mechanism of streptavidin displacement by varying the length of the oligonucleotide substrate. We found that a monomeric form of Dda catalyzed streptavidin displacement; however, the activity increased when multiple helicase molecules bound to the biotinylated oligonucleotide. The activity does not result from cooperative binding of Dda to the oligonucleotide. Rather, the increase in activity is a consequence of the directional bias in translocation of individual helicase monomers. Such a bias leads to protein-protein interactions when the lead monomer stalls owing to the presence of the streptavidin block.     

Estimation of Reversible Substitution Matrices from Multiple Pairs of Sequences

We present a method for estimating the most general reversible substitution matrix corresponding to a given collection of pairwise aligned DNA sequences. This matrix can then be used to calculate evolutionary distances between pairs of sequences in the collection. If only two sequences are considered, our method is equivalent to that of Lanave et al. (1984). The main novelty of our approach is in combining data from different sequence pairs. We describe a weighting method for pairs of taxa related by a known tree that results in uniform weights for all branches. Our method for estimating the rate matrix results in fast execution times, even on large data sets, and does not require knowledge of the phylogenetic relationships among sequences. In a test case on a primate pseudogene, the matrix we arrived at resembles one obtained using maximum likelihood, and the resulting distance measure is shown to have better linearity than is obtained in a less general model.     

Recombination as a Point Process along Sequences

Histories of sequences in the coalescent model with recombination can be simulated using an algorithm that takes as input a sample of extant sequences. The algorithm traces the history of the sequences going back in time, encountering recombinations and coalescence (duplications) until the ancestral material is located on one sequence for homologous positions in the present sequences. Here an alternative algorithm is formulated not as going back in time and operating on sequences, but by moving spatially along the sequences, updating the history of the sequences as recombination points are encountered. This algorithm focuses on spatial aspects of the coalescent with recombination rather than on temporal aspects as is the case of familiar algorithms. Mathematical results related to spatial aspects of the coalescent with recombination are derived.     

Correlations in Protein Sequences and Property Codes

Correlation functions in large sets of non-homologous protein sequences are analysed. Finite size corrections are applied and fluctuations are estimated. As symbol sequences have to be mapped to sequences of numbers to calculate correlation functions, several property codes are tested as such mappings. We found hydrophobicity autocorrelation functions to be strongly oscillating. Another strong signal is the monotonously decaying α-helix propensity autocorrelation function. Furthermore, we detected signals corresponding to an alteration of positively and negatively charged residues at a distance of 3–4 amino acids.To look beyond the property codes gained by the methods of physical chemistry, mappings yielding a strong correlation signal are sought for using a Monte Carlo simulation. The mappings leading to strong signals are found to be related to hydrophobicity of α-helix propensity. A cluster analysis of the top scoring mappings leads to two novel property codes. These two property codes are gained from sequence data only. They turn out to be similar to known property codes for hydrophobicity or polarity.     

Heterogeneity of pumpkin ribosomal DNA

The ribosomal DNA (rDNA) of Cucurbita pepo L. has been found to consist of tandemly arrayed repeat units, most of which are 10 kilobases in length. Thirty-six repeat units, cloned into the HindIII site of pACYC 177, fall into seven classes which differ from each other in length and/or nucleotide sequence. Most of the heterogeneity occurs in noncoding portions of the repeat unit although there is some nucleotide sequence variation in the coding portion as well. Heterogeneity of base modification was observed in genomic rDNA of which two examples are: (a) all of the repeat units have three BamHI sites, one of which is unavailable for restriction in about half of the units and (b) all of the CCGG sites except one are methylated at the internal cytidine in many of the units; a second site is unmethylated in some of the units and in a very few units a third site remains unmethylated.     

CodonTest: Modeling Amino Acid Substitution Preferences in Coding Sequences

Codon models of evolution have facilitated the interpretation of selective forces operating on genomes. These models, however, assume a single rate of non-synonymous substitution irrespective of the nature of amino acids being exchanged. Recent developments have shown that models which allow for amino acid pairs to have independent rates of substitution offer improved fit over single rate models. However, these approaches have been limited by the necessity for large alignments in their estimation. An alternative approach is to assume that substitution rates between amino acid pairs can be subdivided into rate classes, dependent on the information content of the alignment. However, given the combinatorially large number of such models, an efficient model search strategy is needed. Here we develop a Genetic Algorithm (GA) method for the estimation of such models. A GA is used to assign amino acid substitution pairs to a series of rate classes, where is estimated from the alignment. Other parameters of the phylogenetic Markov model, including substitution rates, character frequencies and branch lengths are estimated using standard maximum likelihood optimization procedures. We apply the GA to empirical alignments and show improved model fit over existing models of codon evolution. Our results suggest that current models are poor approximations of protein evolution and thus gene and organism specific multi-rate models that incorporate amino acid substitution biases are preferred. We further anticipate that the clustering of amino acid substitution rates into classes will be biologically informative, such that genes with similar functions exhibit similar clustering, and hence this clustering will be useful for the evolutionary fingerprinting of genes.     

On Graphical and Numerical Representation of Protein Sequences

Abstract

A new approach using a 3-D Cartesian coordinate system to represent protein sequences has been derived. By the 3-D Graphical representation we make a comparison of sequences belonging to nine different proteins.