Wavelet transforms for the characterization and detection of repeating motifs

The role of repeating motifs in protein structures is thought to be as modular building blocks which allow an economic way of constructing complex proteins. In this work novel wavelet transform analysis techniques are used to detect and characterize repeating motifs in protein sequence and structure data, where the Kyte-Doolittle hydrophobicity scale (Eta Phi) and relative accessible surface area (rASA) data provide residue information about the protein sequence and structure, respectively. We analyze a variety of repeating protein motifs, TIM barrels, propellor blades, coiled coils and leucine-rich repeat structures. Detection and characterization of these motifs is performed using techniques based on the continuous wavelet transform (CWT). Results indicate that the wavelet transform techniques developed herein are a promising approach for the detection and characterization of repeating motifs for both structural and in some instances sequence data.     

A Fourier analysis of symmetry in protein structure

The score matrix from a structure comparison program (SAP) was used to search for repeated structures using a Fourier analysis. When tested with artificial data, a simple Fourier transform of the smoothed matrix provided a clear signal of the repeat periodicity that could be used to extract the repeating units with the SAP program. The strength of the Fourier signal was calibrated against the signal from model proteins. The most useful of these was the novel random-walk approach employed to generate realistic 'fake' structures. On the basis of these it was possible to conclude that only a small proportion of protein structures have an unexpected degree of symmetry. Artificially generated 'ideal' folds provided an upper limit on the strength of signal that could be expected from a 'perfectly' repeating compact structure. Unexpectedly, some of the very regular beta-propellor folds attained the same strength but the majority of symmetric structures lay below this region. When native proteins were ranked by the power of their spectrum a wide variety of fold types were seen to score highly. In the betaalpha class, these included the globular betaalpha proteins and the more repetitive leucine-rich betaalpha folds. In the all-beta class; beta-propellors, beta-prisms and beta-helices were found as well as the more globular gamma-crystalin domains. When this ranked list was filtered to remove proteins that contained detectable internal sequence similarity (using the program REPRO), the list became exclusively composed of just globular betaalpha class proteins and in the top 50 re-ranked proteins, only a single 4-fold propellor structure remained.     

Long-range side-chain-main-chain interactions play crucial roles in stabilizing the (betaalpha)8 barrel motif of the alpha subunit of tryptophan synthase

The role of hither-to-fore unrecognized long-range hydrogen bonds between main-chain amide hydrogens and polar side chains on the stability of a well-studied (betaalpha)8, TIM barrel protein, the alpha subunit of tryptophan synthase (alphaTS), was probed by mutational analysis. The F19-D46 and I97-D124 hydrogen bonds link the N terminus of a beta-strand with the C terminus of the succeeding antiparallel alpha-helix, and the A103-D130 hydrogen bond links the N terminus of an alpha-helix with the C terminus of the succeeding antiparallel beta-strand, forming clamps for the respective betaalpha or alphabeta hairpins. The individual replacement of these aspartic acid side chains with alanine leads to what appear to be closely related partially folded structures with significantly reduced far-UV CD ellipticity and thermodynamic stability. Comparisons with the effects of eliminating another main-chain-side-chain hydrogen bond, G26-S33, and two electrostatic side-chain-side-chain hydrogen bonds, D38-H92 and D112-H146, all in the same N-terminal folding unit of alphaTS, demonstrated a unique role for the clamp interactions in stabilizing the native barrel conformation. Because neither the asparagine nor glutamic acid variant at position 46 can completely reproduce the spectroscopic, thermodynamic, or kinetic folding properties of aspartic acid, both size and charge are crucial to its unique role in the clamp hydrogen bond. Kinetic studies suggest that the three clamp hydrogen bonds act in concert to stabilize the transition state leading to the fully folded TIM barrel motif.     

Barrel structures in proteins: automatic identification and classification including a sequence analysis of TIM barrels.

