Electronic Properties of the Amino Acid Side Chains Contribute to the Structural Preferences in Protein Folding

Abstract

A database of 118 non-redundant proteins was examined to determine the preferences of amino acids for secondary structures: α-helix, β-strand and coil conformations. To better understand how the physicochemical properties of amino acid side chains might influence protein folding, several new scales have been suggested for quantifying the electronic effects of amino acids. These include the pKa at the amino group, localized effect substituent constants (eσ), and a composite of these two scales (ε). Amino acids were also classified into 5 categories on the basis of their electronic properties: O (strong electron donor), U (weak donor), Z (ambivalent), B (weak electron acceptor), and X (strong acceptor). Certain categories of amino acid appeared to be critical for particular conformations, e.g., O and U-type residues for α-helix formation. Pairwise analysis of the database according to these categories revealed significant context effects in the structural preferences. In general, the propensity of an amino acid for a particular conformation was related to the electronic features of the side chain. Linear regression analyses revealed that the electronic properties of amino acids contributed about as much to the folding preferences as hydrophobicity, which is a well-established determinant of protein folding. A theoretical model has been proposed to explain how the electronic properties of the side chain groups might influence folding along the peptide backbone.     

Exceptional pairs of amino acid neighbors in alpha-helices

Amino acids seem to have specific preferences for various locations in alpha-helices. These specific preferences, called singlet local propensity (SLP), have been determined by calculating the preference of occurrence of each amino acid in different positions of the alpha-helix. We have studied the occurrence of amino acids, single or pairs, in different positions, singlet or doublet, of alpha-helices in a database of 343 non-homologous proteins representing a unique superfamily from the SCOP database with a resolution better than 2.5 A from the Protein Data Bank. The preference of single amino acids for various locations of the helix was shown by the relative entropy of each amino acid with respect to the background. Based on the total relative entropy of all amino acids occurring in a single position, the N(cap) position was found to be the most selective position in the alpha-helix. A rigorous statistical analysis of amino acid pair occurrences showed that there are exceptional pairs for which, the observed frequency of occurrence in various doublet positions of the alpha-helix is significantly different from the expected frequency of occurrence in that position. The doublet local propensity (DLP) was defined as the preference of occurrences of amino acid pairs in different doublet positions of the alpha-helix. For most amino acid pairs, the observed DLP (DLP(O)) was nearly equal to the expected DLP (DLP(E)), which is the product of the related SLPs. However, for exceptional pairs of amino acids identified above, the DLP(O) and DLP(E) values were significantly different. Based on the relative values of DLP(O) and DLP(E), exceptional amino acid pairs were divided into two categories. Those, for which the DLP(O) values are higher than DLP(E), should have a strong tendency to pair together in the specified position. For those pairs which the DLP(O) values are less than DLP(E), there exists a hindrance in neighboring of the two amino acids in that specific position of the alpha-helix. These cases have been identified and listed in various tables in this paper. The amount of mutual information carried by the exceptional pairs of amino acids was significantly higher than the average mutual information carried by other amino acid pairs. The average mutual information conveyed by amino acid pairs in each doublet position was found to be very small but non-zero.     

Peptide design using alpha,beta-dehydro amino acids: from beta-turns to helical hairpins

Incorporation of alpha,beta-dehydrophenylalanine (DeltaPhe) residue in peptides induces folded conformations: beta-turns in short peptides and 3(10)-helices in larger ones. A few exceptions-namely, alpha-helix or flat beta-bend ribbon structures-have also been reported in a few cases. The most favorable conformation of DeltaPhe residues are (phi,psi) approximately (-60 degrees, -30 degrees ), (-60 degrees, 150 degrees ), (80 degrees, 0 degrees ) or their enantiomers. DeltaPhe is an achiral and planar residue. These features have been exploited in designing DeltaPhe zippers and helix-turn-helix motifs. DeltaPhe can be incorporated in both right and left-handed helices. In fact, consecutive occurrence of three or more DeltaPhe amino acids induce left-handed screw sense in peptides containing L-amino acids. Weak interactions involving the DeltaPhe residue play an important role in molecular association. The C--H.O==C hydrogen bond between the DeltaPhe side-chain and backbone carboxyl moiety, pi-pi stacking interactions between DeltaPhe side chains belonging to enantiomeric helices have shown to stabilize folding. The unusual capability of a DeltaPhe ring to form the hub of multicentered interactions namely, a donor in aromatic C--H.pi and C--H.O==C and an acceptor in a CH(3).pi interaction suggests its exploitation in introducing long-range interactions in the folding of supersecondary structures.     

