Protein folds in the all-beta and all-alpha classes

Interaction of negative regulator OsWD40-193 with OseEF1A1 inhibits Oryza sativa resistance to Hirschmanniella mucronata infection

Rice is a crucial food crop worldwide, but it is highly susceptible to Hirschmanniella mucronata, a migratory parasitic nematode. No rice variety has been identified that could resist H. mucronata infection. Therefore, it is very important to study the interaction between rice and H. mucronata to breed resistant rice varieties. Here, we demonstrated that protein OsWD40-193 interacted with the extension factor OseEF1A1 and both were negative regulators inhibiting rice resistance to H. mucronata infection. Overexpression of either OsWD40-193 or OseEF1A1 led to enhance susceptibility to H. mucronata, whereas the absence of OsWD40-193 or OseEF1A1 led to resistance. Further transcriptomic analysis showed that OseEF1A1 deletion altered the expression of genes association with salicylic acid, jasmonic acid and abolic acid signaling pathways and increased the accumulation of secondary metabolites to enhance resistance in rice. Our study showed that H. mucronata infection affected the expression of negative regulators in rice and inhibited rice resistance, which was conducive to the infection of nematode. Together, our data showed that H. mucronata affected the expression of negative regulators to facilitate its infection and provided potential target genes to engineering resistance germplasm via gene editing of the negative regulators.

Computational design of a cyclic peptide that inhibits the CTLA4 immune checkpoint

Proteins involved in immune checkpoint pathways, such as CTLA4, PD1, and PD-L1, have become important targets for cancer immunotherapy; however, development of small molecule drugs targeting these pathways has proven difficult due to the nature of their protein-protein interfaces. Here, using a hierarchy of computational techniques, we design a cyclic peptide that binds CTLA4 and follow this with experimental verification of binding and biological activity, using bio-layer interferometry, cell culture, and a mouse tumor model. Beginning from a template excised from the X-ray structure of the CTLA4:B7-2 complex, we generate several peptide sequences using flexible docking and modeling steps. These peptides are cyclized head-to-tail to improve structural and proteolytic stability and screened using molecular dynamics simulation and MM-GBSA calculation. The standard binding free energies for shortlisted peptides are then calculated in explicit-solvent simulation using a rigorous multistep technique. The most promising peptide, cyc(EIDTVLTPTGWVAKRYS), yields the standard free energy -6.6 ± 3.5 kcal mol-1, which corresponds to a dissociation constant of ∼15 μmol L-1. The binding affinity of this peptide for CTLA4 is measured experimentally (31 ± 4 μmol L-1) using bio-layer interferometry. Treatment with this peptide inhibited tumor growth in a co-culture of Lewis lung carcinoma (LLC) cells and antigen primed T cells, as well as in mice with an orthotropic Lewis lung carcinoma allograft model.

SPACA6 ectodomain structure reveals a conserved superfamily of gamete fusion-associated proteins

SPACA6 is a sperm-expressed surface protein that is critical for gamete fusion during mammalian sexual reproduction. Despite this fundamental role, little is known about how SPACA6 specifically functions. We elucidated the crystal structure of the SPACA6 ectodomain at 2.2-Å resolution, revealing a two-domain protein containing a four-helix bundle and Ig-like β-sandwich connected via a quasi-flexible linker. This structure is reminiscent of IZUMO1, another gamete fusion-associated protein, making SPACA6 and IZUMO1 founding members of a superfamily of fertilization-associated proteins, herein dubbed the IST superfamily. The IST superfamily is defined structurally by its distorted four-helix bundle and a pair of disulfide-bonded CXXC motifs. A structure-based search of the AlphaFold human proteome identified more protein members to this superfamily; remarkably, many of these proteins are linked to gamete fusion. The SPACA6 structure and its connection to other IST-superfamily members provide a missing link in our knowledge of mammalian gamete fusion.

Genome-Wide Analysis of Anthocyanin Biosynthesis Regulatory WD40 Gene FcTTG1 and Related Family in Ficus carica L

WD40 proteins serve as crucial regulators in a broad spectrum of plant developmental and physiological processes, including anthocyanin biosynthesis. However, in fig (Ficus carica L.), neither the WD40 family nor any member involved in anthocyanin biosynthesis has been elucidated. In the present study, 204 WD40 genes were identified from the fig genome and phylogenetically classified into 5 clusters and 12 subfamilies. Bioinformatics analysis prediction localized 109, 69, and 26 FcWD40 proteins to the cytoplasm, nucleus and other cellular compartments, respectively. RNA-seq data mining revealed 127 FcWD40s expressed at FPKM > 10 in fig fruit. Most of these genes demonstrated higher expression in the early stages of fruit development. FcWD40-97 was recruited according to three criteria: high expression in fig fruit, predicted nuclear localization, and closest clustering with TTG1s identified in other plants. FcWD40-97, encoding 339 amino acids including 5 WD-repeat motifs, showed 88.01 and 87.94% amino acid sequence similarity to apple and peach TTG1, respectively. The gene is located on fig chromosome 4, and is composed of 1 intron and 2 exons. Promoter analysis revealed multiple light-responsive elements, one salicylic acid-responsive element, three methyl jasmonate-responsive elements, and one MYB-binding site involved in flavonoid biosynthesis gene regulation. FcWD40-97 was in the FPKM > 100 expression level group in fig fruit, and higher expression was consistently found in the peel compared to the flesh at the same development stages. Expression level did not change significantly under light deprivation, whereas in leaves and roots, its expression was relatively low. Transient expression verified FcWD40-97's localization to the nucleus. Yeast two-hybrid (Y2H) and biomolecular fluorescence complementation (BiFC) assays revealed that FcWD40-97 interacts with FcMYB114, FcMYB123, and FcbHLH42 proteins in vitro and in vivo, showing that FcWD40-97 functions as a member of the MYB-bHLH-WD40 (MBW) complex in anthocyanin-biosynthesis regulation in fig. We therefore renamed FcWD40-97 as FcTTG1. Our results provide the first systematic analysis of the FcWD40 family and identification of FcTTG1 in fig pigmentation.

Genome Wide Identification and Characterization of Apple WD40 Proteins and Expression Analysis in Response to ABA, Drought, and Low Temperature

