Jumber Kereselidze*, George Mikuchadze and Lia Bobokhidze
Department of Chemistry, Ivane Javkhishvili Tbilisi State University, Tbilisi 0179, Georgia
Received date: January 16, 2017; Accepted date: February 01, 2017; Published date: February 03, 2017
Citation: Kereselidze J, Mikuchadze G, Bobokhidze L (2017) Quantum-Chemical Study of the Propensity of the Amino Acid Pairs for the Peptide Bond Formation. J Comput Sci Syst Biol 10: 014-016. doi:10.4172/jcsb.1000241
Copyright: © 2017 Kereselidze J, et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Visit for more related articles at Journal of Computer Science & Systems Biology
The Modern method of quantum chemistry - the density functional theory (DFT) was used to quantitative describe the formation of the peptide bond between amino acid pairs. In particular the formula of parameter of propensity (Kp), which is a function of six variables: the length and order of the CO and NH bonds, (RCO, RNH, PCO PNH), the activation energy for the formation of the peptide bond (ΔE#), and difference between charges of the carbon atom of the carbonyl group and the amino nitrogen atom (Δq) was constructed. By means of the proposed formula the Kp parameter for 400 amino acid pairs was calculated. Among them only 26 amino acid pairs are most likely to take part in the synthesis of proteins that have been selected based on the value of the parameter Kp. This approach may have important meaning for quantitative description of the amino acid sequences in proteins.
Amino acids; Peptide bonds; Parameter of propensity; DFT calculations
The theoretical description of biochemical processes is the main focus of modern natural science - Biophysical Chemistry. In recent years, for the quantitative description of complex biochemical processes are widely used modern quantum chemistry methods based on density functional theory (DFT). Including for research of peptide bond formation mechanism . It is assumed that the inductive and field effects of amino acid side chains have an essential effect on peptide bond formation . The quantum-mechanical study of different possible mechanisms of peptide synthesis in the ribosome has been carried out using density functional also . Analysis of a database of protein sequences for all possible binary patterns of polar and non-polar amino acid residues revealed that alternating patterns occur significantly less often than others with similar composition . To facilitate understanding of the information available for protein structures, has been constructed the structural classification of proteins (scop) database. This database provides a detailed and comprehensive description of the structural and evolutionary relationships of the proteins of known structure . Analysis of extant proteomes has the potential of revealing how amino acid frequencies within proteins have evolved over biological time. Evidence presented here indicates, that cysteine, tyrosine, and phenylalanine residues have substantially increased in frequency . To understand more fully how amino acid composition of proteins has changed over the course of evolution, a method has been developed for estimating the composition of proteins in an ancestral genome. The method was used to infer the amino acid composition of a large protein set in the Last Universal Ancestor (LUA) of all extant species. It is proposed that the inferred amino acid composition of proteins in the LUA probably reflects historical events in the establishment of the genetic code .
Protein sequences contain many local regions of low compositional complexity. These include different types of residue clusters, some of which contain homopolymers, short period repeats or aperiodic mosaics of a few residue types. Several different formal definitions of local complexity and probability are presented and are compared for their utility in algorithms for localization of such regions in amino acid sequences and sequence databases . The occurrence of all diand tripeptide segments of proteins was counted in a large data base containing about 119 000 residues . Systematic conformational analysis study of the tripeptide units (Gly-X-Pro) and (Gly-Pro-X), with X=Pro, Ala, Ser, Val, Leu, Ile, and Phe it has been reported. The lowenergy conformers obtained by quantum computations are discussed with respect to other theoretical investigations and experimental protein structures . Pairing of Amino acids is only possible on a parallel β ribbon and involves both the polypeptide backbones and the side chains. Model building revealed that of the 210 possible amino acid pairs of the standard 20 amino acids, no more than 26 could be built to meet standard criteria for bonding. Of these 26, 14 were found to be genetically encoded when the codons are read as if they paired in a parallel manner .
Density functional theory (DFT) is a computational quantum mechanical method used in physics, chemistry and biology for investigate the electronic structure in particular atoms, and molecules . The properties of a many-electron system can be determined by using functionals, which in this case is the spatially dependent electron density. Hence the name density functional theory comes from the use of functionals of the electron density. DFT is among the most popular and versatile methods available in computational biology. Unlike the wave function, which is not a physical reality, electron density is a physical characteristic of all molecules. The electron density is a function with three variables - x−, y−, and z− position of the electrons. Hybrid methods, as the name suggests, attempt to incorporate some of the more useful features from ab initio methods (specifically Hartree-Fock methods) with some of the improvements of DFT mathematics. Hybrid methods, such as B3LYP [13-15], tend to be the most commonly used methods for computational chemistry and Biology. Calculations were performed using software, “Priroda” in regime of the reaction coordinate .
The aim of this work is a quantitative description of the propensity of 400 pairs of 20 amino acids for the formation of peptide bond using quantum - chemical method of density functional theory (DFT). For the study of the conformation of 20 amino acids by Kolaskar and Ramabraham 400 options were selected also . As know, at the stage of activation of the formation of the peptide bond between two amino acids the C-O bond of the first and N-H bond of the second one are ruptured, what leads to C-N bond formation with release of water molecule (Figure 1). With the aim of quantitative description of peptide bond formation by the quantum-chemical method DFT the bond orders of carbonyl and amino groups (PCO, PNH) and corresponding bond lengths (RCO, RNH), difference values of charge between C and N atoms (Δq), as well as the activation energy (ΔE#) and the reaction energy (ΔE) have been calculated. In addition, it must satisfy the condition of exothermic formation of the peptide bond (ΔE<0). The increase of bond lengths RCO and RNH, as well as the difference of charges between carbon atom of carbonyl group and nitrogen atom the amino group (Δq) causes an increase in the value of the parameter of the propensity of amino acids of formation of peptide bond (KP). Hence, these values in the formula Kp placed in the numerator. On the other hand, the decrease of bond orders PCO and PNH, as well as the activation energy ΔE# the parameter of propensity of the peptide bond formation (KP) also increases. Hence, these values are in denominator of the formula. Consequently, the formula, for the quantitative description of the effect of R-groups of the formation of peptide bond was built.
From the analysis of the data of Table 1, one can make a general conclusion relatively of the parameter of propensity of amino acids for the formation of peptide bond (KP). In Particular, KP inversely proportional of the activation energy ΔE# of peptide bond formation. In this row the pairs of the amino acids are listed in order of increase of the values of KP and we expect that, they must be presented in proteins more often (for example 21-26 amino acid pairs, from Table 1). Such conclusion, at first glance may seem a less convincing, however from the point of view of the R-group’s influence on the ability of the peptide bond formation of amino acids Kp parameter can be successfully applied to the quantitative description of this process. In addition, the value of the reaction energy is always negative, which is common for exothermic reactions. In other cases, values of Kp is low, which indicates on difficulty of peptide bond formation. It is known that the pairing of amino acids in proteins are encoded by nucleotide base sequences of RNA. On the other hand, the formation of peptide bonds is much dependent on the influence of inductive and steric effects of R-groups. One can consider that, the energy (ΔE#, ΔE), electronic (qN, qC, PCO and PNH) and structural (RCO and RNH) values contained in the Table 1, has completely reasonable values, from the viewpoint of chemical reactions. Such choice of amino acid pairs is probably caused by different direction of influence of inductive and steric effect of R-groups with respect to formation of peptide bonds. It means, that one R-group can promote, but other one prevent the formation of peptide bond. The proposed quantum-chemical formula of amino acid pairs for the peptide bond formation, in our opinion may contribute to the prediction of the direction of peptide bond formation of a quantitative description of the propensity.
|N||Amino acid pairs||∆E#, kJ/mol||∆q||PCO||PNH||RCO, Å||RNH, Å||KP . 10-2||∆E, kJ/mol|
Table 1: The values of activation energy (ΔE#), difference of charge between carbon atom of carbonyl group and nitrogen atom of amino group (Δq), bond orders and bond lengths of CO and NH bonds (PCO, PNH, RCO, (Å) RNH, (Å) and the parameter of propensity of peptide bond formation (KP).
It is unknown to what extent the condensation of mixtures of amino acids occurs preferably randomly or between certain amino acids. In this case, the terminal amino acid residue of the polypeptide could influence the choice of the next amino acid. Such a mechanism could the determine the organization of the first polypeptides. The proposed quantum-chemical formula, which includes the calculated values of energetic, electronic and structural characteristics of amino acids can be applied for quantitative description of the propensity of amino acids for the peptide bond formation as these values quantitatively describe inductive and steric influence of R-groups of amino acids on the reaction center during the peptide bond formation. These factors together encoding nucleotide bases can have an important role in the study of the amino acid sequences in proteins. According to Atkins the orderliness can turn into chaos, and in certain conditions, from chaos it is possible the formation of the orderliness also . On the other hand, F. Chapeville  wrote that in nature, there are many different amino acids in the free state, but only 20 of them are used for protein synthesis. The reason for this limitation is not known, but it is apparently related to the fact that during the evolution of selection showed severe requirements for the physicochemical properties of these molecules, resulting in all the other amino acids were discarded. Based on the above we can assume that our results quantitatively describe prediction F. Chapeville. The proposed theoretical formula for the probability of the sequence of amino acids in the polypeptides can be used for building an adequate correlation with data of the morphological study of the structure of proteins.
This work was supported by the Shota Rustaveli National Science Foundation Grant No: 217732.
Make the best use of Scientific Research and information from our 700 + peer reviewed, Open Access Journals