Amino Acid Frequencies Research Articles

Optimal computational comparison of mass spectrometric peptide profiles of alternative hydrolysates from the same starting material

Food protein hydrolysates are complex mixtures that are increasingly being analysed by tandem mass spectrometry. Given a single starting material, many alternative peptide profiles are achievable under varying hydrolysis conditions. To date, characterisation of the relative similarities and differences between such peptide profiles remains largely unstudied. Here, we investigate optimal computational methods for grouping peptide profiles of hydrolysates derived from the same starting material. Using an experimental bovine milk dataset, we evaluated how these methods grouped either technical replicates, or distinct samples with known cleavage profiles. Analyses performed using only the presence and abundance of peptides were found to be suboptimal for achieving effective sample grouping. Using the amino acid distribution at both termini of peptides was more efficient at grouping replicate samples; however, this approach lacked suitable discrimination between distinct samples. By extending the termini approach to incorporate the abundance associated with terminal amino acids, optimal grouping was achieved. We therefore suggest that grouping peptide profiles of hydrolysates from the same starting material should rely on a combination of N and C terminal amino acid frequency and abundance. Importantly, this approach requires no a priori knowledge of enzyme specificities, making it generally applicable to diverse sets of food matrices.

ScaffoldSeq: Software for characterization of directed evolution populations.

ScaffoldSeq is software designed for the numerous applications-including directed evolution analysis-in which a user generates a population of DNA sequences encoding for partially diverse proteins with related functions and would like to characterize the single site and pairwise amino acid frequencies across the population. A common scenario for enzyme maturation, antibody screening, and alternative scaffold engineering involves naïve and evolved populations that contain diversified regions, varying in both sequence and length, within a conserved framework. Analyzing the diversified regions of such populations is facilitated by high-throughput sequencing platforms; however, length variability within these regions (e.g., antibody CDRs) encumbers the alignment process. To overcome this challenge, the ScaffoldSeq algorithm takes advantage of conserved framework sequences to quickly identify diverse regions. Beyond this, unintended biases in sequence frequency are generated throughout the experimental workflow required to evolve and isolate clones of interest prior to DNA sequencing. ScaffoldSeq software uniquely handles this issue by providing tools to quantify and remove background sequences, cluster similar protein families, and dampen the impact of dominant clones. The software produces graphical and tabular summaries for each region of interest, allowing users to evaluate diversity in a site-specific manner as well as identify epistatic pairwise interactions. The code and detailed information are freely available at http://research.cems.umn.edu/hackel. Proteins 2016; 84:869-874. © 2016 Wiley Periodicals, Inc.

Resolving the Phylogenetic Position of Coelacanth: The Closest Relative Is Not Always the Most Appropriate Outgroup.

Determining the phylogenetic relationship of two extant lineages of lobe-finned fish, coelacanths and lungfishes, and tetrapods is important for understanding the origin of tetrapods. We analyzed data sets from two previous studies along with a newly collected data set, each of which had varying numbers of species and genes and varying extent of missing sites. We found that in all the data sets the sister relationship of lungfish and tetrapods was constructed with the use of cartilaginous fish as the outgroup with a high degree of statistical support. In contrast, when ray-finned fish were used as the outgroup, which is taxonomically an immediate outgroup of lobe-finned fish and tetrapods, the sister relationship of coelacanth and tetrapods was supported most strongly, although the statistical support was weaker. Even though it is generally accepted that the closest relative is an appropriate outgroup, our analysis suggested that the large divergence of the ray-finned fish as indicated by their long branch lengths and different amino acid frequencies made them less suitable as an outgroup than cartilaginous fish.

Immunological evasion of immediate-early varicella zoster virus proteins.

The varicella zoster virus (VZV) causes the childhood disease commonly known as chickenpox and can later in life reactivate as herpes zoster. The adaptive immune system is known to play an important role in suppressing VZV reactivation. A central aspect of this system is the presentation of VZV-derived peptides by the major histocompatibility complex (MHC) proteins. Here, we investigate if key VZV proteins have evolved their amino acid sequence to avoid presentation by MHC based on predictive models of MHC-peptide affinity. This study shows that the immediate-early proteins of all characterized VZV strains are profoundly depleted for high-affinity MHC-I-restricted epitopes. The same depletion can be found in its closest animal analog, the simian varicella virus. Further orthology analysis towards other herpes viruses suggests that the protein amino acid frequency is one of the primary drivers of targeted epitope depletion.

