Resonant Recognition Model Defines the Secondary Structure of Bioactive Proteins

The Resonant Recognition Model (RRM) of protein bioactivity is applied to the protein secondary structure prediction. The method is based on the physical and mathematical model of the electron-ion interaction pseudopotential (EIIP) and uses signal analysis to interpret linear information contained in a macromolecular sequence. The method of analysis is based on a two-step procedure. Protein sequence is first transformed into a numerical series by means of the individual EIIP amino acid values. The second step of the model involves the Fourier spectral analysis of the obtained numerical series. Ćosić et al.1-8 have shown that single frequency peaks of the spectrum define characteristic positions of the amino acids, i.e., hot spots, correlated to the biological function of the protein. We have analysed the secondary protein structure by comparing the patterns of 20 most prominent frequency peaks of the single-series Fourier RRM periodogram. The patterns within 140 nonhomologous a- and p-protein folds obtained from the Jpred and SCOP databases were analysed by means of the classification tree in order to obtain the algorithm for the α- and β-fold classification. This quick and simple procedure of the secondary fold prediction showed high accuracy of 98.55%. The stability of the tree algorithm solution was confirmed by jack-knife testing of the tree algorithm (mean error 2.6). This method of the secondary structure prediction is presented in more detail on a subset of 30 different cytokines, hormones, enzymes and viral proteins. Our results indicate that resonant spectral analysis of the protein primary amino acid sequence may be used to extract information about its secondary structure

Understanding Hydrogen-Bond Patterns in Proteins using a Novel Statistical Model

Proteins are built from basic structural elements and their systematic characterization is of interest. Searching for recurring patterns in protein contact maps, we found several network motifs, patterns that occur more frequently in experimentally determined protein contact maps than in randomized contact maps with the same properties. Some of these network motifs correspond to sub-structures of alpha helices, including topologies not previously recognized in this context. Other motifs characterize beta-sheets, again some of which appear to be novel. This topological characterization of patterns serves as a tool to characterize proteins, and to reveal a high detailed differences map for comparing protein structures solved by X-ray crystallography, NMR and molecular dynamics (MD) simulations. Both NMR and MD show small but consistent differences from the crystal structures of the same proteins, possibly due to the pair-wise energy functions used. Network motifs analysis can serve as a base for many-body energy statistical energy potential, and suggests a dictionary of basic elements of which protein secondary structure is made

A solvable model of the genesis of amino-acid sequences via coupled dynamics of folding and slow genetic variation

We study the coupled dynamics of primary and secondary structure formation (i.e. slow genetic sequence selection and fast folding) in the context of a solvable microscopic model that includes both short-range steric forces and and long-range polarity-driven forces. Our solution is based on the diagonalization of replicated transfer matrices, and leads in the thermodynamic limit to explicit predictions regarding phase transitions and phase diagrams at genetic equilibrium. The predicted phenomenology allows for natural physical interpretations, and finds satisfactory support in numerical simulations.Comment: 51 pages, 13 figures, submitted to J. Phys.

A topological approach for protein classification

Protein function and dynamics are closely related to its sequence and structure. However prediction of protein function and dynamics from its sequence and structure is still a fundamental challenge in molecular biology. Protein classification, which is typically done through measuring the similarity be- tween proteins based on protein sequence or physical information, serves as a crucial step toward the understanding of protein function and dynamics. Persistent homology is a new branch of algebraic topology that has found its success in the topological data analysis in a variety of disciplines, including molecular biology. The present work explores the potential of using persistent homology as an indepen- dent tool for protein classification. To this end, we propose a molecular topological fingerprint based support vector machine (MTF-SVM) classifier. Specifically, we construct machine learning feature vectors solely from protein topological fingerprints, which are topological invariants generated during the filtration process. To validate the present MTF-SVM approach, we consider four types of problems. First, we study protein-drug binding by using the M2 channel protein of influenza A virus. We achieve 96% accuracy in discriminating drug bound and unbound M2 channels. Additionally, we examine the use of MTF-SVM for the classification of hemoglobin molecules in their relaxed and taut forms and obtain about 80% accuracy. The identification of all alpha, all beta, and alpha-beta protein domains is carried out in our next study using 900 proteins. We have found a 85% success in this identifica- tion. Finally, we apply the present technique to 55 classification tasks of protein superfamilies over 1357 samples. An average accuracy of 82% is attained. The present study establishes computational topology as an independent and effective alternative for protein classification

