scholarly journals Residue Geometry Networks: A Rigidity-Based Approach to the Amino Acid Network and Evolutionary Rate Analysis

2016 ◽  
Vol 6 (1) ◽  
Author(s):  
Alexander S. Fokas ◽  
Daniel J. Cole ◽  
Sebastian E. Ahnert ◽  
Alex W. Chin

Abstract Amino acid networks (AANs) abstract the protein structure by recording the amino acid contacts and can provide insight into protein function. Herein, we describe a novel AAN construction technique that employs the rigidity analysis tool, FIRST, to build the AAN, which we refer to as the residue geometry network (RGN). We show that this new construction can be combined with network theory methods to include the effects of allowed conformal motions and local chemical environments. Importantly, this is done without costly molecular dynamics simulations required by other AAN-related methods, which allows us to analyse large proteins and/or data sets. We have calculated the centrality of the residues belonging to 795 proteins. The results display a strong, negative correlation between residue centrality and the evolutionary rate. Furthermore, among residues with high closeness, those with low degree were particularly strongly conserved. Random walk simulations using the RGN were also successful in identifying allosteric residues in proteins involved in GPCR signalling. The dynamic function of these residues largely remain hidden in the traditional distance-cutoff construction technique. Despite being constructed from only the crystal structure, the results in this paper suggests that the RGN can identify residues that fulfil a dynamical function.

GigaScience ◽  
2019 ◽  
Vol 8 (10) ◽  
Author(s):  
Xiao Hu ◽  
Iddo Friedberg

Abstract Background Gene homology type classification is required for many types of genome analyses, including comparative genomics, phylogenetics, and protein function annotation. Consequently, a large variety of tools have been developed to perform homology classification across genomes of different species. However, when applied to large genomic data sets, these tools require high memory and CPU usage, typically available only in computational clusters. Findings Here we present a new graph-based orthology analysis tool, SwiftOrtho, which is optimized for speed and memory usage when applied to large-scale data. SwiftOrtho uses long k-mers to speed up homology search, while using a reduced amino acid alphabet and spaced seeds to compensate for the loss of sensitivity due to long k-mers. In addition, it uses an affinity propagation algorithm to reduce the memory usage when clustering large-scale orthology relationships into orthologous groups. In our tests, SwiftOrtho was the only tool that completed orthology analysis of proteins from 1,760 bacterial genomes on a computer with only 4 GB RAM. Using various standard orthology data sets, we also show that SwiftOrtho has a high accuracy. Conclusions SwiftOrtho enables the accurate comparative genomic analyses of thousands of genomes using low-memory computers. SwiftOrtho is available at https://github.com/Rinoahu/SwiftOrtho


2017 ◽  
Vol 7 (1) ◽  
Author(s):  
Alexander S. Fokas ◽  
Daniel J. Cole ◽  
Sebastian E. Ahnert ◽  
Alex W. Chin

2020 ◽  
Vol 16 (4) ◽  
pp. 451-459 ◽  
Author(s):  
Fortunatus C. Ezebuo ◽  
Ikemefuna C. Uzochukwu

Background: Sulfotransferase family comprises key enzymes involved in drug metabolism. Oxamniquine is a pro-drug converted into its active form by schistosomal sulfotransferase. The conformational dynamics of side-chain amino acid residues at the binding site of schistosomal sulfotransferase towards activation of oxamniquine has not received attention. Objective: The study investigated the conformational dynamics of binding site residues in free and oxamniquine bound schistosomal sulfotransferase systems and their contribution to the mechanism of oxamniquine activation by schistosomal sulfotransferase using molecular dynamics simulations and binding energy calculations. Methods: Schistosomal sulfotransferase was obtained from Protein Data Bank and both the free and oxamniquine bound forms were subjected to molecular dynamics simulations using GROMACS-4.5.5 after modeling it’s missing amino acid residues with SWISS-MODEL. Amino acid residues at its binding site for oxamniquine was determined and used for Principal Component Analysis and calculations of side-chain dihedrals. In addition, binding energy of the oxamniquine bound system was calculated using g_MMPBSA. Results: The results showed that binding site amino acid residues in free and oxamniquine bound sulfotransferase sampled different conformational space involving several rotameric states. Importantly, Phe45, Ile145 and Leu241 generated newly induced conformations, whereas Phe41 exhibited shift in equilibrium of its conformational distribution. In addition, the result showed binding energy of -130.091 ± 8.800 KJ/mol and Phe45 contributed -9.8576 KJ/mol. Conclusion: The results showed that schistosomal sulfotransferase binds oxamniquine by relying on hybrid mechanism of induced fit and conformational selection models. The findings offer new insight into sulfotransferase engineering and design of new drugs that target sulfotransferase.