Expression-Linked Patterns of Codon Usage, Amino Acid Frequency, and Protein Length in the Basally Branching Arthropod Parasteatoda tepidariorum.

Spiders belong to the Chelicerata, the most basally branching arthropod subphylum. The common house spider, Parasteatoda tepidariorum, is an emerging model and provides a valuable system to address key questions in molecular evolution in an arthropod system that is distinct from traditionally studied insects. Here, we provide evidence suggesting that codon usage, amino acid frequency, and protein lengths are each influenced by expression-mediated selection in P. tepidariorum. First, highly expressed genes exhibited preferential usage of T3 codons in this spider, suggestive of selection. Second, genes with elevated transcription favored amino acids with low or intermediate size/complexity (S/C) scores (glycine and alanine) and disfavored those with large S/C scores (such as cysteine), consistent with the minimization of biosynthesis costs of abundant proteins. Third, we observed a negative correlation between expression level and coding sequence length. Together, we conclude that protein-coding genes exhibit signals of expression-related selection in this emerging, noninsect, arthropod model.

High Throughput Sequencing of T Cell Antigen Receptors Reveals a Conserved TCR Repertoire.

The T-cell receptor (TCR) repertoire is a mirror of the human immune system that reflects processes caused by infections, cancer, autoimmunity, and aging. Next-generation sequencing has become a powerful tool for deep TCR profiling. Herein, we used this technology to study the repertoire features of TCR beta chain in the blood of healthy individuals.Peripheral blood samples were collected from 10 healthy donors. T cells were isolated with anti-human CD3 magnetic beads according to the manufacturer's protocol. We then combined multiplex-PCR, Illumina sequencing, and IMGT/High V-QUEST to analyze the characteristics and polymorphisms of the TCR.Most of the individual T cell clones were present at very low frequencies, suggesting that they had not undergone clonal expansion. The usage frequencies of the TCR beta variable, beta joining, and beta diversity gene segments were similar among T cells from different individuals. Notably, the usage frequency of individual nucleotides and amino acids within complementarity-determining region (CDR3) intervals was remarkably consistent between individuals. Moreover, our data show that terminal deoxynucleotidyl transferase activity was biased toward the insertion of G (31.92%) and C (27.14%) over A (21.82%) and T (19.12%) nucleotides.Some conserved features could be observed in the composition of CDR3, which may inform future studies of human TCR gene recombination.

Statistical correlation of nonconservative substitutions of HIV gp41 variable amino acid residues with the R5X4 HIV-1 phenotype.

BackgroundThe interaction of the envelope glycoprotein of HIV-1 (gp120/gp41) with coreceptor molecules has important implications for specific cellular targeting and pathogenesis. Experimental and theoretical evidences have shown a role for gp41 in coreceptor tropism, although there is no consensus about the positions involved. Here we analyze the association of physicochemical properties of gp41 amino acid residues with viral tropism (X4, R5, and R5X4) using a large set of HIV-1 sequences. Under the assumption that conserved regions define the complex structural features essential for protein function, we focused our search only on amino acids in the gp41 variable regions.MethodsGp41 amino acid sequences of 2823 HIV-1 strains from all clades with known coreceptor tropism were retrieved from Los Alamos HIV Database. Consensus sequences were constructed for homologous sequences (those obtained from the same patient and having the same tropism) in order to avoid bias due to sequence overrepresentation, and the variability (entropy) per site was determined. Comparisons of hydropathy index (HI) and charge (Q) of amino acid residues at highly variable positions between coreceptor groups were performed using two non-parametrical tests and Benjamini-Hochberg correction. Pearson’s correlation analysis was performed to determine covariance of HI and Q values.ResultsCalculation of variability per site rendered 58 highly variable amino acid positions. Of these, statistical analysis rendered significantly different HI or Q only for the R5 vs. R5X4 comparison at twelve positions: 535, 602, 619, 636, 640, 641, 658, 662, 667, 723, 756 and 841. The largest differences in particular amino acid frequencies between coreceptor groups were found at 619, 636, 640, 641, 662, 723 and 756. A hydrophobic tendency of residues 619, 640, 641, 723 and 756, along with a hydrophilic/charged tendency at residues 636 and 662 was observed in R5X4 with respect to R5 sequences. HI of position 640 covariated with that of 602, 619, 636, 662, and 756.ConclusionsVariability and significant correlations of physicochemical properties with viral phenotype suggest that substitutions at residues in the loop (602 and 619), the HR2 (636, 640, 641, 662), and the C-terminal tail (723, 756) of gp41 may contribute to phenotype of R5X4 strains.Electronic supplementary materialThe online version of this article (doi:10.1186/s12985-016-0486-6) contains supplementary material, which is available to authorized users.