Automated methods for identifying and characterizing regular beta-barrels from coordinate data have been developed to analyze and classify various kinds of barrel structures based on geometric parameters such as the barrel strand number (n) and shear number (S). In total, we find 1,316 barrels in the January 1998 release of Protein Data Bank. Of 1,316 barrels, 1,277 barrels had an even shear number, corresponding to 50 nonhomologous families. The (beta alpha)8 triose phosphate isomerase (TIM) barrel (n = 8, S = 8) fold has the largest number of apparently nonhomologous entries, 16, although the trypsin like antiparallel (n = 6, S = 8) barrels (representing only three families) are the most common with 527 barrels. Of all the protein families that exhibit barrel structures, 68% are found to be various kinds of enzymes, the remainder being binding proteins and transport membrane proteins. In addition, the layers of side chains, which form the cores of barrels with S = n and S = 2n, are also analyzed. More sophisticated methods were developed for detecting TIM barrels specifically, including consideration of the amino acid propensities for the side chains that form the layers. We found that the residues on the outside of the eight stranded parallel beta-barrel, buried by the alpha-helices, are much more hydrophobic than the residues inside the barrel.     

Consistency analysis of similarity between multiple alignments: prediction of protein function and fold structure from analysis of local sequence motifs

A new method to analyze the similarity between multiply aligned protein motifs (blocks) was developed. It identifies sets of consistently aligned blocks. These are found to be protein regions of similar function and structure that appear in different contexts. For example, the Rossmann fold ligand-binding region is found similar to TIM barrel and methylase regions, various protein families are predicted to have a TIM-barrel fold and the structural relation between the ClpP protease and crotonase folds is identified from their sequence. Besides identifying local structure features, sequence similarity across short sequence-regions (less than 20 amino acid regions) also predicts structure similarity of whole domains (folds) a few hundred amino acid residues long. Most of these relations could not be identified by other advanced sequence-to-sequence or sequence-to-multiple alignments comparisons. We describe the method (termed CYRCA), present examples of our findings, and discuss their implications.     

Structural principles of alpha/beta barrel proteins: the packing of the interior of the sheet

Alpha/beta barrel structures very similar to that first observed in triose phosphate isomerase are now known to occur in 14 enzymes. To understand the origin of this fold, we analyzed in three of these proteins the geometry of the eight-stranded beta-sheets and the packing of the residues at the center of the barrel. The packing in this region is seen in its simplest form in glycolate oxidase. It consists of 12 residues arranged in three layers. Each layer contains four side chains. The packing of RubisCO and TIM can be understood in terms of distortions of this simple pattern, caused by residues with small side chains at some of the positions inside the barrel. Two classes of packing are found. In one class, to which RubisCO and TIM belong, the central layer is formed by a residue from the first, third, fifth, and seventh strands; the upper and lower layers are formed by residues from the second, fourth, sixth, and eighth strands. In the second class, to which GAO belongs, this is reversed: it is side chains from the even-numbered strands that form the central layer, and side chains from the odd-numbered strands that form the outer layers. Our results suggest that not all proteins with this fold are related by evolution, but that they represent a common favorable solution to the structural problems involved in the creation of a closed beta barrel.     

Targeted impairment of thymidine kinase 2 expression in cells induces mitochondrial DNA depletion and reveals molecular mechanisms of compensation of mitochondrial respiratory activity

Recently, significant attention has been drawn to the biology of small leucine-rich repeat proteoglycans (SLRPs) due to their multiple functionalities in various cell types and tissues. Here, we characterize a novel SLRP member, “Podocan-like (Podnl) protein” identified by a bioinformatics approach. The Podnl protein has a signal peptide, a unique cysteine-rich N-terminal cluster, 21 leucine-rich repeat (LRR) motifs, and one putative N-glycosylation site. This protein is structurally similar to podocan in SLRPs. The gene was highly expressed in mineralized tissues and in osteoblastic cells and the high expression level was observed at and after matrix mineralization in vitro. Podnl was enriched in newly formed bones based on immunohistochemical analysis. When Podnl was transfected into osteoblastic cells, the protein with N-glycosylation was detected mainly in the cultured medium, indicating that Podnl is a secreted N-glycosylated protein. The endogenous Podnl protein was also present in bone matrix. These data provide a new insight into our understanding of the emerging SLRP functions in bone formation.     