Basic WD40 proteins, which are characterized by a conserved WD40 domain, comprise a superfamily of regulatory proteins in plants and play important roles in plant growth and development. However, WD40 genes have been rarely studied in apple (Malus × domestica Borkh.). In this study, 346 WD40 genes classified in 12 subfamilies, were identified in the apple genome. Evolutionary analysis of WD40 proteins in apple and Arabidopsis revealed that the genes were classifiable into 14 groups, and the exon/intron structure of each group showed a similar structure. Analysis of collinearity showed that the large-scale amplification of WD40 genes in apple was largely attributable to recent whole-genome replication events. Nineteen candidate stress-related genes, selected by GO annotation and comparison with Arabidopsis homologs, showed different expression profiles in six organs at different developmental stages in response to exogenous abscisic acid (ABA), drought, and low temperature. Eight genes (MdWD40-17, 24, 70, 74, 219, 256, 283, and 307) showed a distinct response to one or more treatments (ABA, drought, and low temperature) as indicated by quantitative real-time PCR analysis. Taken together, these data provide rich resources for further study of MdWD40 genes and their potential roles in stress responses in apple.

Genome-Wide Identification and Characterization of Polygalacturonase Gene Family in Maize (Zea mays L.)

Polygalacturonase (PG, EC 3.2.1.15) is a crucial enzyme for pectin degradation and is involved in various developmental processes such as fruit ripening, pollen development, cell expansion, and organ abscission. However, information on the PG gene family in the maize (Zea mays L.) genome and the specific members involved in maize anther development are still lacking. In this study, we identified 55 PG family genes from the maize genome and further characterized their evolutionary relationship and expression patterns. Phylogenetic analysis revealed that ZmPGs are grouped into six Clades, and gene structures of the same Clade are highly conserved, suggesting their functional conservation. The ZmPGs are randomly distributed across maize chromosomes, and collinearity analysis showed that many ZmPGs might be derived from tandem duplications and segmental duplications, and these genes are under purifying selection. Furthermore, gene expression analysis provided insights into possible functional divergence among ZmPGs. Based on the RNA-seq data analysis, we found that many ZmPGs are expressed in various tissues while 18 ZmPGs are highly expressed in maize anther, and their detailed expression profiles in different anther developmental stages were further investigated by using RT-qPCR analysis. These results provide valuable information for further functional characterization and application of the ZmPGs in maize.

OsTTG1, a WD40 repeat gene, regulates anthocyanin biosynthesis in rice

Anthocyanins play an important role in the growth of plants, and are beneficial to human health. In plants, the MYB-bHLH-WD40 (MBW) complex activates the genes for anthocyanin biosynthesis. However, in rice, the WD40 regulators remain to be conclusively identified. Here, a crucial anthocyanin biosynthesis gene was fine mapped to a 43.4-kb genomic region on chromosome 2, and a WD40 gene OsTTG1 (Oryza sativa TRANSPARENT TESTA GLABRA1) was identified as ideal candidate gene. Subsequently, a homozygous mutant (osttg1) generated by CRISPR/Cas9 showed significantly decreased anthocyanin accumulation in various rice organs. OsTTG1 was highly expressed in various rice tissues after germination, and it was affected by light and temperature. OsTTG1 protein was localized to the nucleus, and can physically interact with Kala4, OsC1, OsDFR and Rc. Furthermore, a total of 59 hub transcription factor genes might affect rice anthocyanin biosynthesis, and LOC_Os01g28680 and LOC_Os02g32430 could have functional redundancy with OsTTG1. Phylogenetic analysis indicated that directional selection has driven the evolutionary divergence of the indica and japonica OsTTG1 alleles. Our results suggest that OsTTG1 is a vital regulator of anthocyanin biosynthesis, and an important gene resource for the genetic engineering of anthocyanin biosynthesis in rice and other plants.

eIF3i regulation of protein synthesis, cell proliferation, cell cycle progression, and tumorigenesis

eIF3i, a 36-kDa protein, is a putative subunit of the eIF3 complex important for translation initiation of mRNAs. It is a WD40 domain-containing protein with seven WD40 repeats that forms a β-propeller structure with an important function in pre-initiation complex formation and mRNA translation initiation. In addition to participating in the eIF3 complex formation for global translational control, eIF3i may bind to specific mRNAs and regulate their translation individually. Furthermore, eIF3i has been shown to bind to TGF-β type II receptor and participate in TGF-β signaling. It may also participate in and regulate other signaling pathways including Wnt/β-catenin pathway via translational regulation of COX-2 synthesis. These multiple canonical and noncanonical functions of eIF3i in translational control and in regulating signal transduction pathways may be responsible for its role in cell differentiation, cell cycle regulation, proliferation, and tumorigenesis. In this review, we will critically evaluate recent progresses and assess future prospects in studying eIF3i.

Essential role of calcium in extending RTX adhesins to their target

•

Elongated beta-sandwich repeats are a major part of bacterial RTX adhesins.

•

The repeats are arranged in tandem to extend away from the bacterial surface.

•

Calcium ions are coordinated in the linkers between repeats to stiffen the protein.

•

Rigidification of the tandem repeats further helps extension of the adhesin.

•

The repeats differ greatly between species, but all have Ca²⁺ in their linkers.

RTX adhesins are long, multi-domain proteins present on the outer membrane of many Gram-negative bacteria. From this vantage point, adhesins use their distal ligand-binding domains for surface attachment leading to biofilm formation. To expand the reach of the ligand-binding domains, RTX adhesins maintain a central extender region of multiple tandem repeats, which makes up most of the proteins’ large molecular weight. Alignments of the 10-15-kDa extender domains show low sequence identity between adhesins. Here we have produced and structurally characterized protein constructs of four tandem repeats (tetra-tandemers) from two different RTX adhesins. In comparing the tetra-tandemers to each other and already solved structures from Marinomonas primoryensis and Salmonella enterica, the extender domains fold as diverse beta-sandwich structures with widely differing calcium contents. However, all the tetra-tandemers have at least one calcium ion coordinated in the linker region between beta-sandwich domains whose role appears to be the rigidification of the extender region to help the adhesin extend its reach.

Systematic analysis of SmWD40s, and responding of SmWD40-170 to drought stress by regulation of ABA- and H2O2-induced stomal movement in Salvia miltiorrhiza bunge

