Search

An analog of the Ramachandran map was drawn, a new representation proposed, and thorough analysis performed using modern recognition and classification methods. Very large maps with a density of more than 50 million dots were created based on the data sets derived from the latest releases of globular protein- structure data banks. A, B, B', C, and D regions that correspond to strongly disallowed conformations were defined and found to occupy 25% of the plot area. A region of noncanonical conformations was determined by subtracting strongly disallowed and permitted conformation regions from the total plot area. Arguments are provided to support the new classification of backbone conformations of the protein polypeptide chain.

The interactions between protein and DNA in essence underlie all processes in a living cell. Understanding the principles of specific recognition of DNA sites will open the way to understand how these processes are controlled and to interfere in their operation. In the paper we studied contacts between the protein and DNA at the atomic level in the structures of all the 3 518 protein­DNA complexes available in PDB by the Voronoi–Delaunay tessellation method. The method unambiguously defines contacts between atoms without any parameters, and characterizes each contact by the distance between atoms and the contact area, which is determined by the corresponding face of the Voronoi polyhedron. It was shown that most contacts are formed between the protein atoms and the sugarphosphate backbone of the DNA (72.9 %). The contact with the atoms of the nucleic bases emerging into the grooves of DNA is 17.0 % for a major groove and 10.1 % for all atomic contacts for a minor groove. Totally, the interaction between protein atoms and nucleic base atoms accounts for 27.1 % of all contacts. Analysis of the accessible surface area of atoms in the major and the minor grooves showed a correlation with the number of contacts (coefficient of linear correlation 0.94 and 0.93, respectively), however, nucleic acid atoms forming hydrogen bonds make contacts more often than may be expected from statistical considerations. It was shown that conformationally stable peptides occur sometimes in the binding regions with DNA. Analysis of the residues in a predefined conformation in 3 518 protein­DNA complexes revealed 159 amino acid residues in a predefined β­bend type I conformation, 15 residues in the conformation of β­bend type I’, and 6 residues in the conformation of β­bend type II. No residues in the conformation of β­bend type II’ were found. Analysis of contacts showed that such residues virtually do not form contacts with DNA. Contacts with nucleic base atoms are found only in the two homologous structures 3qea and 3qe9, where threonine atoms form contacts with atoms of nucleotide bases of the AT­pair.

A detailed analysis of polypeptide-chain backbone conformations was carried out for polypeptide-chain segments adjacent to β-turn regions, including the sites of disallowed conformations. A cross comparison of conformations was performed for disallowed regions of the Ramachandran plot and main types of β-turns and adjacent secondary structures. Based on the results, disallowed region 2 (II, II') in the Ramachandran plot was shown to coincide mainly with β-hairpins and, more exactly, twisted β-hairpins. The frequency of residues with angles ϕi, ψi that fall in region 2 (II, II') in the latter is 140 times higher than in common β-hairpins.

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.

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.

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 Mg2+, Ca2+, and Zn2+ 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 Mg2+ and Ca2+ than of Zn2+. 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.

On the basis of the theory of specific long-range interaction between long molecules, an approach has been elaborated for «fast reading» of nucleotide sequences in one DNA molecule. First, a stretching force is applied to the molecule that causes its unwinding from B-form to S-form. Then, the molecule is placed in the stretched state on a support. After this, the electrostatic potential is estimated in a space along the DNA filament. The information obtained is sufficient for deducing the nucleotide sequence. Another approach to the «reading» of information reduces to measurement of the deformation of filament elements induced by the electric field from the electrode that stretches the filament by an alternating current applied.

The distribution of secondary structure elements along the polypeptide chains of spider silk proteins spidroins 1 and 2 and their recombinant analogs has been studied by statistical methods. It was found that these proteins as monomers contain only traces of β-structure, while the Ala-rich and the Gly-rich regions are predicted as α-helices and as left-handed helices of polyproline II type. Analysis of literature and our CD data shows that the major polypeptide chain conformation of spidroins 1 and 2 and their recombinant analogs in aqueous solutions is the polyproline II helix, with some α-helices and a very small share of β-structures. The transition to the state with extended conformations, which are characteristic of mature silk fibers, requires dehydration of the polypeptide backbone. Thus, the genesis of β-structure in spider web proteins is determined by the conditions of transitions between the main regular backbone conformations.

