Amino Acid Preferences Research Articles

PO-406 Investigation of the repertoire of peptides bound to MHC class I molecules in tasmanian devil transmissible cancers for the development of a peptide vaccine

IntroductionThe marsupial species Tasmanian devil harbours two contagious cancers, Devil Facial Tumour 1 and 2 (DFT1 and DFT2); both are passed between individuals as an allogeneic graft, initiating tumours around the neck and face. DFT1 emerged 20 years ago, causes 100% mortality and has decimated the devil population. In contrast, DFT2 emerged only recently and is still being characterised. As allografts, both tumours should be recognised by T cells interacting with non-self Major Histocompatibility Complex (MHC) molecules and associated peptides. We aim to isolate tumour-specific MHC-bound peptides in order to determine MHC/peptide complexes breaking or inducing tolerance and design a vaccine.Material and methodsWe have characterised the immunopeptidomes of DFT1-IFNγ, DFT2 and devil fibroblast cell lines. The DFT1 cell line was stimulated with devil interferon-gamma to up-regulate MHC class I expression. MHC class I molecule-peptide complexes were isolated by immunoaffinity purification using an antibody against Tasmanian devil β2-microblobulin. Eluted peptide/MHC complexes were fractionated by reversed phase HPLC and analysed by Liquid Chromatography tandem mass spectrometry (LC-MS/MS). All experiments were performed in triplicate. Peptides were identified by searching spectra against custom Tasmanian devil databases using the PEAKS software and data visualised in R studio and GraphPad Prism.Results and discussionsBetween 6373 and 2243 peptides were identified for each cell line, but only peptides found in all replicates for each cell line analysed further. 1806 peptides were identified in all three replicates of the fibroblast cell line, while 455 and 379 peptides were common to all three replicates for the DFT1-IFNγ and DFT2 cell lines respectively. Analysis of peptide length shows a preference for 8mers and 9mers in devil fibroblasts and DFT2 cells and a preference for 9mers in DFT1-IFNγ. We next searched for binding motifs for the 8mers and 9mers across all cell lines and found potential anchor residues at position 3 and position 8/9, where there was a preference for hydrophobic amino acids (in particular Leucine). We then identified 61 and 55 peptides unique to DFT1-IFNγ and DFT2 respectively.ConclusionThis is the first study to characterise the repertoire of peptides bound to MHC molecules in contagious cancers and represents a pivotal step for understanding the immunological features of transmissible cancers.

MipLAAO, a new L-amino acid oxidase from the redtail coral snake Micrurus mipartitus.

L-amino acid oxidases (LAAOs) are ubiquitous enzymes in nature. Bioactivities described for these enzymes include apoptosis induction, edema formation, induction or inhibition of platelet aggregation, as well as antiviral, antiparasite, and antibacterial actions. With over 80 species, Micrurus snakes are the representatives of the Elapidae family in the New World. Although LAAOs in Micrurus venoms have been predicted by venom gland transcriptomic studies and detected in proteomic studies, no enzymes of this kind have been previously purified from their venoms. Earlier proteomic studies revealed that the venom of M. mipartitus from Colombia contains ∼4% of LAAO. This enzyme, here named MipLAAO, was isolated and biochemically and functionally characterized. The enzyme is found in monomeric form, with an isotope-averaged molecular mass of 59,100.6 Da, as determined by MALDI-TOF. Its oxidase activity shows substrate preference for hydrophobic amino acids, being optimal at pH 8.0. By nucleotide sequencing of venom gland cDNA of mRNA transcripts obtained from a single snake, six isoforms of MipLAAO with minor variations among them were retrieved. The deduced sequences present a mature chain of 483 amino acids, with a predicted pI of 8.9, and theoretical masses between 55,010.9 and 55,121.0 Da. The difference with experimentally observed mass is likely due to glycosylation, in agreement with the finding of three putative N-glycosylation sites in its amino acid sequence. A phylogenetic analysis of MmipLAAO placed this new enzyme within the clade of homologous proteins from elapid snakes, characterized by the conserved Serine at position 223, in contrast to LAAOs from viperids. MmipLAAO showed a potent bactericidal effect on S. aureus (MIC: 2 µg/mL), but not on E. coli. The former activity could be of interest to future studies assessing its potential as antimicrobial agent.

