S S symmetry
Article Protein Helical Structures: Defining Handedness and Localization Features
Alla E. Sidorova *,†, Ekaterina V. Malyshko *,† , Aleksey O. Lutsenko, Denis K. Shpigun and Olga E. Bagrova
Faculty of Physics, Lomonosov Moscow State University, Leninskie gory 1-2, 119234 Moscow, Russia; [email protected] (A.O.L.); [email protected] (D.K.S.); [email protected] (O.E.B.) * Correspondence: [email protected] (A.E.S.); [email protected] (E.V.M.) † Alla E. Sidorova and Ekaterina V. Malyshko contributed equally to this work as first authors.
Abstract: The quantitative evaluation of the chirality of macromolecule structures remains one of the exciting issues in biophysics. In this paper, we propose methods for quantitative analysis of the chirality of protein helical and superhelical structures. The analysis of the chirality sign of the protein
helical structures (α-helices and 310-helices) is based on determining the mixed product of every three consecutive vectors between neighboring reference points—α-carbons atoms. The method for evaluating the chirality sign of coiled-coil structures is based on determining the direction and value of the angle between the coiled-coil axis and the α-helices axes. The chirality sign of the coiled coil is calculated by averaging the value of the cosine of the corresponding angle for all helices forming the superhelix. Chirality maps of helical and superhelical protein structures are presented. Furthermore, we propose an analysis of the distributions of helical and superhelical structures in polypeptide chains of several protein classes. The features common to all studied classes and typical for each protein class are revealed. The data obtained, in all likelihood, can reflect considerations about molecular machines as chiral formations. Citation: Sidorova, A.E.; Malyshko, E.V.; Lutsenko, A.O.; Shpigun, D.K.; Keywords: chirality; proteins; secondary structure; helix; coiled coil; chirality map Bagrova, O.E. Protein Helical Structures: Defining Handedness and Localization Features. Symmetry 2021, 13, 879. https://doi.org/10.3390/ 1. Introduction sym13050879 Since the time of Louis Pasteur, one of the key categories in the triad of “symmetry— dissymmetry—asymmetry” in living systems has been considered chirality—the property Academic Editor: Victor Borovkov of objects or systems to be non-superimposable with their mirror images in any combination of translations and rotations in three-dimensional space [1]. It is well-known that proteins Received: 14 April 2021 synthesized in ribosomes contain only L-amino acids, and nucleic acids contain D-ribose Accepted: 11 May 2021 Published: 15 May 2021 and D-deoxyribose. For more than 100 years, the attention of researchers has been attracted to the mystery of the origin and evolutionary fixation of chiral dissymmetry in living
Publisher’s Note: MDPI stays neutral nature. There is no generally accepted answer to this question, although there are numerous with regard to jurisdictional claims in hypotheses [2–9]. However, we are interested in what unique or simply beneficial living published maps and institutional affil- nature has acquired on Earth, having taken the path of consistent symmetry breaking in iations. its development. In previously published works, relying on numerous sources, we developed and substantiated the tendency according to which the hierarchies of the structures of proteins and nucleic acids are characterized by a change in the type of symmetry and chirality sign. Starting from the level of asymmetric carbon in the amino acids of the primary structure of Copyright: © 2021 by the authors. Licensee MDPI, Basel, Switzerland. proteins, there is a tendency to alternate the chirality sign of the structural levels of proteins This article is an open access article L-D-L-D [10–12]. distributed under the terms and If we talk about the general trend, then the primary structure (polypeptide chain) is conditions of the Creative Commons characterized by completely left-handed chirality of amino acid residues. Higher levels Attribution (CC BY) license (https:// are represented by helical structures of various types (in this work, we do not consider β- creativecommons.org/licenses/by/ structures). Secondary structures are represented mainly by right-handed α-helices [13–16]. 4.0/). This conformation is energetically more preferred [17,18]. The next most common helical
Symmetry 2021, 13, 879. https://doi.org/10.3390/sym13050879 https://www.mdpi.com/journal/symmetry Symmetry 2021, 13, 879 2 of 17
secondary structure of proteins is 310-helices. They account for up to 15–20% of all protein helices [19]. Although left-handed helices are rare, they also have structural or functional significance [18,20,21]. At the next level of the hierarchical organization of protein structures, there are superhelical structures—coiled coil composed of α-helices [22–24]. In addition to well-known coiled-coil structures with a heptade sequence, there are also nontrivial ones. In the case of nonheptad sequence repeats, the heterochiral structure of L- and D-polypeptides contains large and likely more stable interfaces [25]. In this paper, the authors expect the heterochiral oligomers would be incompatible with α-helical supercoiling. To date, there is a significant number of works that offer various methods for assessing chirality [26–35]. The Ramachandran plot is a significant tool for describing a peptide backbone in the context of the backbone dihedral angles [36,37]. One of the concepts associated with the Ramachandran plot is a metric for backbone handedness [38]. It is based on interpreting a peptide backbone as a helix with axial and angular displacement. Every region handedness of the Ramachandran plot could be characterized by this metric for cis and trans backbones. The most comprehensive theoretical approach to assessing the chirality of helical structures was developed by M. Petitjean. Following this concept, the measure of chirality should be a continuous characteristic and be determined for a space of any dimension, and the chirality index should not depend on the method of choosing the mirror image [39–41]. The above methods are either highly specialized or rather complicated to implement, so it is impossible to determine the sign and measure of chirality of helical protein structures for large data sets. Since the discovery of coiled-coil structures, several approaches have been proposed to describe them as chiral objects. In [42], the authors proposed to relate the pitch and radius of the coiled coil with a change in the twist angle and a shift along the coiled-coil axis per amino acid. In another work [43], the calculation of the pitch of the double-helical coiled coil was made, taking into account the Cartesian vector coordinates of the average atomic positions of the core atoms of the N- and C-terminus. It should be noted that both methods presented do not allow one to reveal the chirality of coiled coils. Furthermore, the possibility of the formation of a right-handed structure is not taken into account. In work [44], it is proposed to determine the sign of the chirality (τ) of the surface curve for a classical coiled coil using the periodicity of the helix, the angle between adjacent residues, and the shift along the axis per amino acid: for right-handed structures, τ is positive, for left-handed ones—negative. However, the authors note that it is not obvious that at negative τ the superhelix is left-handed. In this work, we present our methods for the quantitative evaluation of the protein helical structure chirality. A necessary and sufficient condition is the mutual arrangement of α-carbons, which allows a reduction in the amount of processed information and is a clear advantage when processing large data arrays. To determine the chirality sign of coiled-coil structures, the twisting direction of each individual α-helix relative to the entire superhelix axis is taken into account. Based on the methods, computer programs in Python 3.7 are implemented. Chirality maps of helical and coiled-coil structures of proteins are obtained. The accuracy of the maps is confirmed by the analysis of real structures. Another aspect of this work is the regularity identification of the protein “chiral” structure using the example of helical and superhelical structures and their functions. The protein structure determines its function and role in the biological system. However, despite the abundant data on the functions of various proteins, the functions of some of them remain unclear. One of the methods for determining the correspondence of the protein structure and function and its role in metabolism is the method of comparison with already known proteins with similar structures. In the present work, the regularities in the distribution of helices and coiled coils along the polypeptide chain, which form a dynamic backbone of protein molecules, are considered. Research on patterns and relationships between the chiral protein structure and its functions could be used to develop the concept of a protein as a molecular machine. Symmetry 2021, 13, 879 3 of 18
relationships between the chiral protein structure and its functions could be used to Symmetry 2021 13 , , 879 develop the concept of a protein as a molecular machine. 3 of 17
2. Computational Methodology 2. ComputationalThe object of Methodology research is helical and coiled-coil structures from the PDB, CC+ databasesThe object [45,46]. of research is helical and coiled-coil structures from the PDB, CC + databases [45,46]. 2.1. Method for Quantitative Analysis of the Chirality of Protein Helical Structures 2.1. MethodPreviously, for Quantitative we proposed Analysis a method of the for Chirality quantita of Proteintive analysis Helical of Structures the chirality of protein helicalPreviously, structures, we based proposed on the a method construction for quantitative of consecutive analysis vectors of the between chirality neighboring of protein helicalcarbons structures, in the chain based and on determination the construction of thei ofr consecutive vector product vectors [47]. between The method neighboring is based carbonson the inmutual the chain arrangement and determination of α-carbons of their (C vectorα) atoms product in the [47 amino]. The methodacid chain—the is based onreference the mutual points arrangement in the model. of α -carbonsThe skeleton (Cα) atomsof α-carbons in the amino and the acid helix chain—the axis are reference directed pointsfrom the in theC-terminus model. The to the skeleton N- terminus. of α-carbons As a and result the of helix constructing axis are directed vectors from between the α C-terminusneighboring to atoms the N- (from terminus. the previous As a result to ofthe constructing subsequent) vectors for n betweenatoms (C neighboringα1…C n), the α α atomsdirection (from vector the previous d and the to sum the subsequent) of the vector for products n atoms s (C were1 ... determined.C n), the direction The sign vector of the dcosineand the of sumthe ofangle the vectorbetween products these svectorswere determined. determines Thethe signchirality of the sign: cosine in of the the case angle of betweencos∠(d, s these) < 0—D; vectors in the determines case of cos the∠( chiralityd, s) > 0—L. sign: The in thelength case of of the cos product∠(d, s)
FigureFigure 1.1.Stages Stages of of graphical graphical construction construction of of triplets triplets of of vectors vectors for for calculating calculating a mixed a mixed vector vector product. C—Cproduct.α atoms C—C inα atoms the helix, in the used helix, in used the model in the asmodel reference as reference points; points; v1, v2—vectors v1, v2—vectors between between two neighboringtwo neighboring reference reference points. points. Construction Constructi stageson stages 1 and 21 areand marked 2 are marked in green: in green: (a) vectors (a) vectors v1, v2, v1, v3; (bv2,) vectors v3; (b) v2,vectors v3, v4. v2, v3, v4.
ForFor everyevery threethree consecutiveconsecutive vectors,vectors, theirtheir mixedmixed productproduct isis calculatedcalculated (taking(taking intointo accountaccount theirtheir coordinatescoordinates in in three-dimensional three-dimensional space): space): v , v , v = − ) + − ) + − ) (1) ([v1, v2], v3) = (y 1z2 − y2z1)x3 + (z 1x2 − z2x1)y3 + (x 1y2 − x2y1)z3 (1) The sum of all mixed products can be represented as a characteristic of the helical structureThe sum chirality, of all and mixed the productsequation canfor its be description represented can as be a characteristic presented as follows: of the helical structure chirality, and the equation for its description can be presented as follows:
= n−3([ ] ) χtotal ∑i=1 v i, vi+1 , vi+2 , (2) The sign of the mixed product determines the chirality sign of the protein helical structure.
Symmetry 2021, 13, 879 4 of 18