Orthogonal beta beta motifs in proteins.

A super-secondary structural motif comprising two orthogonally oriented beta-strands connected by short linking segments of less than or equal to 5 residues has been identified from a data set of 65 independent protein crystal structures. Of the 42 examples from 14 proteins, a vast majority have only a single residue as the linking element. Analysis of the conformational angles at the junction reveals that the recently described type VIII beta-turn occurs frequently at the connecting hinge, while the type II beta-turn is also fairly common.     

Recurrent alpha beta loop structures in TIM barrel motifs show a distinct pattern of conserved structural features.

A systematic survey of seven parallel alpha/beta barrel protein domains, based on exhaustive structural comparisons, reveals that a sizable proportion of the alpha beta loops in these proteins--20 out of a total of 49--belong to either one of two loop types previously described by Thornton and co-workers. Six loops are of the alpha beta 1 type, with one residue between the alpha-helix and beta-strand, and 13 are of the alpha beta 3 type, with three residues between the helix and the strand. Protein fragments embedding the identified loops, and termed alpha beta connections since they contain parts of the flanking helix and strand, have been analyzed in detail revealing that each type of connection has a distinct set of conserved structural features. The orientation of the beta-strand relative to the helix and loop portions is different owing to a very localized difference in backbone conformation. In alpha beta 1 connections, the chain enters the beta-strand via a residue adopting an extended conformation, while in alpha beta 3 it does so via a residue in a near alpha-helical conformation. Other conserved structural features include distinct patterns of side chain orientation relative to the beta-sheet surface and of main chain H-bonds in the loop and the beta-strand moieties. Significant differences also occur in packing interactions of conserved hydrophobic residues situated in the last turn of the helix. Yet the alpha-helix surface of both types of connections adopts similar orientations relative to the barrel sheet surface. Our results suggest furthermore that conserved hydrophobic residues along the sequence of the connections, may be correlated more with specific patterns of interactions made with neighboring helices and sheet strands than with helix/strand packing within the connection itself. A number of intriguing observations are also made on the distribution of the identified alpha beta 1 and alpha beta 3 loops within the alpha/beta-barrel motifs. They often occur adjacent to each other; alpha beta 3 loops invariably involve even numbered beta-strands, while alpha beta 1 loops involve preferentially odd beta-strands; all the analyzed proteins contain at least one alpha beta 3 loop in the first half of the eightfold alpha/beta barrel. Possible origins of all these observations, and their relevance to the stability and folding of parallel alpha/beta barrel motifs are discussed.     

Functional analysis of conserved motifs in the mechanosensitive channel homolog MscS-Like2 from Arabidopsis thaliana

The Mechanosensitive channel of Small conductance (MscS) of Escherichia coli has become an excellent model system for the structural, biophysical, and functional study of mechanosensitive ion channels. MscS, a complex channel with multiple states, contributes to protection against lysis upon osmotic downshock. MscS homologs are widely and abundantly dispersed among the bacterial and plant lineages, but are not found in animals. Investigation into the eukaryotic branch of the MscS family is in the beginning stages, and it remains unclear how much MscS homologs from eukaryotes resemble E. coli MscS with respect to structure, function, and regulation. Here we test the effect of mutating three conserved motifs on the function of MscS-Like (MSL)2, a MscS homolog localized to the plastids of Arabidopsis thaliana. We show that 1) a motif at the top of the cytoplasmic domain, referred to here as the PN(X)(9)N motif, is essential for MSL2 function and for its proper intraplastidic localization; 2) substituting polar residues for two large hydrophobic residues located in the predicted pore-lining transmembrane helix of MSL2 produces a likely gain-of-function allele, as previously shown for MscS; and 3) mis-expression of this allele causes severe defects in leaf growth, loss of chloroplast integrity, and abnormal starch accumulation. Thus, two of the three conserved motifs we analyzed are critical for MSL2 function, consistent with the conservation of structure and function between MscS family members in bacteria and plants. These results underscore the importance of plastidic mechanosensitive channels in the maintenance of normal plastid and leaf morphology.     

Vitellogenin motifs conserved in nematodes and vertebrates