Statistical Mechanics and Kinetics of Amyloid Fibrillation

Amyloid fibrillation is a protein self-assembly phenomenon that is intimately related to well-known human neurodegenerative diseases. During the past few decades, striking advances have been achieved in our understanding of the physical origin of this phenomenon and they constitute the contents of this review. Starting from a minimal model of amyloid fibrils, we explore systematically the equilibrium and kinetic aspects of amyloid fibrillation in both dilute and semi-dilute limits. We then incorporate further molecular mechanisms into the analyses. We also discuss the mathematical foundation of kinetic modeling based on chemical mass-action equations, the quantitative linkage with experimental measurements, as well as the procedure to perform global fitting.Comment: 68 pages, 18 figures, 201 reference

Feedback control architecture & the bacterial chemotaxis network

Bacteria move towards favourable and away from toxic environments by changing their swimming pattern. This response is regulated by the chemotaxis signalling pathway, which has an important feature: it uses feedback to ‘reset’ (adapt) the bacterial sensing ability, which allows the bacteria to sense a range of background environmental changes. The role of this feedback has been studied extensively in the simple chemotaxis pathway of Escherichia coli. However it has been recently found that the majority of bacteria have multiple chemotaxis homologues of the E. coli proteins, resulting in more complex pathways. In this paper we investigate the configuration and role of feedback in Rhodobacter sphaeroides, a bacterium containing multiple homologues of the chemotaxis proteins found in E. coli. Multiple proteins could produce different possible feedback configurations, each having different chemotactic performance qualities and levels of robustness to variations and uncertainties in biological parameters and to intracellular noise. We develop four models corresponding to different feedback configurations. Using a series of carefully designed experiments we discriminate between these models and invalidate three of them. When these models are examined in terms of robustness to noise and parametric uncertainties, we find that the non-invalidated model is superior to the others. Moreover, it has a ‘cascade control’ feedback architecture which is used extensively in engineering to improve system performance, including robustness. Given that the majority of bacteria are known to have multiple chemotaxis pathways, in this paper we show that some feedback architectures allow them to have better performance than others. In particular, cascade control may be an important feature in achieving robust functionality in more complex signalling pathways and in improving their performance

Structural alphabets derived from attractors in conformational space

Background: The hierarchical and partially redundant nature of protein structures justifies the definition of frequently occurring conformations of short fragments as 'states'. Collections of selected representatives for these states define Structural Alphabets, describing the most typical local conformations within protein structures. These alphabets form a bridge between the string-oriented methods of sequence analysis and the coordinate-oriented methods of protein structure analysis.Results: A Structural Alphabet has been derived by clustering all four-residue fragments of a high-resolution subset of the protein data bank and extracting the high-density states as representative conformational states. Each fragment is uniquely defined by a set of three independent angles corresponding to its degrees of freedom, capturing in simple and intuitive terms the properties of the conformational space. The fragments of the Structural Alphabet are equivalent to the conformational attractors and therefore yield a most informative encoding of proteins. Proteins can be reconstructed within the experimental uncertainty in structure determination and ensembles of structures can be encoded with accuracy and robustness.Conclusions: The density-based Structural Alphabet provides a novel tool to describe local conformations and it is specifically suitable for application in studies of protein dynamics. © 2010 Pandini et al; licensee BioMed Central Ltd