Conformation biases of amino acids based on tripeptide microenvironment from PDB database

We have constructed a bank (FTTP) of tendentious factors of three states of three-peptide units from PDB database based on conformational dihedral angle library and demonstrated that amino acid biases toward protein secondary structure are present in natural protein sequences. Our research results reveal that 20 standard amino acids fall into three groups: nine residues inclined to alpha-helix with a common character (e.g. direct side chain aliphatic residues or positive/negative charged residues) arrange in three grades, viz EA, QKRLD, and MN, in turn; seven residues are apt to beta-strand with 2'-branched side chain aliphatic residues or benzyl-included residues, namely PV, IYTC, and F, in three ranks; and four residues SHWG show a double tendency to both alpha and beta. Noticeably, proline has the strongest ability to form extended conformation, especially the Re value up to 9.5298 at position 3 (Table 3). Thus, biases of codons show an evident tendency in protein folding, where GC-rich codons are mainly in charge of forming contracted conformation, especially the codon's first letter plays a dominant role in translating the genomic GC signature into protein sequences and structures. So, biases of amino acids will play an important role in protein folding, folding codons, refining domain, structure prediction, and structural genomics/proteomics.     

Purification of gamma-glutamyltransferase by phenyl boronate affinity chromatography. Studies on the acceptor specificity of transpeptidation by rat kidney gamma-glutamyltransferase

gamma-Glutamyltransferase has been purified from rat kidney by a novel procedure using phenyl boronate affinity chromatography. The highly purified enzyme has been studied with respect to acceptor specificity for a number of amino acids, amino acid analogues, dipeptides and tripeptides. The acceptor activity is specific for L-amino acids. The amino acids and the majority of the essential amino acids are poor acceptors while the sulphur-containing amino acids are the best acceptors. The acceptor activity is modulated by the substitution of the amino acid side chain. Substitution of the side chain at the delta, gamma or beta positions results in a proportionally decreasing ability to act as acceptor. The carbonyl moiety of the gamma-carboxy group of the acceptor appears to be essential for acceptor activity, absence of an alpha-carboxy carbonyl group increases the Kappm of the acceptor approximately 100-fold.     

Nearest-neighbor effects and structural preferences in dipeptides are a function of the electronic properties of amino acid side-chains

The electronic properties of amino acid side-chains are emerging as an important factor in the preference for secondary structure in proteins. These properties have not been fully characterized, nor has their role in the behavior of peptides been explored in any detail. The present studies sought to evaluate several possibilities: 1) that hydrophilicity can be expressed solely in electronic terms, 2) that substituent effects of side-chains extend across the peptide bond, and (3) nearest-neighbor effects in dipeptides correlate with secondary structural preferences. Quantum mechanics (QM) calculations were used to define the electronic properties of individual amino acids and dipeptides. It was found that the hydrophilicity of an amino acid side-chain can be accurately represented as a function of the electron densities of its component atoms. In addition, the nature of an amino acid in the second position of a dipeptide affects the electronic properties (Mulliken populations and electron densities) of the main-chain atoms of the first residue. Certain electronic features of the dipeptides strongly correlated with propensity for secondary structure. Specifically, Mulliken population data at the Calpha atom and N atom predicted preference for alpha-helices versus coil and strand conformations, respectively. Analysis of dipeptides arrayed in either helical or extended structures revealed lengthening of main-chain bonds in the alpha-helical conformations. A thorough characterization of the electronic properties of amino acids and short peptide segments may provide a better understanding of the forces that determine secondary structure in proteins.     

Reconciling observations of sequence-specific conformational propensities with the generic polymeric behavior of denatured proteins

Radii of gyration of denatured proteins vary with chain length and are insensitive to details of amino acid sequence. Observations of sequence independence in polymeric properties conflict with results from spectroscopic experiments, which suggest the presence of sequence-specific residual structure in denatured states. Can we reconcile the two apparently conflicting sets of observations? To answer this question, we need knowledge of the ensemble of conformations accessible to proteins in good solvents. The excluded-volume limit provides an ideal mimic of polymers in good solvents. Therefore, we attempt to solve the "reconciliation problem" by simulating conformational ensembles accessible to peptides and proteins in the excluded-volume limit. Analysis of these ensembles for a variety of polypeptide sequences leads to results that are consistent with experimental observations of sequence-specific conformational preferences in short peptides and the scaling behavior of polymeric quantities for denatured proteins. Reconciliation in the excluded-volume limit comes about due to a tug of war between two factors, namely, minimization of steric overlap and the competing effects of conformational entropy. Minimization of steric overlap promotes chain stretching and leads to experimentally observed sequence-dependent preferences for locally extended segments such as polyproline II helices, beta-strands, and very short stretches of alpha-helix. Conformational entropy opposes chain stretching, and the calculated persistence length for sequence-dependent conformational preferences is less than five amino acids. This estimate does not vary with amino acid sequence. The short persistence lengths lead directly to experimental observations of generic sequence-independent behavior of radii of gyration for denatured proteins.     