Summary Caenorhabditis elegans vitellogenins are encoded by a family of six genes, one of which,vit-5, has been previously sequenced and shown to be surprisingly closely related to the vertebrate vitellogenin genes. Here we report an alignment of the amino acid sequences of vitellogenins from frog and chicken with those from threeC. elegans genes:vit-5 and two newly sequenced genes,vit-2 andvit-6. The four introns ofvit-6 are all in different places from the four introns ofvit-5, but three of these eight positions are identical or close to intron locations in the vertebrate vitellogenin genes. The encoded polypeptides have diverged from one another sufficiently to allow us to draw some conclusions about conserved positions. Many cysteine residues have been conserved, suggesting that vitellogenin structure has been maintained over a long evolutionary distance and is dependent upon disulfide bonds. In addition, a 20-residue segment shows conservation between the vertebrate and the nematode vitellogenins. This sequence may play a highly conserved role in vitellogenesis, such as specific recognition by oocytes. On the whole, however, selection may be acting more strongly on amino acid composition and codon usage than on amino acid sequence, as might be expected for abundant storage proteins: The amino acid compositions ofvit-2, vit-5, andvit-6 products are remarkably similar, despite the fact that the sequence of thevit-2 protein is only 22% and 50% identical to the sequences ofvit-6 andvit-5 proteins, respectively.     

Predicting conserved protein motifs with Sub-HMMs

Background  

Profile HMMs (hidden Markov models) provide effective methods for modeling the conserved regions of protein families. A limitation of the resulting domain models is the difficulty to pinpoint their much shorter functional sub-features, such as catalytically relevant sequence motifs in enzymes or ligand binding signatures of receptor proteins.     

Structural assessment of glycyl mutations in invariantly conserved motifs

Motifs that are evolutionarily conserved in proteins are crucial to their structure and function. In one of our earlier studies, we demonstrated that the conserved motifs occurring invariantly across several organisms could act as structural determinants of the proteins. We observed the abundance of glycyl residues in these invariantly conserved motifs. The role of glycyl residues in highly conserved motifs has not been studied extensively. Thus, it would be interesting to examine the structural perturbations induced by mutation in these conserved glycyl sites. In this work, we selected a representative set of invariant signature (IS) peptides for which both the PDB structure and mutation information was available. We thoroughly analyzed the conformational features of the glycyl sites and their local interactions with the surrounding residues. Using Ramachandran angles, we showed that the glycyl residues occurring in these IS peptides, which have undergone mutation, occurred more often in the L-disallowed as compared with the L-allowed region of the Ramachandran plot. Short range contacts around the mutation site were analyzed to study the steric effects. With the results obtained from our analysis, we hypothesize that any change of activity arising because of such mutations must be attributed to the long-range interaction(s) of the new residue if the glycyl residue in the IS peptide occurred in the L-allowed region of the Ramachandran plot. However, the mutation of those conserved glycyl residues that occurred in the L-disallowed region of the Ramachandran plot might lead to an altered activity of the protein as a result of an altered conformation of the backbone in the immediate vicinity of the glycyl residue, in addition to long range effects arising from the long side chains of the new residue. Thus, the loss of activity because of mutation in the conserved glycyl site might either relate to long range interactions or to local perturbations around the site depending upon the conformational preference of the glycyl residue.     

CMDWave: conserved motifs detection using wavelets

CMDWave (Conserved Motif Detection using WAVElets) is a web server that predicts conserved motifs in protein sequences. A set of query protein sequences are first aligned using ClustalW to obtain equal sized sequences. CMDWave then converts the sequences into a numerical representation using electron-ion interaction potential (EIIP). This is followed by a wavelet decomposition and reconstruction. A new similarity metric along with thresholding is then used to identify conserved motifs across all the query sequences. Users need not specify the number of motifs to be identified. For larger groups of sequences, results can be emailed to the users.     

Fast, sensitive discovery of conserved genome-wide motifs

Regulatory sites that control gene expression are essential to the proper functioning of cells, and identifying them is critical for modeling regulatory networks. We have developed Magma (Multiple Aligner of Genomic Multiple Alignments), a software tool for multiple species, multiple gene motif discovery. Magma identifies putative regulatory sites that are conserved across multiple species and occur near multiple genes throughout a reference genome. Magma takes as input multiple alignments that can include gaps. It uses efficient clustering methods that make it about 70 times faster than PhyloNet, a previous program for this task, with slightly greater sensitivity. We ran Magma on all non-coding DNA conserved between Caenorhabditis elegans and five additional species, about 70 Mbp in total, in <4 h. We obtained 2,309 motifs with lengths of 6-20 bp, each occurring at least 10 times throughout the genome, which collectively covered about 566 kbp of the genomes, approximately 0.8% of the input. Predicted sites occurred in all types of non-coding sequence but were especially enriched in the promoter regions. Comparisons to several experimental datasets show that Magma motifs correspond to a variety of known regulatory motifs.     