Identification and Analysis of Key Residues Involved in Folding and Binding of Protein-carbohydrate Complexes.

Protein-carbohydrate interactions play vital roles in several biological processes in living organisms. The comparative analysis of binding site residues along with stabilizing residues in protein-carbohydrate complexes provides ample insights to understand the structure, function and recognition mechanism. The main objective of this study is to identify and analyze the residues, which are involved in both folding and binding of the protein-carbohydrate complexes. We have identified the stabilizing residues using the knowledge of hydrophobicity, longrange interactions and conservation, as well as binding site residues using a distance cutoff of 3.5Å between any heavy atoms in protein and ligand. Residues, which are common in stabilizing and binding, are termed as key residues. These key resides are analyzed with various sequence and structure based parameters such as frequency of occurrence, surrounding hydrophobicity, longrange order and conservation score. In this work, we have identified 2.45% binding site residues in a non-redundant dataset of 1130 complexes using distance-based criteria and 7.07% stabilizing residues using the concepts of hydrophobicity, long-range interactions and conservation of residues. Further, 5.9% of binding and 2.04% of stabilizing residues are common to each other, which are termed as key residues. The key residues have been analysed based on protein classes, carbohydrate types, gene ontology functional classifications, amino acid preference and structure-based parameters. We found that all-β, α+β and α/β have more key residues than other protein classes and most of the KRs are present in β-strands, which shows their importance in stability and binding of complexes. On the ligand side, Lsaccharide has the highest number of key residues and it has a high percentage of KRs in SRs and BRs than other carbohydrate types. Further, polar and charged residues have a high tendency to serve as key residues. Classifications based on gene ontology terms revealed that Lys is preferred in all the three groups: molecular functions, biological processes and cellular components. Key residues have 6 to 9 contacts within the protein and make only one contact with the carbohydrate ligand. These contacts are dominant to form polar-nonpolar contacts followed by the contacts between charged atoms. Further, the influence of sequence and structural parameters such as surrounding hydrophobicity, solvent accessibility, secondary structure, long-range order and conservation score has been discussed. The results obtained in the present work provide deep insights for understanding the interplay between stability and binding in protein-carbohydrate complexes.

Molecular dynamics-derived rotamer libraries for d-amino acids within homochiral and heterochiral polypeptides.

Computational resources have contributed to the design and engineering of novel proteins by integrating genomic, structural and dynamic aspects of proteins. Non-canonical amino acids, such as d-amino acids, expand the available sequence space for designing and engineering proteins; however, the rotamer libraries for d-amino acids are usually constructed as the mirror images of l-amino acid rotamer libraries, an assumption that has not been tested. To this end, we have performed molecular dynamics (MD) simulations of model host-guest peptide systems containing d-amino acids. Our simulations systematically address the applicability of the mirror image convention as well as the effects of neighboring residue chirality. Rotamer libraries derived from these systems provide realistic rotamer distributions suitable for use in both rational and computational design workflows. Our simulations also address the impact of chirality on the intrinsic conformational preferences of amino acids, providing fundamental insights into the relationship between chirality and biomolecular dynamics. While d-amino acids are rare in naturally occurring proteins, they are used in designed proteins to stabilize a desired conformation, increase bioavailability or confer favorable biochemical and physical attributes. Here, we present d-amino acid rotamer libraries derived from MD simulations of alanine-based host-guest pentapeptides and show how certain residues can deviate from mirror image symmetry. Our simulations directly model d-amino acids as guest residues within the chiral l-Ala and d-Ala pentapeptide series to explicitly incorporate any contributions resulting from the chiralities of neighboring residues.

GC bias lead to increased small amino acids and random coils of proteins in cold-water fishes

