Prediction of β-Turns

Kohonen's self-organization model, a neural network model, is applied to predict the β-turns in proteins. There are 455 β-turn tetrapeptides and 3807 non-β-turn tetrapeptides in the training database. The rates of correct prediction for the 110 β-turn tetrapeptides and 30,229 non-β-turn tetrapeptides in the testing database are 81.8% and 90.7%, respectively. The high quality of prediction of neural network model implies that the residue-coupled effect along a polypeptide chain is important for the formation of reversal turns, such as β-turns, during the process of protein folding.     

Classification and Prediction of β-Turn Types

Although a β-turn consists of only four amino acids, it assumes many different types in proteins. Is this basically dependent on the tetrapeptide sequence alone or is it due to a variety of interactions with the other part of a protein? To answer this question, a residue-coupled model is proposed that can reflect the sequence-coupling effect for a tetrapeptide in not only a β-turn or non-β-turn, but also different types of a β-turn. The predicted results by the model for 6022 tetrapeptides indicate that the rates of correct prediction for β-turn types I, I′, II, II′, VI, and VIII and non-β-turns are 68.54%, 93.60%, 85.19%, 97.75%, 100%, 88.75%, and 61.02%, respectively. Each of these seven rates is significantly higher than $\frac{1}{7}$ = 14.29%, the completely randomized rate, implying that the formation of different β-turn types or non-β-turns is considerably correlated with the sequences of a tetrapeptide.     

Conformational analysis of a snake neurotoxin by prediction from sequence, circular dichroism, and raman spectroscopy.

The secondary structure of the major neurotoxin from the sea snake Lapemis hardwickii was investigated by several methods of conformational analysis: structure prediction, circular dichroism, and laser Raman spectroscopy. From the primary structure, secondary structure prediction yielded two regions of β-sheet structure at residues 1–7 and 41–45. β-Turns were predicted at residues 14–17, 20–23, 30–33, 37–40, and 46–49. From the predictions, the toxin appears to be composed of approximately 20% β-sheet and 33% β-turn. The CD spectrum of the native toxin appears to be a hybrid of model spectra for β-sheet and β-turn proteins. The pH perturbation studies on the toxin observed by CD demonstrated that the toxin is a very stable molecule except at extremely high or low pH values. The Raman data indicated that the toxin contains both antiparallel β-sheet and β-turn structure. Using two methods of secondary structure quantitation from Raman spectra the molecule was calculated to contain 35% β-sheet from one method and 27% from the other. Overall, the various methods demonstrate that the toxin is composed of β-sheet and β-turn structure with little or no α-helix present. From the comparison of these different techniques appreciation can be gained for the necessity of several methods when identifying and quantitating secondary structure.     

Prediction of β-Turns

Kohonen's self-organization model, a neural network model, is applied to predict the -turns in proteins. There are 455 -turn tetrapeptides and 3807 non--turn tetrapeptides in the training database. The rates of correct prediction for the 110 -turn tetrapeptides and 30,229 non--turn tetrapeptides in the testing database are 81.8% and 90.7%, respectively. The high quality of prediction of neural network model implies that the residue-coupled effect along a polypeptide chain is important for the formation of reversal turns, such as -turns, during the process of protein folding.     

An effect of ATP on the permeability of transport-competent plasma membrane vesicles from mouse fibroblasts.

Methods for predicting peptide chain conformation have been applied to amino acid sequences adjacant to the carbohydrate attachment sites of glycoproteins containing the $\text{N}$-glycosylamine type of protein-carbohydrate linkage. Of 31 glycosylated residues examined 30 occur in sequences favouring turn or loop structures. Twentytwo of the glycosylated asparagine residues occur in tetrapeptides predicted to have the β-turn conformation. Carbohydrate attachment is therefore associated with peptide sequences which favour the formation of β-turn or other turn or loop structures.     

Peptides containing the sulfonamide junction: Synthesis,structure, and conformation of Z-Tau-Pro-Phe-NHiPr

