Amino acid distributions around O-linked glycosylation sites.
Ontology highlight
ABSTRACT:
SUBMITTER: Wilson IB
PROVIDER: S-EPMC1150083 | biostudies-other | 1991 Apr
REPOSITORIES: biostudies-other
Similar Datasets
Project description:The CpG dinucleotide and its methylation behaviors play vital roles in gene regulation. Previous studies have divided genes into several categories based on the CpG intensity around transcription starting sites and found that housekeeping genes tend to possess high CpG density, whereas tissue-specific genes are generally characterized by low CpG density. In this study, we investigated how the CpG density distribution of a gene affects its transcription and regulation pattern. Based on the CpG density distribution around transcription starting site, by means of a semi-supervised neural network we designed, which took data augmentation into account, we divided the human genes into three categories, and genes within each cluster shared similar CpG density distribution. Not only sequence properties, these different clusters exhibited distinctly different structural features, regulatory mechanisms, correlation patterns between the expression level and CpG/TpG density, and expression and epigenetic mark variations during tumorigenesis. For instance, the activation of cluster 3 genes relies more on 3D genome reorganization, compared with cluster 1 and 2 genes, whereas cluster 2 genes showed the strongest correlation between gene expression and H3K27me3. Genes exhibiting uncoupled correlation between gene regulation and histone modifications are mainly in cluster 3. These results emphasized that the usage of epigenetic marks in gene regulation is partially rooted in the sequence property of genes such as their CpG density distribution and explained to some extent why the relation between epigenetic marks and gene expression is controversial.
| S-EPMC9113945 | biostudies-literature
Project description:BACKGROUND: While occurring enzymatically in biological systems, O-linked glycosylation affects protein folding, localization and trafficking, protein solubility, antigenicity, biological activity, as well as cell-cell interactions on membrane proteins. Catalytic enzymes involve glycotransferases, sugar-transferring enzymes and glycosidases which trim specific monosaccharides from precursors to form intermediate structures. Due to the difficulty of experimental identification, several works have used computational methods to identify glycosylation sites. RESULTS: By investigating glycosylated sites that contain various motifs between Transmembrane (TM) and non-Transmembrane (non-TM) proteins, this work presents a novel method, GlycoRBF, that implements radial basis function (RBF) networks with significant amino acid pairs (SAAPs) for identifying O-linked glycosylated serine and threonine on TM proteins and non-TM proteins. Additionally, a membrane topology is considered for reducing the false positives on glycosylated TM proteins. Based on an evaluation using five-fold cross-validation, the consideration of a membrane topology can reduce 31.4% of the false positives when identifying O-linked glycosylation sites on TM proteins. Via an independent test, GlycoRBF outperforms previous O-linked glycosylation site prediction schemes. CONCLUSION: A case study of Cyclic AMP-dependent transcription factor ATF-6 alpha was presented to demonstrate the effectiveness of GlycoRBF. Web-based GlycoRBF, which can be accessed at http://GlycoRBF.bioinfo.tw, can identify O-linked glycosylated serine and threonine effectively and efficiently. Moreover, the structural topology of Transmembrane (TM) proteins with glycosylation sites is provided to users. The stand-alone version of GlycoRBF is also available for high throughput data analysis.
| S-EPMC2989983 | biostudies-literature
Project description:ObjectiveThe current study aimed to understand epidemiological feature and critical factors associated with pathogenesis of circulating recombinant form (CRF) 01_AE strains in Northeast China.DesignCompared analysis was made between CRF01_AE and non-CRF01_AE samples to understand the pathogenicity features of CRF01_AE. Further analyses between CRF01_AE samples with high or low CD4 cell counts and between samples with different coreceptor usages were done to explore the possible factors correlating to the pathogenesis of CRF01_AE viruses.MethodsThe genotypes of newly identified strains were determined by phylogenetic analyses using Mega 6.06. Coreceptor usage was predicted by Geno2Pheno algorithm. Potential N-linked glycosylation site (PNGS) number was calculated using the online N-glycosite software. The properties of amino acid sequences were analyzed by the online ProtParam tool.ResultsCRF01_AE become the main HIV-1 genotype since 2010. Compared with non-CRF01_AE group, the CRF01_AE group showed a higher proportion of samples with CD4 cell count less than 200 cells/μl. Shorter amino acid length, fewer PNGSs and the presence of a basic motif R/KNXT or NR/KT in V4 correlated to a lower CD4 cell count, and existence or coexistence of Thr12, Arg13, Val21 and Lys33, presence of more than 4 of net charges and lack of the PNGS within V3 favored to the X4/R5X4 coreceptor usage of CRF01_AE viruses.ConclusionCRF01_AE has dominated HIV-1 genotype in Northeast China. Infection with CRF01_AE exhibited a fast disease progression, which may be associated with specific amino acid residues and PNGSs in V3 and V4 regions as well as amino acid length of V4 region.
| S-EPMC6635051 | biostudies-literature
Project description:Disulfide bonds provide an inexhaustible source of information on molecular evolution and biological specificity. In this work, we described the amino acid composition around disulfide bonds in a set of disulfide-rich proteins using appropriate descriptors, based on ANOVA (for all twenty natural amino acids or classes of amino acids clustered according to their chemical similarities) and Scheffé (for the disulfide-rich proteins superfamilies) statistics. We found that weakly hydrophilic and aromatic amino acids are quite abundant in the regions around disulfide bonds, contrary to aliphatic and hydrophobic amino acids. The density distributions (as a function of the distance to the center of the disulfide bonds) for all defined entities presented an overall unimodal behavior: the densities are null at short distances, have maxima at intermediate distances and decrease for long distances. In the end, the amino acid environment around the disulfide bonds was found to be different for different superfamilies, allowing the clustering of proteins in a biologically relevant way, suggesting that this type of chemical information might be used as a tool to assess the relationship between very divergent sets of disulfide-rich proteins.
| S-EPMC3000107 | biostudies-literature
Project description:BackgroundSince thermodynamic stability is a global property of proteins that has to be conserved during evolution, the selective pressure at a given site of a protein sequence depends on the amino acids present at other sites. However, models of molecular evolution that aim at reconstructing the evolutionary history of macromolecules become computationally intractable if such correlations between sites are explicitly taken into account.ResultsWe introduce an evolutionary model with sites evolving independently under a global constraint on the conservation of structural stability. This model consists of a selection process, which depends on two hydrophobicity parameters that can be computed from protein sequences without any fit, and a mutation process for which we consider various models. It reproduces quantitatively the results of Structurally Constrained Neutral (SCN) simulations of protein evolution in which the stability of the native state is explicitly computed and conserved. We then compare the predicted site-specific amino acid distributions with those sampled from the Protein Data Bank (PDB). The parameters of the mutation model, whose number varies between zero and five, are fitted from the data. The mean correlation coefficient between predicted and observed site-specific amino acid distributions is larger than <r> = 0.70 for a mutation model with no free parameters and no genetic code. In contrast, considering only the mutation process with no selection yields a mean correlation coefficient of <r> = 0.56 with three fitted parameters. The mutation model that best fits the data takes into account increased mutation rate at CpG dinucleotides, yielding <r> = 0.90 with five parameters.ConclusionThe effective selection process that we propose reproduces well amino acid distributions as observed in the protein sequences in the PDB. Its simplicity makes it very promising for likelihood calculations in phylogenetic studies. Interestingly, in this approach the mutation process influences the effective selection process, i.e. selection and mutation must be entangled in order to obtain effectively independent sites. This interdependence between mutation and selection reflects the deep influence that mutation has on the evolutionary process: The bias in the mutation influences the thermodynamic properties of the evolving proteins, in agreement with comparative studies of bacterial proteomes, and it also influences the rate of accepted mutations.
| S-EPMC1570368 | biostudies-literature
Project description:1. When pig heart pyruvate dehydrogenase complex was phosphorylated to completion with [gamma-32P]ATP by its intrinsic kinase, three phosphorylation sites were observed. The amino acid sequences around these sites were: sequence 1, Tyr-Gly-Met-Gly-Thr-Ser(P)-Val-Glu-Arg; and sequence 2, Tyr-His-Gly-His-Ser(P)-Met-Ser-Asp-Pro-Gly-Val-Ser(P)-Tyr-Arg. 2. When phosphorylated to inactivation by repetitive additions of limiting quantities of [gamma-32P]ATP, phosphate was incorporated mainly (more than 90%) into Ser-5 of sequence 2. Phosphorylation of this site thus results in activation of pyruvate dehydrogenase. 3. If Ser-5 is phosphorylated with ATP and the enzyme then incubated with [gamma-32P]ATP, phosphorylation of the remaining sites occurred. Ser-12 of sequence 2 is phosphorylated about twice as rapidly as Ser-6 of sequence 1. 4. Incubation of pyruvate dehydrogenase with excess [gamma-32P]ATP with termination of phosphorylation at about 50% complete inactivation showed that Ser-5 of sequence 2 was phosphorylated most rapidly, but also that Ser-12 of sequence 2 was significantly (15% of total) phosphorylated. Ser-6 sequence 1 contained about 1% total P. 5. These results suggest that addition of limiting amounts of ATP produces primarily phosphorylation of Ser-5 of sequence 2 (inactivating site). This also occurs during incubation with excess ATP before complete inactivation occurs, but a greater occupancy of other sites also occurs during this treatment.
| S-EPMC1161173 | biostudies-other
Project description:Mycobacterium tuberculosis (Mtb) causes tuberculosis, one of the leading causes of fatal infectious diseases worldwide. Cell-cell recognition between the pathogen Mtb and its host is mediated in part by glycosylated proteins. So far, glycoproteins in Mtb are understudied and for only very few glycoproteins glycosylation sites have been described, e.g., alanine and proline rich secreted protein apa, superoxide dismutase SODC, lipoprotein lpqH and MPB83/MPT83. In this study, glycosylated proteins in Mtb culture filtrate were investigated using liquid chromatography-mass spectrometry approaches and bioinformatic analyses. To validate the presence of glycoproteins, several strategies were pursued including collision induced dissociation, high energy collision dissociation and electron transfer dissociation techniques, and bioinformatics analyses involving a neutral loss search for glycosylated moieties. After extensive data curation, we report glycosylation sites for thirteen Mtb glycoproteins using a combination of mass spectrometry techniques on a dataset collected from culture filtrate proteins. This is the first glycoproteomics study identifying glycosylation sites on mycobacterial culture filtrate proteins (CFP) on a global scale.Biological significanceIn this study, glycosylation sites in Mtb were characterized by collision-induced dissociation, electron-transfer dissociation and high energy collision dissociation techniques. The identification of glycosylation sites is important for our understanding of the physiology and pathophysiology of Mtb. Glycoproteins are often responsible for protein-protein interactions between host and pathogen and thus represent interesting targets for vaccine development. In addition, our strategy is not limited to Mtb, but could be extended to other organisms. This article is part of a Special Issue entitled: Trends in Microbial Proteomics.
| S-EPMC3806883 | biostudies-literature
Project description:Our manuscript presents a novel approach to protein structure analyses. We have organized an 8-dimensional data cube with protein 3D-structural information from 8706 high-resolution non-redundant protein-chains with the aim of identifying packing rules at the amino acid pair level. The cube contains information about amino acid type, solvent accessibility, spatial and sequence distance, secondary structure and sequence length. We are able to pose structural queries to the data cube using program ProPack. The response is a 1, 2 or 3D graph. Whereas the response is of a statistical nature, the user can obtain an instant list of all PDB-structures where such pair is found. The user may select a particular structure, which is displayed highlighting the pair in question. The user may pose millions of different queries and for each one he will receive the answer in a few seconds. In order to demonstrate the capabilities of the data cube as well as the programs, we have selected well known structural features, disulphide bridges and salt bridges, where we illustrate how the queries are posed, and how answers are given. Motifs involving cysteines such as disulphide bridges, zinc-fingers and iron-sulfur clusters are clearly identified and differentiated. ProPack also reveals that whereas pairs of Lys residues virtually never appear in close spatial proximity, pairs of Arg are abundant and appear at close spatial distance, contrasting the belief that electrostatic repulsion would prevent this juxtaposition and that Arg-Lys is perceived as a conservative mutation. The presented programs can find and visualize novel packing preferences in proteins structures allowing the user to unravel correlations between pairs of amino acids. The new tools allow the user to view statistical information and visualize instantly the structures that underpin the statistical information, which is far from trivial with most other SW tools for protein structure analysis.
| S-EPMC3235099 | biostudies-literature
Project description:In many applications of evolutionary inference, a model of protein evolution needs to be fitted to the amino acid variation at individual sites in a multiple sequence alignment. Most existing models fall into one of two extremes: Either they provide a coarse-grained description that lacks biophysical realism (e.g., dN/dS models), or they require a large number of parameters to be fitted (e.g., mutation-selection models). Here, we ask whether a middle ground is possible: Can we obtain a realistic description of site-specific amino acid frequencies while severely restricting the number of free parameters in the model? We show that a distribution with a single free parameter can accurately capture the variation in amino acid frequency at most sites in an alignment, as long as we are willing to restrict our analysis to predicting amino acid frequencies by rank rather than by amino acid identity. This result holds equally well both in alignments of empirical protein sequences and of sequences evolved under a biophysically realistic all-atom force field. Our analysis reveals a near universal shape of the frequency distributions of amino acids. This insight has the potential to lead to new models of evolution that have both increased realism and a limited number of free parameters.
| S-EPMC7717668 | biostudies-literature
Project description:We present a new method for analyzing ion, or molecule, distributions around helical nucleic acids and illustrate the approach by analyzing data derived from molecular dynamics simulations. The analysis is based on the use of curvilinear helicoidal coordinates and leads to highly localized ion densities compared to those obtained by simply superposing molecular dynamics snapshots in Cartesian space. The results identify highly populated and sequence-dependent regions where ions strongly interact with the nucleic and are coupled to its conformational fluctuations. The data from this approach is presented as ion populations or ion densities (in units of molarity) and can be analyzed in radial, angular and longitudinal coordinates using 1D or 2D graphics. It is also possible to regenerate 3D densities in Cartesian space. This approach makes it easy to understand and compare ion distributions and also allows the calculation of average ion populations in any desired zone surrounding a nucleic acid without requiring references to its constituent atoms. The method is illustrated using microsecond molecular dynamics simulations for two different DNA oligomers in the presence of 0.15 M potassium chloride. We discuss the results in terms of convergence, sequence-specific ion binding and coupling with DNA conformation.
| S-EPMC4081102 | biostudies-literature