BackgroundTemperature adaptation of biological molecules is fundamental in evolutionary studies but remains unsolved. Fishes living in cold water are adapted to low temperatures through adaptive modification of their biological molecules, which enables their functioning in extreme cold. To study nucleotide and amino acid preference in cold-water fishes, we investigated the substitution asymmetry of codons and amino acids in protein-coding DNA sequences between cold-water fishes and tropical fishes., The former includes two Antarctic fishes, Dissostichus mawsoni (Antarctic toothfish), Gymnodraco acuticeps (Antarctic dragonfish), and two temperate fishes, Gadus morhua (Atlantic cod) and Gasterosteus aculeatus (stickleback), and the latter includes three tropical fishes, including Danio rerio (zebrafish), Oreochromis niloticus (Nile tilapia) and Xiphophorus maculatus (Platyfish).ResultsCold-water fishes showed preference for Guanines and cytosines (GCs) in both synonymous and nonsynonymous codon substitution when compared with tropical fishes. Amino acids coded by GC-rich codons are favored in the temperate fishes, while those coded by AT-rich codons are disfavored. Similar trends were discovered in Antarctic fishes but were statistically weaker. The preference of GC rich codons in nonsynonymous substitution tends to increase ratio of small amino acid in proteins, which was demonstrated by biased small amino acid substitutions in the cold-water species when compared with the tropical species, especially in the temperate species. Prediction and comparison of secondary structure of the proteomes showed that frequency of random coils are significantly larger in the cold-water fish proteomes than those of the tropical fishes.ConclusionsOur results suggested that natural selection in cold temperature might favor biased GC content in the coding DNA sequences, which lead to increased frequency of small amino acids and consequently increased random coils in the proteomes of cold-water fishes.

DNA Methylation of T1R1 Gene in the Vegetarian Adaptation of Grass Carp Ctenopharyngodon idella

Although previous studies have indicated importance of taste receptors in food habits formation in mammals, little is known about those in fish. Grass carp is an excellent model for studying vegetarian adaptation, as it shows food habit transition from carnivore to herbivore. In the present study, pseudogenization or frameshift mutations of the umami receptors that hypothesized related to dietary switch in vertebrates, were not found in grass carp, suggesting other mechanisms for vegetarian adaptation in grass carp. T1R1 and T1R3 strongly responded to L-Arg and L-Lys, differing from those of zebrafish and medaka, contributing to high species specificity in amino acid preferences and diet selection of grass carp. After food habit transition of grass carp, DNA methylation levels were higher in CPG1 and CPG3 islands of upstream control region of T1R1 gene. Luciferase activity assay of upstream regulatory region of T1R1 (−2500-0 bp) without CPG1 or CPG3 indicated that CPG1 and CPG3 might be involved in transcriptional regulation of T1R1 gene. Subsequently, high DNA methylation decreased expression of T1R1 in intestinal tract. It could be a new mechanism to explain, at least partially, the vegetarian adaptation of grass carp by regulation of expression of umami receptor via epigenetic modification.

An Unbiased Determination of pMHC Repertoires for Improved Antigen Prediction

Abstract Successful understanding and treatment of immune diseases spanning autoimmunity, cancer, and infectious disease relies upon understanding what is seen by the immune system. For T cell responses, this is controlled by MHCs, specialized polymorphic proteins that present a snapshot of the proteome of their parent cell or surroundings. Importantly, MHC polymorphisms cluster at their peptide-binding surface, where even single substitutions can greatly alter the peptide-MHC (pMHC) repertoire seen by T cells. Due to extreme person-to-person variability in MHC usage, data for peptide MHC specificity is largely limited to a handful of the most prevalent alleles. While algorithms that predict pMHC interactions perform well for these alleles, they can suffer from inaccuracy for less-studied alleles, due in part to biases and limited diversity in their training data. We present a platform that provides large unbiased data sets to improve pMHC repertoire definitions and predictions. Using yeast-displayed MHC molecules displaying a random peptide library, we identified hundreds of thousands of unique peptide binders for both well- and poorly-studied human MHC alleles. Our data demonstrates strong amino acid preferences at defined “anchor” positions for each allele, but also finds preferences at auxiliary sites that are not considered by most prediction software. We also find noncanonical amino acid preferences at defined anchor positions that confound current software, suggesting alternative peptide conformations. These data show the importance of unbiased pMHC repertoires to improve existing antigen prediction software, and suggest our approach can be used to define presently under-studied MHC alleles.

Byssus Structure and Protein Composition in the Highly Invasive Fouling Mussel Limnoperna fortunei.