WD40 proteins play crucial roles in response to abiotic stress. By screening the genome sequences of Salvia miltiorrhiza Bunge, 225 SmWD40 genes were identified and divided into 9 subfamilies (I-IX). Physiological, biochemical, gene structure, conserved protein motif and GO annotation analyses were performed on SmWD40 family members. The SmWD40-170 was found in 110 SmWD40 genes that contain drought response elements, SmWD40-170 was one of these genes whose response in terms of expression under drought was significant. The expression of SmWD40-170 was also up-regulated by ABA and H₂O₂. Through observed the stomatal phenotype of SmWD40-170 transgenic lines, the stomatal closure was abolished under dehydration, ABA and H₂O₂ treatment in SmWD40-170 knockdown lines. Abscisic acid (ABA), as the key phytohormone, elevates reactive oxygen species (ROS) levels under drought stress. The ABA-ROS interaction mediated the generation of H₂O₂ and the activation of anion channel in guard cells. The osmolality alteration of guard cells further accelerated the stomatal closure. As a second messenger, nitric oxide (NO) regulated ABA signaling, the NO stimulated protein kinase activity inhibited the K⁺ influx which result in stomatal closure. These NO-relevant events were essential for ABA-induced stomatal closure. The reduction of NO production was also observed in the guard cells of SmWD40-170 knockdown lines. The abolished of stomatal closure attributed to the SmWD40-170 deficiency induced the reduction of NO content. In general, the SmWD40-170 is a critical drought response gene in SmWD40 gene family and regulates ABA- and H₂O₂-induced stomatal movement by affecting the synthesis of NO.

Quadruply Stranded Metallo-Supramolecular Helicate [Pd2(hextrz)4]4+ Acts as a Molecular Mimic of Cytolytic Peptides

[Pd₂(hextrz)₄]⁴⁺ is a quadruply stranded helicate, a novel bioinorganic complex designed to mimic the structure and function of proteins due to its high stability and supramolecular size. We have previously reported that [Pd₂(hextrz)₄]⁴⁺ exhibited cytotoxicity toward a range of cell lines, with IC₅₀ values ranging from 3 to 10 μM. Here we demonstrate that [Pd₂(hextrz)₄]⁴⁺ kills cells by forming pores within the cell membrane, a mechanism of cell death analogous to the naturally occurring cytolytic peptides. [Pd₂(hextrz)₄]⁴⁺ induced cell death is characterized by an initial influx of Ca²⁺, followed by nuclear condensation and mitochondrial swelling. This is accompanied by progressive cell membrane damage that results in the formation of large blebs at the cell surface. This allows the efflux of molecules from the cell leading to loss of cell viability. These data suggest that it may be possible to design metallo-supramolecular complexes to mimic the cytotoxic action of pore forming proteins and peptides and so provide a new class of drug to treat cancer, autoimmune disorders, and microbial infection.

Simulated Epidemics in 3D Protein Structures to Detect Functional Properties

The outcome of an epidemic is closely related to the network of interactions between individuals. Likewise, protein functions depend on the 3D arrangement of their residues and the underlying energetic interaction network. Borrowing ideas from the theoretical framework that has been developed to address the spreading of real diseases, we study for the first time the diffusion of a fictitious epidemic inside the protein nonbonded interaction network, aiming to study network features and properties. Our approach allows us to probe the overall stability and the capability of propagating information in complex 3D structures, proving to be very efficient in addressing different problems, from the assessment of thermal stability to the identification of functional sites.

An exhaustive survey of regular peptide conformations using a new metric for backbone handedness (h)

The Ramachandran plot is important to structural biology as it describes a peptide backbone in the context of its dominant degrees of freedom—the backbone dihedral angles φ and ψ (Ramachandran, Ramakrishnan & Sasisekharan, 1963). Since its introduction, the Ramachandran plot has been a crucial tool to characterize protein backbone features. However, the conformation or twist of a backbone as a function of φ and ψ has not been completely described for both cis and trans backbones. Additionally, little intuitive understanding is available about a peptide’s conformation simply from knowing the φ and ψ values of a peptide (e.g., is the regular peptide defined by φ = ψ =  − 100°  left-handed or right-handed?). This report provides a new metric for backbone handedness (h) based on interpreting a peptide backbone as a helix with axial displacement d and angular displacement θ, both of which are derived from a peptide backbone’s internal coordinates, especially dihedral angles φ, ψ and ω. In particular, h equals sin(θ)d∕|d|, with range [−1, 1] and negative (or positive) values indicating left(or right)-handedness. The metric h is used to characterize the handedness of every region of the Ramachandran plot for both cis (ω = 0°) and trans (ω = 180°) backbones, which provides the first exhaustive survey of twist handedness in Ramachandran (φ, ψ) space. These maps fill in the ‘dead space’ within the Ramachandran plot, which are regions that are not commonly accessed by structured proteins, but which may be accessible to intrinsically disordered proteins, short peptide fragments, and protein mimics such as peptoids. Finally, building on the work of (Zacharias & Knapp, 2013), this report presents a new plot based on d and θ that serves as a universal and intuitive alternative to the Ramachandran plot. The universality arises from the fact that the co-inhabitants of such a plot include every possible peptide backbone including cis and trans backbones. The intuitiveness arises from the fact that d and θ provide, at a glance, numerous aspects of the backbone including compactness, handedness, and planarity.

A WD40 protein, AtGHS40, negatively modulates abscisic acid degrading and signaling genes during seedling growth under high glucose conditions

The Arabidopsis thaliana T-DNA insertion mutant glucose hypersensitive (ghs) 40-1 exhibited hypersensitivity to glucose (Glc) and abscisic acid (ABA). The ghs40-1 mutant displayed severely impaired cotyledon greening and expansion and showed enhanced reduction in hypocotyl elongation of dark-grown seedlings when grown in Glc concentrations higher than 3 %. The Glc-hypersensitivity of ghs40-1 was correlated with the hyposensitive phenotype of 35S::AtGHS40 seedlings. The phenotypes of ghs40-1 were recovered by complementation with 35S::AtGHS40. The AtGHS40 (At5g11240) gene encodes a WD40 protein localized primarily in the nucleus and nucleolus using transient expression of AtGHS40-mRFP in onion cells and of AtGHS40-EGFP and EGFP-AtGHS40 in Arabidopsis protoplasts. The ABA biosynthesis inhibitor fluridone extensively rescued Glc-mediated growth arrest. Quantitative real time-PCR analysis showed that AtGHS40 was involved in the control of Glc-responsive genes. AtGHS40 acts downstream of HXK1 and is activated by ABI4 while ABI4 expression is negatively modulated by AtGHS40 in the Glc signaling network. However, AtGHS40 may not affect ABI1 and SnRK2.6 gene expression. Given that AtGHS40 inhibited ABA degrading and signaling gene expression levels under high Glc conditions, a new circuit of fine-tuning modulation by which ABA and ABA signaling gene expression are modulated in balance, occurred in plants. Thus, AtGHS40 may play a role in ABA-mediated Glc signaling during early seedling development. The biochemical function of AtGHS40 is also discussed.

