Structural and Functional Analysis of Bacterial Sulfonosphingolipids and Rosette-Inducing Factor 2 (RIF-2) by Mass Spectrometry-Guided Isolation and Total Synthesis.

We have analyzed the abundance of bacterial sulfonosphingolipids, including rosette-inducing factors (RIFs), in seven bacterial prey strains by using high-resolution tandem mass spectrometry (HRMS2 ) and molecular networking (MN) within the Global Natural Product Social Molecular Networking (GNPS) web platform. Six sulfonosphingolipids resembling RIFs were isolated and their structures were elucidated based on comparative MS and NMR studies. Here, we also report the first total synthesis of two RIF-2 diastereomers and one congener in 15 and eight synthetic steps, respectively. For the total synthesis of RIF-2 congeners, we employed a decarboxylative cross-coupling reaction to synthesize the necessary branched α-hydroxy fatty acids, and the Garner-aldehyde approach to generate the capnine base carrying three stereogenic centers. Bioactivity studies in the choanoflagellate Salpingoeca rosetta revealed that the rosette inducing activity of RIFs is inhibited dose dependently by the co-occurring sulfonosphingolipid sulfobacins D and F and that activity of RIFs is specific for isolates obtained from Algoriphagus.

Domain Analysis and Motif Matcher (DAMM): A Program to Predict Selectivity Determinants in Monosiga brevicollis PDZ Domains Using Human PDZ Data.

Choanoflagellates are single-celled eukaryotes with complex signaling pathways. They are considered the closest non-metazoan ancestors to mammals and other metazoans and form multicellular-like states called rosettes. The choanoflagellate Monosiga brevicollis contains over 150 PDZ domains, an important peptide-binding domain in all three domains of life (Archaea, Bacteria, and Eukarya). Therefore, an understanding of PDZ domain signaling pathways in choanoflagellates may provide insight into the origins of multicellularity. PDZ domains recognize the C-terminus of target proteins and regulate signaling and trafficking pathways, as well as cellular adhesion. Here, we developed a computational software suite, Domain Analysis and Motif Matcher (DAMM), that analyzes peptide-binding cleft sequence identity as compared with human PDZ domains and that can be used in combination with literature searches of known human PDZ-interacting sequences to predict target specificity in choanoflagellate PDZ domains. We used this program, protein biochemistry, fluorescence polarization, and structural analyses to characterize the specificity of A9UPE9_MONBE, a M. brevicollis PDZ domain-containing protein with no homology to any metazoan protein, finding that its PDZ domain is most similar to those of the DLG family. We then identified two endogenous sequences that bind A9UPE9 PDZ with <100 µM affinity, a value commonly considered the threshold for cellular PDZ-peptide interactions. Taken together, this approach can be used to predict cellular targets of previously uncharacterized PDZ domains in choanoflagellates and other organisms. Our data contribute to investigations into choanoflagellate signaling and how it informs metazoan evolution.

Structural characterization and computational analysis of PDZ domains in Monosiga brevicollis.

Identification of the molecular networks that facilitated the evolution of multicellular animals from their unicellular ancestors is a fundamental problem in evolutionary cellular biology. Choanoflagellates are recognized as the closest extant nonmetazoan ancestors to animals. These unicellular eukaryotes can adopt a multicellular-like "rosette" state. Therefore, they are compelling models for the study of early multicellularity. Comparative studies revealed that a number of putative human orthologs are present in choanoflagellate genomes, suggesting that a subset of these genes were necessary for the emergence of multicellularity. However, previous work is largely based on sequence alignments alone, which does not confirm structural nor functional similarity. Here, we focus on the PDZ domain, a peptide-binding domain which plays critical roles in myriad cellular signaling networks and which underwent a gene family expansion in metazoan lineages. Using a customized sequence similarity search algorithm, we identified 178 PDZ domains in the Monosiga brevicollis proteome. This includes 11 previously unidentified sequences, which we analyzed using Rosetta and homology modeling. To assess conservation of protein structure, we solved high-resolution crystal structures of representative M. brevicollis PDZ domains that are homologous to human Dlg1 PDZ2, Dlg1 PDZ3, GIPC, and SHANK1 PDZ domains. To assess functional conservation, we calculated binding affinities for mbGIPC, mbSHANK1, mbSNX27, and mbDLG-3 PDZ domains from M. brevicollis. Overall, we find that peptide selectivity is generally conserved between these two disparate organisms, with one possible exception, mbDLG-3. Overall, our results provide novel insight into signaling pathways in a choanoflagellate model of primitive multicellularity.