Biofouling mediated by byssus adhesion in invasive bivalves has become a global environmental problem in aquatic ecosystems, resulting in negative ecological and economic consequences. Previous studies suggested that mechanisms responsible for byssus adhesion largely vary among bivalves, but it is poorly understood in freshwater species. Understanding of byssus structure and protein composition is the prerequisite for revealing these mechanisms. Here, we used multiple methods, including scanning electron microscope, liquid chromatography–tandem mass spectrometry, transcriptome sequencing, real-time quantitative PCR, inductively coupled plasma mass spectrometry, to investigate structure, and protein composition of byssus in the highly invasive freshwater mussel Limnoperna fortunei. The results indicated that the structure characteristics of adhesive plaque, proximal and distal threads were conducive to byssus adhesion, contributing to the high biofouling capacity of this species. The 3,4-dihydroxyphenyl-α-alanine (Dopa) is a major post-transnationally modification in L. fortunei byssus. We identified 16 representative foot proteins with typical repetitive motifs and conserved domains by integrating transcriptomic and proteomic approaches. In these proteins, Lfbp-1, Lffp-2, and Lfbp-3 were specially located in foot tissue and highly expressed in the rapid byssus formation period, suggesting the involvement of these foot proteins in byssus production and adhesion. Multiple metal irons, including Ca2+, Mg2+, Zn2+, Al3+, and Fe3+, were abundant in both foot tissue and byssal thread. The heavy metals in these irons may be directly accumulated by L. fortunei from surrounding environments. Nevertheless, some metal ions (e.g., Ca2+) corresponded well with amino acid preferences of L. fortunei foot proteins, suggesting functional roles of these metal ions by interacting with foot proteins in byssus adhesion. Overall, this study provides structural and molecular bases of adhesive mechanisms of byssus in L. fortunei, and findings here are expected to develop strategies against biofouling by freshwater organisms.

Extended cleavage specificity of human neutrophil cathepsin G: A low activity protease with dual chymase and tryptase-type specificities.

Human neutrophils express at least four active serine proteases, cathepsin G, N-elastase, proteinase 3 and neutrophil serine protease 4 (NSP4). They have all been extensively studied due to their importance in neutrophil biology and immunity. However, their extended cleavage specificities have never been determined in detail. Here we present a detailed cleavage specificity analysis of human cathepsin G (hCG). The specificity was determined by phage display analysis and the importance of individual amino acids in and around the cleavage site was then validated using novel recombinant substrates. To provide a broader context to this serine protease, a comparison was made to the related mast cell protease, human chymase (HC). hCG showed similar characteristics to HC including both the primary and extended specificities. As expected, Phe, Tyr, Trp and Leu were preferred in the P1 position. In addition, both proteases showed a preference for negatively charged amino acids in the P2´ position of substrates and a preference for aliphatic amino acids both upstream and downstream of the cleavage site. However, overall the catalytic activity of hCG was ~10-fold lower than HC. hCG has previously been reported to have a dual specificity consisting of chymase and tryptase-type activities. In our analysis, tryptase activity against substrates with Lys in P1 cleavage position was indeed only 2-fold less efficient as compared to optimal chymase substrates supporting strong dual-type specificity. We hope the information presented here on extended cleavage specificities of hCG and HC will assist in the search for novel in vivo substrates for these proteases as well as aid in the efforts to better understand the role of hCG in immunity and bacterial defence.

Mapping mutational effects along the evolutionary landscape of HIV envelope.

The immediate evolutionary space accessible to HIV is largely determined by how single amino acid mutations affect fitness. These mutational effects can shift as the virus evolves. However, the prevalence of such shifts in mutational effects remains unclear. Here, we quantify the effects on viral growth of all amino acid mutations to two HIV envelope (Env) proteins that differ at [Formula: see text]100 residues. Most mutations similarly affect both Envs, but the amino acid preferences of a minority of sites have clearly shifted. These shifted sites usually prefer a specific amino acid in one Env, but tolerate many amino acids in the other. Surprisingly, shifts are only slightly enriched at sites that have substituted between the Envs-and many occur at residues that do not even contact substitutions. Therefore, long-range epistasis can unpredictably shift Env's mutational tolerance during HIV evolution, although the amino acid preferences of most sites are conserved between moderately diverged viral strains.

Accelerated Estimation of Frequency Classes in Site-Heterogeneous Profile Mixture Models

