A structural census of the current population of protein sequences

We examine the occurrence of the ≈300 known protein folds in different groups of organisms. To do this, we characterize a large fraction of the currently known protein sequences (≈140,000) in structural terms, by matching them to known structures via sequence comparison (or by secondary-structure class prediction for those without structural homologues). Overall, we find that an appreciable fraction of the known folds are present in each of the major groups of organisms (e.g., bacteria and eukaryotes share 156 of 275 folds), and most of the common folds are associated with many families of nonhomologous sequences (i.e., >10 sequence families for each common fold). However, different groups of organisms have characteristically distinct distributions of folds. So, for instance, some of the most common folds in vertebrates, such as globins or zinc fingers, are rare or absent in bacteria. Many of these differences in fold usage are biologically reasonable, such as the folds of metabolic enzymes being common in bacteria and those associated with extracellular transport and communication being common in animals. They also have important implications for database-based methods for fold recognition, suggesting that an unknown sequence from a plant is more likely to have a certain fold (e.g., a TIM barrel) than an unknown sequence from an animal.

Antisense oligonucleotides against rat brain α1E DNA and its atrial homologue decrease T-type calcium current in atrial myocytes

Low voltage-activated, or T-type, calcium currents are important regulators of neuronal and muscle excitability, secretion, and possibly cell growth and differentiation. The gene (or genes) coding for the pore-forming subunit of low voltage-activated channel proteins has not been unequivocally identified. We have used reverse transcription–PCR to identify partial clones from rat atrial myocytes that share high homology with a member of the E class of calcium channel genes. Antisense oligonucleotides targeting one of these partial clones (raE1) specifically block the increase in T-current density that normally results when atrial myocytes are treated with insulin-like growth factor 1 (IGF-1). Antisense oligonucleotides targeting portions of the neuronal rat α1E sequence, which are not part of the clones detected in atrial tissue, also block the IGF-1-induced increase in T-current, suggesting that the high homology to α1E seen in the partial clone may be present in the complete atrial sequence. The basal T-current expressed in these cells is also blocked by antisense oligonucleotides, which is consistent with the notion that IGF-1 up-regulates the same gene that encodes the basal current. These results support the hypothesis that a member of the E class of calcium channel genes encodes a low voltage-activated calcium channel in atrial myocytes.

Asparagine-proline sequence within membrane-spanning segment of SREBP triggers intramembrane cleavage by Site-2 protease

The NH2-terminal domains of membrane-bound sterol regulatory element-binding proteins (SREBPs) are released into the cytosol by regulated intramembrane proteolysis, after which they enter the nucleus to activate genes encoding lipid biosynthetic enzymes. Intramembrane proteolysis is catalyzed by Site-2 protease (S2P), a hydrophobic zinc metalloprotease that cleaves SREBPs at a membrane-embedded leucine-cysteine bond. In the current study, we use domain-swapping methods to localize the residues within the SREBP-2 membrane-spanning segment that are required for cleavage by S2P. The studies reveal a requirement for an asparagine-proline sequence in the middle third of the transmembrane segment. We propose a model in which the asparagine-proline sequence serves as an NH2-terminal cap for a portion of the transmembrane α-helix of SREBP, allowing the remainder of the α-helix to unwind partially to expose the peptide bond for cleavage by S2P.

ACTIVITY: a database on DNA/RNA sites activity adapted to apply sequence-activity relationships from one system to another

ACTIVITY is a database on DNA/RNA site sequences with known activity magnitudes, measurement systems, sequence-activity relationships under fixed experimental conditions and procedures to adapt these relationships from one measurement system to another. This database deposits information on DNA/RNA affinities to proteins and cell nuclear extracts, cutting efficiencies, gene transcription activity, mRNA translation efficiencies, mutability and other biological activities of natural sites occurring within promoters, mRNA leaders, and other regulatory regions in pro- and eukaryotic genomes, their mutant forms and synthetic analogues. Since activity magnitudes are heavily system-dependent, the current version of ACTIVITY is supplemented by three novel sub-databases: (i) SYSTEM, measurement systems; (ii) KNOWLEDGE, sequence-activity relationships under fixed experimental conditions; and (iii) CROSS_TEST, procedures adapting a relationship from one measurement system to another. These databases are useful in molecular biology, pharmacogenetics, metabolic engineering, drug design and biotechnology. The databases can be queried using SRS and are available through the Web, http://wwwmgs.bionet.nsc.ru/systems/Activity/.

Excitation–contraction coupling is unaffected by drastic alteration of the sequence surrounding residues L720–L764 of the α1S II-III loop