Increasing genomic diversity and evidence of constrained lifestyle evolution due to insertion sequences in Aeromonas salmonicida.

BackgroundAeromonads make up a group of Gram-negative bacteria that includes human and fish pathogens. The Aeromonas salmonicida species has the peculiarity of including five known subspecies. However, few studies of the genomes of A. salmonicida subspecies have been reported to date.ResultsWe sequenced the genomes of additional A. salmonicida isolates, including three from India, using next-generation sequencing in order to gain a better understanding of the genomic and phylogenetic links between A. salmonicida subspecies. Their relative phylogenetic positions were confirmed by a core genome phylogeny based on 1645 gene sequences. The Indian isolates, which formed a sub-group together with A. salmonicida subsp. pectinolytica, were able to grow at either at 18 °C and 37 °C, unlike the A. salmonicida psychrophilic isolates that did not grow at 37 °C. Amino acid frequencies, GC content, tRNA composition, loss and gain of genes during evolution, pseudogenes as well as genes under positive selection and the mobilome were studied to explain this intraspecies dichotomy.ConclusionInsertion sequences appeared to be an important driving force that locked the psychrophilic strains into their particular lifestyle in order to conserve their genomic integrity. This observation, based on comparative genomics, is in agreement with previous results showing that insertion sequence mobility induced by heat in A. salmonicida subspecies causes genomic plasticity, resulting in a deleterious effect on the virulence of the bacterium. We provide a proof-of-concept that selfish DNAs play a major role in the evolution of bacterial species by modeling genomes.Electronic supplementary materialThe online version of this article (doi:10.1186/s12864-016-2381-3) contains supplementary material, which is available to authorized users.

Data mining and analysis of bacillus virus

Anthrax- can be found naturally in soil and commonly affects domestic and wild animals around the world. Bacillus anthracis is mostly common viruses of Anthrax. There are some more occurred in similar DNA sequences: 3 Bacillus Viruses: Bacillus anthracis, Bacillus cereus, and Bacillus thuringiensis. This is a report about analyzing the similarities between the Bacillus viruses by investigating the frequency of amino acid and finding the difference between those three viruses based on the gene. Those diseases are infected by parasite and host animals, and cause muscle pain. Therefore, we can conclude that Leucine is a protein that plays a significant role in causing muscle pain. Secondly, In Analysis of decision tree, there are only little differences between each classes. The classes represent positions that include representative protein. The very position that windows mention is their difference.

Analysis of the Repertoire Features of TCR Beta Chain CDR3 in Human by High-Throughput Sequencing

Background/Aims: To ward off a wide variety of pathogens, the human adaptive immune system harbors a vast array of T-cell receptors, collectively referred to as the TCR repertoire. Assessment of the repertoire features of TCR is vital for us to deeper understand of immune behaviour and immune response. Methods: In this study, we used a combination of multiplex-PCR, Illumina sequencing and IMGT (ImMunoGeneTics)/HighV-QUEST for a standardized analysis of the repertoire features of TCR beta chain in the blood of healthy individuals, including the repertoire features of public TCR complementarity-determining regions (CDR3) sequences, highly expanded clones, long TCR CDR3 sequences. Results: We found that public CDR3 sequences and high-frequency sequences had the same characteristics, both of them had fewer nucleotide additions and shorter CDR3 length, which were closer to the germline sequence. Moreover, our studies provided evidence that public amino acid sequences are produced by multiple nucleotide sequences. Notably, there was skewed VDJ segment usage in long CDR3 sequences, the expression levels of 10 TRβV segments, 7 TRβJ segments and 2 TRβD segments were significantly different in the long CDR3 sequences compared to the short CDR3 sequences. Moreover, we identified that extensive N additions and increase of D gene usage contributing to TCR CDR3 length, and observed there was distinct usage frequency of amino acids in long CDR3 sequences compared to the short CDR3 sequences. Conclusions: Some repertoire features could be observed in the public sequences, highly abundance clones, and long TCR CDR3 sequences, which might be helpful for further study of immune behavior and immune response.