The Ramachandran Number: An Order Parameter for Protein Geometry

Three-dimensional protein structures usually contain regions of local order, called secondary structure, such as α-helices and β-sheets. Secondary structure is characterized by the local rotational state of the protein backbone, quantified by two dihedral angles called ϕ and ψ. Particular types of secondary structure can generally be described by a single (diffuse) location on a two-dimensional plot drawn in the space of the angles ϕ and ψ, called a Ramachandran plot. By contrast, a recently-discovered nanomaterial made from peptoids, structural isomers of peptides, displays a secondary-structure motif corresponding to two regions on the Ramachandran plot [Mannige et al., Nature 526, 415 (2015)]. In order to describe such ‘higher-order’ secondary structure in a compact way we introduce here a means of describing regions on the Ramachandran plot in terms of a single Ramachandran number, R, which is a structurally meaningful combination of ϕ and ψ. We show that the potential applications of R are numerous: it can be used to describe the geometric content of protein structures, and can be used to draw diagrams that reveal, at a glance, the frequency of occurrence of regular secondary structures and disordered regions in large protein datasets. We propose that R might be used as an order parameter for protein geometry for a wide range of applications.

Cloning and characterization of a potato StAN11 gene involved in anthocyanin biosynthesis regulation

Anthocyanins are a class of products of plant secondary metabolism and are responsible for tubers color in potato. The biosynthesis of anthocyanins is a complex biological process, in which multiple genes are involved including structural genes and regulatory genes. In this study, StAN11, a WD40-repeat gene, was cloned from potato cultivar Chieftain (Solanum tuberosum L.). StAN11 (HQ599506) contained no intron and its open reading frame (ORF) was 1,029 bp long, encoding a putative protein of 342 amino acids. In order to verify its role in anthocyanin biosynthesis, StAN11 was inserted behind the CaMV-35S promoter of pCMBIA1304 and the recombination vector was introduced into the potato cultivar Désirée plants by Agrobacterium-mediated transformation. The color of transgenic tuber skin was significantly deepened, compared to the wild-type control, which was highly consistent with the accumulation of anthocyanin and expression of StAN11 in transgenic lines tuber skin. Further analysis on the expression of Flavonone-3-hydroxylase (F3H), Dihydroflavonol reductase (DFR), Anthocyanidin synthase (ANS), and Flavonoid 3-O-glucosyl transferase (3GT) in transgenic plants revealed that only DFR was upregulated. This result suggested that StAN11 regulated anthocyanin biosynthesis in potato by controlling DFR expression and accumulation of anthocyanin could be increased through overexpression of StAN11 in the tubers with the genetic background of anthocyanin biosynthesis.

Origination of the Protein Fold Repertoire from Oily Pluripotent Peptides

While the repertoire of protein folds that exists today underlies most of life’s capabilities, our mechanistic picture of protein fold origination is incomplete. This paper discusses a hypothetical mechanism for the emergence of the protein fold repertoire from highly dynamic and collapsed peptides, exemplified by peptides with high oil content or hydrophobicity. These peptides are called pluripotent to emphasize their capacity to evolve into numerous folds transiently available to them. As evidence, the paper will discuss previous simulation work on the superior fold evolvability of oily peptides, trace (“fossil”) evidence within proteomes seen today, and a general relationship between protein dynamism and evolvability. Aside from implications on the origination of protein folds, the hypothesis implies that the vanishing utility of a random peptide in protein origination may be relatively exaggerated, as some random peptides with a certain composition (e.g., oily) may fare better than others. In later sections, the hypothesis is discussed in the context of existing discussions regarding the spontaneous origination of biomolecules.

Dynamic New World: Refining Our View of Protein Structure, Function and Evolution

Proteins are crucial to the functioning of all lifeforms. Traditional understanding posits that a single protein occupies a single structure ("fold"), which performs a single function. This view is radically challenged with the recognition that high structural dynamism-the capacity to be extra "floppy"-is more prevalent in functional proteins than previously assumed. As reviewed here, this dynamic take on proteins affects our understanding of protein "structure", function, and evolution, and even gives us a glimpse into protein origination. Specifically, this review will discuss historical developments concerning protein structure, and important new relationships between dynamism and aspects of protein sequence, structure, binding modes, binding promiscuity, evolvability, and origination. Along the way, suggestions will be provided for how key parts of textbook definitions-that so far have excluded membership to intrinsically disordered proteins (IDPs)-could be modified to accommodate our more dynamic understanding of proteins.

Amino acid distribution rules predict protein fold

In the present article, we provide a brief overview of the main approaches to analysing the sequence-structure relationship of proteins and outline a novel method of structure prediction. The proposed method involves finding a set of rules that describes a correlation between the distribution of residues in a sequence and the essential structural characteristics of a protein structure. The residue distribution rules specify the 'favourable' residues that are required in certain positions of a polypeptide chain in order for it to assume a particular protein fold, and the 'unfavourable' residues incompatible with the given fold. Identification of amino acid distribution rules derives from examination of inter-residue contacts. We describe residue distribution rules for a large group of β-sandwich-like proteins characterized by a specific arrangement of strands in their two β-sheets. It was shown that this method has very high accuracy (approximately 85%). The advantage of the residue rule approach is that it makes possible prediction of protein folding even in polypeptide chains that have very low global sequence similarities, as low as 18%. Another potential benefit is that a better understanding of which residues play essential roles in a given protein fold may facilitate rational protein engineering design.

A method to predict edge strands in beta-sheets from protein sequences

There is a need for rules allowing three-dimensional structure information to be derived from protein sequences. In this work, consideration of an elementary protein folding step allows protein sub-sequences which optimize folding to be derived for any given protein sequence. Classical mechanics applied to this system and the energy conservation law during the elementary folding step yields an equation whose solutions are taken over the field of rational numbers. This formalism is applied to beta-sheets containing two edge strands and at least two central strands. The number of protein sub-sequences optimized for folding per amino acid in beta-strands is shown in particular to predict edge strands from protein sequences. Topological information on beta-strands and loops connecting them is derived for protein sequences with a prediction accuracy of 75%. The statistical significance of the finding is given. Applications in protein structure prediction are envisioned such as for the quality assessment of protein structure models.

ST proteins, a new family of plant tandem repeat proteins with a DUF2775 domain mainly found in Fabaceae and Asteraceae

BACKGROUND