The taurine (Tau) containing tripeptide derivative Z-Tau-Pro-Phe-NHiPr (1) has been synthesized as suitable sulfonamido-pseudopeptide model to investigate formation and conformational properties of folded secondary structures stabilized by intramolecular H bonds directly involving the sulfonamide junction. In the crystal the pseudopeptide 1 adopts a type I β-turn with the Pro and Phe residues located at the (i + 1) and (i + 2) corner positions, respectively. The turn is stabilized by a 4 → 1 H bond engaging one of the SO₂ oxygen atoms and the isopropylamide NH. In CDCl₃ solution the β-turn folding is accompanied by a γ-turn centered at the Pro and involving a 3 → 1 H bond between the SO₂ and the Phe NH. A comparison of the structural and conformational properties found in 1 with those of the already known sulfonamido-pseudopeptides, with particular reference to the models containing the Tau-Pro junction, is also reported. © 1997 John Wiley & Sons, Inc. Biopoly 41: 555–567, 1997.     

Crystallization,molecular replacement solution,and refinement of tetrameric β-amylase from sweet potato

Sweet potato β-amylase is a tetramer of identical subunits, which are arranged to exhibit 222 molecular symmetry. Its subunit consists of 498 amino acid residues (Mr 55,880). It has been crystallized at room temperature using polyethylene glycol 1500 as precipitant. The crystals, growing to dimensions of 0.4 mm × 0.4 mm × 1.0 mm within 2 weeks, belong to the tetragonal space group P4₂2₁2 with unit cell dimensions of a = b = 129.63 Å and c = 68.42 Å. The asymmetric unit contains 1 subunit of β-amylase, with a crystal volume per protein mass (V_M) of 2.57 ų/Da and a solvent content of 52% by volume. The three-dimensional structure of the tetrameric β-amylase from sweet potato has been determined by molecular replacement methods using the monomeric structure of soybean enzyme as the starting model. The refined subunit model contains 3,863 nonhydrogen protein atoms (488 amino acid residues) and 319 water oxygen atoms. The current R-value is 20.3% for data in the resolution range of 8–2.3 Å (with 2 σ cut-off) with good stereochemistry. The subunit structure of sweet potato β-amylase (crystallized in the absence of α-cyclodextrin) is very similar to that of soybean β-amylase (complexed with α-cyclodextrin). The root-mean-square (RMS) difference for 487 equivalent Cα atoms of the two β-amylases is 0.96 Å. Each subunit of sweet potato β-amylase is composed of a large (α/β)₈ core domain, a small one made up of three long loops [L3 (residues 91–150), LA (residues 183–258), and L5 (residues 300–327)], and a long C-terminal loop formed by residues 445–493. Conserved Glu 187, believed to play an important role in catalysis, is located at the cleft between the (α/β)₈ barrel core and a small domain made up of three long loops (L3, L4, and L5). Conserved Cys 96, important in the inactivation of enzyme activity by sulfhydryl reagents, is located at the entrance of the (α/β)₈ barrel. © 1995 Wiley-Liss, Inc.     

β-Turn structure and intramolecular interaction of tetrapeptides containing Asp and Lys

Nmr and CD studies of terminally protected tetrapeptides were carried out in aqueous and DMSO solutions to investigate the formation and stabilization of the β-turn structure. Boc-Gly-Lys-Asp-Gly-OMe and Boc-Asp-Lys-Asp-Gly-OMe appear to have a tendency to adopt a β-turn structure in aqueous solution from the CD spectra and temperature-dependence studies of the amide proton chemical shifts. The side-chain conformation of the Asp residue depends greatly on its ionization state but was not affected by the deprotonation of the neighboring Lys side chain. There is evidence for an intramolecular interaction between the Asp and Lys side chains of Boc-Gly-Lys-Asp-Gly-OMe. Such an interaction can contribute to the stabilization of the β-turn structure.     

A comparative NMR and molecular dynamics study of the conformations of bradykinin B1 and B2, B2, and B1-specific receptor antagonists B-9430, B-9436, and B-9858