MotifAnalyzer-PDZ: A computational program to investigate the evolution of PDZ-binding target specificity.

Recognition of short linear motifs (SLiMs) or peptides by proteins is an important component of many cellular processes. However, due to limited and degenerate binding motifs, prediction of cellular targets is challenging. In addition, many of these interactions are transient and of relatively low affinity. Here, we focus on one of the largest families of SLiM-binding domains in the human proteome, the PDZ domain. These domains bind the extreme C-terminus of target proteins, and are involved in many signaling and trafficking pathways. To predict endogenous targets of PDZ domains, we developed MotifAnalyzer-PDZ, a program that filters and compares all motif-satisfying sequences in any publicly available proteome. This approach enables us to determine possible PDZ binding targets in humans and other organisms. Using this program, we predicted and biochemically tested novel human PDZ targets by looking for strong sequence conservation in evolution. We also identified three C-terminal sequences in choanoflagellates that bind a choanoflagellate PDZ domain, the Monsiga brevicollis SHANK1 PDZ domain (mbSHANK1), with endogenously-relevant affinities, despite a lack of conservation with the targets of a homologous human PDZ domain, SHANK1. All three are predicted to be signaling proteins, with strong sequence homology to cytosolic and receptor tyrosine kinases. Finally, we analyzed and compared the positional amino acid enrichments in PDZ motif-satisfying sequences from over a dozen organisms. Overall, MotifAnalyzer-PDZ is a versatile program to investigate potential PDZ interactions. This proof-of-concept work is poised to enable similar types of analyses for other SLiM-binding domains (e.g., MotifAnalyzer-Kinase). MotifAnalyzer-PDZ is available at http://motifAnalyzerPDZ.cs.wwu.edu.

Evolutionary conservation of complexins: from choanoflagellates to mice.

Complexins are synaptic SNARE complex-binding proteins that cooperate with synaptotagmins in activating Ca(2+)-stimulated, synaptotagmin-dependent synaptic vesicle exocytosis and in clamping spontaneous, synaptotagmin-independent synaptic vesicle exocytosis. Here, we show that complexin sequences are conserved in some non-metazoan unicellular organisms and in all metazoans, suggesting that complexins are a universal feature of metazoans that predate metazoan evolution. We show that complexin from Nematostella vectensis, a cnidarian sea anemone far separated from mammals in metazoan evolution, functionally replaces mouse complexins in activating Ca(2+)-triggered exocytosis, but is unable to clamp spontaneous exocytosis. Thus, the activating function of complexins is likely conserved throughout metazoan evolution.

Novel leucine rich repeat domains in proteins from unicellular eukaryotes and bacteria.

Leucine rich repeats (LRRs) are present in over 20,000 proteins from viruses to eukaryotes. Two to sixty-two LRRs occur in tandem. Each repeat is typically 20-30 residues long and can be divided into an HCS (Highly conserved segment) and a VS (Variable segment). The HCS part consists of an eleven or a twelve residue stretch, LxxLxLxxNx(x/-)L, in which "L" is Leu, Ile, Val, or Phe, "N" is Asn, Thr, Ser, or Cys, "x" is a non-conserved residue, and "-" is a possible deletion site. Eight classes have been recognized. However, there are many unclassified or unrecognized LRRs. Here we performed to search novel LRRs using protein sequence database. The novel LRR domains are present over three hundred proteins, which include fungal ECM33 protein and Monosiga brevicollis LRR receptor kinase, from unicellular eukaryotes and bacteria. The HCS part is clearly different from that of the known LRRs and consists of a twelve or a thirteen residue stretch, VxGx(L/F)x(L/C)xxNx(x/-)L, that is characterized by the addition of Gly between the first conserved Val and the second conserved Leu. The novel LRRs identified here form a new family. The novel LRR domains were classified into four classes. The VS parts of the two classes are consistent with those of known, normal "SDS22-like" and "IRREKO" classes, while the other two classes have unique VS parts. The structures, functions, and evolution of the novel LRR domains and their proteins are described. The present results should stimulate various experimental studies.

In silico resurrection of the major vault protein suggests it is ancestral in modern eukaryotes.