Improvement of alignment accuracy utilizing sequentially conserved motifs

Background  

Multiple sequence alignment algorithms are very important tools in molecular biology today. Accurate alignment of proteins is central to several areas such as homology modelling, docking studies, understanding evolutionary trends and study of structure-function relationships. In recent times, improvement of existing progressing programs and implementation of new iterative algorithms have made a significant change in this field.     

Identification of four conserved motifs among the RNA-dependent polymerase encoding elements.

Four consensus sequences are conserved with the same linear arrangement in RNA-dependent DNA polymerases encoded by retroid elements and in RNA-dependent RNA polymerases encoded by plus-, minus- and double-strand RNA viruses. One of these motifs corresponds to the YGDD span previously described by Kamer and Argos (1984). These consensus sequences altogether lead to 4 strictly and 18 conservatively maintained amino acids embedded in a large domain of 120 to 210 amino acids. As judged from secondary structure predictions, each of the 4 motifs, which may cooperate to form a well-ordered domain, places one invariant amino acid in or proximal to turn structures that may be crucial for their correct positioning in a catalytic process. We suggest that this domain may constitute a prerequisite 'polymerase module' implicated in template seating and polymerase activity. At the evolutionary level, the sequence similarities, gap distribution and distances between each motif strongly suggest that the ancestral polymerase module was encoded by an individual genetic element which was most closely related to the plus-strand RNA viruses and the non-viral retroposons. This polymerase module gene may have subsequently propagated in the viral kingdom by distinct gene set recombination events leading to the wide viral variety observed today.     

Functional significance of conserved amino acid residues.

A systemic study of single amino acid substitutions in bacteriophage T4 lysozyme permitted a test of the concept that conserved amino acid residues are more functionally important than nonconserved residues. Substitutions of amino acid residues that are conserved among five bacteriophage-encoded lysozymes were found to lead more frequently to loss of function than substitutions of nonconserved residues. Of 163 residues tested, only 74 (45%) are sensitive to at least one substitution; however, all 14 residues that are fully conserved are sensitive to substitutions.     

Recognition of distantly related protein sequences using conserved motifs and neural networks.

A sensitive technique for protein sequence motif recognition based on neural networks has been developed. It involves three major steps. (1) At each appropriate alignment position of a set of N matched sequences, a set of N aligned oligopeptides is specified with preselected window length. N neural nets are subsequently and successively trained on N-1 amino acid spans after eliminating each ith oligopeptide. A test for recognition of each of the ith spans is performed. The average neural net recognition over N such trials is used as a measure of conservation for the particular windowed region of the multiple alignment. This process is repeated for all possible spans of given length in the multiple alignment. (2) The M most conserved regions are regarded as motifs and the oligopeptides within each are used to train intensively M individual neural networks. (3) The M networks are then applied in a search for related primary structures in a databank of known protein sequences. The oligopeptide spans in the database sequence with strongest neural net output for each of the M networks are saved and then scored according to the output signals and the proper combination that follows the expected N- to C-terminal sequence order. The motifs from the database with highest similarity scores can then be used to retrain the M neural nets, which can be subsequently utilized for further searches in the databank, thus providing even greater sensitivity to recognize distant familial proteins. This technique was successfully applied to the integrase, DNA-polymerase and immunoglobulin families.     

Functional consequences of mutating conserved SF2 helicase motifs in the Type III restriction endonuclease EcoP15I translocase domain

For efficient DNA hydrolysis, Type III restriction endonuclease EcoP15I interacts with two inversely oriented recognition sites in an ATP-dependent process. EcoP15I consists of two methylation (Mod) subunits and a single restriction (Res) subunit yielding a multifunctional enzyme complex able to methylate or to hydrolyse DNA. Comprehensive sequence alignments, limited proteolysis and mass spectroscopy suggested that the Res subunit is a fusion of a motor or translocase (Tr) domain of superfamily II helicases and an endonuclease domain with a catalytic PD…EXK motif. In the Tr domain, seven predicted helicase motifs (I, Ia, II–VI), a recently discovered Q-tip motif and three additional regions (IIIa, IVa, Va) conserved among Type III restriction enzymes have been identified that are predicted to be involved in DNA binding and ATP hydrolysis. Because DNA unwinding activity for EcoP15I (as for bona fide helicases) has never been found and EcoP15I ATPase rates are only low, the functional importance of the helicase motifs and regions was questionable and has never been probed systematically. Therefore, we mutated all helicase motifs and conserved regions predicted in Type III restriction enzyme EcoP15I and examined the functional consequences on EcoP15I enzyme activity and the structural integrity of the variants by CD spectroscopy. The resulting eleven enzyme variants all, except variant IVa, are properly folded showing the same secondary structure distribution as the wild-type enzyme. Classical helicase motifs I–VI are important for ATP and DNA cleavage by EcoP15I and mutations therein led to complete loss of ATPase and cleavage activity. Among the catalytically inactive enzyme variants three preserved the ability to bind ATP. In contrast, newly assigned motifs Q-tip, Ia and Va are not essential for EcoP15I activity and the corresponding enzyme variants were still catalytically active. DNA binding was only marginally reduced (2–7 fold) in all enzyme variants tested.     