2021 ◽  
Vol 11 (1) ◽  
Author(s):  
Kiyoto Kamagata ◽  
Rika Chiba ◽  
Ichiro Kawahata ◽  
Nanako Iwaki ◽  
Saori Kanbayashi ◽  
...  

AbstractLiquid droplets of aggregation-prone proteins, which become hydrogels or form amyloid fibrils, are a potential target for drug discovery. In this study, we proposed an experiment-guided protocol for characterizing the design grammar of peptides that can regulate droplet formation and aggregation. The protocol essentially involves investigation of 19 amino acid additives and polymerization of the identified amino acids. As a proof of concept, we applied this protocol to fused in sarcoma (FUS). First, we evaluated 19 amino acid additives for an FUS solution and identified Arg and Tyr as suppressors of droplet formation. Molecular dynamics simulations suggested that the Arg additive interacts with specific residues of FUS, thereby inhibiting the cation–π and electrostatic interactions between the FUS molecules. Second, we observed that Arg polymers promote FUS droplet formation, unlike Arg monomers, by bridging the FUS molecules. Third, we found that the Arg additive suppressed solid aggregate formation of FUS, while Arg polymer enhanced it. Finally, we observed that amyloid-forming peptides induced the conversion of FUS droplets to solid aggregates of FUS. The developed protocol could be used for the primary design of peptides controlling liquid droplets and aggregates of proteins.


2021 ◽  
Vol 0 (0) ◽  
Author(s):  
Pablo Mier ◽  
Miguel A. Andrade-Navarro

Abstract According to the amino acid composition of natural proteins, it could be expected that all possible sequences of three or four amino acids will occur at least once in large protein datasets purely by chance. However, in some species or cellular context, specific short amino acid motifs are missing due to unknown reasons. We describe these as Avoided Motifs, short amino acid combinations missing from biological sequences. Here we identify 209 human and 154 bacterial Avoided Motifs of length four amino acids, and discuss their possible functionality according to their presence in other species. Furthermore, we determine two Avoided Motifs of length three amino acids in human proteins specifically located in the cytoplasm, and two more in secreted proteins. Our results support the hypothesis that the characterization of Avoided Motifs in particular contexts can provide us with information about functional motifs, pointing to a new approach in the use of molecular sequences for the discovery of protein function.


2021 ◽  
Vol 12 (1) ◽  
Author(s):  
Erna Davydova ◽  
Tadahiro Shimazu ◽  
Maren Kirstin Schuhmacher ◽  
Magnus E. Jakobsson ◽  
Hanneke L. D. M. Willemen ◽  
...  

AbstractPost-translational methylation plays a crucial role in regulating and optimizing protein function. Protein histidine methylation, occurring as the two isomers 1- and 3-methylhistidine (1MH and 3MH), was first reported five decades ago, but remains largely unexplored. Here we report that METTL9 is a broad-specificity methyltransferase that mediates the formation of the majority of 1MH present in mouse and human proteomes. METTL9-catalyzed methylation requires a His-x-His (HxH) motif, where “x” is preferably a small amino acid, allowing METTL9 to methylate a number of HxH-containing proteins, including the immunomodulatory protein S100A9 and the NDUFB3 subunit of mitochondrial respiratory Complex I. Notably, METTL9-mediated methylation enhances respiration via Complex I, and the presence of 1MH in an HxH-containing peptide reduced its zinc binding affinity. Our results establish METTL9-mediated 1MH as a pervasive protein modification, thus setting the stage for further functional studies on protein histidine methylation.


Genetics ◽  
2000 ◽  
Vol 155 (1) ◽  
pp. 431-449 ◽  
Author(s):  
Ziheng Yang ◽  
Rasmus Nielsen ◽  
Nick Goldman ◽  
Anne-Mette Krabbe Pedersen