Prediction of protein subcellular multisite localization using a new feature extraction method.

A basic problem of proteomics is identifying the subcellular locations of a protein. One factor making the problem more complicated is that some proteins may simultaneously exist in two or more than two subcellular locations. To improve multisite prediction quality, it is necessary to use effective feature extraction methods. Here, we developed a new feature extraction method based on the pK value and frequencies of amino acids to represent a protein as a real values vector. Using this novel feature extraction method, the multi-label k-nearest neighbors (ML-KNN) algorithm and setting different weights into different attributes' ML-KNN, known as wML-KNN, were employed to predict multiplex protein subcellular locations. The best overall accuracy rate on dataset S1 from the predictor of Virus-mPLoc was 59.92 and 86.04% on dataset S2 from Gpos-mPLoc, respectively.

Genome-Wide Relative Analysis of Codon Usage Bias and Codon Context Pattern in the Bacteria Salinibacter Ruber, Chromohalobacter Salexigens and Rhizobium Etli

Codon is the basic unit for biological message transmission during synthesis of proteins in an organism. Codon Usage Bias is preferential usage among synonymous codons, in an organisms. This preferential use of a synonymous codon was found not only among species but also occurs among genes within the same genome of a species. This variation of codon usage patterns are controlled by natural processes such as mutation, drift and pressure. In this study, we have used computational as well as statistical techniques for finding codon usage bias and codon context pattern of Salinibacter ruber (extreme halophilic), Chromohalobacter salexigens (moderate halophilic) and Rhizobium etli (nonhalophilic). In addition to this, compositional variation in translated amino acid frequency, effective number of codons and optimal codons were also studied. A plot of ENc versus GC3s suggests that both mutation bias and translational selection contribute to these differences of codon bias. However, mutation bias is the driving force of the synonymous codon usage patterns in halophilic bacteria (Salinibacter ruber and Chromohalobacter salexigens) and translational selection seems to affect codon usage pattern in non-halophilic bacteria (Rhizobium etli). Correspondence analysis of Relative Synonymous Codon Usage revealed different clusters of genes varying in numbers in the bacteria under study. Moreover, codon context pattern was also seen variable in these bacteria. These results clearly indicate the variation in the codon usage pattern in these bacterial genomes.

Codon and Amino Acid Usage Are Shaped by Selection Across Divergent Model Organisms of the Pancrustacea.

In protein-coding genes, synonymous codon usage and amino acid composition correlate to expression in some eukaryotes, and may result from translational selection. Here, we studied large-scale RNA-seq data from three divergent arthropod models, including cricket (Gryllus bimaculatus), milkweed bug (Oncopeltus fasciatus), and the amphipod crustacean Parhyale hawaiensis, and tested for optimization of codon and amino acid usage relative to expression level. We report strong signals of AT3 optimal codons (those favored in highly expressed genes) in G. bimaculatus and O. fasciatus, whereas weaker signs of GC3 optimal codons were found in P. hawaiensis, suggesting selection on codon usage in all three organisms. Further, in G. bimaculatus and O. fasciatus, high expression was associated with lowered frequency of amino acids with large size/complexity (S/C) scores in favor of those with intermediate S/C values; thus, selection may favor smaller amino acids while retaining those of moderate size for protein stability or conformation. In P. hawaiensis, highly transcribed genes had elevated frequency of amino acids with large and small S/C scores, suggesting a complex dynamic in this crustacean. In all species, the highly transcribed genes appeared to favor short proteins, high optimal codon usage, specific amino acids, and were preferentially involved in cell-cycling and protein synthesis. Together, based on examination of 1,680,067, 1,667,783, and 1,326,896 codon sites in G. bimaculatus, O. fasciatus, and P. hawaiensis, respectively, we conclude that translational selection shapes codon and amino acid usage in these three Pancrustacean arthropods.

Composition of Overlapping Protein-Protein and Protein-Ligand Interfaces.