Many proteins with tandem repeats in their sequence have been described and classified according to the length of the repeats: I) Repeats of short oligopeptides (from 2 to 20 amino acids), including structural cell wall proteins and arabinogalactan proteins. II) Repeats that range in length from 20 to 40 residues, including proteins with a well-established three-dimensional structure often involved in mediating protein-protein interactions. (III) Longer repeats in the order of 100 amino acids that constitute structurally and functionally independent units. Here we analyse ShooT specific (ST) proteins, a family of proteins with tandem repeats of unknown function that were first found in Leguminosae, and their possible similarities to other proteins with tandem repeats.

RESULTS

ST protein sequences were only found in dicotyledonous plants, limited to several plant families, mainly the Fabaceae and the Asteraceae. ST mRNAs accumulate mainly in the roots and under biotic interactions. Most ST proteins have one or several Domain(s) of Unknown Function 2775 (DUF2775). All deduced ST proteins have a signal peptide, indicating that these proteins enter the secretory pathway, and the mature proteins have tandem repeat oligopeptides that share a hexapeptide (E/D)FEPRP followed by 4 partially conserved amino acids, which could determine a putative N-glycosylation signal, and a fully conserved tyrosine. In a phylogenetic tree, the sequences clade according to taxonomic group. A possible involvement in symbiosis and abiotic stress as well as in plant cell elongation is suggested, although different STs could play different roles in plant development.

CONCLUSIONS

We describe a new family of proteins called ST whose presence is limited to the plant kingdom, specifically to a few families of dicotyledonous plants. They present 20 to 40 amino acid tandem repeat sequences with different characteristics (signal peptide, DUF2775 domain, conservative repeat regions) from the described group of 20 to 40 amino acid tandem repeat proteins and also from known cell wall proteins with repeat sequences. Several putative roles in plant physiology can be inferred from the characteristics found.

SHEEP: A TOOL FOR DESCRIPTION OF β-SHEETS IN PROTEIN 3D STRUCTURES

The description of a protein fold is a hard problem due to significant variability of main structural units, β-sheets and α-helixes, and their mutual arrangements. An adequate description of the structural units is an important step in objective protein structure classification, which to date is based on expert judgment in a number of cases. Explicit determination and description of structural units is more complicated for β-sheets than for α-helixes due to β-sheets variability both in composition and geometry. We have developed an algorithm that can significantly modify β-sheets detected by commonly used DSSP and Stride algorithms and represent the result as a “β-sheet map,” a table describing certain β-sheet features. In our approach, β-sheets (rather than β-strands) are considered as holistic objects. Both hydrogen bonds and geometrical restrains are explored for the determination of β-sheets. The algorithm is implemented in SheeP program. It was tested for prediction architectures of domains from 93 well-defined all-β and α/β SCOP protein domain families, and showed 93% of correct results. The Web-service http://mouse.belozersky.msu.ru/sheep allows to detect β-sheets in a given protein structure, visualize β-sheet maps, as well as input three-dimensional structures with highlighted β-sheets and their structural features.

Genomic survey, expression profile and co-expression network analysis of OsWD40 family in rice

Background

WD40 proteins represent a large family in eukaryotes, which have been involved in a broad spectrum of crucial functions. Systematic characterization and co-expression analysis of OsWD40 genes enable us to understand the networks of the WD40 proteins and their biological processes and gene functions in rice.

Results

In this study, we identify and analyze 200 potential OsWD40 genes in rice, describing their gene structures, genome localizations, and evolutionary relationship of each member. Expression profiles covering the whole life cycle in rice has revealed that transcripts of OsWD40 were accumulated differentially during vegetative and reproductive development and preferentially up or down-regulated in different tissues. Under phytohormone treatments, 25 OsWD40 genes were differentially expressed with treatments of one or more of the phytohormone NAA, KT, or GA3 in rice seedlings. We also used a combined analysis of expression correlation and Gene Ontology annotation to infer the biological role of the OsWD40 genes in rice. The results suggested that OsWD40 genes may perform their diverse functions by complex network, thus were predictive for understanding their biological pathways. The analysis also revealed that OsWD40 genes might interact with each other to take part in metabolic pathways, suggesting a more complex feedback network.

Conclusions

All of these analyses suggest that the functions of OsWD40 genes are diversified, which provide useful references for selecting candidate genes for further functional studies.

Structure, stability, and mobility of a lyophilized IgG1 monoclonal antibody as determined using second-derivative infrared spectroscopy

There are many aspects of stabilization of lyophilized proteins. Of these various factors, retention of native structure, having sufficient amount of stabilizer to embed the protein within an amorphous matrix, and dampening β-relaxations have been shown to be critical in optimizing protein stability during storage. In this study, an IgG1 was lyophilized with varying amounts of sucrose. In some formulations, a small amount of sorbitol was added as a plasticizer. The structure of the protein in dried state was monitored using infrared (IR) spectroscopy. The IR spectra indicated increasing retention of the native structure, which correlated with stability as indicated by size-exclusion chromatography as well as micro-flow imaging. Maximal stability was achieved with a 2:1 mass ratio of sucrose to protein, which is more than that would be expected based on earlier studies. Analysis of both high and low frequency bands associated with intramolecular β-sheet structure provides additional information on the structure of antibodies in the solid state. Finally, there is a correlation between the bandwidth of the β-sheet bands and the enthalpy of relaxation, suggesting that amide I bands can provide some indication of the degree of coupling to the sugar matrix, as well as structural heterogeneity of the protein.

Structural classification of proteins and structural genomics: new insights into protein folding and evolution

During the past decade, the Protein Structure Initiative (PSI) centres have become major contributors of new families, superfamilies and folds to the Structural Classification of Proteins (SCOP) database. The PSI results have increased the diversity of protein structural space and accelerated our understanding of it. This review article surveys a selection of protein structures determined by the Joint Center for Structural Genomics (JCSG). It presents previously undescribed β-sheet architectures such as the double barrel and spiral β-roll and discusses new examples of unusual topologies and peculiar structural features observed in proteins characterized by the JCSG and other Structural Genomics centres.

A surface loop directs conformational switching of a lipoyl domain between a folded and a novel misfolded structure

A prominent surface loop links the first two beta strands of the lipoyl domain (E2plip) from the pyruvate dehydrogenase multienzyme complex of Escherichia coli. We show here that shortening this loop by two residues generates a protein that populates two structurally distinct stable conformers: an active, native-like monomer (HM) and a functionally compromised misfolded dimer (LM). Conversion of LM to HM was observed after exposure to temperatures above 50 degrees C. Removal of two additional residues from the loop caused the protein to adopt exclusively the misfolded conformation. Detailed NMR structural studies of the misfolded dimer reveal that the N-terminal half of the domain was unfolded and dynamic, whereas the C-terminal halves of two monomers had associated to form a structure with two-fold symmetry and a topology mimicking that of the folded monomer. The surface loop is therefore a hitherto unsuspected determinant in the folding process that leads to a functional protein.