Extensive proton magnetic resonance experiments were carried out on three bradykinin peptide antagonists B-9430, B-9436, and B-9858 in aqueous solutions as well as in sodium dodecylsulphate micelles (B-9430 and B-9436) and CD₃OH/H₂O (60%/40%) mixtures for B-9858. All three peptides showed no observable secondary structure in aqueous solution. However, in their respective structure-inducing solvents, B-9430 (B₁ and B₂ receptor antagonist) and B-9436 (a B₂ receptor antagonist) exhibit a type II β-turn involving residues 2–5, and B-9430 also exhibits a type II′ β-turn involving residues 6–9 (in sodium dodecylsulfate micellar solutions), whereas B-9858, a B₁-specific receptor antagonist, exhibits only a type II β-turn involving residues 2–5 (in CD₃OH/H₂O solutions). Simulated annealing calculations on B-9858 confirm the experimental conclusions based on the nmr data. In addition, simulated annealing of the (2S, 3aS, 7aS)-octahydroindole-2-carboxylic acid (Oic residue), which is present in two of the three decapeptides studied, show that the one-chair conformation of the six-membered ring predominates, in agreement with the experimental data. The activities of these peptides are compared with their secondary structures and the specific receptor activity appears to depend on the presence of specific amino acid residues, such as N-(2-indanyl)glycine (Nig) and D[α-(2-indanyl)glycine] (D-Igl) as well as on elements of secondary structure. © 1997 John Wiley & Sons, Inc. Biopoly 42: 521–535, 1997     

β-Turn conformations in crystal structures of model peptides containing α,α-Di-n-propylglycine and α,α-Di-n-butylglycine

The crystal state conformations of three peptides containing the α,α-dialkylated residues. α,α-di-n-propylglycine (Dpg) and α,α-di-n-butylglycine (Dbg), have been established by x-ray diffraction. Boc-Ala-Dpg-Alu-OMe (I) and Boc-Ala-Dbg-Ala-OMe (III) adopt distorted type II β-turn conformations with Ala (1) and Dpg/Dbg (2) as the corner residues. In both peptides the conformational angles at the Dxg residue (I: ? = 66.2°, ψ = 19.3°; III: ? = 66.5°. ψ = 21.1°) deviate appreciably from ideal values for the i + 2 residue in a type II β-turn. In both peptides the observed (N…O) distances between the Boc CO and Ala (3) NH groups are far too long (1: 3.44 Å: III: 3.63 Å) for an intramolecular 4 → 1 hydrogen bond. Boc-Ala-Dpg-Ata-NHMe (II) crystallizes with two independent molecules in the asymmetric unit. Both molecules HA and HB adopt consecutive β-turn (type III-III in HA and type III-I in IIB) or incipient 3₁₀-helical structures, stabilized by two intramolecular 4 → 1 hydrogen bonds. In all four molecules the bond angle N-C^α-C′ (τ) at the Dxg residues are ≥ 110°. The observation of conformational angles in the helical region of ?,ψ space at these residues is consistent with theoretical predictions. © 1995 John Wiley & Sons, Inc.     

Conformational interconversions in peptide β-turns: Discrimination between enantiomeric conformations by chiral perturbation

The crystal structure of the model tripeptide Boc-Aib-Gly-Leu-OMe ( 1 ) reveals two independent molecules in the asymmetric unit that adopt “enantiomeric” type I and type I′ β-turn conformations with the Aib and Gly residues occupying the corner (i + 1 and i + 2) positions. ¹³C cross polarization and magic angle sample spinning spectra in the solid state also support the coexistence of two conformational species. ¹³C-nmr in CDCl₃ establishes the presence of a single species or rapid exchange between conformations. 400 MHz ¹H-nmr provides evidence for conformational exchange involving a major and minor species, with β-turn conformations supported by the low solvent exposure of Leu(3) NH and the observation of N_iH ↔ N_i+1H nuclear Overhauser effects. CD bands in the region 190–230 nm are positive, supporting a major population of type I′ β-turns. The isomeric peptide, Boc-Gly-Leu-Aib-OMe ( 2 ), adopts an “open” type II′ β-turn conformation in crystals. Solid state and solution nmr support population of a single conformational species. Chiral perturbation introduced outside the folded region of peptides may provide a means of modulating screw sense in achiral sequences. © 1998 John Wiley & Sons, Inc. Biopoly 45: 191–202, 1998     

β-Turn preference of tetrapeptide sequences as analyzed by CD spectra of their Dnp-pNA derivatives

The Dnp-pNA derivatives of more than 30 tetrapeptides were synthesized and subjected to CD measurements in MeOH. The CD spectra were analyzed on the basis of exciton chirality theory and were found to be useful for estimating β-turn preference of the parent tetrapeptides based on the following facts: (1) when a compound of the series takes the βturn conformation, a pair of strong CD bands with opposite sign are observed at 305 and 350 nm, and (2) when it is in a random or nonfolded conformation, no significant bands are observed at the wavelengths. The intensity of the 350-nm band reflects especially well the βturn preference. The CD spectra can be classified into two groups according to the sign of the band at 350 nm, and the groups seem to correspond to the type of βturn. By summarizing the results obtained, some generalizations concerning the sequence–conformation relationship are proposed.     