As a consequence of structural and functional constraints, proteins tend to have site-specific preferences for particular amino acids. Failing to adjust for heterogeneity of frequencies over sites can lead to artifacts in phylogenetic estimation. Site-heterogeneous mixture-models have been developed to address this problem. However, due to prohibitive computational times, maximum likelihood implementations utilize fixed component frequency vectors inferred from sequences in a database that are external to the alignment under analysis. Here, we propose a composite likelihood approach to estimation of component frequencies for a mixture model that directly uses the data from the alignment of interest. In the common case that the number of taxa under study is not large, several adjustments to the default composite likelihood are shown to be necessary. In simulations, the approach is shown to provide large improvements over hierarchical clustering. For empirical data, substantial improvements in likelihoods are found over mixtures using fixed components.

A conserved leucine zipper-like motif accounts for strong tetramerization capabilities of SEPALLATA-like MADS-domain transcription factors.

The development of angiosperm flowers is regulated by homeotic MIKC-type MADS-domain transcription factors that activate or repress target genes via the formation of DNA-bound, organ-specific tetrameric complexes. The protein-protein interaction (PPI) capabilities differ considerably between different MIKC-type proteins. In Arabidopsis thaliana the floral homeotic protein SEPALLATA3 (SEP3) acts as a hub that incorporates numerous other MADS-domain proteins into tetrameric complexes that would otherwise not form. However, the molecular mechanisms that underlie these promiscuous interactions remain largely unknown. In this study, we created a collection of amino acid substitution mutants of SEP3 to quantify the contribution of individual residues on protein tetramerization during DNA-binding, employing methods of molecular biophysics. We show that leucine residues at certain key positions form a leucine-zipper structure that is essential for tetramerization of SEP3, whereas the introduction of physicochemically very similar residues at respective sites impedes the formation of DNA-bound tetramers. Comprehensive molecular evolutionary analyses of MADS-domain proteins from a diverse set of flowering plants revealed exceedingly high conservation of the identified leucine residues within SEP3-subfamily proteins throughout angiosperm evolution. In contrast, MADS-domain proteins that are unable to tetramerize among themselves exhibit preferences for other amino acids at homologous sites. Our findings indicate that the subfamily-specific conservation of amino acid residues at just a few key positions accounts for subfamily-specific interaction capabilities of MADS-domain transcription factors and this has shaped the present-day structure of the PPI network controlling flower development.

High throughput protease profiling comprehensively defines active site specificity for thrombin and ADAMTS13

We have combined random 6 amino acid substrate phage display with high throughput sequencing to comprehensively define the active site specificity of the serine protease thrombin and the metalloprotease ADAMTS13. The substrate motif for thrombin was determined by >6,700 cleaved peptides, and was highly concordant with previous studies. In contrast, ADAMTS13 cleaved only 96 peptides (out of >107 sequences), with no apparent consensus motif. However, when the hexapeptide library was substituted into the P3-P3′ interval of VWF73, an exosite-engaging substrate of ADAMTS13, 1670 unique peptides were cleaved. ADAMTS13 exhibited a general preference for aliphatic amino acids throughout the P3-P3′ interval, except at P2 where Arg was tolerated. The cleaved peptides assembled into a motif dominated by P3 Leu, and bulky aliphatic residues at P1 and P1′. Overall, the P3-P2′ amino acid sequence of von Willebrand Factor appears optimally evolved for ADAMTS13 recognition. These data confirm the critical role of exosite engagement for substrates to gain access to the active site of ADAMTS13, and define the substrate recognition motif for ADAMTS13. Combining substrate phage display with high throughput sequencing is a powerful approach for comprehensively defining the active site specificity of proteases.

Biochemical Properties of Human D-Amino Acid Oxidase.

D-amino acid oxidase catalyzes the oxidative deamination of D-amino acids. In the brain, the NMDA receptor coagonist D-serine has been proposed as its physiological substrate. In order to shed light on the mechanisms regulating D-serine concentration at the cellular level, we biochemically characterized human DAAO (hDAAO) in greater depth. In addition to clarify the physical-chemical properties of the enzyme, we demonstrated that divalent ions and nucleotides do not affect flavoenzyme function. Moreover, the definition of hDAAO substrate specificity demonstrated that D-cysteine is the best substrate, which made it possible to propose it as a putative physiological substrate in selected tissues. Indeed, the flavoenzyme shows a preference for hydrophobic amino acids, some of which are molecules relevant in neurotransmission, i.e., D-kynurenine, D-DOPA, and D-tryptophan. hDAAO shows a very low affinity for the flavin cofactor. The apoprotein form exists in solution in equilibrium between two alternative conformations: the one at higher affinity for FAD is favored in the presence of an active site ligand. This may represent a mechanism to finely modulate hDAAO activity by substrate/inhibitor presence. Taken together, the peculiar properties of hDAAO seem to have evolved in order to use this flavoenzyme in different tissues to meet different physiological needs related to D-amino acids.

