- Research article
- Open Access
A nested leucine rich repeat (LRR) domain: The precursor of LRRs is a ten or eleven residue motif
BMC Microbiology volume 10, Article number: 235 (2010)
Leucine rich repeats (LRRs) are present in over 60,000 proteins that have been identified in viruses, bacteria, archae, and eukaryotes. All known structures of repeated LRRs adopt an arc shape. Most LRRs are 20-30 residues long. All LRRs contain LxxLxLxxNxL, in which "L" is Leu, Ile, Val, or Phe and "N" is Asn, Thr, Ser, or Cys and "x" is any amino acid. Seven classes of LRRs have been identified. However, other LRR classes remains to be characterized. The evolution of LRRs is not well understood.
Here we describe a novel LRR domain, or nested repeat observed in 134 proteins from 54 bacterial species. This novel LRR domain has 21 residues with the consensus sequence of LxxLxLxxNxLxxLDLxx(N/L/Q/x)xx or LxxLxCxxNxLxxLDLxx(N/L/x)xx. This LRR domain is characterized by a nested periodicity; it consists of alternating 10- and 11- residues units of LxxLxLxxNx(x/-). We call it "IRREKO" LRR, since the Japanese word for "nested" is "IRREKO". The first unit of the "IRREKO" LRR domain is frequently occupied by an "SDS22-like" LRR with the consensus of LxxLxLxxNxLxxLxxLxxLxx or a "Bacterial" LRR with the consensus of LxxLxLxxNxLxxLPxLPxx. In some proteins an "SDS22-like" LRR intervenes between "IRREKO" LRRs.
Proteins having "IRREKO" LRR domain are almost exclusively found in bacteria. It is suggested that IRREKO@LRR evolved from a common ancestor with "SDS22-like" and "Bacterial" classes and that the ancestor of IRREKO@LRR is 10 or 11 residues of LxxLxLxxNx(x/-). The "IRREKO" LRR is predicted to adopt an arc shape with smaller curvature in which β-strands are formed on both concave and convex surfaces.
LRR (leucine rich repeat) domains are present in over 60, 000 proteins listed in PFAM, PRINTS, SMART, InterPro and PANTHER databases . LRR-containing proteins have been identified in viruses, bacteria, archae, and eukaryotes. Most LRR proteins are involved in protein, ligand and in protein, protein interactions; these include plant immune response and the mammalian innate immune response [2–6].
All LRR units can be divided into a HCS (highly conserved segment) and a VS (variable segment). The HCS part consists of an eleven residue stretch, LxxLxLxxNxL, or a twelve residue stretch, LxxLxLxxCxxL, in which "L" is Leu, Ile, Val, or Phe, "N" is Asn, Thr, Ser, or Cys, and "C" is Cys, Ser or Asn. Three residues at positions 3 to 5 in the highly conserved segments form a short β-strand. The β-strands stack parallel and the multiple LRRs then form an arc. The concave face consists of a parallel β-sheet and the convex face is made of a variety of secondary structures including the a-helix, 310-helix, polyproline II helix, and an extended structure or a tandem arrangement of β-turns. In most LRR proteins the β-strands on the concave surface and (mostly) helical elements on the convex surface are connected by short loops or β-turns. Seven classes of LRRs have been recognized, characterized by different lengths and consensus sequences of the VS part of the repeats [7, 8]. They are "RI-like", "CC", "Bacterial", "SDS22-like", "plant specific", "typical", and "TpLRR". The seven classes of LRR domains adopt a variety of structures.
"Typical" LRRs are the most abundant LRR class. The consensus sequence is L xxL xL xxN xL xxL pxxoFxxL xx. The repeat length is 20-27 residues. Bold uppercase letters indicate more than 70% occurrence of a given residue in a certain position; normal letters indicate 40-70% occurrence and lowercase letters indicate 30-40% occurrence; "o" indicates a non-polar residue, and "x" indicates nonconserved residues. Their variable segments adopt mainly polyproline II plus β-turn, consecutive β-turns or β-turn plus polyproline II in the convex faces; the structural units may be represented by β - (βt + PPII). "RI-like" LRRs are contained in proteins such as ribonuclease inhibitor and Ran GTPase activating protein. The consensus sequence is L xxL xL xxN x(L/C)xxxgoxxLxxoLxxxxx. The repeat length is 28-29. Their VSs mainly adopt α-helix (β - α structural units). Cysteine-containing (CC) LRR proteins include GRR1 proteins from Saccharomyces cerevisiae. The consensus sequence is L xxL xL xxC xxI TD xxoxxL(a/g)xx(C/L)xx. The repeat length is 25-27. Their VSs mainly adopt α-helix (β - α structural units). A GALA-LRR is a subclass of CC-LRR; its consensus sequence is L xxL xL xxN xIgd x(g/a)axxLax(n/s/d)xx of 24 residues . Plant-specific (PS) LRR proteins include PGIP and Cf-2.1. The consensus sequence is L xxL xL xxN xL( t/s)G xIP xxLGxLxx. The repeat length is 23-25. The VSs mainly adopt 310 - helix. Also in individual LRRs the β-strand on the concave face at the N-terminus and the 310 - helix on the convex face at the C-terminus is connected by a β-turn; the structural units are β - (βt + 310). "SDS22-like" LRRs are included in SDS22 and internalins. The consensus sequence is L xxL xL xxN( r/k)I( r/k)(r/k)IE(N/G)LExL xx. The repeat length is 21-23. The structural units of individual repeats are β - 310. "Bacterial" LRRs are found in YopM from Yersinia pestis, and IpaH from Shigella flexneri. The consensus sequence is L xxL xV xxN xL xxLP( D/E)LP xx. The repeat length is 20-22. The structural units are β - pII. "TpLRR" are found in Treponema pallidum LRR protein and in Bacteroides forsythus surface antigen. The consensus sequence is L xxL xLxxxL xxIgxxAFxx(C/N)xx. The repeat length is 23-25. The dominant feature is a highly conserved segment of ten residues, differing from the corresponding eleven residues of other LRRs. The structure of this class remains unknown.
Most of the known LRR structures have a cap, which shields the hydrophobic core of the first unit of LRR domain at the N-terminus and/or the last unit at the C-terminus. In extracellular proteins or extracellular regions, these caps frequently consist of Cys clusters including two or four Cys residues; the Cys clusters on the N- and C-terminal sides of the LRR arcs are called LRRNT and LRRCT, respectively [4–6]. Non-LRR, island regions interrupting LRRs are widely distributed. Island regions are observed in many LRR proteins including plant LRR-RLKs, plant LRR-RLPs, insect Toll and Toll-related proteins, Slit proteins, fungi adenylate cyclases, and Leishmania proteophosphoglycans [10–14].
The evolution of LRRs is not well understood. It is not even known whether all LRR's share a common ancestor. Kobe and Deisenhofer  pointed out the possibility of their having been at least a few independent occurrences of LRRs. Kajava  also suggested separate origins for several different classes of LRRs based on the high levels of conservation within each LRR class. In contrast, Andrade et al.,  found that searches by a homology-based method, REP, could not absolutely partition LRRs into these separate classes and thus they suggested that these proteins have a common origin, rather than separate origins as proposed by Kajava.
Duplication and recombination as a mechanism of the evolution of the disease resistance gene (R-gene) from various plant species has been proposed by many investigators [16–24]. Distinct higher-order repeating units of LRR's occur in a group of LRR proteins including ribonuclease inhibitor, the subfamily of small leucine-rich repeat proteoglycan (SLRP), and the subfamily of Toll-like receptors (TLR7, TLR8 and TLR9) [4, 25–28]. An evolutionary model has been proposed that involves duplication of the higher-order LRR repeating units [26, 28]. Moreover, the possibility of horizontal gene transfer (HGT) has been discussed .
Escherichia coli yddk is 318 residues long and contains 13 tandem repeats of LRRs; six of the 13 repeats have the consensus of LxxLxLxxNxLxxLxLxxxxx with 21 residues (Figure 1A). The variable segment differs significantly from those of the above seven classes. The purpose of this paper is to investigate the occurrence of this novel domains. We identified many LRR proteins having the novel domain (called IRREKO@LRR) and analyzed their sequences. We discuss the evolution and structure of "IRREKO" LRR.
Proteins having IRREKO@LRRs
We identified a total of 134 IRREKO@LRR proteins from 54 bacterial species including Escherichia, Shigella, Vibrio, Shewanella, Photobacterium, Bifidobacterium, Porphyromonas, Treponema, Listeria, Alistipes, Bacteroides, Clostridium, Cytophaga, and Flavobacterium (Additional file 1, Table 1). A group of these proteins contain a signal peptide (but have no transmembrane helix), indicating that they are extracellular. The others lack both a signal peptide and a transmembrane helix, indicating that they are intracellular.
Some extracellular IRREKO@LRR proteins contain Cys clusters on the N-terminal side of the IRREKO@LRR domain (LRRNT); while LRRCT is not observed. For examples, IRREKO@LRR proteins from Vibrio, Shewanella, and Photobacterium have an LRRNT with the pattern of Cx 16 C (Additional file 1, Table 1). Three Vibrio IRREKO@LRR proteins (VV2_1682, CPS_3882 and VVA0501) have an LRRNT of Cx 20 C. Cysteine in the first LRR sometimes participates in LRRNT (Figure 1).
Some IRREKO@LRR proteins have non-LRR, island regions interrupting LRRs (Figure 1 and Additional files 1 and 2: Table 1 and Figure S1, respectively). They include KAOT1_04155 from Kordia algicida OT-1, CPS_3882 from Vibrio psychroerythus, Fjoh_1188/FjohDRAFT_4748 and Fjoh_1189/FjohDRAFT_4747 from Flavobacterium johnsoniae, Cbac1_010100006401 from Clostridiale bacterium 1_7_47_FAA, SCB49_05520 from unidentified eubacterium, ALIPUT_01468 from Alistipes putredinis, and FAEPRAM212_00794 from Faecalibacterium prausnitzii M21/2.
There is a single example of an "IRREKO" domain from a eukaryote and a single example from a virus. The eukaryote protein is TVAG_084780 from Trichomonas vaginalis G3 (Figure 1Q and Additional file 2, Figure S1). TVAG_084780 contains 10 LRRs. Two of the 10 repeats are clearly "IRREKO" domains. The virus protein is MSV251 from Melanoplus sanguinipes entomopoxvirus [Q9YVJ1]. This protein contains 11 LRRs with the consensus of L kyL dC sNN xL xnL xiN(n/d)n (Additional file 1, Table 1). The repeating unit length is 19 residues and thus shorter than that of typical "IRREKO" LRR.
Two subtypes of IRREKO@LRR domains
IRREKO@LRRs that are 21 residues long may be classified into two subtypes (Figure 1). The first subtype has the consensus of LxxLxLxxNxLxxLDLxx(N/L/Q/x)xx, while the second has the consensus of LxxLxCxxNxLxxLDLxx(N/L/x)xx, where "L" is Leu, Val, Ile, Phe, Met or Ala, "N " is Asn, Thr or Ser, "D" is Asp or Asn, "Q" is Gln, and "x" is nonconserved residues. As well as the other seven classes, "x" is generally hydrophilic or neutral residues (Figure 1 and Additional files 1 and 2: Table 1 and Figure S1, respectively).
In these two subgroups, "L" at positions 1, 4, 14 and 16 is predominantly Leu, while "L" or "C" at position 6 is not only Leu or Cys but also Val or Ile, and frequently Ala and Phe. "N" at position 9 is predominantly Asn and often Thr, Ser or Cys. "D" at position 15 is predominantly occupied by Asp and frequently by Asn. Position 19 is often occupied by Leu, Asn, or Gln. Some IRREKO@LRR proteins such as Listeria internalin-J homologs and four Bacteroides proteins include LRRs in which the HCS part consists of a twelve residue stretch, LxxLxLxx(N/C)xxL As LRRs with 20 or 22 residues sometimes keep the most conserved segments of Lx(L/C) in both HCS and VS parts, we regard those as IRREKO@LRR.
IRREKO@LRR domains that mainly consist of the first subtype are observed in 61 proteins (Additional file 1, Table 1). Some proteins have the consensus of LxxLxLxxNxLxxLDLxxNxx. These include BIFLAC_05879 and BLA_0865 from Bifidobacterium animalis, A1Q_3393, VAS14_09189, VAS14_14509, and CPS_2313 from Vibrio species, SwooDRAFT_0647, SwooDRAFT_0647, and Shal_3481 from Shewanella species, and SKA34_06710 and SKA34_09358 from Photobacterium sp. SKA34 (Figures 1B, C and 1D, and Additional file 2, Figure S1). Also, the consensus of LxxLxLxxNxLxxLDLxxLxx is observed in a few proteins including SCB49_09905 from unidentified eubacterium SCB49 (Figure 1E). The pattern of LxxLxLxxNxLxxLDLxxQxx is observed in only CPS_3882 from Vibrio psychroerythus (Figure 1F).
IRREKO@LRR domains that consist mainly of the second subtype are observed in 57 proteins (Additional file 1, Table 1). The consensus of LxxLxCxxNxLxxLDLxxNxx in which "L" at position 16 is more frequently occupied by Val or Ile than by Leu is observed in some proteins. They include Listeria lmo0331 homologs, CHU_0515 from Cytophaga hutchinsonii and PORUE0001_1723 from Porphyromonas uenonis 60-3 (Figure 1G). Also, the pattern of LxxLxCxxNxLxxLDLxxLxx is observed in TDE_0593, TDE_2231, and TDE_2003 from Treponema denticola (Figure 1H, and Additional file 2, Figure S1). Moreover, the pattern of LxxLxCxxNxLxxLDLxxVxx is observed in Pnap_3264 from Polaromonas naphthalenivorans and MldDRAFT_4836 from Delta proteobacterium MLMS-1 (Figures 1I and 1J, and Additional file 2, Figure S1).
The coexistence of the first and the second subtypes is observed in the LRR domains in at least six IRREKO@LRR proteins. They include KAOT1_04155 from Kordia algicida OT-1, COPEUT_03021 from Coprococcus eutactus ATCC 27759, Fjoh_1188/FjohDRAFT_4748 and Fjoh_1189/FjohDRAFT_4747 from Flavobacterium johnsoniae, RUMGNA_03120 from Ruminococcus gnavus ATCC 29149, DORFOR_03338 from Dorea formicigenerans ATCC 27755, and internain-J homologs from eleven Listeria monocytogenes strains (Figures 1K and 1L, and Additional file 2, Figure S1).
Nested periodicity of IRREKO@LRRs
IRREKO@LRRs show a characteristic, nested periodicity; the domains consist of alternating 10- and 11- residue units of LxxLxLxxNx(x/-). To confirm this periodic nesting we performed detailed sequence analysis of IRREKO@LRR proteins using dot plots analysis and a radar chart analysis.
Self dot plots were performed for four IRRECO@LRR proteins - BIFLAC_05879 from Bifidobacterium animalis, A1Q_3393 from Vibrio harveyi HY01, lmo0331 protein from Listeria monocytogenes and an internalin-related protein, TDE_0593, from Treponema denticola - (Additional file 3, Figure S2). The self dot plots indicate that these proteins demonstrate tandem repeats of short residues that is ~10-11 residues long, in addition to tandem repeats of IRRECO@LRR with 21 residues.
Radar charts were drawn for three families of IRREKO@LRRs proteins, in which the occurrence frequency of amino acids is compared between positions 1-10 and positions 11-21. Figure 2A shows a radar chart of Vibrio proteins. Seven Vibrio species encode twelve IRREKO@LRR proteins which are potential homologs (Additional file 1, Table 1). The IRREKO@LRRs domains in their proteins contain 158 LRR repeats. One hundred thirty-seven of the 158 repeats are complete "IRREKO" domains with 21 residues. The radar chart of the 137 LRRs is shown in Figure 2. As expected, "L" at positions 1, 4, and 6 is highly conserved with positions 11, 14 and 16, respectively. In addition, a significant, weak conservation is observed between positions 10 and 21 but not 20, because amino acid distribution of positions 10 and 21 is very similar and are relatively rich in Lys, Asn and Gln. Alsor, positions 3 and 13 show a conservation in which the amino acids are relatively rich Ser, Thr, Asp and Glu. Moreover, positions 7 and 17 show a weak conservation at which those are relatively rich in Ser and Thr.
Similarly, in addition to high conservation of positions of 1-11, 4-14, and 6-16, a weak conservation among even "x" positions occupied by non-conserved residues is also observed in IRREKO@LRRs within nine, potential homologs from four Shewanella species; positions 2-12, 3-13, and 7-17 are relatively rich in Thr and Ser, and in those within four, potential homologs from two Photobacterium species; positions 3-13 are relatively rich in Thr, Ser, Asp and Glu, and positions 7-17 are relatively rich in Ser and Thr, and positions 10-21 are relatively rich in Gln and Lys (Figures 2B and 2C).
The analyses of both dot plots analysis and radar chart demonstrate that IRREKO@LRRs show a nested periodicity consisting of alternating 10- and 11- residue units with the consensus of LxxLxLxxNx(x/-).
Secondary structure prediction
The protein secondary structure prediction of IRREKO@LRR proteins was performed (Additional file 4, Figure S3). E. coli yddk contains 13 LRRs (Figure 1A). Proteus and SSpro4.0 [30, 31] predict that 12 of the 13 LRRs prefer β-strands at positions 3 through 5 and/or its neighboring positions in the HCS part; although only the eighth LRR does not prefer β-strand, its HCS part - V TYF SA AHN QL- is clearly a canonical LRR. Similarly, all or most LRRs in other proteins prefer β-strands at the corresponding positions in the HCS part.
Both methods of secondary structure prediction indicate that residues at positions 13 through 15 and/or its neighboring positions prefer coil conformations in most LRRs of E. coli yddk, Listeria lmo0331 protein, and Treponema TDE_0593. On the other hand, in most LRRs of Bifidobacterium BIFLAC_05879, Vibrio A1Q_3393 and Shewanella SwooDRAFT_0647, residues at the corresponding positions prefer β-strands. It is concluded that individual three residues at positions 3 to 5 and 13 to 15 could form a short β-strand.
Occurrence of "SDS22-like" and "Bacterial" LRR domains within IRREKO@LRR domains
The first LRR of LRR domain in a large number of IRREKO@LRR proteins are an "SDS22-like" domain, LxxLxLxxNxLxxLxxLxxLxx; even though "N" at position 9 is sometimes occupied by Lys, Gln or Leu (which is frequently seen in the first LRR of LRR domains consisting of only other LRR classes) (Additional file 1, Table 1). Their proteins include eleven proteins from seven Vibrio species, eight proteins from five Shewanella species, eleven internalin-J homologs from eleven Listeria monocytogenes strains, nine lmo0331 homologs from eight L. monocytogenes strains and L. innocua, and nine proteins from three Flavobacterium species.
"SDS22-like" LRR occurs even in the middle position in the IRREKO@LRR domains in some proteins. Cbac1_010100006401 from Clostridiale bacterium 1_7_47_FAA with 1,002 residues contains 16 tandem repeats of LRRs; one non-LRR, island region is observed between the seventh and eighth LRRs (Figure 1M, and Additional file 2, Figure S1). Twelve of the 16 repeats are "IRREKO" domain with 20-22 residues. On the other hand, the remaining (LRRs 3, 5, 10 and 11) belong to "SDS22-like" class with the consensus is LxxLxCxxNxLxxLxxLxxLxx.
The three Listeria lin1204 homologs - LMOf6854_0364, LMOh7858_0369, and LMOf2365_0349 - have 993-1,099 residues and contain 25 tandem repeats of LRRs (Figure 1N and Additional file 2, Figure S1). Six of the 25 repeats are "IRREKO" domain, while eight repeats are "SDS22-like" class.
Other examples include FB2170_11006 from Flavobacteriale bacterium HTCC2170 and three proteins - BACOVA_03150 from Bacteroides ovatus, BACCAC_03004 from Bacteroides caccae ATCC 43185, and BACFIN_03505 from Bacteroides finegoldii DSM 17565 - that are homologous to each other (Additional file 1, Table 1). The former contains nine tandem repeats of LRRs and the third LRR of L VLV EI LAN EL HTI KGL SKM TQ is an "SDS22-like" class. The latter three proteins contains eight tandem repeats of LRRs. The fifth LRR is I AI LI GC AF QSL DI LC CPS and thus appears to be a "SDS22-like" domain.
Five ECUMM_1703 homologs from three Escherichia coli strains and two Shigella species contain 11-15 tandem repeats of LRRs (Figure 1O and Additional file 1, Table 1). Three ECs2075/Z2240 homologs from several Escherichia coli strains and two Shigella strains contain four or five tandem repeats of LRRs (Figure 1P and Additional file 1, Table 1). The first LRR are all M ASL DL SYL DL SELP PIP ST and thus belongs to "Bacterial" class with the consensus of LxxLxLxxNxLxxLPxLPxx (although "N" at position 9 is often occupied by Leu) . Three ECUMM_1723 homologs occur in three E. coli strains with 11 repeats of IRREKO@LRR. The first LRR is QNDI DL SGL NL (T/S)TQP PGL QN. It may belong to "Bacterial" LRR.
IRREKO@LRR as new class of LRR
The present observations indicate that IRREKO@LRR is a new class of LRR. This is supported by several additional observations. The identification of LRRs by PFAM or SMART occurs in a large number of IRREKO@LRR proteins including E. coli yddK; this results from the significant similarity of their HCSs with those of the other LRR classes. There are many LRR proteins that contain the LRR domain consisting mainly of "SDS22-like" domain. The "SDS22-like" LRRs in Listeria lin1204/LMOf6854_0364 and Microcoleus chthonoplastes PCC 7420 MC7420_1958 [B4VM60] also have some IRREKO@LRR domains.
The IRREKO@LRRs show a nested periodicity consisting of alternating 10- and 11- residue units with the consensus of Lxx(L/C)xLxxNx(x/-). The IRREKO@LRR domains in many proteins contain a mixture of both subtypes. The first LRR of the LRR domains is frequently "SDS22-like" or "Bacterial" classes. In addition, among the IRREKO@LRR domain "SDS22-like" class occurs in some proteins. The two subtypes of IRREKO@LRR appear to have evolved from a common precursor. Further, the "IRREKO" domain evolved from a precursor common to "SDS22-like" and "Bacterial" classes. The precursor of IRREKO@LRR is shorter sequence - LxxLxLxxNx(x/-) -. This parsimonious evolutionary scenario for three LRR classes, "IRREKO", "SDS22-like", and "Bacterial" LRRs is shown in Figure 3.
Previous studies revealed that LRR domains in many LRR proteins contain tandem repeats of a super-domain of STT, where "T" is "typical" LRR and "S" is "Bacterial" LRR; they include the SLRP subfamily (biglycan, decorin, asporin, lumican, fibromodulin, PRELP, keratocan, osteoadherin, epiphycan, osteoglycin, opticin, and podocan), the TLR7 family (TLR7, TLR8 and TLR9), the FLRT family (FLRT1, FLRT2, and FLRT3), and OMGP [4, 25–27]. The combination of the previous and the present observations suggest that the four LRR classes of "Bacterial", "typical", "SDS22-like" and "IRREKO" might evolve from a common precursor.
The known LRR structures reveal that conserved hydrophobic residues in the consensus contribute to the hydrophobic cores in the LRR arcs [2–6]. As noted, the consensus of IRREKO@LRR is LxxLxLxxNxLxxLDLxx(N/L/Q/x)xx or LxxLxCxxNxLxxLDLxx(N/L/x)xx. It is likely that the conserved hydrophobic residues at the six (or seven) positions of 1, 4, 6 and 11, 14 and 16 (and 19) participate in the hydrophobic core (Figure 4).
The LRR structures with α-helices in their convex faces have more pronounced curvature than structures with 310 or polyproline II helices [4, 32]. This difference in curvature is attributed to the differences in diameter of the different secondary structure elements on the convex face, α-helices being wider than 310-helices, polyproline II helices or tandem β-turns. IRREKO@LRR is predicted to adopt β-β structural units, because individual three residues at positions 3 to 5 and 13 to 15 could form a short β-strand (Figure 4). β-strands have the smallest diameter. Moreover, the loops that link the C-terminal ends of the β-strands in the HCS to the N termini of those in the VS appear to be different from the loops that link the C-terminal ends of those in the VS to the N termini of the following β-strands, as the HCS is one residue longer than the VS. Thus, an inferred arc structure of IRREKO@LRR has a smaller curvature.
Position 2 in the i-th and the (i+1)-th repeats of IRREKO@LRRs is alternatively occupied by positive and negative charged amino acids in some proteins. Examples include CdifQCD-2_010100017965 and CdifQ_04001775 from Clostridium difficile and CHU_1860 from Cytophaga hutchinsonii, as well as FjohDRAFT_1094 and Fjoh_0631 from Flavobacterium johnsoniae (Additional file 1, Table 1). The inferred arc structure of IRREKO@LRRs will enable them to form polar hydrogen bond interactions which lead to its structural stability.
It is possible that the β-solenoid structure of IRREKO@LRRs is related to β-helix proteins [33–35]. A β-β structural unit that is responsible for tandem repeats of GGxGxD is also observed in serralysin . The β-solenoids with β-β structural units in IRREKO@LRR protein and serralysin represent an example of convergent evolution. Future studies should resolve this question.
IRREKO@LRR is a new, unique class of LRR. IRREKO@LRR with the consensus of LxxLx(L/C) xxNxLxxLxLxx(L/Q/x)xx is a nested sequence consisting of alternating 10 - and 11-residue units of LxxLxLxxNx(x/-). The IRREKO@LRR domains frequently coexist with "SDS22-like" or "Bacterial" LRR. These findings suggest that the ancestor of IRREKO@LRR is shorter residues of LxxLxLxxNx(x/-) and that IRREKO@LRR evolved from a common ancestor with "SDS22-like" and "Bacterial" classes. IRREKO@LRRs are predicted to adopt an arc shape with smaller curvature in which individual repeats adopt β-β structural units.
The putative uncharacterized protein yddK from Escherichia coli (strain K12) with 318 residues [YDDK_ECOLI] is an LRR protein. It is identified in the data bases of InterPro, PFAM, PRINTS and SMART. The InterPro data base indicates that the LRR domain contains nine repeats. The PFAM program predicts that yddK contain one significant LRR (residues 216-238) and seven insignificant LRRs (12-30; 33-53; 109-131; 153-175; 196-213; 260-282; 284-306).
We recently developed a new method that utilizes known LRR structures to recognize and align new LRR domains and incorporate multiple sequence alignments and secondary structure predictions . This method predicts correctly the number of LRRs, their lengths and their boundaries. Its usefulness was confirmed by crystal structures of TLR1, TLR2, and TLR4 [37, 38].
This new method was used for multiple sequence alignments of LRRs in the yddK protein. This analysis predicted not nine repeats of the LRRs but 13 repeats and also revealed that their "phasing" differ significantly. We noticed that LRRs, 1, 5 7, 8, 9, and 10 contain a unique domain whose consensus is LxxLxLxxNxLxxLxLxxxxx with 21 residues. The variable segment offers a characteristic hydrophobic pattern unidentified previously (Figure 1A). Each LRR domain is a nested sequence and consists of repeats alternating 10- and 11- residue units of LxxLxLxxNx(x/-).
LRR proteins having the IRREKO@LRR domains were identified in three steps:
Step 1: Detection of LRR proteins containing the six, novel LRRs in E-coli yddk by using FASTA
Step 2: Identification of the IRREKO@LRRs in individual LRR proteins by a new
Step 3: Iteration of these two steps using novel LRRs in newly identified LRR proteins
In step 1, we performed similarity search using the six, novel LRRs as probes by FASTA at the Bioinformatic Center, Institute for Chemical Research, Kyoto University on April 27, 2009 http://www.genome.ad.jp/. This procedure detected many yddK homologs from Escherichia coli strains and Shigella flexneri [Q0T447 and Q83R94] with significant similarity (E-values < 6.5 × 10-29). In addition, two other proteins were detected with significant similarity (E-value < 3.3 × 10-9). One is SSON_1653 that is 387 residues long [Q3Z1L5]. The other is SD1012_2081 with 163 residues [B3WXZ7]. In step 2, we performed multiple sequence alignment among their LRR domains of SSON_1653 and Sd1012_2081. SSON_1653 contains 14 LRRs and 9 of the 12 repeats consist of LxxLxLxxNxLxxL(D/N)(L/F)xxxxx where "L" is Leu, Val, or Ile. Sd1012_2081 contains 4.5 LRRs; 3.5 of these repeats consist of LxxLxLxxNxLxxIx(I/A/F)xxaxx In step 3, the above procedures were iterated to identify other LRR proteins having this IRREKO@LRR domain.
The dot-matrix comparisons were performed using the BLOSUM62 scoring matrix and a window size of 21 residues http://emboss.bioinformatics.nl/cgi-bin/emboss/dotmatcher. A radar chart is a graphical method displaying multivariate data in the form of a two-dimensional chart of three or more quantitative variables represented on axes starting from the same point http://en.wikipedia.org/wiki/Radar_chart. For a given observation, the length of each ray is the occurrence frequency of each amino acid at two positions of "IRREKO" LRR with 21 residues. Multiple sequence alignments were performed by CLUSTALW at the Bioinformatic Center. The protein secondary structure prediction was performed by SSpro4.0 http://contact.ics.uci.edu/sspro4.html and Proteus http://22.214.171.124/proteus/#. Signal sequence analysis was carried out using the program SignalP .
Mistry J, Finn R: Pfam: a domain-centric method for analyzing proteins and proteomes. Methods Mol Biol. 2007, 396: 43-58. full_text.
Kobe B, Deisenhofer J: The leucine-rich repeat: a versatile binding motif. Trends Biochem Sci. 1994, 19 (10): 415-421. 10.1016/0968-0004(94)90090-6.
Kobe B, Kajava AV: The leucine-rich repeat as a protein recognition motif. Curr Opin Struct Biol. 2001, 11 (6): 725-732. 10.1016/S0959-440X(01)00266-4.
Matsushima N, Enkhbayar P, Kamiya M, Osaki M, Kretsinger R: Leucine-Rich Repeats (LRRs): Structure, Function, Evolution and Interaction with Ligands. Drug Design Reviews. 2005, 2 (4): 305-322. 10.2174/1567269054087613.
Matsushima N, Tachi N, Kuroki Y, Enkhbayar P, Osaki M, Kamiya M, Kretsinger RH: Structural analysis of leucine-rich-repeat variants in proteins associated with human diseases. Cell Mol Life Sci. 2005, 62 (23): 2771-2791. 10.1007/s00018-005-5187-z.
Bella J, Hindle KL, McEwan PA, Lovell SC: The leucine-rich repeat structure. Cell Mol Life Sci. 2008, 65 (15): 2307-2333. 10.1007/s00018-008-8019-0.
Kajava AV: Structural diversity of leucine-rich repeat proteins. J Mol Biol. 1998, 277 (3): 519-527. 10.1006/jmbi.1998.1643.
Ohyanagi T, Matsushima N: Classification of tandem leucine-rich repeats within a great variety of proteins. FASEB J. 1997, 11: A949-
Kajava AV, Anisimova M, Peeters N: Origin and evolution of GALA-LRR, a new member of the CC-LRR subfamily: from plants to bacteria?. PLoS One. 2008, 3 (2): e1694-10.1371/journal.pone.0001694.
Torii KU: Leucine-rich repeat receptor kinases in plants: structure, function, and signal transduction pathways. Int Rev Cytol. 2004, 234: 1-46. full_text.
van der Hoorn RA, Wulff BB, Rivas S, Durrant MC, van der Ploeg A, de Wit PJ, Jones JD: Structure-function analysis of cf-9, a receptor-like protein with extracytoplasmic leucine-rich repeats. Plant Cell. 2005, 17 (3): 1000-1015. 10.1105/tpc.104.028118.
Fritz-Laylin LK, Krishnamurthy N, Tor M, Sjolander KV, Jones JD: Phylogenomic analysis of the receptor-like proteins of rice and Arabidopsis. Plant Physiol. 2005, 138 (2): 611-623. 10.1104/pp.104.054452.
Wang G, Ellendorff U, Kemp B, Mansfield JW, Forsyth A, Mitchell K, Bastas K, Liu CM, Woods-Tor A, Zipfel C: A genome-wide functional investigation into the roles of receptor-like proteins in Arabidopsis. Plant Physiol. 2008, 147 (2): 503-517. 10.1104/pp.108.119487.
Matsushima N, Mikami T, Tanaka T, Miyashita H, Yamada K, Kuroki Y: Analyses of non-leucine-rich repeat (non-LRR) regions intervening between LRRs in proteins. Biochim Biophys Acta. 2009, 1790 (10): 1217-1237.
Andrade MA, Ponting CP, Gibson TJ, Bork P: Homology-based method for identification of protein repeats using statistical significance estimates. J Mol Biol. 2000, 298 (3): 521-537. 10.1006/jmbi.2000.3684.
Lehmann P: Structure and evolution of plant disease resistance genes. J Appl Genet. 2002, 43 (4): 403-414.
Leister D: Tandem and segmental gene duplication and recombination in the evolution of plant disease resistance gene. Trends Genet. 2004, 20 (3): 116-122. 10.1016/j.tig.2004.01.007.
Hulbert SH, Webb CA, Smith SM, Sun Q: Resistance gene complexes: evolution and utilization. Annu Rev Phytopathol. 2001, 39: 285-312. 10.1146/annurev.phyto.39.1.285.
Young ND: The genetic architecture of resistance. Curr Opin Plant Biol. 2000, 3 (4): 285-290. 10.1016/S1369-5266(00)00081-9.
Ellis J, Dodds P, Pryor T: Structure, function and evolution of plant disease resistance genes. Curr Opin Plant Biol. 2000, 3 (4): 278-284. 10.1016/S1369-5266(00)00080-7.
Richter TE, Ronald PC: The evolution of disease resistance genes. Plant Mol Biol. 2000, 42 (1): 195-204. 10.1023/A:1006388223475.
Ronald PC: Resistance gene evolution. Curr Opin Plant Biol. 1998, 1 (4): 294-298. 10.1016/1369-5266(88)80049-9.
Michelmore RW, Meyers BC: Clusters of resistance genes in plants evolve by divergent selection and a birth-and-death process. Genome Res. 1998, 8 (11): 1113-1130.
Couch BC, Spangler R, Ramos C, May G: Pervasive purifying selection characterizes the evolution of I2 homologs. Mol Plant Microbe Interact. 2006, 19 (3): 288-303. 10.1094/MPMI-19-0288.
Matsushima N, Kamiya M, Suzuki N, Tanaka T: Super-motifs of leucine-rich repeats (LRRs) proteins. Genome informatics. 2000, 11: 343-345.
Matsushima N, Ohyanagi T, Tanaka T, Kretsinger RH: Super-motifs and evolution of tandem leucine-rich repeats within the small proteoglycans--biglycan, decorin, lumican, fibromodulin, PRELP, keratocan, osteoadherin, epiphycan, and osteoglycin. Proteins. 2000, 38 (2): 210-225. 10.1002/(SICI)1097-0134(20000201)38:2<210::AID-PROT9>3.0.CO;2-1.
Matsushima N, Tanaka T, Enkhbayar P, Mikami T, Taga M, Yamada K, Kuroki Y: Comparative sequence analysis of leucine-rich repeats (LRRs) within vertebrate toll-like receptors. BMC Genomics. 2007, 8: 124-10.1186/1471-2164-8-124.
Eugster M, Roten CA, Greub G: Analyses of six homologous proteins of Protochlamydia amoebophila UWE25 encoded by large GC-rich genes (lgr): a model of evolution and concatenation of leucine-rich repeats. BMC Evol Biol. 2007, 7: 231-10.1186/1471-2148-7-231.
Hirt RP, Harriman N, Kajava AV, Embley TM: A novel potential surface protein in Trichomonas vaginalis contains a leucine-rich repeat shared by micro-organisms from all three domains of life. Mol Biochem Parasitol. 2002, 125 (1-2): 195-199. 10.1016/S0166-6851(02)00211-6.
Cheng J, Randall AZ, Sweredoski MJ, Baldi P: SCRATCH: a protein structure and structural feature prediction server. Nucleic Acids Res. 2005, W72-76. 10.1093/nar/gki396. 33 Web Server
Montgomerie S, Cruz JA, Shrivastava S, Arndt D, Berjanskii M, Wishart DS: PROTEUS2: a web server for comprehensive protein structure prediction and structure-based annotation. Nucleic Acids Res. 2008, W202-209. 10.1093/nar/gkn255. 36 Web Server
Enkhbayar P, Kamiya M, Osaki M, Matsumoto T, Matsushima N: Structural principles of leucine-rich repeat (LRR) proteins. Proteins. 2004, 54 (3): 394-403. 10.1002/prot.10605.
Jenkins J, Mayans O, Pickersgill R: Structure and evolution of parallel beta-helix proteins. J Struct Biol. 1998, 122 (1-2): 236-246. 10.1006/jsbi.1998.3985.
Jenkins J, Pickersgill R: The architecture of parallel beta-helices and related folds. Prog Biophys Mol Biol. 2001, 77 (2): 111-175. 10.1016/S0079-6107(01)00013-X.
Kobe B, Kajava AV: When protein folding is simplified to protein coiling: the continuum of solenoid protein structures. Trends Biochem Sci. 2000, 25 (10): 509-515. 10.1016/S0968-0004(00)01667-4.
Baumann U: Crystal structure of the 50 kDa metallo protease from Serratia marcescens. J Mol Biol. 1994, 242 (3): 244-251. 10.1006/jmbi.1994.1576.
Kim HM, Park BS, Kim JI, Kim SE, Lee J, Oh SC, Enkhbayar P, Matsushima N, Lee H, Yoo OJ: Crystal structure of the TLR4-MD-2 complex with bound endotoxin antagonist Eritoran. Cell. 2007, 130 (5): 906-917. 10.1016/j.cell.2007.08.002.
Jin MS, Kim SE, Heo JY, Lee ME, Kim HM, Paik SG, Lee H, Lee JO: Crystal structure of the TLR1-TLR2 heterodimer induced by binding of a tri-acylated lipopeptide. Cell. 2007, 130 (6): 1071-1082. 10.1016/j.cell.2007.09.008.
Bendtsen JD, Nielsen H, von Heijne G, Brunak S: Improved prediction of signal peptides: SignalP 3.0. J Mol Biol. 2004, 340 (4): 783-795. 10.1016/j.jmb.2004.05.028.
We thank Dr. Robert H. Kretsinger of the University of Virginia for his valuable suggestion and comments.
NM (corresponding author) carried out the molecular genetic studies, participated in the sequence alignment and drafted the manuscript. HM performed dot plot analysis and radar chart analysis. TM contributed to the data analysis including the sequence alignment. KY conceived of the study, and participated in its design and coordination. All authors read and approved the final manuscript.