Structural mapping of nucleotide binding sites on chloroplast coupling factor

Fluorescence resonance energy transfer was used to measure the distances between three nucleotide binding sites on solubilized chloroplast coupling factor from spinach and between each nucleotide site and two tyrosine residues which are important for catalytic activity. The nucleotide energy donor was 1,N6-ethenoadenosine di- or triphosphate, and the nucleotide energy acceptor was 2'(3')-(trinitrophenyl)adenosine diphosphate. The tyrosine residues were specifically labeled with 7-chloro-4-nitro-2,1,3-benzoxadiazole, which served as an energy acceptor. The results obtained indicate the three nucleotide binding sites form a triangle with sides of 44, 48, and 36 A. (The assumption has been made in calculating these distances that the energy donor and acceptor rotate rapidly relative to the fluorescence lifetime.) Two of the nucleotide sites are approximately equidistant from each of the two tyrosines: one of the nucleotide sites is about 37 A and the other about 41 A from each tyrosine. The third nucleotide site is about 41 A from one of the tyrosines and greater than or equal to 41 A from the other tyrosine.     

New elements in the structure of binding sites of glucocorticoid-receptor complex in hormone-regulated genes

An analysis of the structure of DNA sites responsible for binding to glucocorticoid-receptor complex (GlRC) was carried out. The use of the frequency matrices and of a variant of the perception method made it possible to establish that in the GlRC binding site on both sides of the known conservative nucleotide sequence (nucleus) there were additional conservative elements which seemed to be able to modulate the efficiency of GlRC binding. A criterion is worked out for detecting the potential GlRC binding sites in given sequences. It is based on the simultaneous use of several perceptron matrices. The efficiency of detection of GlRC binding sites by means of the proposed criterion is by an order higher than that performed according to the GlRC binding site consensus (Beato et al. [2]).     

Prediction of human mRNA donor and acceptor sites from the DNA sequence

Artificial neural networks have been applied to the prediction of splice site location in human pre-mRNA. A joint prediction scheme where prediction of transition regions between introns and exons regulates a cutoff level for splice site assignment was able to predict splice site locations with confidence levels far better than previously reported in the literature. The problem of predicting donor and acceptor sites in human genes is hampered by the presence of numerous amounts of false positives: here, the distribution of these false splice sites is examined and linked to a possible scenario for the splicing mechanism in vivo. When the presented method detects 95% of the true donor and acceptor sites, it makes less than 0.1% false donor site assignments and less than 0.4% false acceptor site assignments. For the large data set used in this study, this means that on average there are one and a half false donor sites per true donor site and six false acceptor sites per true acceptor site. With the joint assignment method, more than a fifth of the true donor sites and around one fourth of the true acceptor sites could be detected without accompaniment of any false positive predictions. Highly confident splice sites could not be isolated with a widely used weight matrix method or by separate splice site networks. A complementary relation between the confidence levels of the coding/non-coding and the separate splice site networks was observed, with many weak splice sites having sharp transitions in the coding/non-coding signal and many stronger splice sites having more ill-defined transitions between coding and non-coding.     

One parameter to describe the mechanism of splice sites competition

The choice of a splice site is not only related to its own intrinsic strength, but also is influenced by its flanking competitors. Splice site competition is an important mechanism for splice site prediction, especially, it is a new insight for alternative splice site prediction. In this paper, the position weight matrix scoring function is used to represent splice site strength, and the mechanism of splice site competition is described by only one parameter: scoring function subtraction. While applying on the alternative splice site prediction, based on the only one parameter, 68.22% of donor sites and 70.86% of acceptor sites are correctly classified into alternative and constitutive. The prediction abilities are approximately equal to the recent method which is based on the mechanism of splice site competition. The results reveal that the scoring function subtraction is the best parameter to describe the mechanism of splice sites competition.     

TRII: A Probabilistic Scoring of Drosophila melanogaster Translation Initiation Sites

