Structural coordinates: A novel approach to predict protein backbone conformation.

Motivation: Local protein structure is usually described via classifying each peptide to a unique class from a set of pre-defined structures. These classifications may differ in the number of structural classes, the length of peptides, or class attribution criteria. Most methods that predict the local structure of a protein from its sequence first rely on some classification and only then proceed to the 3D conformation assessment.

The ranging of amino acids substitution matrices of various types in accordance with the alignment accuracy criterion.

Background: The alignment of character sequences is important in bioinformatics. The quality of this procedure is determined by the substitution matrix and parameters of the insertion-deletion penalty function. These matrices are derived from sequence alignment and thus reflect the evolutionary process.

Numeric analysis of reversibility of classic movement equations and constructive criteria of estimating quality of molecular dynamic simulations.

The fundamental criteria of the quality of molecular dynamics (MD) simulation represent a pivotal challenge, especially in the case of MD simulations of large systems (in particular, proteins).This work presents a simple theoretical analysis of time reversibility in classical mechanics that has allowed us to formulate a number of constructive criteria for evaluating the quality of the trajectories, generated in MD simulations. The results of testing the criteria on the structures of eight small proteins are presented.

Role of structural water for prediction of cation binding sites in apoproteins.

Structures of many metal-binding proteins are often obtained without structural cations in their apoprotein forms. Missing cation coordinates are usually updated from structural templates constructed from many holoprotein structures. Such templates usually do not include structural water, the important contributor to the ion binding energy.

Omnipresence of the polyproline II helix in fibrous and globular proteins.

Left-handed helical conformation of a polypeptide chain (PPII) is the third type of the protein backbone structure. This conformation universally exists in fibrous, globular proteins, and biologically active peptides. It has unique physical and chemical properties determining a wide range of biological functions, from the protein folding to the tissue differentiation.

Alternatingly twisted β-hairpins and nonglycine residues in the disallowed II' region of the Ramachandran plot.

The structure of the SH3 domain of α-spectrin (PDB code 1SHG) features Asn47 in the II' area of the Ramachandran plot, which as a rule admits only glycine residues, and this phenomenon still awaits its explanation. Here, we undertook a computational study of this particular case by means of molecular dynamics and bioinformatics approaches. We found that the region of the SH3 domain in the vicinity of Asn47 remains relatively stable during denaturing molecular dynamics simulations of the entire domain and of its parts.

Structural attributes of nucleotide sequences in promoter regions of supercoiling-sensitive genes: how to relate microarray expression data with genomic sequences.

The level of supercoiling in the chromosome can affect gene expression. To clarify the basis of supercoiling sensitivity, we analyzed the structural features of nucleotide sequences in the vicinity of promoters for the genes with expression enhanced and decreased in response to loss of chromosomal supercoiling in Escherichia coli. Fourier analysis of promoter sequences for supercoiling-sensitive genes reveals the tendency in selection of sequences with helical periodicities close to 10nt for relaxation-induced genes and to 11nt for relaxation-repressed genes.

Comparative analysis of the quality of a global algorithm and a local algorithm for alignment of two sequences.

Background: Algorithms of sequence alignment are the key instruments for computer-assisted studies of biopolymers. Obviously, it is important to take into account the "quality" of the obtained alignments, i.e.

Coexistence of different base periodicities in prokaryotic genomes as related to DNA curvature, supercoiling, and transcription.

We analyzed the periodic patterns in E. coli promoters and compared the distributions of the corresponding patterns in promoters and in the complete genome to elucidate their function. Except the three-base periodicity, coincident with that in the coding regions and growing stronger in the region downstream from the transcriptions start (TS), all other salient periodicities are peaked upstream of TS.

The Intramolecular Impact to the Sequence Specificity of B→A Transition: Low Energy Conformational Variations in AA/TT and GG/CC Steps.

Abstract It is well known, that local B→A transformation in DNA is involved in several biological processes. In vitro B↔A transition is sequence-specific. The physical basis of this specificity is not known yet.

The intramolecular impact to the sequence specificity of B-->A transition: low energy conformational variations in AA/TT and GG/CC steps.

It is well known, that local B--> A transformation in DNA is involved in several biological processes. In vitro B<--> A transition is sequence-specific. The physical basis of this specificity is not known yet.

A novel model system for design of biomaterials based on recombinant analogs of spider silk proteins.

Spider dragline silk possesses impressive mechanical and biochemical properties. It is synthesized by a couple of major ampullate glands in spiders and comprises of two major structural proteins--spidroins 1 and 2. The relationship between structure and mechanical properties of spider silk is not well understood.