The structures of D-period segments of collagen (234 amino residues or ∼1/4 of whole length) are established by methods of molecular mechanics and geometry analysis. Each D-period segment proves to have a unique spatial structure. The distributions of local helical parameters along the molecule are calculated. It is found that a second hydrogen bond is formed in every case when the second residue in the tripeptide G-X-Y is an amino acid. With such a combined H-bond network, all the peptide CO groups of glycines and of the third residues in tripeptides have quasi-equivalent positions on the surface of the collagen molecule. The local deformations of the polyproline II helix in the triple complex give rise to the observed modulation of structure at the macromolecular level, which may be important for the mutual orientation of collagen molecules during fibrillogenesis.

In the previous paper in the series a full set of possible structures of double-helical DNA was obtained on the basis of complete scanning throughout conformational space of dihedral angles. In this work, these conformations have been used as starting points for energy optimization. All structures with Watson-Crick pairing coinciding with local minima of energy are deduced. These are right helices, B (a classical Wilkins B form), B1, B2 (the Watson and Crick form), B3, left helices, L1 (the Sasisekharan form), L2 (the Tumanyan et al. form), L3, and some helices without stacking of nucleotide pairs, K1-K6. Conformations B2, L2, and L1 are close to the B form with respect to energy (L3, B1, and B3 have rather poor energy). Incidentally, the similarity expands to values of helical parameters. For the structures described, the various energetic contributions in total energy have been estimated. The types of sugar conformations are analyzed. The possible biological significance of the elaborated optimal structures is discussed.

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).

We suggest a new simple approach for comparing the primary structure of proteins and their spatial structure. It relies on the one-to-one correspondence between each residue of the polypeptide chain and the energy of van der Waals interactions between the regions of the native globule flanking this residue. The method obviates the sophisticated geometrical criteria for estimating similarity between spatial structures. Besides, it permits one to analyze structural units of different scale.

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.

