Search

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. The criteria can be useful for solving different MD problems, such as: choosing the most appropriate thermostats for a MD system under study, the methods for sampling conformations, etc.Communicated by Ramaswamy H. Sarma.

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. However, most classification methods rely on homologous proteins' existence, unavoidably lose information by attributing a peptide to a single class or suffer from a suboptimal choice of the representative classes., Results: To alleviate the above challenges, we propose a method that constructs a peptide's structural representation from the sequence, reflecting its similarity to several basic representative structures. For 5-mer peptides and 16 representative structures, we achieved the Q16 classification accuracy of 67.9%, which is higher than what is currently reported in the literature. Our prediction method does not utilize information about protein homologues but relies only on the amino acids' physicochemical properties and the resolved structures' statistics. We also show that the 3D coordinates of a peptide can be uniquely recovered from its structural coordinates, and show the required conditions under various geometric constraints., Competing Interests: The authors have declared that no competing interests exist.

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. Structural templates are also inconvenient for taking into account structural modifications around the binding site at apo-/holo- transitions. An approach based upon statistical potentials readily takes into account structural modifications associated with binding as well as contribution of structural water molecules. Here, we construct a set of statistical potentials for Mg 2+ , Ca 2+ , and Zn 2+ contacting with protein atoms of a different type or structural water oxygens. Each type of the cations tends to form tight contacts with protein atoms of specific types. Structural water contributes relatively more into the binding pseudo-energy of Mg 2+ and Ca 2+ than of Zn 2+ . We have developed PIONCA (Protein-Ion Calculator), a fast CUDA GPGPU-based algorithm that predicts ion-binding sites in apoproteins. Comparative tests demonstrate that PIONCA outperforms most of the tools based on structural templates or docking. Our software can be also used for locating bound cations in holoprotein structures with missing cation heteroatoms. PIONCA is equipped with an interactive web interface based upon JSmol.

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. New examples of the structure have been appearing in spite of difficulties in their detection and investigation. The annotation and prediction of the PPII was also a challenging task. Recently, many PPII motifs with new and/or unexpected functions are being accumulated in databases. In this review we describe the major structural and dynamic forms of PPII, the diversity of its functions, and the role in different biological processes., (Copyright © 2016 Elsevier Ltd. All rights reserved.)

The analysis of conformationally stable (conformational conservative) tetrapeptides selected from protein structures deposited in PDBSelect data bank has been fulfilled. The subset contained 943 tetrapeptide amino acid sequences and there were merely five 3D protein segment representatives for each sequence. As a result, the conclusion has been drawn on the basis of DSSP annotation analysis that in the majority of cases (900 of 943) alpha-helical conformation is obvious. Different than alpha-helix, in particular, the left-handed polyproline II helical conformation was observed in 43 sequences. The physical and chemical properties of conformationally stable peptides taken from the appropriate sample were estimated by the average hydrophobicity/hydrophilicity of tetrapeptides. The results of calculations show that the "neutrality" towards hydrophobicity/hydrophilicity is representative of conformationally stable oligopeptides. It should be noted, that dispersion of hydrophobicity/hydrophilicity distribution is sufficiently lower than for the test subsets. Thus, the conformationally stable oligopeptides present a distinct group of local protein structures which are very close with respect to conformational and physicochemical properties. In accordance with our developed theory of specific long range interactions these peptides are the objects being quite useful for effective mutual molecular recognition.

Conservative regions varying in length from 30 to 70 nucleotides (nt) are revealed as a result of interspecies comparison of homologous nucleotide sequences of the loci, genes controlling early development in Drosophila species. These regions are observed at the intermediate level between transcription factor binding sites (usually about 7-10 nt) and the nucleosome repeat unit (165-210 nt). Segments found are located mainly in the area of known functional elements of the locus (cis-regulatory modules: enhancers, a proximal promoter, a coding segment). Conservative domains occupy less than a half of the full enhancer length and contain practically all annotated transcription factor binding sites. A quasi periodical pattern of conservative region disposition is in accord with the experimental nucleosome localization pattern. The distance between neighboring domains is about 84 nt, i.e. a pitch of nucleosome DNA super helix. A repeatable accuracy of the distance permits of treating the delineated conservative regions as co-phased blocks and stating that these domains are close in the space being located at neighboring coils of the nucleosome DNA super helix.