Insights into the fold organization of TIM barrel from interaction energy based structure networks

There are many well-known examples of proteins with low sequence similarity, adopting the same structural fold. This aspect of sequence-structure relationship has been extensively studied both experimentally and theoretically, however with limited success. Most of the studies consider remote homology or "sequence conservation" as the basis for their understanding. Recently "interaction energy" based network formalism (Protein Energy Networks (PENs)) was developed to understand the determinants of protein structures. In this paper we have used these PENs to investigate the common non-covalent interactions and their collective features which stabilize the TIM barrel fold. We have also developed a method of aligning PENs in order to understand the spatial conservation of interactions in the fold. We have identified key common interactions responsible for the conservation of the TIM fold, despite high sequence dissimilarity. For instance, the central beta barrel of the TIM fold is stabilized by long-range high energy electrostatic interactions and low-energy contiguous vdW interactions in certain families. The other interfaces like the helix-sheet or the helix-helix seem to be devoid of any high energy conserved interactions. Conserved interactions in the loop regions around the catalytic site of the TIM fold have also been identified, pointing out their significance in both structural and functional evolution. Based on these investigations, we have developed a novel network based phylogenetic analysis for remote homologues, which can perform better than sequence based phylogeny. Such an analysis is more meaningful from both structural and functional evolutionary perspective. We believe that the information obtained through the "interaction conservation" viewpoint and the subsequently developed method of structure network alignment, can shed new light in the fields of fold organization and de novo computational protein design.     

Partial NMR assignments and secondary structure mapping of the isolated alpha subunit of Escherichia coli tryptophan synthase,a 29-kD TIM barrel protein

The alpha subunit of tryptophan synthase (alphaTS) from S. typhimurium belongs to the triosephosphate isomerase (TIM) or the (beta/alpha)(8) barrel fold, one of the most common structures in biology. To test the conservation of the global fold in the isolated Escherichia coli homolog, we have obtained a majority of the backbone assignments for the 29-kD alphaTS by using standard heteronuclear multidimensional NMR methods on uniformly (15)N- and (15)N/(13)C-labeled protein and on protein selectively (15)N-labeled at key hydrophobic residues. The secondary structure mapped by chemical shift index, nuclear Overhauser enhancements (NOEs), and hydrogen-deuterium (H-D) exchange, and several abnormal chemical shifts are consistent with the conservation of the global TIM barrel fold of the isolated E. coli alphaTS. Because most of the amide protons that are slow to exchange with solvent correspond to the beta-sheet residues, the beta-barrel is likely to play an important role in stabilizing the previously detected folding intermediates for E. coli alphaTS. A similar combination of uniform and selective labeling can be extended to other TIM barrel proteins to obtain insight into the role of the motif in stabilizing what appear to be common partially folded forms.     

Extension of a de novo TIM barrel with a rationally designed secondary structure element

The ability to construct novel enzymes is a major aim in de novo protein design. A popular enzyme fold for design attempts is the TIM barrel. This fold is a common topology for enzymes and can harbor many diverse reactions. The recent de novo design of a four‐fold symmetric TIM barrel provides a well understood minimal scaffold for potential enzyme designs. Here we explore opportunities to extend and diversify this scaffold by adding a short de novo helix on top of the barrel. Due to the size of the protein, we developed a design pipeline based on computational ab initio folding that solves a less complex sub‐problem focused around the helix and its vicinity and adapt it to the entire protein. We provide biochemical characterization and a high‐resolution X‐ray structure for one variant and compare it to our design model. The successful extension of this robust TIM‐barrel scaffold opens opportunities to diversify it towards more pocket like arrangements and as such can be considered a building block for future design of binding or catalytic sites.     