Conformational preferences of amino acid side chains in collagen

Conformational energy computations were carried out on collagenlike triple-stranded conformations of several poly(tripeptide)s with the general structure CH₃CO? (Gly? X? Y)₃? NHCH₃. The sequences considered had various amino acid residues in position X or Y of the central tripeptide, with either Pro or Ala as a neighbor, i.e., Gly-X-Pro, Gly-X-Ala, Gly-Pro-Y, and Gly-Ala-Y. Minimum-energy conformations were computed for the side chains, and their distributions were compared for the four sequences. The residues used were Abu (= α-aminobutyric acid), Leu, Phe, Ser, Asp, Asn, Val, Ile, and Thr. The conformational energy of a ? Ch₂? CH₃ side chain in Abu was mapped as a function of the dihedral angle χ¹. Intrastrand interactions with neighboring residues do not affect the conformations of a side chain in position Y, and they have a minor effect on it in the X-Ala sequence, but they strongly restrict the conformational freedom of the side chain in the X-Pro sequence. Conversely, interstrand interactions do not affect side chains in position X, but they strongly restrict the conformational freedom of a side chain in position Y if there is a nearby Pro residue in a neighboring strand. Hydrogen bonds with the backbone can be formed in some conformations of long polar side chains, such as Asp, Asn, or Gln. All amino acid residues can be accommodated in collagen. Because of the interactions mentioned above, steric and energetic constraints can be correlated with observed preferences of certain amino acids for positions X or Y in collagen. Hence, these preferences may be explained, in part, in terms of differences in the conformational freedom of the side chains in the triple-stranded structure.     

Comparison of scales of amino acid side chain properties by conservation during evolution of four proteins

As the amino acid sequence of a given protein changes along the phylogenetic tree, enough of the overall folding pattern must be conserved to ensure that the protein still fulfils its biological function. Eighteen published scales which tabulate various side chain properties are compared here by computing the variance of each scale when applied to each of several protein families. The conservation of each scale of side chain properties is examined for the 20,627 residues in 60 mammalian myoglobins, 31 mammalian ribonucleases, insulin A and B chains (29 sequences each), 29 vertebrate and 28 plant cytochrome c's. Those scales which are the most highly conserved through the evolution of each protein family may well be the best predictors of protein folding patterns. The mean-area-buried scale and the optimized matching hydrophobicities scale are more conserved than other scales. An additional result is the relatively poor conservation across evolution of the Chou-Fasman secondary structure predictors.     

Using quantum mechanics to improve estimates of amino acid side chain rotamer energies

Amino acid side chains adopt a discrete set of favorable conformations typically referred to as rotamers. The relative energies of rotamers partially determine which side chain conformations are more often observed in protein structures and accurate estimates of these energies are important for predicting protein structure and designing new proteins. Protein modelers typically calculate side chain rotamer energies by using molecular mechanics (MM) potentials or by converting rotamer probabilities from the protein database (PDB) into relative free energies. One limitation of the knowledge‐based energies is that rotamer preferences observed in the PDB can reflect internal side chain energies as well as longer‐range interactions with the rest of the protein. Here, we test an alternative approach for calculating rotamer energies. We use three different quantum mechanics (QM) methods (second order Møller‐Plesset (MP2), density functional theory (DFT) energy calculation using the B3LYP functional, and Hartree‐Fock) to calculate the energy of amino acid rotamers in a dipeptide model system, and then use these pre‐calculated values in side chain placement simulations. Energies were calculated for over 36,000 different conformations of leucine, isoleucine, and valine dipeptides with backbone torsion angles from the helical and strand regions of the Ramachandran plot. In a subset of cases these energies differ significantly from those calculated with standard molecular mechanics potentials or those derived from PDB statistics. We find that in these cases the energies from the QM methods result in more accurate placement of amino acid side chains in structure prediction tests. Proteins 2008. © 2007 Wiley‐Liss, Inc.     

A computer model to dynamically simulate protein folding: studies with crambin