Local and global anatomy of antibody-protein antigen recognition.

Deciphering antibody-protein antigen recognition is of fundamental and practical significance. We constructed an antibody structural dataset, partitioned it into human and murine subgroups, and compared it with nonantibody protein-protein complexes. We investigated the physicochemical properties of regions on and away from the antibody-antigen interfaces, including net charge, overall antibody charge distributions, and their potential role in antigen interaction. We observed that amino acid preference in antibody-protein antigen recognition is entropy driven, with residues having low side-chain entropy appearing to compensate for the high backbone entropy in interaction with protein antigens. Antibodies prefer charged and polar antigen residues and bridging water molecules. They also prefer positive net charge, presumably to promote interaction with negatively charged protein antigens, which are common in proteomes. Antibody-antigen interfaces have large percentages of Tyr, Ser, and Asp, but little Lys. Electrostatic and hydrophobic interactions in the Ag binding sites might be coupled with Fab domains through organized charge and residue distributions away from the binding interfaces. Here we describe some features of antibody-antigen interfaces and of Fab domains as compared with nonantibody protein-protein interactions. The distributions of interface residues in human and murine antibodies do not differ significantly. Overall, our results provide not only a local but also a global anatomy of antibody structures.

Bioisostere Identification by Determining the Amino Acid Binding Preferences of Common Chemical Fragments.

To assist in the structural optimization of hit/lead compounds during drug discovery, various computational approaches to identify potentially useful bioisosteric conversions have been reported. Here, the preference of chemical fragments to hydrogen bonds with specific amino acid residues was used to identify potential bioisosteric conversions. We first compiled a data set of chemical fragments frequently occurring in complex structures contained in the Protein Data Bank. We then used a computational approach to determine the amino acids to which these chemical fragments most frequently hydrogen bonded. The results of the frequency analysis were used to hierarchically cluster chemical fragments according to their amino acid preferences. The Euclid distance between amino acid preferences of chemical fragments for hydrogen bonding was then compared to MMP information in the ChEMBL database. To demonstrate the applicability of the approach for compound optimization, the similarity of amino acid preferences was used to identify known bioisosteric conversions of the epidermal growth factor receptor inhibitor gefitinib. The amino acid preference distance successfully detected bioisosteric fragments corresponding to the morpholine ring in gefitinib with a higher ROC score compared to those based on topological similarity of substituents and frequency of MMP in the ChEMBL database.

Quality Characteristics of Mealworm Seasoning Paste prepared by Chungukjang Recipe

In this study, two kinds of mealworm seasoning pastes (naturally fermented or inoculation of Bacillus subtilis) were prepared by applying fermentation method of chungkukjang. The samples were labeled into soybean control (SC; no fermentation), mealworm control (MC; no fermentation), mealworm natural fermentation (MR) and mealworm inoculated fermentation of Bacillus (MB). Water content, color, pH, titratable acidity, salinity, amino acid nitrogen, reducing sugar and preference test were examined on the pastes. Moisture contents were 46.9 with soybean control (MC) and 54.5～73.3 with mealworm pastes. The L-value was not significantly different between the samples. In the case of the a-value and the b-value, SC showed the highest value whereas mealworm pastes showed lower values than SC. pH and salinity were higher in mealworm pastes. Amino nitrogen and reducing sugar showed the highest values of MR 326,67 and 63.33 ㎎%, respectively. In the preference test, the appearance and flavor were the highest values with SC, but texture, taste and overall acceptability were with MR, Therefore, soybean fermentation method with Bacillus species could be applied to use protein of mealworm and possibly edible insects as well.

Improved Modeling of Compositional Heterogeneity Supports Sponges as Sister to All Other Animals