β-Alanine containing cyclic peptides with predetermined turned structure. V

In the present paper we describe the synthesis, purification, single crystal x-ray analysis, and solution structural characterization by nmr spectroscopy, combined with restrained molecular dynamic simulations, of the cyclic hexapeptide cyclo-(Pro-Phe-β-Ala-Phe-Phe-β-Ala). The peptide was synthesized by classical solution methods and the cyclization of the free hexapeptide was accomplished in good yields in diluted methylenechloride solution using N, N-dicyclohexyl-carbodiimide. The compound crystallizes in the monoclinic space group P2₁ from methanol/ethyl acetate. The molecule adopts in the solid state a conformation characterized by cis β-Ala⁶-Pro¹ peptide bond. The α-amino acid residues are at the corner positions of turned structures. The Pro¹-Phe² segment is incorporated in a pseudo type I β-turn, while Phe⁴-Phe⁵ is in a typical type I β-turn. Assignment of all ¹H and ¹³C resonances was achieved by homo- and heteronuclear two-dimensional techniques in dimethylsulfoxide (DMSO) solutions. The conformational analysis was based on inter-proton distances derived from rotating frame nuclear Overhauser effect spectroscopy spectra and homonuclear coupling constants. Restrained molecular dynamic simulation in vacuo was also performed to built refined molecular models. The molecule is present in DMSO solution as two slowly interconverting conformers, characterized by a cis-tran isomerism around the β-Ala⁶-Pro¹ peptide bond. This work confirms our expectations on the low propensity of β-alanyl residues to be positioned at the corners of turned structure. © 1994 John Wiley & Sons, Inc.     

Type I and II β-turns prediction using NMR chemical shifts

A method for predicting type I and II β-turns using nuclear magnetic resonance (NMR) chemical shifts is proposed. Isolated β-turn chemical-shift data were collected from 1,798 protein chains. One-dimensional statistical analyses on chemical-shift data of three classes β-turn (type I, II, and VIII) showed different distributions at four positions, (i) to (i + 3). Considering the central two residues of type I β-turns, the mean values of C_ο, C_α, H_N, and N_H chemical shifts were generally (i + 1) > (i + 2). The mean values of C_β and H_α chemical shifts were (i + 1) < (i + 2). The distributions of the central two residues in type II and VIII β-turns were also distinguishable by trends of chemical shift values. Two-dimensional cluster analyses on chemical-shift data show positional distributions more clearly. Based on these propensities of chemical shift classified as a function of position, rules were derived using scoring matrices for four consecutive residues to predict type I and II β-turns. The proposed method achieves an overall prediction accuracy of 83.2 and 84.2 % with the Matthews correlation coefficient values of 0.317 and 0.632 for type I and II β-turns, indicating that its higher accuracy for type II turn prediction. The results show that it is feasible to use NMR chemical shifts to predict the β-turn types in proteins. The proposed method can be incorporated into other chemical-shift based protein secondary structure prediction methods.     

Type II β-turn conformation of pivaloyl-L-prolyl-α-aminoisobutyryl-N-methylamide: Theoretical,spectroscopic, and X-ray studies

Pivaloyl-L -Pro-Aib-N-methylamide has been shown to possess one intramolecular hydrogen bond in (CD₃)₂SO solution, by ¹H-nmr methods, suggesting the existence of β-turns, with Pro-Aib as the corner residues. Theoretical conformational analysis suggests that Type II β-turn conformations are about 2 kcal mol^?1 more stable than Type III structures. A crystallographic study has established the Type II β-turn in the solid state. The molecule crystallizes in the space group P2₁ with a = 5.865 Å, b = 11.421 Å, c = 12.966 Å, β = 97.55°, and Z = 2. The structure has been refined to a final R value of 0.061. The Type II β-turn conformation is stabilized by an intramolecular 4 → 1 hydrogen bond between the methylamide NH and the pivaloyl CO group. The conformational angles are ?_Pro = ?57.8°, ψ_Pro = 139.3°, ?_Aib = 61.4°, and ψ_Aib = 25.1°. The Type II β-turn conformation for Pro-Aib in this peptide is compared with the Type III structures observed for the same segment in larger peptides.     