AbstractComparison of relative fixation rates of synonymous (silent) and nonsynonymous (amino acid-altering) mutations provides a means for understanding the mechanisms of molecular sequence evolution. The nonsynonymous/synonymous rate ratio (ω = dN/dS) is an important indicator of selective pressure at the protein level, with ω = 1 meaning neutral mutations, ω < 1 purifying selection, and ω > 1 diversifying positive selection. Amino acid sites in a protein are expected to be under different selective pressures and have different underlying ω ratios. We develop models that account for heterogeneous ω ratios among amino acid sites and apply them to phylogenetic analyses of protein-coding DNA sequences. These models are useful for testing for adaptive molecular evolution and identifying amino acid sites under diversifying selection. Ten data sets of genes from nuclear, mitochondrial, and viral genomes are analyzed to estimate the distributions of ω among sites. In all data sets analyzed, the selective pressure indicated by the ω ratio is found to be highly heterogeneous among sites. Previously unsuspected Darwinian selection is detected in several genes in which the average ω ratio across sites is <1, but in which some sites are clearly under diversifying selection with ω > 1. Genes undergoing positive selection include the β-globin gene from vertebrates, mitochondrial protein-coding genes from hominoids, the hemagglutinin (HA) gene from human influenza virus A, and HIV-1 env, vif, and pol genes. Tests for the presence of positively selected sites and their subsequent identification appear quite robust to the specific distributional form assumed for ω and can be achieved using any of several models we implement. However, we encountered difficulties in estimating the precise distribution of ω among sites from real data sets.


2020 ◽  
Vol 21 (12) ◽  
pp. 4447
Author(s):  
Pedro A. Lazo ◽  
Juan L. García ◽  
Paulino Gómez-Puertas ◽  
Íñigo Marcos-Alcalde ◽  
Cesar Arjona ◽  
...  

Complex neurodevelopmental syndromes frequently have an unknown etiology, in which genetic factors play a pathogenic role. This study utilizes whole-exome sequencing (WES) to examine four members of a family with a son presenting, since birth, with epileptic-like crises, combined with cerebral palsy, severe neuromotor and developmental delay, dystonic tetraparexia, axonal motor affectation, and hyper-excitability of unknown origin. The WES study detected within the patient a de novo heterozygous in-frame duplication of thirty-six nucleotides within exon 7 of the human KCNQ2 gene. This insertion duplicates the first twelve amino acids of the calmodulin binding site I. Molecular dynamics simulations of this KCNQ2 peptide duplication, modelled on the 3D structure of the KCNQ2 protein, suggest that the duplication may lead to the dysregulation of calcium inhibition of this protein function.


2015 ◽  
Vol 24 (4) ◽  
pp. 197-205
Author(s):  
Dwi Wulandari ◽  
Lisnawati Rachmadi ◽  
Tjahjani M. Sudiro

Background: E6 and E7 are oncoproteins of HPV16. Natural amino acid variation in HPV16 E6 can alter its carcinogenic potential. The aim of this study was to analyze phylogenetically E6 and E7 genes and proteins of HPV16 from Indonesia and predict the effects of single amino acid substitution on protein function. This analysis could be used to reduce time, effort, and research cost as initial screening in selection of protein or isolates to be tested in vitro or in vivo.Methods: In this study, E6 and E7 gene sequences were obtained from 12 samples of  Indonesian isolates, which  were compared with HPV16R (prototype) and 6 standard isolates in the category of European (E), Asian (As), Asian-American (AA), African-1 (Af-1), African-2 (Af-2), and North American (NA) branch from Genbank. Bioedit v.7.0.0 was used to analyze the composition and substitution of single amino acids. Phylogenetic analysis of E6 and E7 genes and proteins was performed using Clustal X (1.81) and NJPLOT softwares. Effects of single amino acid substitutions on protein function of E6 and E7 were analysed by SNAP.Results: Java variants and isolate ui66* belonged to European branch, while the others belonged to Asian and African branches. Twelve changes of amino acids were found in E6 and one in E7 proteins. SNAP analysis showed two non neutral mutations, i.e. R10I and C63G in E6 proteins. R10I mutations were found in Af-2 genotype (AF472509) and Indonesian isolates (Af2*), while C63G mutation was found only in Af2*.Conclusion: E6 proteins of HPV16 variants were more variable than E7. SNAP analysis showed that only E6 protein of African-2 branch had functional differences compared to HPV16R.


Sign in / Sign up

Export Citation Format

Share Document