The II-III loop of the skeletal muscle dihydropyridine receptor (DHPR) α1S subunit is responsible for bidirectional-signaling interactions with the ryanodine receptor (RyR1): transmitting an orthograde, excitation–contraction (EC) coupling signal to RyR1 and receiving a retrograde, current-enhancing signal from RyR1. Previously, several reports argued for the importance of two distinct regions of the skeletal II-III loop (residues R681–L690 and residues L720–Q765, respectively), claiming for each a key function in DHPR–RyR1 communication. To address whether residues 720–765 of the II-III loop are sufficient to enable skeletal-type (Ca2+ entry-independent) EC coupling and retrograde interaction with RyR1, we constructed a green fluorescent protein (GFP)-tagged chimera (GFP-SkLM) having rabbit skeletal (Sk) DHPR sequence except for a II-III loop (L) from the DHPR of the house fly, Musca domestica (M). The Musca II-III loop (75% dissimilarity to α1S) has no similarity to α1S in the regions R681–L690 and L720–Q765. GFP-SkLM expressed in dysgenic myotubes (which lack endogenous α1S subunits) was unable to restore EC coupling and displayed strongly reduced Ca2+ current densities despite normal surface expression levels and correct triad targeting (colocalization with RyR1). Introducing rabbit α1S residues L720–L764 into the Musca II-III loop of GFP-SkLM (substitution for Musca DHPR residues E724–T755) completely restored bidirectional coupling, indicating its dependence on α1S loop residues 720–764 but its independence from other regions of the loop. Thus, 45 α1S-residues embedded in a very dissimilar background are sufficient to restore bidirectional coupling, indicating that these residues may be a site of a protein–protein interaction required for bidirectional coupling.

DNA sequence insertion and evolutionary variation in gene regulation.

Current evidence on the long-term evolutionary effect of insertion of sequence elements into gene regions is reviewed, restricted to cases where a sequence derived from a past insertion participates in the regulation of expression of a useful gene. Ten such examples in eukaryotes demonstrate that segments of repetitive DNA or mobile elements have been inserted in the past in gene regions, have been preserved, sometimes modified by selection, and now affect control of transcription of the adjacent gene. Included are only examples in which transcription control was modified by the insert. Several cases in which merely transcription initiation occurred in the insert were set aside. Two of the examples involved the long terminal repeats of mammalian endogenous retroviruses. Another two examples were control of transcription by repeated sequence inserts in sea urchin genomes. There are now six published examples in which Alu sequences were inserted long ago into human gene regions, were modified, and now are central in control/enhancement of transcription. The number of published examples of Alu sequences affecting gene control has grown threefold in the last year and is likely to continue growing. Taken together, all of these examples show that the insertion of sequence elements in the genome has been a significant source of regulatory variation in evolution.

Sequence similarity analysis of Escherichia coli proteins: functional and evolutionary implications.

A computer analysis of 2328 protein sequences comprising about 60% of the Escherichia coli gene products was performed using methods for database screening with individual sequences and alignment blocks. A high fraction of E. coli proteins--86%--shows significant sequence similarity to other proteins in current databases; about 70% show conservation at least at the level of distantly related bacteria, and about 40% contain ancient conserved regions (ACRs) shared with eukaryotic or Archaeal proteins. For > 90% of the E. coli proteins, either functional information or sequence similarity, or both, are available. Forty-six percent of the E. coli proteins belong to 299 clusters of paralogs (intraspecies homologs) defined on the basis of pairwise similarity. Another 10% could be included in 70 superclusters using motif detection methods. The majority of the clusters contain only two to four members. In contrast, nearly 25% of all E. coli proteins belong to the four largest superclusters--namely, permeases, ATPases and GTPases with the conserved "Walker-type" motif, helix-turn-helix regulatory proteins, and NAD(FAD)-binding proteins. We conclude that bacterial protein sequences generally are highly conserved in evolution, with about 50% of all ACR-containing protein families represented among the E. coli gene products. With the current sequence databases and methods of their screening, computer analysis yields useful information on the functions and evolutionary relationships of the vast majority of genes in a bacterial genome. Sequence similarity with E. coli proteins allows the prediction of functions for a number of important eukaryotic genes, including several whose products are implicated in human diseases.

Mapping the mouse genome: current status and future prospects.

The mouse is the best model system for the study of mammalian genetics and physiology. Because of the feasibility and importance of studying genetic crosses, the mouse genetic map has received tremendous attention in recent years. It currently contains over 14,000 genetically mapped markers, including 700 mutant loci, 3500 genes, and 6500 simple sequence length polymorphisms (SSLPs). The mutant loci and genes allow insights and correlations concerning physiology and development. The SSLPs provide highly polymorphic anchor points that allow inheritance to be traced in any cross and provide a scaffold for assembling physical maps. Adequate physical mapping resources--notably large-insert yeast artificial chromosome (YAC) libraries--are available to support positional cloning projects based on the genetic map, but a comprehensive physical map is still a few years away. Large-scale sequencing efforts have not yet begun in mouse, but comparative sequence analysis between mouse and human is likely to provide tremendous information about gene structure and regulation.