Dbf4 motifs: conserved motifs in activation subunits for Cdc7 kinases essential for S-phase

Dbf4 and its related molecules were originally identified as cyclin-like partners for Cdc7 kinases, essential for S-phase. Recent reports and database search indicate the presence of multiple Dbf4-related molecules with distinct functions. We have identified three stretches of amino acids which are conserved in various Dbf4-related molecules and possibly play distinct functions in binding to and activation of the catalytic subunits as well as in interactions with other proteins. Discovery of conserved motifs for this possible new protein family would serve as a useful framework for future identification of new members of this family as well as for probing their functions.     

Function assignment to conserved residues in mammalian alkaline phosphatases

We have probed the structural/functional relationship of key residues in human placental alkaline phosphatase (PLAP) and compared their properties with those of the corresponding residues in Escherichia coli alkaline phosphatase (ECAP). Mutations were introduced in wild-type PLAP, i.e. [E429]PLAP, and in some instances also in [G429]PLAP, which displays properties characteristic of the human germ cell alkaline phosphatase isozyme. All active site metal ligands, as well as residues in their vicinity, were substituted to alanines or to the homologous residues present in ECAP. We found that mutations at Zn2 or Mg sites had similar effects in PLAP and ECAP but that the environment of the Zn1 ion in PLAP is less affected by substitutions than that in ECAP. Substitutions of the Mg and Zn1 neighboring residues His-317 and His-153 increased k(cat) and increased K(m) when compared with wild-type PLAP, contrary to what was predicted by the reciprocal substitutions in ECAP. All mammalian alkaline phosphatases (APs) have five cysteine residues (Cys-101, Cys-121, Cys-183, Cys-467, and Cys-474) per subunit, not homologous to any of the four cysteines in ECAP. By substituting each PLAP Cys by Ser, we found that disrupting the disulfide bond between Cys-121 and Cys-183 completely prevents the formation of the active enzyme, whereas the carboxyl-terminally located Cys-467-Cys-474 bond plays a lesser structural role. The substitution of the free Cys-101 did not significantly affect the properties of the enzyme. A distinguishing feature found in all mammalian APs, but not in ECAP, is the Tyr-367 residue involved in subunit contact and located close to the active site of the opposite subunit. We studied the A367 and F367 mutants of PLAP, as well as the corresponding double mutants containing G429. The mutations led to a 2-fold decrease in k(cat), a significant decrease in heat stability, and a significant disruption of inhibition by the uncompetitive inhibitors l-Phe and l-Leu. Our mutagenesis data, computer modeling, and docking predictions indicate that this residue contributes to the formation of the hydrophobic pocket that accommodates and stabilizes the side chain of the inhibitor during uncompetitive inhibition of mammalian APs.     

Two conserved essential motifs of the murine immunoglobulin lambda enhancers bind B-cell-specific factors.

Two highly homologous enhancers associated with the two murine immunoglobulin lambda constant-region clusters were recently identified. In order to better understand the molecular basis for the developmental stage- and cell-type-restricted expression of lambda genes, we have undertaken an analysis of the putative regulatory domains of these enhancers. By using a combination of DNase I footprinting, electrophoretic mobility shift assay, and site-specific mutations, four candidate protein binding sites have been identified at analogous positions in both enhancers. A mutation of any of these sites decreases enhancer activity. Two of the sites, lambda A and lambda B, are essential for enhancer function, and both of these sites appear to bind both B-cell-specific and general factors. Nevertheless, isolated lambda A and lambda B sites show no evidence of inherent transactivating potential, alone or together, even when present in up to three copies. We suggest that the generation of transactivating signals from these enhancers may require the complex interaction of multiple B-cell-specific and nonspecific DNA-binding factors.