Monte Carlo simulations [(N, V, T)-ensemble] were performed for the hydration shell of poly(dA-dT).poly(dA-dT) in canonical B form and for the hydration shell of poly(dA).poly(dT) in canonical B conformation and in a conformation with narrow minor groove, highly inclined bases, but with a nearly zero-inclined base pair plane (B' conformation). We introduced helical periodic boundary conditions with a rather small unit cell and a limited number of water molecules to reduce the dimensionality of the configuration space. The coordinates of local maxima of water density and the properties of one- and two-membered water bridges between polar groups of the DNA were obtained. The AT-alternating duplex hydration mirrors the dyad symmetry of polar group distribution. At the dApdT step, a water bridge between the two carbonyl oxygens O2 of thymines is formed as in the central base-pair step of Dickerson's dodecamer. In the major groove, 5-membered water chains along the tetranucleotide pattern d(TATA).d(TATA) are observed. The hydration geometry of poly(dA).poly(dT) in canonical B conformation is distinguished by autonomous primary hydration of the base-pair edges in both grooves. When this polymer adopts a conformation with highly inclined bases and narrow minor groove, the water density distribution in the minor groove is in excellent agreement with Dickerson's spine model. One local maximum per base pair of the first layer is located near the dyad axis between adjacent base pairs, and one local maximum per base pair in the second shell lies near the dyad axis of the base pair itself. The water bridge between the two strands formed within the first layer was observed with high probability. But the water molecules of the second layer do not have a statistically favored orientation necessary for bridging first layer waters. In the major groove, the hydration geometry of the (A.T) base-pair edge resembles the main features of the AT-pair hydration derived from other sequences for the canonical B form. The preference of the B' conformation for oligo(dA).oligo(dT) tracts may express the tendency to common hydration of base-pair edges of successive base pairs in the grooves of B-type DNA. The mean potential energy of hydration of canonical B-DNA was estimated to be -60 to -80 kJ/mole nucleotides in dependence on the (G.C) contents. Because of the small system size, this estimation is preliminary.

We calculated interchain contacts on the atomic level for nonredundant set of 4602 pro- tein-protein interfaces using an unbiased Voronoi- Delaune tessellation method, and made 20320 resi- due contact matrixes both for homodimers and heterocomplexes. The area of contacts and the dis- tance distribution for these contacts were calculated on both the residue and the atomic levels. We ana- lyzed residue area distribution and showed the exis- tence of two types of interresidue contacts: stochas- tic and specific. We also derived formulas describing the distribution of contact area for stochastic and specific interactions in parametric form. Maximum pairing preference index was found for Cys-Cys con- tacts and for oppositely charged interactions. A sig- nificant difference in residue contacts was observed between homodimers and heterocomplexes. Interfa- ces in homodimers were enriched with contacts between residues of the same type due to the effects of structure symmetry. Proteins 2007;67:1060-1077. V C 2007 Wiley-Liss, Inc.

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 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.

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. Various levels of hierarchy in the protein spatial structure are discerned by analysis of the energy interaction between structural units of different scales. Thus the level of energy hierarchy plays role of sole parameter, and the method obviates the use of complicated geometrical criteria with numerous fitting parameters. The algorithm readily and accurately locates domains formed by continuous segments of the protein chain as well as those comprising non-sequential segments, sets no limit to the number of segments in a domain. We have analyzed 309 protein structures. Among 277 structures for which our results could be compared with the domain definitions made in other works, 243 showed complete or partial coincidence, and only in 34 cases the domain structures proved substantially different. The domains delineated with our approach may coincide with reference definition at different levels of the globule hierarchy. Along with defining the domain structure, our approach allows one to consider the protein spatial structure in terms of the spatial distribution of the interaction energy in order to establish the correspondence between the hierarchy of energy distribution and the hierarchy of structural elements.

Linker segments assuming the polyproline II type conformation within DNA-protein complexes were sought in protein and linker databases. Seventy-three linker-DNA complexes were found. The mean length of polyproline II type segments was six residues, and prolines were not predominant there. It was shown that the symmetrical position of prolines in these segments prevented the formation of the cooperative water network involving amide groups. An example of specific proline location in some motility apparatus proteins is presented.

The effect of imino acids on thermodynamic characteristics of the collagen-type structures in collagens from different sources is analyzed. It is demonstrated that the main mechanism of the entropy increase in the system water-protein with an increase in the imino acid content in a polypeptide chain is a change in the number of cooperative regions, detectable during the melting of a fibrous macromolecule. The variation of physical characteristics of the cooperative units determines, in particular, the variation of the hydrogen bond parameters, which appears as a broadening of the bands of NH valence vibrations and the half-widths of transitions in postdenaturation structures. Thus, the main mechanism of how imino acids act on the thermodynamic characteristics of collagens is connected with a complex nature of the denaturation process, when any substitutions of imino acids for amino acids significantly change the dehydration-hydration of the native and denatured states.

A complete classification of types of the protein secondary structure is developed on the basis of computer analysis of the crystallographic structural data deposited in the protein Data Bank. The majority of amino acid residues fall into five conformation types. A conclusion is drawn that the number of sequence variants of torsion angles phi, psi in globular proteins is limited and is essentially less than the number of possible amino acid sequences for this chain length. Along with alpha-helix and beta-structure, the distribution analysis assigning every maximum of distribution of amino acid conformations on Ramachandran map to a certain type of the secondary structure exposed a third type of the secondary structure that was previously neglected. This type of the structure is extended left-handed helical conformation, designated as mobile (M-) conformation. A full set of M-conformation fragments that seems to play a major role in protein globule dynamics has been obtained, a small radius of correlation for the polypeptide chain in M-conformation is demonstrated. It explains a prevalence of short segments of mobile conformation revealed in globular proteins. For secondary structure types, the frequency of occurrence of amino acid residues has been computed.

On the basis of complete scanning through conformational space of dihedral angles, twelve structural genera were obtained. Subsequent energy minimization within these genera yielded a limited set of duplexes with stacking: right-handed B-form (Wilkins type), B2-form (Watson and Crick type) and left-handed Ll-form (Sasisekharan type) and the new L2-form. In the polymeric DNA only right-handed double-helices are possible, the left-handed helices are forbidden due to poor 1–5 interchain contacts. In contrast, for short fragments the left- and right-handed helicek have practically the same energies providing some physical ground for side-by-side form, which biologically is possible as recombination form and may be as replication form.

The problem of deducing the DNA structure with correct base pairing and specific symmetry is formulated in the form of algebraic equations. In this way the number of independent variables determining double-helix conformation can be reduced from six to four. The effectiveness of the method is illustrated by the computations for A-DNA and B-DNA. The method allows one to predict all possible conformations of the complementary nucleic acids.

A computational method is elaborated for studying the water environment around regular polynucleotide duplexes; it allows rigorous structural information on the hydration shell of DNA to be obtained. The crucial aspect of this Monte Carlo simulation is the use of periodical boundary conditions. The output data consists of local maxima of water density in the space near the DNA molecule and the properties of one- and two-membered water bridges as function of pairs of polar groups of DNA. In the present paper the results for poly(dG) · poly(dC) and poly(dG-dC) · poly(dG-dC) are presented. The differences in their hydration shells are of a purely structural nature and are caused by the symmetry of the polar groups of the polymers under study, the symmetry being reflected by the hydration shell. The homopolymer duplex hydration shell mirrors the mononucleotide repeat. The water molecules contacting the polynucleotide in the minor groove are located nearly in the plane midway between the planes of successive base pairs. One water molecule per base pair forms a water bridge facing two polar groups of bases from adjacent base pairs and on different strands making a “spine”-like structure. In contrast, the major groove hydration is stabilized exclusively by two-membered water bridges; the water molecules deepest in the groove are concentrated near the plane of the corresponding base pair. The alternating polymer is characterized by a marked dyad symmetry of the hydration shell corresponding to the axis between two successive base pairs. The minor groove hydration of the dCpdG step resembles the characteristic features of the homopolymer, but the bridge between the O2 oxygens of the other base-stacking type is formed by two water molecules. The major groove hydration is characterized by high probability of one-membered water bridges and by localization of a water molecule on the dyad axis of the dGpdC step. The found structural elements are discussed as reasonable invariants of a dynamic hydration shell.

Conformational analysis of triple helics of a type of collagen was performed with typical collagen tripeptide sequences based on Gly-Pro-Ala, Gly-Ala-Hyp, and Gly-Ala-Ala. During energy minimization, the possibility of continual deformation of the pyrrolidine cycle was taken into account in order to achieve better accuracy in the resulting structure. The (Gly-Pro-Ala)n structure is almost isomorphic to the (Gly-Pro-Hyp)n structure obtained in the previous work [Tumanyan, V. G. & Esipova, N.G. (1982) Biopolymers21, 475–497]. For a collagen-type structure, the optimal conformation of (Gly-Ala-Hyp)n tends to have a decreased unit twist (t = 15°), although the energy advantage with respect to the conformation with t = 45° is not so significant. A similar situation is observed for (Gly-Ala-Ala)n. In this case, the energy decrease during unwinding to t = 15° from t = 45° is quite small. The conformations of (Gly-Ala-Hyp)n and (Gly-Ala-Ala)n with t = 15° exhibit a similarity with a triple complex of polyproline II helices—a noncoiled coil such as (Gly-Pro-Hyp)n and (Gly-Pro-Ala)n. A similar structure may be postulated for subcomponent cq1 of the first component of a human complement containing substantial Gly-X-Pro and Gly-X-Y tripeptide derivatives in the primary structure (X, Y = any amino acid). The results suggest that the observed helical symmetry of collagen (t = 36°) is a consequence of superposition of diffraction patterns (for sufficiently long segments) from various helices (t varies from ∼15° for Gly-X-Hyp and Gly-X-Y to ∼56° for Gly-Pro-Ala). For short alternating segments, some unification of different helical structures is possible.

The structure of collagen is investigated by means of an energy minimization procedure. There is a rather broad trough of potential energy which involves a variety of structures (all collagen II type). It is seen that the unit height and the unit twist can strongly deviate from 2.86 A, 36° values. Some of the good conformations are presented.

Analysis of 68 proteins from Protein Data Bank disclosed a new widely spread type of the secondary structure that is designated as mobile (M-) conformation. Helical parameters of M-conformation are close to the poly-L-proline II type helix. Its occurrence in globular proteins approximates that of the beta-sheet. The angles corresponding to the position of the M-conformation maximum in distribution of amino acid residues on a conformational map are phi: -65 degrees, psi: 140 degrees. Unique features and high occurrence in proteins make it possible to distinguish the M-conformation as an independent third type of the secondary structure in globular proteins, that should be included in the present classification.

A complete analysis of all possible conformations with correct hydrogen bonds of the collagen II type was performed on the basis of developed simultaneous equations. Using a unimodal search (by varying Ψ3), the energetically favorable structure was obtained. No other energetically satisfactory structural solutions are possible. The next aim was to obtain a precise model of the molecule. The program used includes a subroutine for continual deformation of the pyrrolidine rings. The set of parameters determining the structure consists of 14 independent variables (8 dihedral and 6 bond angles). As starting points for the energy optimization, conformations produced by scanning and some structures from previous work were used. The final structures (practically the same for both polymers) have helix parameters h = 0.285 nm and t = 52°, which are in excellent agreement with the 7/2 symmetry of diffraction data. The conformations of the pyrrolidine rings are of the B type, i.e., C2-Cβ-exo-Cγ-endo. For both polypeptides, the conformations of imino acids in position 3 of the triplet are the same; in position 2, however, they are slightly different. The difference in diffraction patterns for the 7/2 and 10/3 helices is discussed.

In the work the arguments are presented in favor of the idea on the role of conformationally stable oligopeptides 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 participating 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 acid sequence in collagen helix formation is shown.