Vaults are very large oligomeric ribonucleoproteins conserved among a variety of species. The rat vault 3D structure shows an ovoid oligomeric particle, consisting of 78 major vault protein monomers, each of approximately 861 amino acids. Vaults are probably the largest ribonucleoprotein structures in eukaryote cells, being approximately 70 nm in length with a diameter of 40 nm--the size of three ribosomes and with a lumen capacity of 50 million Å(3). We use both protein sequences and inferred ancestral sequences for in silico virtual resurrection of tertiary and quaternary structures to search for vaults in a wide variety of eukaryotes. We find that the vault's phylogenetic distribution is widespread in eukaryotes, but is apparently absent in some notable model organisms. Our conclusion from the distribution of vaults is that they were present in the last eukaryote common ancestor but they have apparently been lost from a number of groups including fungi, insects, and probably plants. Our approach of inferring ancestral 3D and quaternary structures is expected to be useful generally.

Cooperatively generated stresslet flows supply fresh fluid to multicellular choanoflagellate colonies.

The flagellated protozoan Salpingoeca rosetta is one of the closest relatives of multicellular animals. Unicellular S. rosetta can be induced to form multicellular colonies, but colonies swim more slowly than individual cells so the advantages conferred by colony formation are uncertain. Here we use theoretical models to show that hydrodynamic cooperation between cells can increase the fluid supply to the colony, an important predictor of feeding rate. Our results suggest that hydrodynamic benefits may have been an important selective factor in the evolution of early multicellular animals.

Ancient complexity, opisthokont plasticity, and discovery of the 11th subfamily of Arf GAP proteins.

The organelle paralogy hypothesis is one model for the acquisition of nonendosymbiotic organelles, generated from molecular evolutionary analyses of proteins encoding specificity in the membrane traffic system. GTPase activating proteins (GAPs) for the ADP-ribosylation factor (Arfs) GTPases are additional regulators of the kinetics and fidelity of membrane traffic. Here we describe molecular evolutionary analyses of the Arf GAP protein family. Of the 10 subfamilies previously defined in humans, we find that 5 were likely present in the last eukaryotic common ancestor. Of the 3 most recently derived subfamilies, 1 was likely present in the ancestor of opisthokonts (animals and fungi) and apusomonads (flagellates classified as the sister lineage to opisthokonts), while 2 arose in the holozoan lineage. We also propose to have identified a novel ancient subfamily (ArfGAPC2), present in diverse eukaryotes but which is lost frequently, including in the opisthokonts. Surprisingly few ancient domains accompanying the ArfGAP domain were identified, in marked contrast to the extensively decorated human Arf GAPs. Phylogenetic analyses of the subfamilies reveal patterns of single and multiple gene duplications specific to the Holozoa, to some degree mirroring evolution of Arf GAP targets, the Arfs. Conservation, and lack thereof, of various residues in the ArfGAP structure provide contextualization of previously identified functional amino acids and their application to Arf GAP biology in general. Overall, our results yield insights into current Arf GAP biology, reveal complexity in the ancient eukaryotic ancestor and integrate the Arf GAP family into a proposed mechanism for the evolution of nonendosymbiotic organelles.

Comment on "Positive selection of tyrosine loss in metazoan evolution".

Tan et al. (Reports, 25 September 2009, p. 1686) argued that loss of tyrosine residues from proteins in metazoans was driven by positive selection to remove potentially deleterious phosphorylation sites. We challenge this hypothesis, providing evidence that the high guanine-cytosine (GC) content of metazoan genomes was the primary driver in the loss of tyrosine residues.

Biosilicification of loricate choanoflagellate: organic composition of the nanotubular siliceous costal strips of Stephanoeca diplocostata.

Loricate choanoflagellates (unicellular, eukaryotic flagellates; phylum Choanozoa) synthesize a basket-like siliceous lorica reinforced by costal strips (diameter of approximately 100 nm and length of 3 µm). In the present study, the composition of these siliceous costal strips is described, using Stephanoeca diplocostata as a model. Analyses by energy-dispersive X-ray spectroscopy (EDX), coupled with transmission electron microscopy (TEM), indicate that the costal strips comprise inorganic and organic components. The organic, proteinaceous scaffold contained one major polypeptide of mass 14 kDa that reacted with wheat germ agglutinin. Polyclonal antibodies were raised that allowed mapping of the proteinaceous scaffold, the (glyco)proteins, within the costal strips. Subsequent in vitro studies revealed that the organic scaffold of the costal strips stimulates polycondensation of ortho-silicic acid in a concentration- and pH-dependent way. Taken together, the data gathered indicate that the siliceous costal strips are formed around a proteinaceous scaffold that supports and maintains biosilicification. A scheme is given that outlines that the organic template guides both the axial and the lateral growth of the strips.