Protein-protein interactions (PPIs) play a major role in many biological processes and they represent an important class of targets for therapeutic intervention. However, targeting PPIs is challenging because often no convenient natural substrates are available as starting point for small-molecule design. Here, we explored the characteristics of protein interfaces in five non-redundant datasets of 174 protein-protein (PP) complexes, and 161 protein-ligand (PL) complexes from the ABC database, 436 PP complexes, and 196 PL complexes from the PIBASE database and a dataset of 89 PL complexes from the Timbal database. In all cases, the small molecule ligands must bind at the respective PP interface. We observed similar amino acid frequencies in all three datasets. Remarkably, also the characteristics of PP contacts and overlapping PL contacts are highly similar.

Conformational Parameters for Amino Acids in Helical, β-Sheet, and Random Coil Regions Calculated from Proteins: After 40 Years

Forty years ago, Peter Y. Chou and Gerald D. Fasman (1974a), relying on the information from fifteen proteins calculated α helix, β-sheet, and coil conformational parameters, Pα,Pβ,and Pc,for the 20 naturally occurring amino acids from the frequency of occurrence of each amino acid residue in the α, β, and coil conformations. Secondary structure of these 15 proteins had been determined by X-ray crystallography. Although the accuracy could not go over to the level of 60% too much, these values utilized for a long time to provide a simple procedure, devoid of complex computer calculations, to predict the secondary structure of proteins from their known amino acid sequences. In the same article of Peter Y. Chou and Gerald D. Fasman, a detailed analysis of the helix and β-sheet boundary residues in proteins provided amino acid frequencies at the N-and C-terminal ends which were used to delineate helical and β regions. Charged residues are found with the greatest frequency at both helical ends, but they were mostly absent in β-sheet regions. In the same article a mechanism of protein folding was proposed, whereby helix nucleation starts at the centers of the helix where the Pα values are highest, and propagates in both directions, until strong helix breakers where Pα values are lowest, terminate the growth at both ends. Similarly, residues with the highest Pβ values will initiate β regions and residues with the lowest Pβ values will terminate β regions. The helical region with the largest Pα was proposed as the site of the first fold during protein renaturation. The mechanism whereby proteins fold into their native conformation, capable of biological activity, has been a long sought after goal. With the elucidation of the three-dimensional structure of many proteins through X-ray crystallography, a new momentum has been given to understanding the factors governing this complex assembly of polypeptide chains. In this paper, using similar statistics from 20 347 proteins, the level of reliability of formerly found results is discussed.

Selection of specific interactors from phage display library based on sea lamprey variable lymphocyte receptor sequences

Variable lymphocyte receptors (VLRs) are non-immunoglobulin components of adaptive immunity in jawless vertebrates. These proteins composed of leucine-rich repeat modules offer some advantages over antibodies in target binding and therefore are attractive candidates for biotechnological applications. In this paper we report the design and characterization of a phage display library based on a previously proposed dVLR scaffold containing six LRR modules [Wezner-Ptasinska et al., 2011]. Our library was designed based on a consensus approach in which the randomization scheme reflects the frequencies of amino acids naturally occurring in respective positions responsible for antigen recognition. We demonstrate general applicability of the scaffold by selecting dVLRs specific for lysozyme and S100A7 protein with KD values in the micromolar range. The dVLR library could be used as a convenient alternative to antibodies for effective isolation of high affinity binders.

Cardiolipin Interactions with Proteins

Cardiolipins (CL) represent unique phospholipids of bacteria and eukaryotic mitochondria with four acyl chains and two phosphate groups that have been implicated in numerous functions from energy metabolism to apoptosis. Many proteins are known to interact with CL, and several cocrystal structures of protein-CL complexes exist. In this work, we describe the collection of the first systematic and, to the best of our knowledge, the comprehensive gold standard data set of all known CL-binding proteins. There are 62 proteins in this data set, 21 of which have nonredundant crystal structures with bound CL molecules available. Using binding patch analysis of amino acid frequencies, secondary structures and loop supersecondary structures considering phosphate and acyl chain binding regions together and separately, we gained a detailed understanding of the general structural and dynamic features involved in CL binding to proteins. Exhaustive docking of CL to all known structures of proteins experimentally shown to interact with CL demonstrated the validity of the docking approach, and provides a rich source of information for experimentalists who may wish to validate predictions.