Relative individual information is a measurement that scores the quality of DNA- and RNA-binding sites for biological machines. The development of analytical approaches to increase the power of this scoring method will improve its utility in evaluating the functions of motifs. In this study, the scoring method was applied to potential translation initiation sites in Drosophila to compute Translation Relative Individual Information (TRII) scores. The weight matrix at the core of the scoring method was optimized based on high-confidence translation initiation sites identified by using a progressive partitioning approach. Comparing the distributions of TRII scores for sites of interest with those for high-confidence translation initiation sites and random sequences provides a new methodology for assessing the quality of translation initiation sites. The optimized weight matrices can also be used to describe the consensus at translation initiation sites, providing a quantitative measure of preferred and avoided nucleotides at each position.     

Mutating protein kinase cAMP-binding sites into cGMP-binding sites. Mechanism of cGMP selectivity.

The cAMP-dependent protein kinase contains two different cAMP-binding sites referred to as the slow and fast sites. Mutation of Ala-334 to a threonine in the slow site of the bovine type I regulatory subunit created a site with marked increase in cGMP affinity without changing cAMP affinity (Shabb, J. B., Ng. L., Corbin, J. D. (1990) J. Biol. Chem. 265, 16031-16034). The corresponding fast site residue (Ala-210) was changed to a threonine by oligonucleotide-directed mutagenesis, and a double mutant containing a threonine in each site was also made. Holoenzymes were formed from native catalytic subunit and each recombinant regulatory subunit. The fast site mutant holoenzyme exhibited an improved cGMP activation constant and an impaired cAMP activation constant. The double mutant cGMP/cAMP selectivity was 200-fold greater than that of wild-type holoenzyme, making it as responsive to cGMP as native cGMP-dependent protein kinase. The increased intrinsic binding energies of mutated sites for cGMP were 2.7-3.0 kcal mol-1, consistent with the presence of an extra hydrogen bond. Cyclic nucleotide analog studies implied that this hydrogen bond was between the threonine hydroxyl and the 2-amino of cGMP. Comparisons of amino acid sequences and cyclic nucleotide specificities suggested that the Ala/Thr difference may also impart cAMP/cGMP binding selectivity to related proteins such as cyclic nucleotide-gated ion channels.     

Evolutionary change of restriction cleavage sites and phylogenetic inference for man and apes

A mathematical theory for the evolutionary change of restriction endonuclease cleavage sites is developed, and the probabilities of various types of restriction-site changes are evaluated. A computer simulation is also conducted to study properties of the evolutionary change of restriction sites. These studies indicate that parsimony methods of constructing phylogenetic trees often make erroneous inferences about evolutionary changes of restriction sites unless the number of nucleotide substitutions per site is less than 0.01 for all branches of the tree. This introduces a systematic error in estimating the number of mutational changes for each branch and, consequently, in constructing phylogenetic trees. Therefore, parsimony methods should be used only in cases where nucleotide sequences are closely related. Reexamination of Ferris et al.'s data on restriction-site differences of mitochondrial DNAs does not support Templeton's conclusions regarding the phylogenetic tree for man and apes and the molecular clock hypothesis. Templeton's claim that Nei and Li's method of estimating the number of nucleotide substitutions per site is seriously affected by parallel losses and loss-gains of restriction sites is also unsupported.      

Analysis of substrate binding in individual active sites of bifunctional human ATIC

Aminoimidazolecarboxamide ribonucleotide formyl transferase (AICARFT): Inosine monophosphate cyclohydrolase (IMPCH, collectively called ATIC) is a bifunctional enzyme that catalyses the penultimate and final steps in the purine de novo biosynthesis pathway. The bifunctional protein is dimeric and each monomer contains two different active sites both of which are capable of binding nucleotide substrates, this means to a potential total of four distinct binding events might be observed. Within this work we used a combination of site-directed and truncation mutants of ATIC to independently investigate the binding at these two sites using calorimetry. A single S10W mutation is sufficient to block the IMPCH active site allowing investigation of the effects of mutation on ligand binding in the AICARFT active site. The majority of nucleotide ligands bind selectively at one of the two active sites with the exception of xanthosine monophosphate, XMP, which, in addition to binding in both AICARFT and IMPCH active sites, shows evidence for cooperative binding with communication between symmetrically-related active sites in the two IMPCH domains. The AICARFT site is capable of independently binding both nucleotide and folate substrates with high affinity however no evidence for positive cooperativity in binding could be detected using the model ligands employed in this study.     