The work presents the results of an exhaustive conformational analysis of β-turns involving amino acid residues with disallowed backbone conformation of the polypeptide chain. It is known that the first residue of the β-turn (Asn47) of the distal β-hairpin in the α-spectrin SH3-domain is characterized by sterically disallowed main chain conformation (values of the dihedral angles (φ and ψ are in the right bottom quadrant of the Ramachandran plot). All α-spectrin structures with the anomalous elements deposited in the PDB were analysed. We hypothesized that the formation of disallowed conformation may occur through the fixation (due to the SH3 domain structure) of the adjacent to the β-turn amino acid residues with the β-structure. These residues are disposed in such a manner that β-turn conformation of the residues contributes just to the disallowed local conformation of this residue whereas any other β-turn conformations (with allowed local conformation) are impossible. To test this hypothesis an exhaustive conformational analysis of the β-bend has been performed by altering internal coordinates (two pairs of φ and ψ angles and two Ω angles). The conformations were selected as a result of grid search procedure with. 1 degrees step that corresponded to stereochemically allowed local deformations of the polypeptide chain segment forming the β-turn. In all conformations obtained the local conformation of Asn47 rests in the disallowed region. The conformations found include conformations coinciding with experimentally determined structures from the PDB as well as an additional variant that differs from X-ray structure in values of a pair of φ and ψ angles of the second residue belonging to the β-bend. Values of these angles fall in the region of the Ramachandran plot near the line φ = 0 (and negative values of ψ) i.e. in strongly disallowed region without experimental points. Therefore the additional variants of the β-turn local deformation are impossible to observe in experiment. Thus, the idea that disallowed conformation is intruded to the β-bend by fixation of adjacent residues receives confirmation in this work. The topological limitations in a context of the structure in such kind of β-hairpins exclude the allowed local conformations.

We propose a method for magnetic sorting of cell suspensions able to differentiate not only cells with single specific antigen on their surface, but also cells with a group of pre-defined antigens. Individually, each antigen of this group may be present on surfaces of non-selected cells. However, only the simultaneous presence of all given antigens on the cell surface means that such a cell should be separated. The method is of interest, for purification of cell suspensions from malignant cells, in particular, for purification of bone marrow material for autologous transplantation in leukemia.

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. This increased stability may be connected with the dynamic hydrogen bonding that is susceptible to targeted in silico mutations of Arg49. Bioinformatics analysis indicated that Asn47 is in the β-turn of a distinctive structural fragment we called 'alternatingly twisted β-hairpin.' Fragments of similar conformation are quite abundant in a nonredundant set of PDB chains and are distinguished from ordinary β-hairpins by some surplus of glycine in their β-turns, lack of certain interpeptide hydrogen bonds, and an increased chirality index. Thus, the disallowed conformation of residues other than glycine is realized in the β-turns of alternatingly twisted β-hairpins.

The classification of amino acid residues based on the events of contact formation between distinct amino acid and selected nucleotides was constructed. Thus, the most integral properties, that characterize interactions in organization of DNA-protein complexes, were used. We applied the Voronoi-Delaunay tessellation to draw statistics of contacts and area of contacts for the set included 1937 DNA-protein complexes. Similarities of amino acid residues have been searched for based on the comparison of corresponded rows and matrixes of contacts and areas of contacts. Nine measures of distance were used for estimation of rows similarity degree. The procedure of clustering amino acids in groups included three hierarchical and two nonhierarchical methods. A total tree was built using nine techniques of estimating distance with three hierarchical clustering methods. It was shown that clustering centers in the main groups are always constant while other relationships between objects vary. Clustering of binary associations was found for the most amino acids. Major classes of up to six amino acids correspond to the certain local structures of the polypeptide chain in the context of amino acid composition. These data should be taken into account when designing DNA-protein ligands.

In the work the arguments are presented in favor of the idea on the role of conformationally stable oligo-peptides in specific long-distance interactions in phenomena of molecular recognition during various biological processes. Original authors' and literature data are taken into account. The examples of conformationally stable short oligopeptides participation in alpha-helix and collagen type structures formation are given simultaneously with theoretical approaches. The conformationally stable oligopeptides obtained in the course of PDB bank analysis are discussed. The role of amino acids sequence in collagen helix formation is shown.

The correlation was observed between DNA segments characterized by periodic patterns of nucleotides and nucleosome DNA binding sites predicted by the profiles of local DNA flexibility. Exons of human type I and VII collagen genes loci have been chosen as a sample set of periodic DNA segments. It was pointed that short periodic exons coding fibrillar part of these proteins are localized near borders of the DNA sequence that binds the histone core. One of the borders of exon and of this sequence coincide as a rule. In non-fibrillar parts of collagen genes exons (longer on average) dispose, in contrast, in any place of the DNA sequence bound to the nucleosome core.

It was shown that selective interactions between helical segments of macromolecules can realize in globular proteins in the segments characterized by the same periodicities of charge distribution i.e. between conformationally conservative oligopeptides. It was found that in the macromolecules of alpha-helical proteins conformationally conservative oligopeptides are disposed at a distance being characteristic of direct interactions. For representatives of many structural families of alpha-type proteins specific disposition of conformationally conservative segments is observed. This disposition is inherent to a particular structural family. Disposition of conformationally conservative segments is not related to homology of the amino acid sequence but reflects peculiarities of native 3D-architectures of protein globules.

On the basis of the theory of specific long-range interaction between long molecules the approach has been elaborated for the "fast reading" of nucleotide sequences in sole DNA molecule. Firstly, a stretching force applies to the molecule that causes its unwinding from B-form to so called S-form. Then, a molecule disposes in a stretched state on the background. After this, the electrostatic potential is estimated in the space along the DNA filament. Information obtained is sufficient for the deduction of the nucleotide sequence in DNA. Another approach to the "reading" of information reduces to the measurement of the filament element deformation induced by electrostatic field from the electrode that stretches the filament by alternating current applied.

Motivation: Modern experimental methods provide substantial information on protein-DNA recognition. Studying arrangements of transcription factor binding sites (TFBSs) of interacting transcription factors (TFs) advances understanding of the transcription regulatory code., Results: We constructed binding motifs for TFs forming a complex with HIF-1α at the erythropoietin 3(')-enhancer. Corresponding TFBSs were predicted in the segments around transcription start sites (TSSs) of all human genes. Using the genome-wide set of regulatory regions, we observed several strongly preferred distances between hypoxia-responsive element (HRE) and binding sites of a particular cofactor protein. The set of preferred distances was called as a preferred pair distance template (PPDT). PPDT dramatically depended on the TF and orientation of its binding sites relative to HRE. PPDT evaluated from the genome-wide set of regulatory sequences was used to detect significant PPDT-consistent binding site pairs in regulatory regions of hypoxia-responsive genes. We believe PPDT can help to reveal the layout of eukaryotic regulatory segments., Contact: ivan.kulakovskiy@gmail.com, Supplementary Information: Supplementary data are available at Bioinformatics online.

The process of globular structure formation from a long molecular chain has been examined. In the course of this process, various regions of the chain interact with one another. We classify the bonds formed during this process as "correct" and "erroneous" ones. The term "correct" bonds implies the bonds characteristic for a completely formed native globular structure. All other bonds can be treated as "erroneous". It was demonstrated that the process of globule formation may proceed actually without the formation and the following decay of "erroneous" bonds. Our model permits one to avoid the examination of numerous "erroneous" variants since, between the regions of the chain that form "correct" bonds, long-distance interactions characterized simultaneously by high selectivity take place. The existence of interactions of this kind facilitates the drawing together and subsequent interaction of just these regions of the chain that yield "correct" bonds. Based on the data bank analysis, it was demonstrated that the model elaborated is valid not only for abstract structures but also for real polypeptide chains capable of forming protein globules and superhelical fibrils.

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. Here we analyze the effect of intramolecular interactions on the structural behavior of the GG/CC and AA/TT steps. These steps exemplify sequence specific bias to the B- or A-form structure. Optimization of potential energy of the molecular systems composed of an octanucle-otide, neutralized by Na(+) and solvated with TIP3P water molecules in rectangular box with periodic boundary conditions gives the statistically representative sets of low energy structures for GG/CC and AA/TT steps in the middle of the diverse flanking sequences. Permissible 3D variations of GG/CC and AA/TT, and correlation of the relative motion of base pairs in these steps were analyzed. AA/TT step permits high variability for low energy conformers in the B-form DNA and small variability for low energy conformers in the A-form DNA. In contrast GG/CC step permits high variability for low energy conformers in the A-form DNA and small variability for low energy conformers in the B-form DNA. The relative motion of base pairs in GG/CC step is high correlated, while in AA/TT step this correlation is notably less. Atom-atom interactions inside-the-step always favors the B-form and their component - stacking interactions (atomatom interactions between nucleic bases) is crucial for the duplex stabilization. Formation of the A-form for both steps is a result of interactions with the flanking sequences and water-cation environment in the box. The average energy difference between conformations presenting B-form and A-form for the GG/CC step is high, while for the AA/TT step it is rather low. Thus, intramolecular interactions in GG/CC and AA/TT steps affect the possible conformational diversity ("conformational entropy") of the A- and B- type structures of DNA step. This determines the known bias of the A-form DNA depending on the enrichment of sequences with GG/CC. If structural tuning during the process of protein-DNA complex formation lead to the local B→A transformation of DNA, it is largely directed by high conformational diversity of GG/CC step in the A-form. In such a case the presence in the target site of both kinds of examined steps ensures the reversible character of ligand binding.

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. Here we analyze the effect of intramolecular interactions on the structural behavior of the GG/CC and AA/TT steps. These steps exemplify sequence specific bias to the B- or A-form structure. Optimization of potential energy of the molecular systems composed of an octanucleotide, neutralized by Na(+) and solvated with TIP3P water molecules in rectangular box with periodic boundary conditions gives the statistically representative sets of low energy structures for GG/CC and AA/TT steps in the middle of the diverse flanking sequences. Permissible 3D variations of GG/CC and AA/TT, and correlation of the relative motion of base pairs in these steps were analyzed. AA/TT step permits high variability for low energy conformers in the B-form DNA and small variability for low energy conformers in the A-form DNA. In contrast GG/CC step permits high variability for low energy conformers in the A-form DNA and small variability for low energy conformers in the B-form DNA. The relative motion of base pairs in GG/CC step is high correlated, while in AA/TT step this correlation is notably less. Atom-atom interactions inside-the-step always favors the B-form and their component - stacking interactions (atom-atom interactions between nucleic bases) is crucial for the duplex stabilization. Formation of the A-form for both steps is a result of interactions with the flanking sequences and water-cation environment in the box. The average energy difference between conformations presenting B-form and A-form for the GG/CC step is high, while for the AA/TT step it is rather low. Thus, intramolecular interactions in GG/CC and AA/TT steps affect the possible conformational diversity ("conformational entropy") of the A- and B- type structures of DNA step. This determines the known bias of the A-form DNA depending on the enrichment of sequences with GG/CC. If structural tuning during the process of protein-DNA complex formation lead to the local B--> A transformation of DNA, it is largely directed by high conformational diversity of GG/CC step in the A-form. In such a case the presence in the target site of both kinds of examined steps ensures the reversible character of ligand binding.

Fourier spectra of nucleotide sequences from chromosome 15 of Gallus gallus have been obtained. The periodical structures of several Ig locus fragments have been characterized. The uniqueness of Fourier spectra of the locus of the Ig light chain and its fragments (such as a variable region, a constant region, an enhancer, and several other fragments) among Fourier spectra of other fragments of chromosome 15 is demonstrated.

The distribution of secondary structures along the polypeptide chains of spider proteins spidroins 1 and 2 and their recombinant analogs has been studied by statistical methods. It was found that these proteins in the monomolecular form contain only trace amounts of beta-structures. At the same time, the regions of the sequence including Ala and Gly are predicted as helical-containing (with alpha-helices and left-helices of polyproline II type). An analysis of literature data and our data obtained in this study shows that the main conformation of the polypeptide chain solutions of spidroins 1 and 2 and their recombinant analogs in water solutions is the left-helix of polyproline II type with some contaminations of alpha-helices and a very small share of beta-structures. The transition to the state with extended conformations, which are peculiar to mature filaments of spider webs, requires the dehydration of the polypeptide chain backbone. Thus, the genesis of beta-structure in spider web proteins is determined by the conditions of conformation transitions between the main regular conformations of the polypeptide chain backbone.

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. Here, we modeled the complete process of the spider silk assembly using two new recombinant analogs of spidroins 1 and 2. The artificial genes sequence of the hydrophobic core regions of spidroin 1 and 2 have been designed using computer analysis of existing databases and mathematical modeling. Both proteins were expressed in Pichia pastoris and purified using a cation exchange chromatography. Despite the absence of hydrophilic N- and C-termini, both purified proteins spontaneously formed the nanofibrils and round micelles of about 1 microm in aqueous solutions. The electron microscopy study has revealed the helical structure of a nanofibril with a repeating motif of 40 nm. Using the electrospinning, the thin films with an antiparallel beta-sheet structure were produced. In summary, we were able to obtain artificial structures with characteristics that are perspective for further biomedical applications, such as producing three-dimensional matrices for tissue engineering and drug delivery.

The influence of imino acids on the thermodynamic characteristics of collagen type structures in various collagens has been analyzed. It was shown that the basic mechanism of entropy increase in the protein-water system consists in the alteration in the number of cooperative segments accompanying the increase in imino acids content, which can be observed during the melting of the fibrous macromolecule. The range of variation in the physical characteristics of cooperative units is determined, in particular, by the variability of hydrogen bond parameters. This is displayed in a broadening of the bands of NH-valence vibrations and the half-widths of transitions in post-denatured structures. Thus, the basic mechanism of the influence of imino acids on thermodynamic characteristics of collagens is related to the complex nature of the melting process. The dehydration-hydration mechanisms of native and denatured states become significantly different upon replacement of any amino acid by an imino acid.

A new method for the rigorous analysis of DNA-protein contacts has been developed on the basis of Voronoi tessellation. This method permits one to determine close neighbors on the atomic level and compute the area of contact by edges of Voronoi polyhedra. Based on the results of the study of 1109 protein-DNA complexes from PDB, it was demonstrated that about one third of the contacts are the contacts with positively charged Arg and Lys. There is distinct amino acid prevalence for nucleotides: for A - Pro; for T - His; for G - Asp; for C - Trp, Asp, Glu. Therefore, GC pairs prefer to interact with negatively charged residues, and alanine and methionine, whereas AT pairs prefer to interact with histidine and unpolar residues.

Regions with a quasiperiodical location of exon--intron sites have been found in the loci of genes of I and VII type collagens (with a total length of exons more than 15% of the entire size of the locus). The periods observed are similar to periods typical for the nucleosomal level of the organization of chromatin. It was shown that the sites consisting of successively arranged exons and introns form groups involving two to five such regions of the same length. The groups encoding the fibrillar regions of the gene product contain more than 50% of exons. The regions are on the average 165 nt long, which is close to the minimal nucleosomal repeat length observed in some regions of the eukaryotic genome. In the nonfibrillar region of the gene of VII type collagen, groups of several exon-intron pairs with an average length of 227 nt were identified. The change in the length of exon-introns sites on going from the nonfibrillar to the fibrillar moiety occurs in a jump, which is clearly seen on a periodogram of the locus.

Promoters recognized by deltaA-RNA and deltaH-RNA polymerase have different periodic patterns their nucleotide disposition. By using the special method of Fourier analysis for symbolic sequences. Fourier spectra of the primary structure of promoters transcribed bv deltaA-RNA and deltaH-RNA polymerase were obtained. For two data sets, a small one and a big one, a stepwise discriminant analysis with jackknife test based on the spectral characteristics of each of the promoters was performed. It was shown that it is possible to classify the data with the accuracy of 100% sets into classes: promoters that are recognized by the deltaA-RNA polymerise and promoters that are recognized by the deltaH-RNA polymerase. Significant correlations between the promoter strength and the characteristics of their Fourier spectra were obtained. Thus, the periodicity in nucleotide distribution along the polynucleotide chain is the attribute sufficient for promoter recognition by RNA polymerase holoenzyme.

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. Another aspect of our approach is that the results of multiple sequence alignment, or the other kinds of sequence alignment, are not used in the frame of the method. Using our approach, we predicted the formation of three main secondary structure types (alpha-helix, beta-structure and coil) with high accuracy, that is Q(3) = 76%. Predicting the formation of alpha-helix and non-alpha-helix states we reached the accuracy which was measured as Q(2) = 86%. Also we have identified certain protein chain properties that affect the formation of membrane protein secondary structure. These protein properties include hydrophobic properties of amino acid residues, presence of Gly, Ala and Val amino acids, and the location of protein chain end.

A new method for determining bound water was developed, which is based on precise measurements of the enthalpy of water evaporation from a sample using differential scanning calorimetry.

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. We also attempted to predict alpha-helices and beta-strands using the same approach. We find that the application of our method reveals interrelation between sequence and structure even for very short oligopeptides (tetrapeptides)., (Proteins 2005. 2005 Wiley-Liss, Inc.)

The formation of beta-turns in globular proteins has been studied by the method of molecular mechanics. Statistical method of discriminant analysis was applied to calculate energy components and sequences of oligopeptide segments, and after this prediction of I type beta-turns has been drawn. The accuracy of true positive prediction is 65%. Components of conformational energy considerably affecting beta-turn formation were delineated. There are torsional energy, energy of hydrogen bonds, and van der Waals energy.

Step-by-step computations of octapeptide structures was performed for human collagens I and III. It was shown that computational results (coordinates of atoms) practically coincide with X-ray data for the collagen fragment.

Sequences of spidroin 1 and 2 from spiders belonging to various species were analyzed by Fourier analysis. Specific periodical patterns were found in various segments of the sequences. These characteristic periodicities vary within the same proteins as well as between spidroin 1 and spidroin 2 sequences. It is possible that alterations in periodicity help to recognize contact sites between the molecules. Spidroins of 2 type have similar sequence structure consisting of four parts with a particular periodical pattern. These parts are a constant C-terminal part, a long-periodical part, a short-periodical part, and a constant N-terminal part.

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%. The model considers only local interactions in the polypeptide chain without the influence of the tertiary structure. It was shown that the local interaction defines the alpha-helical conformation for 85% of the native alpha-helical regions. The relative energy contributions to the energy of the model were analyzed with the finding that the van der Waals component determines the formation of alpha-helices. Hydrogen bonds remain at constant energy independently whether alpha-helix or non-alpha-helix occurs in the native protein, and do not determine the location of helical regions. In contrast to existing methods, this approach additionally permits the prediction of conformations of side chains. The model suggests the correct values for ~60% of all chi-angles of alpha-helical residues.

The quantitative characteristics of a spectral portrait of the nucleotide sequence of E. coli oriC were determined by the method of matrix Fourier analysis.

A new approach to the analysis of regular structures in proteins that is based on the method of molecular mechanics is proposed. The method uses only the information about the amino acid sequence. The alpha-helical conformation was simulated using the ICM program of molecular mechanics. Energy profiles of the sequences in the alpha-helical conformation, spanning the entire polypeptide chain, were plotted for eight proteins from the Protein Data Bank. The regions of each profile that exhibit energy minima were found to correspond to the alpha-helical regions of the real spatial structure of the protein. Twenty-four out of 25 helices were distinctly pronounced, which indicates a rather high accuracy of the prediction. The energy profiles also help reveal the short regions that correspond to 3/10-helices and the turns that include local alpha-helical conformations. Unlike the known statistical methods of prediction, this method makes it possible to establish the physical principles of the formation of alpha-helical conformations. The English version of the paper: Russian Journal of Bioorganic Chemistry, 2002, vol. 28, no. 6; see also http://www.maik.ru.

The regions of initiation of replication of some bacterial genomes were studied by the method of Fourier matrix analysis. A generalized spectral portrait of the primary structures of E. coli-like regions of initiation of replication in bacteria was obtained, which reflects the features of their structural and functional organization. It contains well-pronounced peaks that correspond to the periods T = 2, 11, 17, 27, 86-105 of nucleotides. The peaks corresponding to T = 9, 13, 14, 18, 19, 33-35, 45-47, 74-85, 106-110 are less pronounced. The uniqueness of the Fourier spectrum corresponding to the region of initiation of replication of E. coli oriC was considered by the example of the complete genome of E. coli. Some regions of the E. coli genome were identified that differ from oriC in the primary structure but have Fourier spectra resembling the spectrum of oriC. A number of these regions are alternative points of initiation of replication in sdrA(rnh) mutants of E. coli, the others are localized in yet unidentified regions of the E. coli genome but are capable, in our opinion, to participate in the initiation of replication. Thus, from the similarity of spectral portraits of different regions of the genome, it was possible to reveal several regions that have similar functions, i.e., are involved in initiation of replication.

Sequence alignment shows that residue Arg 284 (according to the numbering of the residues in formate dehydrogenase, FDH, from the methylotrophic bacterium Pseudomonas sp. 101) is conserved in NAD-dependent FDHs and D-specific 2-hydroxyacid dehydrogenases. Mutation of Arg 284 to glutamine and alanine results in a change of the catalytic, thermodynamic and spectral properties of FDH. In comparison to wild-type, the affinity of the mutants for the substrate (K(formate)m) or the transition state analogue (K(azide)i) decreases and correlates with the ability of the side chain of residue 284 to form H-bonds. In contrast, the affinity for the coenzyme (K(NAD)d or K(NAD)m) is either not affected or increases and correlates inversely with the partial positive charge of the side chain. The temperature dependence of circular dichroism (CD) spectra of the wild-type FDH and its Ala mutant has been studied over the 5-90 degrees C temperature range. Both proteins reveal regions of enhanced conformational mobility at the predenaturing temperatures (40-55 degrees C) associated with a change of enzyme kinetic parameters and a co-operative transition around 55-70 degrees C which is followed by the loss of enzyme activity. CD spectra of the wild-type and mutant proteins were deconvoluted and contributions from various types of secondary structure estimated. It is shown that the co-operative transition at 55-70 degrees C in the FDH protein globule is triggered by a loss of alpha-helical secondary structure. The results confirm the conclusion, from the crystal structures, that Arg 284 is directly involved in substrate binding. In addition this residue seems to exert a major structural role by supporting the catalytic conformation of the enzyme active centre.