The pili of Pseudomonas aeruginosa strains PAK and PAO bind specifically to the carbohydrate sequence βGalNAc(1–4)βGal found in glycosphingolipids asialo-GM1 and asialo-GM2

Pseudomonas aeruginosa employs pili to mediate adherence to epithelial cell surfaces. The pilus adhesin of P. aeruginosa strains PAK and PAO has been shown to bind to the glycolipid asialo-GM₁ (Lee et al., 1994 —accompanying article). PAK and PAO pili were examined for their abilities to bind to the synthetic βGalNAc(1–4)βGal (a minimal structural carbohydrate receptor sequence of asialo-GM₁ and asialo-GM₂ proposed by Krivan et al., 1988a) using solid-phase binding assays. Both pill specifically bound to βGalNAc(1–4)βGal. The binding of βGal-NAc(1–4)βGal-Biotin to the Immobilized PAK and PAO pili was inhibited by corresponding free pili. The receptor binding domain of the PAK pilus resides in the C-terminal disulphide-looped region (residues 128–144) of the pilin structural subunit (Irvin et al., 1989). Biotinylated synthetic peptides corresponding the C-terminal residues 128–144 of P. aeruginosa PAK and PAO pilin molecules were shown to bind to the βGalNAc(1–4)βGal-(bovine serum albumin (BSA)). The binding of biotinylated peptides to βGalNAc-(1–4)βGal-BSA was inhibited by PAK pili, Ac-KCTSDQDEOFIPKGCSK-OH (AcPAK(128–144)_ox-OH) and Ac-ACKSTQDPMFTPKGCDN-OH (AcPAO(128–144)_ox-OH) peptides. (In these peptides Ac denotes N^α -acetylation of the N-terminus, -OH means a peptide with a free a-carboxyl group at the C-terminus and the‘ox’denotes the oxidation of the sulphhydryl groups of Cys–129 and Cys–142.) Both acetylated peptides were also able to inhibit the binding of βGalNAc(1–4)βGal-biotin to the corresponding BSA-Peptide(128–144)_ox-OH conjugates. The βGlcNAc(1–3)βGal(1–4)βGlc-biotin conjugate was unable to specifically bind to either Immobilized PAK and PAO pili or the respective C-termlnal peptides. The data above demonstrated that the P. aeruginosa pili recognize asialo-GM₁ receptor analogue and that βGalNAc(1–4)βGal disaccharlde is sufficient for binding. Furthermore, the binding to βGalNAc(1–4)βGal was mediated by residues 128–144 of the pilin subunit.     

Sequential polypeptides containing arginine as histone models: synthesis and conformational studies

The sequential polypeptides (L -Arg-X-Gly)_n, where X represents amino acid residues Ala, Val, and Leu, were prepared as models of arginine-rich histones to be used in studying their structure and their interactions with DNA. The polymerization was carried out on the pentachlorophenyl active esters of the appropriate tripeptides, while the toluene-4-sulfonyl group was used for protecting the arginine guanido group. CD was employed to investigate the conformation of (L -Arg-X-Gly)_n polymers in aqueous solutions, at different pH, as well as in trifluoroenthanol and hexafluoroisopropyl alcohol solutions. In aqueous solutions (at pH 7 and 12) the prepared sequential polymers behaved as a random coil. The CD spectra in various trifluoroethanol–water or hexafluoroisopropyl alcohol–water mixtures indicated that the degree of helical conformation of the studied polytripeptides increased in the order of Ala → Val → Leu. The opposite was true for the β-structure. Characteristics of β-turn are excluded from the poly(L -Arg-L -Leu-Gly), which assumed the most pronounced helical conformation. The poly(L -Arg-L -Val-Gly) exerts a significant preference to the β-turn structure compared to that of poly(L -Arg-L -Ala-Gly). Thus the probability for helical, β-structure or β-turn conformations of the polymers was analyzed in relation to the bulkiness and length, and to the special features of the X-residue side chain (β-branching). We concluded that the prepared sequential arginine-containing polypeptides are plausible models for histone fractions, f₃ and f₂α₁.     