Structural principles of leucine-rich repeat (LRR) proteins

LRR-containing proteins are present in over 2000 proteins from viruses to eukaryotes. Most LRRs are 20-30 amino acids long, and the repeat number ranges from 2 to 42. The known structures of 14 LRR proteins, each containing 4-17 repeats, have revealed that the LRR domains fold into a horseshoe (or arc) shape with a parallel beta-sheet on the concave face and with various secondary structures, including alpha-helix, 3(10)-helix, and pII helix on the convex face. We developed simple methods to charactere quantitatively the arc shape of LRR and then applied them to all known LRR proteins. A quantity of 2Rsin(phi/2), in which R and phi are the radii of the LRR arc and the rotation angle about the central axis per repeating unit, respectively, is highly conserved in all the LRR proteins regardless of a large variety of repeat number and the radius of the LRR arc. The radii of the LRR arc with beta-alpha structural units are smaller than those with beta-3(10) or beta-pII units. The concave face of the LRR beta-sheet forms a surface analogous to a part of a M?bius strip.     

Dipole moment in TIM alpha/beta fold proteins

TIM proteins of alpha/beta barrel fold from alpha/beta class as given in SCOP database were taken for dipole moment analysis. In all, 32 structures were analyzed for their dipole moment contributions. Representative structures from 20 super families in the alpha/beta fold, with different enzyme functions and 12 protein domains of TIM family in TIM super family were considered. The active sites of these proteins are located on the C-terminal side of the beta-strands. The molecules of same alpha/beta fold, but differing in their functionality also showed a common electrostatic field pattern along the barrel axis and had the dipole moment along the barrel axis and towards C-terminal end of the beta-strands. However, it is observed from our calculations that the dipole moment direction is possibly a consequence of the structural fold, with distribution of charges playing a modulatory role, and does not contribute to the location of active site. We show here that apart from the commonly held view as proposed by Hol et al [Hol W G L, van Duijnen PT and Berendsen H J C (1978) Nature (London), 273, 443-446] of the role of the alpha helical dipole moment, the beta-sheets in the barrel can also have a considerable dipole moment contribution. Taken together with our dipole moment analysis on integral membrane proteins [Vasanthi G and Krishnaswamy S (2002) Indian J Biochem Biophys 39, 93-100], this suggests the need to examine the role of dipole moment in the case of especially beta sheets forming barrels.     

PIRLs: a novel class of plant intracellular leucine-rich repeat proteins

Leucine-rich repeat (LRR) proteins feature tandem leucine-rich motifs that form a protein-protein interaction domain. Plants contain diverse classes of LRR proteins, many of which take part in signal transduction. We have identified a novel family of nine Arabidopsis LRR proteins that, based on predicted intracellular location and LRR motif consensus sequence, are related to Ras-binding LRR proteins found in signaling complexes in animals and yeast. This new class has been named plant intracellular Ras group-related LRR proteins (PIRLs). We have characterized PIRL cDNAs, rigorously defined gene and protein annotations, investigated gene family evolution and surveyed mRNA expression. While LRR regions suggested a relationship to Ras group LRR proteins, outside of their LRR domains PIRLs differed from Ras group proteins, exhibiting N- and C-terminal regions containing low complexity stretches and clusters of charged amino acids. PIRL genes grouped into three subfamilies based on sequence relationships and gene structures. Related gene pairs and dispersed chromosomal locations suggested family expansion by ancestral genomic or segmental duplications. Expression surveys revealed that all PIRL mRNAs are actively transcribed, with three expressed differentially in leaves, roots or flowers. These results define PIRLs as a distinct, plant-specific class of intracellular LRR proteins that probably mediate protein interactions, possibly in the context of signal transduction. T-DNA knock-out mutants have been isolated as a starting point for systematic functional analysis of this intriguing family.     