Protein Structure Hierarchy

The complexity of protein structures requires a description of their structural components. This chapter describes the elements of protein secondary structure — regular local structural patterns — such as helices, sheets, turns, and loops. Helices and sheets tend to fall into specific regions in the {ϕ, ψ} space of the Ramachandran plot (see Figures 28 and 29). The corresponding width and shape of each region reflects the spread of that motif as found in proteins.

Study on β-sheet packing, stabilized by aromatic interactions in a tri-peptide crystal

The role of unconventional aromatic interactions in the β-sheet packing, inside the crystal, has been highlighted. Crystal structure of terminally blocked tri-peptide, Z-L-Ala-L-Ala-L-Leu-pNA has been determined. There are four different conformers of tri-peptide inside the unit-cell, in space group P1. All the four independent molecules are described by semi-extended backbone conformations which form an anti-parallel β-sheet. The inter β-sheet packing is predominantly stabilized by C—H…π and π…π interactions. In π…π interactions, the center-to-center distance between aromatic rings varies from 3.8 Å to 4.6 Å while the closest distance of approach ranges from 3.4 Å to 3.8 Å. The associations of aromatic-aromatic rings are described by either face-to-face or inclined arrangements.

A free-energy approach for all-atom protein simulation

All-atom free-energy methods offer a promising alternative to kinetic molecular mechanics simulations of protein folding and association. Here we report an accurate, transferable all-atom biophysical force field (PFF02) that stabilizes the native conformation of a wide range of proteins as the global optimum of the free-energy landscape. For 32 proteins of the ROSETTA decoy set and six proteins that we have previously folded with PFF01, we find near-native conformations with an average backbone RMSD of 2.14 A to the native conformation and an average Z-score of -3.46 to the corresponding decoy set. We used nonequilibrium sampling techniques starting from completely extended conformations to exhaustively sample the energy surface of three nonhomologous hairpin-peptides, a three-stranded beta-sheet, the all-helical 40 amino-acid HIV accessory protein, and a zinc-finger beta beta alpha motif, and find near-native conformations for the minimal energy for each protein. Using a massively parallel evolutionary algorithm, we also obtain a near-native low-energy conformation for the 54 amino-acid engrailed homeodomain. Our force field thus stabilized near-native conformations for a total of 20 proteins of all structure classes with an average RMSD of only 3.06 A to their respective experimental conformations.

SRWD: a novel WD40 protein subfamily regulated by salt stress in rice (OryzasativaL.)

By analysis with microarray data, we found that a gene encoding a novel protein containing five WD40 repeats, was regulated by salt stress in rice and named as SRWD1 (Salt responsive WD40 protein 1). By database searching, additional four SRWD1-like genes (SRWD2-SRWD5) were found in rice genome, and these five SRWD genes formed a novel WD40 subfamily. Phylogenetic analysis showed that plant SRWD proteins divided into four groups. The significant functional divergences during SRWD evolution were found. The tissue-specific and salt responsive expression profiling for SRWD genes was investigated based on microarray data. It was found that all five SRWD genes in rice were regulated by salt stress. Further, we found that SRWD1 was regulated with different patterns by salt stress in two rice cultivars responding differently to salt stress. Our study correlates WD40 proteins with salt stress in plants and provides fundamental information for the further investigation of plant SRWD proteins.

De novo design of beta-helical polypeptides

Many proteins, including several growth factor receptors such as the IGF-1R and EGFR family, contain variants of the beta-helix fold. Inspection of the irregular protein beta-helices suggested that different families of regular beta-helical polypeptides can be designed using a series of hinged vectors and the constraints imposed by the geometry of a peptide backbone. We have conceived beta-helices with five and six beta-strands per turn and designed, in detail, a series of regular beta-helices with rhomboidal or triangular cross-sections. Each beta-helix was modeled by threading C(alpha) atoms to follow the vectorial beta-helix and then creating the H-bonded polypeptide backbone and appropriate side-chain orientations. The conformational stability of these regular beta-helices was assessed using molecular dynamics simulations. Several potential repeat amino acid sequences were identified for different geometries of beta-helix. Regular beta-helices offer new possibilities for the study of protein folding, the production of nanofibers, catalysts, inhibitors of growth factor receptors and drug carriers.

Molecular dynamics studies of hexamers of amyloid-beta peptide (16-35) and its mutants: influence of charge states on amyloid formation

To study the early stage of amyloid-beta peptide (Abeta) aggregation, hexamers of the wild-type (WT) Abeta(16-35) and its mutants with amyloid-like conformations have been studied by molecular dynamics simulations in explicit water for a total time of 1.7 micros. We found that the amyloid-like structures in the WT oligomers are destabilized by the solvation of ionic D23/K28 residues, which are buried in the fibrils. This means that the desolvation of D23/K28 residues may contribute to the kinetic barrier of aggregation in the early stage. In the E22Q/D23N, D23N/K28Q, and E22Q/D23N/K28Q mutants, hydration becomes much less significant because the mutated residues have neutral amide side-chains. These amide side-chains can form linear cross-strand hydrogen bond chains, or "polar zippers", if dehydrated. These "polar zippers" increase the stability of the amyloid-like conformation, reducing the barrier for the early-stage oligomerization. This is in accord with experimental observations that both the D23/K28 lactamization and the E22Q/D23N mutation promote aggregation. We also found that the E22Q/D23N mutant prefers an amyloid-like conformation that differs from the one found for WT Abeta. This suggests that different amyloid structures may be formed under different conditions.

Evolution of protein structural classes and protein sequence families

In protein structure space, protein structures cluster into four elongated regions when mapped based solely on similarity among the 3D structures. These four regions correspond to the four major classes of present-day proteins defined by the contents of secondary structure types and their topological arrangement. Evolution of and restriction to these four classes suggest that, in most cases, the evolution of genes may have been constrained or selected to those genetic changes that results in structurally stable proteins occupying one of the four "allowed" regions of the protein structure space, "structural selection," an important component of natural selection in gene evolution. Our studies on tracing the "common structural ancestor" for each protein sequence family of known structure suggest that: (i) recently emerged proteins belong mostly to three classes; (ii) the proteins that emerged earlier evolved to gain a new class; and (iii) the proteins that emerged earliest evolved to become the present-day proteins in the four major classes, with the fourth-class proteins becoming the most dominant population. Furthermore, our studies also show that not all present-day proteins evolved from one single set of proteins in the last common ancestral organism, but new common ancestral proteins were "born" at different evolutionary times, not traceable to one or two ancestral proteins: "the multiple birth model" for the evolution of protein sequence families.