Adenine nucleotide-binding sites on beef heart F1-ATPase. Conditions that affect occupancy of catalytic and noncatalytic sites

Beef heart mitochondrial F1 contains a total of six adenine nucleotide-binding sites including at least two different types of sites. Three "exchangeable" sites exchange rapidly during hydrolysis of MgATP, whereas three "nonexchangeable" sites do not (Cross, R. L. and Nalin, C. M. (1982) J. Biol. Chem. 257, 2874-2881). When F1 that has been stored as a suspension in (NH4)2SO4/ATP/EDTA/sucrose/Tris, pH 8.0, is pelleted, rinsed with (NH4)2SO4, dissolved, and desalted, it retains three bound adenine nucleotides. We find that two of these endogenous nucleotides are bound at nonexchangeable sites and one at an exchangeable site. The vacant nonexchangeable site is highly specific for adenine nucleotide and is rapidly filled by ADP upon addition of ADP or during ATP hydrolysis. ADP bound at this site can be removed by reprecipitating the enzyme with (NH4)2SO4. The single nucleotide retained by desalted F1 at an exchangeable site is displaced during hydrolysis of ATP, GTP, or ITP. The binding of PPi at two sites on the enzyme also promotes its dissociation. Neither procedure affects retention of nucleotide at the nonexchangeable sites. These observations, combined with the finding that PPi is much more easily removed from exchangeable sites than ADP, have led to the development of a procedure for preparing F1 with uniform and well-defined nucleotide site occupancy. This involves sequential exposure to MgATP, PPi, and high concentrations of Pi. Unbound ligand is removed between each step. The resulting enzyme, F1[3,0], has three occupied nonexchangeable sites and three vacant exchangeable sites. Evidence that nonexchangeable and exchangeable sites represent noncatalytic and catalytic sites, respectively, suggest that this form of the enzyme will prove useful in numerous applications, including transient kinetic measurements and affinity labeling of active site residues.     

Distances between the functional sites of the (Ca2+ + Mg2+)-ATPase of sarcoplasmic reticulum

Luminescence energy transfer measurements have been used to determine the distances between the two high affinity Ca2+ binding-transport sites of the (Ca2+ + Mg2+)-ATPase of skeletal muscle sarcoplasmic reticulum. The lanthanide Tb3+ situated at one high affinity Ca2+ site was used as the transfer donor, and acceptors at the other Ca2+ site were the lanthanides Nd3+, Pr3+, Ho3+, or Er3+. Terbium bound to the enzyme was excited directly with a pulsed dye laser. Analysis of the changes in the terbium luminescence lifetime due to the presence of the acceptor indicates that the distance between the Ca2+ sites is 10.7 A. The distance between the Ca2+ sites and the nucleotide-binding catalytic site was determined using Tb3+ at the Ca2+ sites and either trinitrophenyl nucleotides (TNP-N) or fluorescein 5-isothiocyanate (FITC) in the catalytic site as energy acceptors. The R0 values for the Tb-acceptor pairs are approximately 30 and approximately 40 A for TNP-N and FITC, respectively. The distance between Tb3+ at the Ca2+ sites and TNP-ATP at the nucleotide site is approximately 35 A and that between the Ca2+ sites and the FITC labeling site is approximately 47 A. Considerations of the molecular dimensions of the ATPase polypeptide indicate that while the two Ca2+ sites are close to each other, the Ca2+ sites and the nucleotide site are quite remote in the three-dimensional structure of the enzyme.     

Substitution rate variation among sites in hypervariable region 1 of human mitochondrial DNA