Molecular modeling of 2-nitropropane dioxygenase domain of Mycobacterium tuberculosis H37Rv and docking of herbal ligands

The 3D structure of enoyl reductase (ER) domain generated by the SWISS MODEL server contains the 2-nitropropane dioxygenase (2NPD) structure displaying the TIM barrel fold. Though TIM barrel fold is made up of both main and inserted domains, in our study, we could only predict the structure of the main domain, which had central barrel of eight beta-strands surrounded by eight alpha-helices. Superimposition of the 2NPD region of ER domain of Mycobacterium tuberculosis H37Rv on to the corresponding region of 2UVA_G revealed a good structural alignment between the two, suggesting this template to be a good structural homologue. Among various herbal ligands that were screened as inhibitors, daucosterol was found to bind in closest proximity to the flavin mono nucleotide (FMN) binding site with the lowest docking energy.     

Crystal structure of a methyltetrahydrofolate- and corrinoid-dependent methyltransferase

BACKGROUND: Methyltetrahydrofolate, corrinoid iron-sulfur protein methyltransferase (MeTr), catalyzes a key step in the Wood-Ljungdahl pathway of carbon dioxide fixation. It transfers the N5-methyl group from methyltetrahydrofolate (CH3-H4folate) to a cob(I)amide center in another protein, the corrinoid iron-sulfur protein. MeTr is a member of a family of proteins that includes methionine synthase and methanogenic enzymes that activate the methyl group of methyltetra-hydromethano(or -sarcino)pterin. We report the first structure of a protein in this family. RESULTS: We determined the crystal structure of MeTr from Clostridium thermoaceticum at 2.2 A resolution using multiwavelength anomalous diffraction methods. The overall architecture presents a new functional class of the versatile triose phosphate isomerase (TIM) barrel fold. The MeTr tertiary structure is surprisingly similar to the crystal structures of dihydropteroate synthetases despite sharing less than 20% sequence identity. This homology permitted the methyl-H4folate binding site to be modeled. The model suggests extensive conservation of the pterin ring binding residues in the polar active sites of the methyltransferases and dihydropteroate synthetases. The most significant structural difference between these enzymes is in a loop structure above the active site. It is quite open in MeTr, where it can be modeled as the cobalamin binding site. CONCLUSIONS: The MeTr structure consists of a TIM barrel that embeds methyl-H4folate and cobamide. All related methyltransferases are predicted to fold into a similar TIM barrel pattern and have a similar pterin and cobamide binding site. The observed structure is consistent with either a 'front' (N5) or 'back' (C8a) side protonation of CH3-H4folate, a key step that enhances the electrophilic character of the methyl group, activating it for nucleophilic attack by Co(I).     

Alternative Splice Variants in TIM Barrel Proteins from Human Genome Correlate with the Structural and Evolutionary Modularity of this Versatile Protein Fold

After the surprisingly low number of genes identified in the human genome, alternative splicing emerged as a major mechanism to generate protein diversity in higher eukaryotes. However, it is still not known if its prevalence along the genome evolution has contributed to the overall functional protein diversity or if it simply reflects splicing noise. The (βα)₈ barrel or TIM barrel is one of the most frequent, versatile, and ancient fold encountered among enzymes. Here, we analyze the structural modifications present in TIM barrel proteins from the human genome product of alternative splicing events. We found that 87% of all splicing events involved deletions; most of these events resulted in protein fragments that corresponded to the (βα)₂, (βα)₄, (βα)₅, (βα)₆, and (βα)₇ subdomains of TIM barrels. Because approximately 7% of all the splicing events involved internal β-strand substitutions, we decided, based on the genomic data, to design β-strand and α-helix substitutions in a well-studied TIM barrel enzyme. The biochemical characterization of one of the chimeric variants suggests that some of the splice variants in the human genome with β-strand substitutions may be evolving novel functions via either the oligomeric state or substrate specificity. We provide results of how the splice variants represent subdomains that correlate with the independently folding and evolving structural units previously reported. This work is the first to observe a link between the structural features of the barrel and a recurrent genetic mechanism. Our results suggest that it is reasonable to expect that a sizeable fraction of splice variants found in the human genome represent structurally viable functional proteins. Our data provide additional support for the hypothesis of the origin of the TIM barrel fold through the assembly of smaller subdomains. We suggest a model of how nature explores new proteins through alternative splicing as a mechanism to diversify the proteins encoded in the human genome.     