NetTurnP--neural network prediction of beta-turns by use of evolutionary information and predicted protein sequence features

β-turns are the most common type of non-repetitive structures, and constitute on average 25% of the amino acids in proteins. The formation of β-turns plays an important role in protein folding, protein stability and molecular recognition processes. In this work we present the neural network method NetTurnP, for prediction of two-class β-turns and prediction of the individual β-turn types, by use of evolutionary information and predicted protein sequence features. It has been evaluated against a commonly used dataset BT426, and achieves a Matthews correlation coefficient of 0.50, which is the highest reported performance on a two-class prediction of β-turn and not-β-turn. Furthermore NetTurnP shows improved performance on some of the specific β-turn types. In the present work, neural network methods have been trained to predict β-turn or not and individual β-turn types from the primary amino acid sequence. The individual β-turn types I, I', II, II', VIII, VIa1, VIa2, VIba and IV have been predicted based on classifications by PROMOTIF, and the two-class prediction of β-turn or not is a superset comprised of all β-turn types. The performance is evaluated using a golden set of non-homologous sequences known as BT426. Our two-class prediction method achieves a performance of: MCC=0.50, Qtotal=82.1%, sensitivity=75.6%, PPV=68.8% and AUC=0.864. We have compared our performance to eleven other prediction methods that obtain Matthews correlation coefficients in the range of 0.17-0.47. For the type specific β-turn predictions, only type I and II can be predicted with reasonable Matthews correlation coefficients, where we obtain performance values of 0.36 and 0.31, respectively. CONCLUSION: The NetTurnP method has been implemented as a webserver, which is freely available at http://www.cbs.dtu.dk/services/NetTurnP/. NetTurnP is the only available webserver that allows submission of multiple sequences.     

Molecular Dynamics Simulations of β-turn Forming Tetra- and Hexapeptides

Abstract

It was previously shown that the structural ensemble of model peptides DDKG and GKDG (H. Ishii et al. Biopolymers 24, 2045–2056, 1985), DEKS (A. Otter et al. J. Biomol. Struct. Dyn. 7, 455–476, 1989) NPGQ (F. R. Carbone et al. Int. J. Pept. Protein. Res. 26, 498–508, 1985), SALN (H. Santa et al. J. Biomol. Struct. Dyn. 16, 1033–1041, 1999), SYPFDV and SYPYDV (J. Yao et al. J. Mol. Biol. 243, 736–753, 1994), VP^DAH and VP^DSH (B. Imperiali et al. J. Am. Chem. Soc. 114, 3182–3188, 1992) in solution contains a significant—or in some cases dominant—proportion of β-turn conformation. In this study, a protein database was searched for the above, unprotected sequences which incorporate only L-amino acid residues. Simulated annealing and 25 ns MD simulations of structures were also performed. The DSSP and STRIDE secondary structure-assigning algorithms and clustering were used to analyze trajectories and i, i+3 hydrogen bonds were also sought. The DSSP analysis showed a fluctuation between β-turn and random meander structure, although bend structures were not detected because of the insufficient length of peptide chains. This alternating trend was confirmed when the STRIDE algorithm was used to analyze trajectories, but STRIDE assigned more turn structures. The population of the strongest clusters was above 40% and the middle structures adopted β-turn structure for most sequences. These results are in good agreement with previous experimental results and support the idea of the ultra-marginal stability of turns in the absence of stabilizing long-range interactions of the neighboring segments of a polypeptide chain. However, interactions between the side-chains in tetrapeptides could also contribute to turn stability and result in unusual stability in some cases. Our observations suggest that such interactions are the consequence rather than the driving force of turn formation.     

Proline-containing β-turns in peptides and proteins. II. Physicochemical studies on tripeptides with the Pro-Gly sequence

Amino acids are known to differ in their individual preferences for each of the four positions of the β-turn conformation formed by tetrapeptide segments. Proline and glycine show relatively high preferences for positions 2 and 3, respectively, of the β-turn. Using tripeptides of the type N-acetyl-Pro-Gly-X-OH, where X = Gly, Ala, Leu, Ile, and Phe, we have sought to study the influence of the 4th residue X on the stability of the β-turn conformation in these tripeptides. Our nmr and CD results show that the β-turn stability is quite significantly governed by the nature of the amino acid residue at this position in the following order: Leu > Ala > Ile, Gly > Phe.