More than an order of magnitude difference in substitution rate exists among sites within hypervariable region 1 of the control region of human mitochondrial DNA. A two-rate Poisson mixture and a negative binomial distribution are used to describe the distribution of the inferred number of changes per nucleotide site in this region. When three data sets are pooled, however, the two-rate model cannot explain the data. The negative binomial distribution always fits, suggesting that substitution rates are approximately gamma distributed among sites. Simulations presented here provide support for the use of a biased, yet commonly employed, method of examining rate variation. The use of parsimony in the method to infer the number of changes at each site introduces systematic errors into the analysis. These errors preclude an unbiased quantification of variation in substitution rate but make the method conservative overall. The method can be used to distinguish sites with highly elevated rates, and 29 such sites are identified in hypervariable region 1. Variation does not appear to be clustered within this region. Simulations show that biases in rates of substitution among nucleotides and non-uniform base composition can mimic the effects of variation in rate among sites. However, these factors contribute little to the levels of rate variation observed in hypervariable region 1.     

Detecting cis-regulatory binding sites for cooperatively binding proteins

Several methods are available to predict cis-regulatory modules in DNA based on position weight matrices. However, the performance of these methods generally depends on a number of additional parameters that cannot be derived from sequences and are difficult to estimate because they have no physical meaning. As the best way to detect cis-regulatory modules is the way in which the proteins recognize them, we developed a new scoring method that utilizes the underlying physical binding model. This method requires no additional parameter to account for multiple binding sites; and the only necessary parameters to model homotypic cooperative interactions are the distances between adjacent protein binding sites in basepairs, and the corresponding cooperative binding constants. The heterotypic cooperative binding model requires one more parameter per cooperatively binding protein, which is the concentration multiplied by the partition function of this protein. In a case study on the bacterial ferric uptake regulator, we show that our scoring method for homotypic cooperatively binding proteins significantly outperforms other PWM-based methods where biophysical cooperativity is not taken into account.     

Recognition of functional sites in protein structures

Recognition of regions on the surface of one protein, that are similar to a binding site of another is crucial for the prediction of molecular interactions and for functional classifications. We first describe a novel method, SiteEngine, that assumes no sequence or fold similarities and is able to recognize proteins that have similar binding sites and may perform similar functions. We achieve high efficiency and speed by introducing a low-resolution surface representation via chemically important surface points, by hashing triangles of physico-chemical properties and by application of hierarchical scoring schemes for a thorough exploration of global and local similarities. We proceed to rigorously apply this method to functional site recognition in three possible ways: first, we search a given functional site on a large set of complete protein structures. Second, a potential functional site on a protein of interest is compared with known binding sites, to recognize similar features. Third, a complete protein structure is searched for the presence of an a priori unknown functional site, similar to known sites. Our method is robust and efficient enough to allow computationally demanding applications such as the first and the third. From the biological standpoint, the first application may identify secondary binding sites of drugs that may lead to side-effects. The third application finds new potential sites on the protein that may provide targets for drug design. Each of the three applications may aid in assigning a function and in classification of binding patterns. We highlight the advantages and disadvantages of each type of search, provide examples of large-scale searches of the entire Protein Data Base and make functional predictions.     

PSSM-based prediction of DNA binding sites in proteins

Background  

Detection of DNA-binding sites in proteins is of enormous interest for technologies targeting gene regulation and manipulation. We have previously shown that a residue and its sequence neighbor information can be used to predict DNA-binding candidates in a protein sequence. This sequence-based prediction method is applicable even if no sequence homology with a previously known DNA-binding protein is observed. Here we implement a neural network based algorithm to utilize evolutionary information of amino acid sequences in terms of their position specific scoring matrices (PSSMs) for a better prediction of DNA-binding sites.     

RNA splice junctions of different classes of eukaryotes: sequence statistics and functional implications in gene expression

A systematic analysis of the RNA splice junction sequences of eukaryotic protein coding genes was carried out using the GENBANK databank. Nucleotide frequencies obtained for the highly conserved regions around the splice sites for different categories of organisms closely agree with each other. A striking similarity among the rare splice junctions which do not contain AG at the 3' splice site or GT at the 5' splice site indicates the existence of special mechanisms to recognize them, and that these unique signals may be involved in crucial gene-regulation events and in differentiation. A method was developed to predict potential exons in a bare sequence, using a scoring and ranking scheme based on nucleotide weight tables. This method was used to find a majority of the exons in selected known genes, and also predicted potential new exons which may be used in alternative splicing situations.