An allosteric pathway explains beneficial fitness in yeast for long‐range mutations in an essential TIM barrel enzyme

Protein evolution proceeds by a complex response of organismal fitness to mutations that can simultaneously affect protein stability, structure, and enzymatic activity. To probe the relationship between genotype and phenotype, we chose a fundamental paradigm for protein evolution, folding, and design, the (βα)₈ TIM barrel fold. Here, we demonstrate the role of long‐range allosteric interactions in the adaptation of an essential hyperthermophilic TIM barrel enzyme to mesophilic conditions in a yeast host. Beneficial fitness effects observed with single and double mutations of the canonical βα‐hairpin clamps and the α‐helical shell distal to the active site revealed an underlying energy network between opposite faces of the cylindrical β‐barrel. We experimentally determined the fitness of multiple mutants in the energetic phase plane, contrasting the energy barrier of the chemical reaction and the folding free energy of the protein. For the system studied, the reaction energy barrier was the primary determinant of organism fitness. Our observations of long‐range epistatic interactions uncovered an allosteric pathway in an ancient and ubiquitous enzyme that may provide a novel way of designing proteins with a desired activity and stability profile.     

Structural basis for the function of pyridoxine 5'-phosphate synthase

BACKGROUND: Pyridoxal 5'-phosphate is the active form of vitamin B(6) that acts as an essential, ubiquitous coenzyme in amino acid metabolism. In Escherichia coli, the pathway of the de novo biosynthesis of vitamin B(6) results in the formation of pyridoxine 5'-phosphate (PNP), which can be regarded as the first synthesized B(6) vitamer. PNP synthase (commonly referred to as PdxJ) is a homooctameric enzyme that catalyzes the final step in this pathway, a complex intramolecular condensation reaction between 1-deoxy-D-xylulose-5'-phosphate and 1-amino-acetone-3-phosphate. RESULTS: The crystal structure of E. coli PNP synthase was solved by single isomorphous replacement with anomalous scattering and refined at a resolution of 2.0 A. The monomer of PNP synthase consists of one compact domain that adopts the abundant TIM barrel fold. Intersubunit contacts are mediated by three additional helices, respective to the classical TIM barrel helices, generating a tetramer of symmetric dimers with 422 symmetry. In the shared active sites of the active dimers, Arg20 is directly involved in substrate binding of the partner monomer. Furthermore, the structure of PNP synthase with its physiological products, PNP and P(i), was determined at 2.3 A resolution, which provides insight into the dynamic action of the enzyme and allows us to identify amino acids critical for enzymatic function. CONCLUSION: The high-resolution structures of the free enzyme and the enzyme-product complex of E. coli PNP synthase suggest essentials of the enzymatic mechanism. The main catalytic features are active site closure upon substrate binding by rearrangement of one C-terminal loop of the TIM barrel, charge-charge stabilization of the protonated Schiff-base intermediate, the presence of two phosphate binding sites, and a water channel that penetrates the beta barrel and allows the release of water molecules in the closed state. All related PNP synthases are predicted to fold into a similar TIM barrel pattern and have comparable active site architecture. Thus, a common mechanism can be anticipated.