Reconstruction of genuine pair-wise sequence alignment.

In many applications, the algorithmically obtained alignment ideally should restore the "golden standard" (GS) alignment, which superimposes positions originating from the same position of the common ancestor of the compared sequences. The average similarity between the algorithmically obtained and GS alignments ("the quality") is an important characteristic of an alignment algorithm. We proposed to determine the quality of an algorithm, using sequences that were artificially generated in accordance with an appropriate evolution model; the approach was applied to the global version of the Smith-Waterman algorithm (SWA).

Comprehensive statistical analysis of residues interaction specificity at protein-protein interfaces.

We calculated interchain contacts on the atomic level for nonredundant set of 4602 protein-protein interfaces using an unbiased Voronoi-Delaune tessellation method, and made 20x20 residue contact matrixes both for homodimers and heterocomplexes. The area of contacts and the distance distribution for these contacts were calculated on both the residue and the atomic levels. We analyzed residue area distribution and showed the existence of two types of interresidue contacts: stochastic and specific.

Study and prediction of secondary structure for membrane proteins.

In this paper we present a novel approach to membrane protein secondary structure prediction based on the statistical stepwise discriminant analysis method. A new aspect of our approach is the possibility to derive physical-chemical properties that may affect the formation of membrane protein secondary structure. The certain physical-chemical properties of protein chains can be used to clarify the formation of the secondary structure types under consideration.

A tetrapeptide-based method for polyproline II-type secondary structure prediction.

We describe a new method for polyproline II-type (PPII) secondary structure prediction based on tetrapeptide conformation properties using data obtained from all globular proteins in the Protein Data Bank (PDB). This is the first method for PPII prediction with a relatively high level of accuracy (approximately 60%). Our method uses only frequencies of different conformations among oligopeptides without any additional parameters.

Analysis of forces that determine helix formation in alpha-proteins.

A model for prediction of alpha-helical regions in amino acid sequences has been tested on the mainly-alpha protein structure class. The modeling represents the construction of a continuous hypothetical alpha-helical conformation for the whole protein chain, and was performed using molecular mechanics tools. The positive prediction of alpha-helical and non-alpha-helical pentapeptide fragments of the proteins is 79%.

Molecular modelling of disease-causing single-nucleotide polymorphisms in collagen.

The purpose of the work was to investigate at the molecular structural and energy levels the consequence of amino acid substitutions in collagen that cause systemic diseases. The data have been systematized on defects in human collagen III, and the patterns of single-nucleotide polymorphisms collected. Then molecular mechanics calculations were performed for native and mutant collagen molecule fragments.

Use of molecular mechanics for secondary structure prediction. Is it possible to reveal alpha-helix?

A new approach to predicting protein standard conformations is suggested. The idea consists in modeling by molecular mechanics tools a continuous alpha-helical conformation for the whole protein. The profile of energy along the model alpha-helix reveals minima corresponding to real alpha-helical segments in the native protein.

Segmentation of long genomic sequences into domains with homogeneous composition with BASIO software.

Unlabelled: We present a software system BASIO that allows one to segment a sequence into regions with homogeneous nucleotide composition at a desired length scale. The system can work with arbitrary alphabet and therefore can be applied to various (e.g.

DNA segmentation through the Bayesian approach.

We present a new approach to DNA segmentation into compositionally homogeneous blocks. The Bayesian estimator, which is applicable for both short and long segments, is used to obtain the measure of homogeneity. An exact optimal segmentation is found via the dynamic programming technique.

Hierarchy of regions of amino acid sequence with respect to their role in the protein spatial structure.

The method of the representation of amino acid sequence by graph of the interactions energy between parts of spatial structure has been elaborated. Our method provides the possibility to establish the compatibility between each point of a polypeptide chain and the Van der Waals interactions energy of regions of a native globule adjacent to this amino acid residue. We have undertaken an exhaustive analysis of a set of proteins.

A new approach for the calculation of the energy of van der Waals interactions in macromolecules of globular proteins.

Van der Waals interaction energy in globular proteins is presented by the interaction energies between regions of protein spatial structure with homogenous medium density distribution. We introduce a notion of the local medium permittivity as a function of absorptance of molecular groups with particular conformation. Proposed theory avoids shortcomings which are typical for the calculations on the basis of the pairwise additive approximation.

Hierarchy of the interaction energy distribution in the spatial structure of globular proteins and the problem of domain definition.

An algorithm for determining of protein domain structure is proposed. Domain structures resulted from the algorithm application have been obtained and compared with available data. The method is based on entirely physical model of van der Waals interactions that reflects as illustrated in this work the distribution of electron density.