The current work describes a simplified representation of protein structure with uses in the simulation of protein folding. The model assumes that a protein can be represented by a freely rotating rigid chain with a single atom approximating the effect of each side chain. Potentials describing the attraction or repulsion between different types of amino acids are determined directly from the distribution of amino acids in the database of known protein structures. The optimization technique of simulated annealing has been used to dynamically sample the conformations available to this simple model, allowing the protein to evolve from an extended, random coil into a compact globular structure. Many characteristics expected of true proteins, such as the sequence-dependent formation of secondary structure, the partitioning of hydrophobic residues, and specific disulfide pairing, are reproduced by the simulation, suggesting the model may accurately simulate the folding process.     

Molecular orbital studies of the conformation around the 5 position in angiotensin II

Analogs of angiotensin are more potent when the side chain at position 5 is branched rather than unbranched. We have performed molecular orbital calculations of conformational preferences of l-valine (branched side chain) and l-a-aminobutyric acid (unbranched side chain) as amino acid residues. The results of these calculations illustrate the differing conformational preferences of the two classes of amino acids and support a stereochemical role for position 5 of angiotensin.     

Structure, stability and folding of the alpha-helix

Pauling first described the alpha-helix nearly 50 years ago, yet new features of its structure continue to be discovered, using peptide model systems, site-directed mutagenesis, advances in theory, the expansion of the Protein Data Bank and new experimental techniques. Helical peptides in solution form a vast number of structures, including fully helical, fully coiled and partly helical. To interpret peptide results quantitatively it is essential to use a helix/coil model that includes the stabilities of all these conformations. Our models now include terms for helix interiors, capping, side-chain interactions, N-termini and 3(10)-helices. The first three amino acids in a helix (N1, N2 and N3) and the preceding N-cap are unique, as their amide NH groups do not participate in backbone hydrogen bonding. We surveyed their structures in proteins and measured their amino acid preferences. The results are predominantly rationalized by hydrogen bonding to the free NH groups. Stabilizing side-chain-side-chain energies, including hydrophobic interactions, hydrogen bonding and polar/non-polar interactions, were measured accurately in helical peptides. Helices in proteins show a preference for having approximately an integral number of turns so that their N- and C-caps lie on the same side. There are also strong periodic trends in the likelihood of terminating a helix with a Schellman or alpha L C-cap motif. The kinetics of alpha-helix folding have been studied with stopped-flow deep ultraviolet circular dichroism using synchrotron radiation as the light source; this gives a far superior signal-to-noise ratio than a conventional instrument. We find that poly(Glu), poly(Lys) and alanine-based peptides fold in milliseconds, with longer peptides showing a transient overshoot in helix content.     

Molecular insights into protein synthesis with proline residues

Proline is an amino acid with a unique cyclic structure that facilitates the folding of many proteins, but also impedes the rate of peptide bond formation by the ribosome. As a ribosome substrate, proline reacts markedly slower when compared with other amino acids both as a donor and as an acceptor of the nascent peptide. Furthermore, synthesis of peptides with consecutive proline residues triggers ribosome stalling. Here, we report crystal structures of the eukaryotic ribosome bound to analogs of mono‐ and diprolyl‐tRNAs. These structures provide a high‐resolution insight into unique properties of proline as a ribosome substrate. They show that the cyclic structure of proline residue prevents proline positioning in the amino acid binding pocket and affects the nascent peptide chain position in the ribosomal peptide exit tunnel. These observations extend current knowledge of the protein synthesis mechanism. They also revise an old dogma that amino acids bind the ribosomal active site in a uniform way by showing that proline has a binding mode distinct from other amino acids.     

From peptide to protein: comparative analysis of the substrate specificity of N-linked glycosylation in C. jejuni

The gram-negative bacterium Campylobacter jejuni was recently discovered to contain a general N-linked protein glycosylation pathway. Central to this pathway is PglB, a homologue of the Stt3p subunit of the eukaryotic oligosaccharyl transferase (OT), which is involved in the transfer of an oligosaccharide from a polyisoprenyl pyrophosphate carrier to the asparagine side chain of proteins within the conserved glycosylation sites D/E-X1-N-X2-S/T, where X1 and X2 can be any amino acids except proline. Using a library of peptide substrates and a quantitative radioactivity-based in vitro assay, we assessed the amino acids at each position of the consensus glycosylation sequence for their impact on glycosylation efficiency, whereby the sequence DQNAT was found to be the optimal acceptor substrate. In the context of a full-length folded protein, the differences between variations of the glycosylation sequences were found to be consistent with the trends observed from their peptidyl counterparts, though less dramatic because of additional influences. In addition to characterizing the acceptor preferences of PglB, we also assessed the selectivity toward the glycan donor. Interestingly, despite recent reports of relaxed selectivity toward the glycan donor, PglB was not found to be capable of utilizing glycosyl donors such as dolichyl-pyrophosphate-chitobiose, which is the minimum substrate for the eukaryotic OT process.     