The cloning and sequencing of a cDNA encoding a WD repeat protein in cotton (Gossypium hirsutum L.)

In this research, one 1156 bp cDNA containing full open reading frame and encoding a novel 24-kDa protein with four tandem WD repeat motifs was cloned from cotton, therefore was named GhWDR and the GenBank accession number is AY870657. By search of GhWDR cDNA and amino acid sequences in the database, we found that GhWDR and OSJNBa0003G23.2 from Oryza sativa show 90% sequence identity and 84% identity to WD-repeat protein from Arabidopsis thaliana, and also has high sequence identity to other WD repeat proteins, most of which are similar to Pop3 from fission yeast (accession number T39922) and Lst8p from Saccharomyces cerevisiae (accession number NP014392). Therefore, we proposed that GhWDR could act in some cellular processes as pop3 or LST8 does. In addition, the expression of GhWDR in various tissues was studied by RT-PCR, and it is expressed in all of the studied tissues, but the level of expression is low in the leaves when compared to that of other tissues.

A geometric construction determines all permissible strand arrangements of sandwich proteins

For a large class of proteins called sandwich-like proteins (SPs), the secondary structures consist of two beta-sheets packed face-to-face, with each beta-sheet consisting typically of three to five beta-strands. An important step in the prediction of the three-dimensional structure of a SP is the prediction of its supersecondary structure, namely the prediction of the arrangement of the beta-strands in the two beta-sheets. Recently, significant progress in this direction was made, where it was shown that 91% of observed SPs form what we here call "canonical motifs." Here, we show that all canonical motifs can be constructed in a simple manner that is based on thermodynamic considerations and uses certain geometric structures. The number of these structures is much smaller than the number of possible strand arrangements. For instance, whereas for SPs consisting of six strands there exist a priori 900 possible strand arrangements, there exist only five geometric structures. Furthermore, the few motifs that are noncanonial can be constructed from canonical motifs by a simple procedure.

Valois, a component of the nuage and pole plasm, is involved in assembly of these structures, and binds to Tudor and the methyltransferase Capsuléen

Using the Capsuleen (Csul) methyltransferase as bait in the yeast two-hybrid system, we have identified a novel Drosophila protein containing multiple WD repeats and encoded by the valois (vsl) gene, which acts in pole plasm function. Vls is homologous to human MEP50, which forms a complex with the PRMT5 methyltransferase--the human homologue of Csul. We found that Vls localizes to the nuage in the nurse cells and to the pole plasm in the oocyte. Moreover vls is required for the synthesis and/or stability of Oskar and the localization of Tudor (Tud) in both the nuage and at the posterior pole of the oocyte. Furthermore, we show that Vls and a fragment of Tud interact directly in binding assay. As the PMRT5/MEP50 complex is involved in ribonucleoprotein complex assembly, we hypothesize that the Vls complex may play a similar function in assembling the nuage in nurse cells and the polar granules in the oocyte.

Apo-azurin folds via an intermediate that resembles the molten-globule

The folding of Pseudomonas aeruginosa apo-azurin was investigated with the intent of identifying putative intermediates. Two apo-mutants were constructed by replacing the main metal-binding ligand C112 with a serine (C112S) and an alanine (C112A). The guanidinium-induced unfolding free energies (DeltaG(U-N)(H2O)) of the C112S and C112A mutants were measured to 36.8 +/- 1 kJ mole(-1) and 26.1 +/- 1 kJ mole(-1), respectively, and the m-value of the transition to 23.5 +/- 0.7 kJ mole(-1) M(-1). The difference in folding free energy (DeltaDeltaG(U-N)(H2O)) is largely attributed to the intramolecular hydrogen bonding properties of the serine Ogamma in the C112S mutant, which is lacking in the C112A structure. Furthermore, only the unfolding rates differ between the two mutants, thus pointing to the energy of the native state as the source of the observed Delta DeltaG(U-N)(H2O). This also indicates that the formation of the hydrogen bonds present in C112S but absent in C112A is a late event in the folding of the apo-protein, thus suggesting that formation of the metal-binding site occurs after the rate-limiting formation of the transition state. In both mutants we also noted a burst-phase intermediate. Because this intermediate was capable of binding 1-anilinonaphtalene-8-sulfonate (ANS), as were an acid-induced species at pH 2.6, we ascribe it molten globule-like status. However, despite the presence of an intermediate, the folding of apo-azurin C112S is well approximated by a two-state kinetic mechanism.

Evolutionary optimization of a nonbiological ATP binding protein for improved folding stability

Structural comparison of in vitro evolved proteins with biological proteins will help determine the extent to which biological proteins sample the structural diversity available in protein sequence space. We have previously isolated a family of nonbiological ATP binding proteins from an unconstrained random sequence library. One of these proteins was further optimized for high-affinity binding to ATP, but biophysical characterization proved impossible due to poor solubility. To determine if such nonbiological proteins can be optimized for improved folding stability, we performed multiple rounds of mRNA-display selection under increasingly denaturing conditions. Starting from a pool of protein variants, we evolved a population of proteins capable of binding ATP in 3 M guanidine hydrochloride. One protein was chosen for further characterization. Circular dichroism, tryptophan fluorescence, and (1)H-(15)N correlation NMR studies show that this protein has a unique folded structure.

Candidates for novel RNA topologies