The relationships at the root of the animal tree have proven difficult to resolve, with the current debate focusing on whether sponges (phylum Porifera) or comb jellies (phylum Ctenophora) are the sister group of all other animals [1-5]. The choice of evolutionary models seems to be at the core of the problem because Porifera tends to emerge as the sister group of all other animals ("Porifera-sister") when site-specific amino acid differences are modeled (e.g., [6, 7]), whereas Ctenophora emerges as the sister group of all other animals ("Ctenophora-sister") when they are ignored (e.g., [8-11]). We show that two key phylogenomic datasets that previously supported Ctenophora-sister [10, 12] display strong heterogeneity in amino acid composition across sites and taxa and that no routinely used evolutionary model can adequately describe both forms of heterogeneity. We show that data-recoding methods [13-15] reduce compositional heterogeneity in these datasets and that models accommodating site-specific amino acid preferences can better describe the recoded datasets. Increased model adequacy is associated with significant topological changes insupport of Porifera-sister. Because adequate modeling of the evolutionary process that generated the data is fundamental to recovering an accurate phylogeny [16-20], our results strongly support sponges as the sister group of all other animals and provide further evidence that Ctenophora-sister represents a tree reconstruction artifact. VIDEO ABSTRACT.

Crotalus durissus terrificus crotapotin naturally displays preferred positions for amino acid substitutions

BackgroundClassically, Crotalus durissus terrificus (Cdt) venom can be described, according to chromatographic criteria, as a simple venom, composed of four major toxins, namely: gyroxin, crotamine, crotoxin and convulxin. Crotoxin is a non-covalent heterodimeric neurotoxin constituted of two subunits: an active phospholipase A2 and a chaperone protein, termed crotapotin. This molecule is composed of three peptide chains connected by seven disulfide bridges. Naturally occurring variants/isoforms of either crotoxin or crotapotin itself have already been reported.MethodsThe crude Cdt venom was separated by using RP-HPLC and the toxins were identified by mass spectrometry (MS). Crotapotin was purified, reduced and alkylated in order to separate the peptide chains that were further analyzed by mass spectrometry and de novo peptide sequencing.ResultsThe RP-HPLC profile of the isolated crotapotin chains already indicated that the α chain would present isoforms, which was corroborated by the MS and tandem mass spectrometry analyses.ConclusionIt was possible to observe that the Cdt crotapotin displays a preferred amino acid substitution pattern present in the α chain, at positions 31 and 40. Moreover, substitutions could also be observed in β and γ chains (one for each). The combinations of these four different peptides, with the already described chains, would produce ten different crotapotins, which is compatible to our previous observations for the Cdt venom.

Biophysical and functional characterization of the N-terminal domain of the cat T1R1 umami taste receptor expressed in Escherichia coli.

Umami taste perception is mediated by the heterodimeric G-protein coupled receptors (GPCRs), formed by the assembly of T1R1 and T1R3 subunits. T1R1 and T1R3 subunits are class C GPCRs whose members share common structural homologies including a long N-terminal domain (NTD) linked to a seven transmembrane domain by a short cysteine-rich region. The NTD of the T1R1 subunit contains the primary binding site for umami stimuli, such as L-glutamate (L-Glu) for humans. Inosine-5’-monophosphate (IMP) binds at a location close to the opening of the T1R1-NTD “flytrap”, thus creating the observed synergistic response between L-Glu and IMP. T1R1/T1R3 binding studies have revealed species-dependent differences. While human T1R1/T1R3 is activated specifically by L-Glu, the T1R1/T1R3 in other species is a broadly tuned receptor, sensitive to a range of L-amino acids. Because domestic cats are obligate carnivores, they display strong preferences for some specific amino acids. To better understand the structural basis of umami stimuli recognition by non-human taste receptors, we measured the binding of selected amino acids to cat T1R1/T1R3 (cT1R1/cT1R3) umami taste receptor. For this purpose, we expressed cT1R1-NTD in bacteria as inclusion bodies. After purification, refolding of the protein was achieved. Circular dichroism spectroscopic studies revealed that cT1R1-NTD was well renatured with evidence of secondary structures. Using size-exclusion chromatography coupled to light scattering, we found that the cT1R1-NTD behaves as a monomer. Ligand binding quantified by intrinsic tryptophan fluorescence showed that cT1R1-NTD is capable of binding L-amino acids with Kd values in the micromolar range. We demonstrated that IMP potentiates L-amino acid binding onto renatured cT1R1-NTD. Interestingly, our results revealed that IMP binds the extracellular domain in the absence of L-amino acids. Thus, this study demonstrates that the feasibility to produce milligram quantities of cT1R1-NTD for functional and structural studies.