Analysis on folding of misgurin using two-dimensional HP model

Misgurin is an antimicrobial peptide from the loach, while the hydrophobic-polar (HP) model is a way to study the folding conformations and native states in peptide and protein although several amino acids cannot be classified either hydrophobic or polar. Practically, the HP model requires extremely intensive computations, thus it has yet to be used widely. In this study, we use the two-dimensional HP model to analyze all possible folding conformations and native states of misgurin with conversion of natural amino acids according to the normalized amino acid hydrophobicity index as well as the shortest benchmark HP sequence. The results show that the conversion of misgurin into HP sequence with glycine as hydrophobic amino acid at pH 2 has 1212 folding conformations with the same native state of minimal energy -6; the conversion of glycine as polar amino acid at pH 2 has 13,386 folding conformations with three native states of minimal energy -5; the conversion of glycine as hydrophobic amino acid at pH 7 has 2538 folding conformations with three native states of minimal energy -5; and the conversion of glycine as polar amino acid at pH 7 has 12,852 folding conformations with three native states of minimal energy -4. Those native states can be ranked according to the normalized amino acid hydrophobicity index. The detailed discussions suggest two ways to modify misgurin.     

Substrate specificity of Escherichia coli peptidyltransferase at the donor site

The substrate specificity of $\text{E.}\text{coli}$ peptidyltransferase at the donor site was investigated by the “50S reaction”. Seventeen N-acetylated or unacetylated aminoacyl-tRNAs and dipeptidyl-tRNAs were used as the donor substrates and puromycin as the acceptor. Results indicated that the nature of amino acid side chain of the donor tRNA has a predominant effect on the reaction rate of peptidyltransferase. Amino acids or dipeptides with high hydrophobicity were transferred faster than those with low hydrophobicity. Amino acids with alkyl side chains are better donors than those with aromatic side chains. Substrates with C-terminal proline were transferred extremely slowly which can probably be attributed to its unusual α-imino structure in addition to its low hydrophobicity.     

Left-handed topology of super-secondary structure formed by aligned alpha-helix and beta-hairpin.

A novel super-secondary structure common for many non-homological proteins is considered. This folding pattern, consisting of adjacent along the chain alpha-helix and beta-hairpin, has an aligned packing. It is found that one of the two possible 'mirror-symmetrical' topologies is observed in proteins. The alpha-helix + beta-hairpin structures have a similar pattern of hydrophobic residues in their amino acid sequences. The remaining part of a molecule or a domain is almost always located on the same side of the considered folding pattern. These results can be used in the prediction of three-dimensional protein structure and protein design.     

Design, synthesis and biological evaluation of selective boron-containing thrombin inhibitors.

Based on the structural comparison of the S-1 pocket in different trypsin-like serine proteases, a series of Boc-D-trimethylsilylalanine-proline-boro-X pinanediol derivatives, with boro-X being different amino boronic acids, have been synthesised as inhibitors of thrombin. The influence of hydrogen donor/acceptor properties of different residues in the P-1 side chain of these inhibitors on the selectivity profile has been investigated. This study confirmed the structure-based working hypothesis: The hydrophobic/hydrophilic character of amino acid residues 190 and 213 in the neighbourhood of Asp 189 in the S-1 pocket of thrombin (Ala/Val), trypsin (Ser/Val) and plasmin (Ser/Thr) define the specificity for the interaction with different P-1 residues of the inhibitors. Many of the synthesised compounds demonstrate potent antithrombin activity with Boc-D-trimethylsilylalanine-proline-boro-methoxypropylglycine++ + pinanediol (9) being the most selective thrombin inhibitor of this series.     

Coumarin as a structural component of substrates and probes for serine and cysteine proteases

Coumarins represent well-established structures to introduce fluorescence into tool compounds for biochemical investigations. They are valued for their small size, chemical stability and accessibility as well as their tunable photochemical properties. As components of fluorophore/quencher pairs or FRET donor/acceptor pairs, coumarins have frequently been applied in substrate mapping approaches for serine and cysteine proteases. This review also focuses on the incorporation of coumarins into the side chain of amino acids and the exploitation of the resulting fluorescent amino acids for the positional profiling of protease substrates. The protease-inhibiting properties of certain coumarin derivatives and the utilization of coumarin moieties to assemble activity-based probes for serine and cysteine proteases are discussed as well.