Because the functional repertiore of RNA molecules, like proteins, is closely linked to the diversity of their shapes, uncovering RNA's structural repertoire is vital for identifying novel RNAs, especially in genomic sequences. To help expand the limited number of known RNA families, we use graphical representation and clustering analysis of RNA secondary structures to predict novel RNA topologies and their abundance as a function of size. Representing the essential topological properties of RNA secondary structures as graphs enables enumeration, generation, and prediction of novel RNA motifs. We apply a probabilistic graph-growing method to construct the RNA structure space encompassing the topologies of existing and hypothetical RNAs and cluster all RNA topologies into two groups using topological descriptors and a standard clustering algorithm. Significantly, we find that nearly all existing RNAs fall into one group, which we refer to as "RNA-like"; we consider the other group "non-RNA-like". Our method predicts many candidates for novel RNA secondary topologies, some of which are remarkably similar to existing structures; interestingly, the centroid of the RNA-like group is the tmRNA fold, a pseudoknot having both tRNA-like and mRNA-like functions. Additionally, our approach allows estimation of the relative abundance of pseudoknot and other (e.g. tree) motifs using the "edge-cut" property of RNA graphs. This analysis suggests that pseudoknots dominate the RNA structure universe, representing more than 90% when the sequence length exceeds 120 nt; the predicted trend for <100 nt agrees with data for existing RNAs. Together with our predictions for novel "RNA-like" topologies, our analysis can help direct the design of functional RNAs and identification of novel RNA folds in genomes through an efficient topology-directed search, which grows much more slowly in complexity with RNA size compared to the traditional sequence-based search.

Conformational transitions and glycation of serum albumin in patients with minimal-change glomerulopathy

BACKGROUND

There has been a lack of study on the structural changes of serum albumin in patients with minimal change disease (MCD). To determine whether glycation and/or conformational transitions of albumin are involved in the pathogenesis of albuminuria, nine patients with MCD were enrolled in a prospective follow-up study for comparison of these parameters in serum albumin during the remission and relapse of nephrotic syndrome.

METHODS

Circular dichroism measurements were made with purified albumin. Ellipticities at each wavelength were transformed to mean residue ellipticity. Monosaccharide composition was analyzed by high-pH anion-exchange chromatography with pulsed amperometric detection.

RESULTS

There was no difference in the proportions of alpha-helix, beta-conformation, and beta-turn of albumin between the sera of control patients and those with nephrotic syndrome. However, the proportion of the random configuration was slightly higher in the plasma albumin of patients in relapse than in those in remission. The proportion of the random configuration was lower in the albumin of the serum than in the urine of patients with nephrotic syndrome, but there was no difference in the proportions of alpha-helix, beta-conformation, and beta-tum of albumin between their plasma and urine.

CONCLUSION

Our results suggest that conformational changes in albumin are involved in albuminuria in patients with MCD.

Regions of minimal structural variation among members of protein domain superfamilies: application to remote homology detection and modelling using distant relationships

Structurally conserved regions or structural templates have been identified and examined for features such as amino acid content, solvent accessibility, secondary structures, non-polar interaction, residue packing and extent of structural deviations in 179 aligned members of superfamilies involving 1208 pairs of protein domains. An analysis of these structural features shows that the retention of secondary structural conservation and similar hydrogen bonding pattern within the templates is 2.5 and 1.8 times higher, respectively, than full-length alignments suggesting that they form the minimum structural requirement of a superfamily. The identification and availability of structural templates find value in different areas of protein structure prediction and modelling such as in sensitive sequence searches, accurate sequence alignment and three-dimensional modelling on the basis of distant relationships.

The amino terminal domain of a novel WD repeat protein from Trypanosoma cruzi contains a non-canonical mitochondrial targeting signal

WD (tryptophan/aspartic acid) repeat proteins perform a wide variety of functions in eukaryotic cells. They are characterised by the presence of a number of conserved repeat motifs that contribute to the beta-propeller structures which are the common feature of this large group of proteins. We report here the properties of the first characterised member of this family in the American trypanosome, Trypanosoma cruzi (TcBPP1). In the CL Brener clone the protein is 482 amino acids long and is predicted to contain four WD repeat motifs, flanked by amino and carboxyl terminal extensions. TcBPP1 is a single copy gene present on a 1.0/1.6 Mb pair of homologous chromosomes in a locus that is syntenic with the corresponding regions of Trypanosoma brucei and Leishmania major chromosomes. Consistent with the proposed hybrid nature of the CL Brener clone, the proteins encoded by the two different alleles share only 97% identity at the amino acid level. To determine subcellular location, we examined transfected parasites for the distribution of green fluorescent protein (GFP) fused with different regions of TcBPP1. These studies demonstrated that a 115 amino acid peptide derived from the amino terminal domain of TcBPP1 is able to target GFP to the mitochondrion. Interestingly this region lacks a typical amino terminal presequence suggesting that mitochondrial import is mediated by an alternative targeting signal.

VARICOSE, a WD-domain protein, is required for leaf blade development

To gain insight into the processes controlling leaf development, we characterized an Arabidopsis mutant, varicose (vcs), with leaf and shoot apical meristem defects. The vcs phenotype is temperature dependent; low temperature growth largely suppressed defects, whereas high growth temperatures resulted in severe leaf and meristem defects. VCS encodes a putative WD-domain containing protein, suggesting a function involving protein-protein interactions. Temperature shift experiments indicated that VCS is required throughout leaf development, but normal secondary vein patterning required low temperature early in leaf development. The low-temperature vcs phenotype is enhanced in axr1-3 vcs double mutants and in vcs mutants grown in the presence of polar auxin transport inhibitors, however, vcs has apparently normal auxin responses. Taken together, these observations suggest a role for VCS in leaf blade formation.

The WD-repeat protein superfamily in Arabidopsis: conservation and divergence in structure and function

BACKGROUND

The WD motif (also known as the Trp-Asp or WD40 motif) is found in a multitude of eukaryotic proteins involved in a variety of cellular processes. Where studied, repeated WD motifs act as a site for protein-protein interaction, and proteins containing WD repeats (WDRs) are known to serve as platforms for the assembly of protein complexes or mediators of transient interplay among other proteins. In the model plant Arabidopsis thaliana, members of this superfamily are increasingly being recognized as key regulators of plant-specific developmental events.

RESULTS

We analyzed the predicted complement of WDR proteins from Arabidopsis, and compared this to those from budding yeast, fruit fly and human to illustrate both conservation and divergence in structure and function. This analysis identified 237 potential Arabidopsis proteins containing four or more recognizable copies of the motif. These were classified into 143 distinct families, 49 of which contained more than one Arabidopsis member. Approximately 113 of these families or individual proteins showed clear homology with WDR proteins from the other eukaryotes analyzed. Where conservation was found, it often extended across all of these organisms, suggesting that many of these proteins are linked to basic cellular mechanisms. The functional characterization of conserved WDR proteins in Arabidopsis reveals that these proteins help adapt basic mechanisms for plant-specific processes.

CONCLUSIONS

Our results show that most Arabidopsis WDR proteins are strongly conserved across eukaryotes, including those that have been found to play key roles in plant-specific processes, with diversity in function conferred at least in part by divergence in upstream signaling pathways, downstream regulatory targets and /or structure outside of the WDR regions.