In Silico Study of Variable Surface Proteins in Plasmodium Species: Perspectives in Drug Design.

The variable surface proteins expressed by P. falciparum and P. vivax are transported to the surface of infected erythrocyte and are exposed to the host immune system. The possibility of using variable surface proteins as a common drug target has been analyzed in both the Plasmodium species. Sequence analysis of variable surface proteins showed a low-level conservation within as well as between the species. Amino acid composition analysis revealed higher frequency of hydrophilic amino acids as compared with that of hydrophobic residues. In order to gain more insight into their diverse functional role, the three-dimensional structure was predicted using comparative modeling approach. These models were evaluated and validated by checking stereochemistry of underlying amino acids. Structural alignment of variable surface proteins by superimposing them shows less conservation. Due to differences at sequence as well as structural level, the variable surface proteins are expected to show difference in their degree of invasiveness. These differences were also cross-examined by evolutionary study, and the results obtained were in accordance with the aforesaid study. The existence of structural differences noticed in the present study showed that the variable surface proteins could not be used as a common drug target in both the malarial species. Therefore, species-specific strategy may be followed for drug targeting against variable surface proteins of P. falciparum and P. vivax.

Massively parallel sampling of lattice proteins reveals foundations of thermal adaptation.

Evolution of proteins in bacteria and archaea living in different conditions leads to significant correlations between amino acid usage and environmental temperature. The origins of these correlations are poorly understood, and an important question of protein theory, physics-based prediction of types of amino acids overrepresented in highly thermostable proteins, remains largely unsolved. Here, we extend the random energy model of protein folding by weighting the interaction energies of amino acids by their frequencies in protein sequences and predict the energy gap of proteins designed to fold well at elevated temperatures. To test the model, we present a novel scalable algorithm for simultaneous energy calculation for many sequences in many structures, targeting massively parallel computing architectures such as graphics processing unit. The energy calculation is performed by multiplying two matrices, one representing the complete set of sequences, and the other describing the contact maps of all structural templates. An implementation of the algorithm for the CUDA platform is available at http://www.github.com/kzeldovich/galeprot and calculates protein folding energies over 250 times faster than a single central processing unit. Analysis of amino acid usage in 64-mer cubic lattice proteins designed to fold well at different temperatures demonstrates an excellent agreement between theoretical and simulated values of energy gap. The theoretical predictions of temperature trends of amino acid frequencies are significantly correlated with bioinformatics data on 191 bacteria and archaea, and highlight protein folding constraints as a fundamental selection pressure during thermal adaptation in biological evolution.

NATURE’S STARSHIPS. II. SIMULATING THE SYNTHESIS OF AMINO ACIDS IN METEORITE PARENT BODIES

Carbonaceous chondrite meteorites are known for having high water and organic material contents, including amino acids. Here we address the origin of amino acids in the warm interiors of their parent bodies (planetesimals) within a few million years of their formation, and connect this with the astrochemistry of their natal protostellar disks. We compute both the total amino acid abundance pattern as well as the relative frequencies of amino acids within the CM2 (e.g. Murchison) and CR2 chondrite subclasses based on Strecker reactions within these bodies. We match the relative frequencies to well within an order of magnitude among both CM2 and CR2 meteorites for parent body temperatures $<$ 200$^{\circ}$C. These temperatures agree with 3D models of young planetesimal interiors. We find theoretical abundances of approximately 7x10$^5$ parts-per-billion (ppb), which is in agreement with the average observed abundance in CR2 meteorites of 4$\pm$7x10$^5$, but an order of magnitude higher than the average observed abundance in CM2 meteorites of 2$\pm$2x10$^4$. We find that the production of hydroxy acids could be favoured over the production of amino acids within certain meteorite parent bodies (e.g. CI1, CM2) but not others (e.g. CR2). This could be due to the relatively lower NH$_3$ abundances within CI1 and CM2 meteorite parent bodies, which leads to less amino acid synthesis. We also find that the water content in planetesimals is likely to be the main cause of variance between carbonaceous chondrites of the same subclass. We propose that amino acid abundances are primarily dependent on the ammonia and water content of planetesimals that are formed in chemically distinct regions within their natal protostellar disks.