Stabilizing nonpolar/polar side-chain interactions in the alpha-helix.

A simplistic, yet often used, view of protein stability is that amino acids attract other amino acids with similar polarity, whereas nonpolar and polar side chains repel. Here we show that nonpolar/polar interactions, namely Val or Ile bonding to Lys or Arg in alpha-helices, can in fact be stabilizing. Residues spaced i, i + 4 in alpha-helices are on the same face of the helix, with potential to favorably interact and stabilize the structure. We observe that the nonpolar/polar pairs Ile-Lys, Ile-Arg, and Val-Lys occur in protein helices more often than expected when spaced i, i + 4. Partially helical peptides containing pairs of nonpolar/polar residues were synthesized. Controls with i, i + 5 spacing have the residues on opposite faces of the helix and are less helical than the test peptides with the i, i + 4 interactions. Experimental circular dichroism results were analyzed with helix-coil theory to calculate the free energy for the interactions. All three stabilize the helix with DeltaG between -0.14 and -0.32 kcal x mol(-1). The interactions are hydrophobic with contacts between Val or Ile and the alkyl groups in Arg or Lys. Side chains such as Lys and Arg can thus interact favorably with both polar and nonpolar residues.     

Relationship between ion pair geometries and electrostatic strengths in proteins

The electrostatic free energy contribution of an ion pair in a protein depends on two factors, geometrical orientation of the side-chain charged groups with respect to each other and the structural context of the ion pair in the protein. Conformers in NMR ensembles enable studies of the relationship between geometry and electrostatic strengths of ion pairs, because the protein structural contexts are highly similar across different conformers. We have studied this relationship using a dataset of 22 unique ion pairs in 14 NMR conformer ensembles for 11 nonhomologous proteins. In different NMR conformers, the ion pairs are classified as salt bridges, nitrogen-oxygen (N-O) bridges and longer-range ion pairs on the basis of geometrical criteria. In salt bridges, centroids of the side-chain charged groups and at least a pair of side-chain nitrogen and oxygen atoms of the ion-pairing residues are within a 4 A distance. In N-O bridges, at least a pair of the side-chain nitrogen and oxygen atoms of the ion-pairing residues are within 4 A distance, but the distance between the side-chain charged group centroids is greater than 4 A. In the longer-range ion pairs, the side-chain charged group centroids as well as the side-chain nitrogen and oxygen atoms are more than 4 A apart. Continuum electrostatic calculations indicate that most of the ion pairs have stabilizing electrostatic contributions when their side-chain charged group centroids are within 5 A distance. Hence, most (approximately 92%) of the salt bridges and a majority (68%) of the N-O bridges are stabilizing. Most (approximately 89%) of the destabilizing ion pairs are the longer-range ion pairs. In the NMR conformer ensembles, the electrostatic interaction between side-chain charged groups of the ion-pairing residues is the strongest for salt bridges, considerably weaker for N-O bridges, and the weakest for longer-range ion pairs. These results suggest empirical rules for stabilizing electrostatic interactions in proteins.     

ICBS: a database of interactions between protein chains mediated by beta-sheet formation

MOTIVATION: Interchain beta-sheet (ICBS) interactions occur widely in protein quaternary structures, interactions between proteins and protein aggregation. These interactions play a central role in many biological processes and in diseases ranging from AIDS and cancer to anthrax and Alzheimer's. RESULTS: We have created a comprehensive database of ICBS interactions that is updated on a weekly basis and allows entries to be sorted and searched by relevance and other criteria through a simple Web interface. We derive a simple ICBS index to quantify the relative contributions of the beta-ladders in the overall interchain interaction and compute first- and second-order statistics regarding amino acid composition and pairing at different relative positions in the beta-strands. Analysis of the database reveals a 15.8% prevalence of significant ICBS interactions, the majority of which involve the formation of antiparallel beta-sheets and many of which involve the formation of dimers and oligomers. The frequencies of amino acids in ICBS interfaces are similar to those in intrachain beta-sheet interfaces. A full range of non-covalent interactions between side chains complement the hydrogen-bonding interactions between the main chains. Polar amino acids pair preferentially with polar amino acids and non-polar amino acids pair preferentially with non-polar amino acids among antiparallel (i, j) pairs. We anticipate that the statistics and insights gained from the database will guide the development of agents that control interchain beta-sheet interactions and that the database will help identify new protein interactions and targets for these agents. AVAILABILITY: The database is available at: http://www.igb.uci.edu/servers/icbs/     

Sequence specificity in the dimerization of transmembrane alpha-helices.

While several reports have suggested a role for helix-helix interactions in membrane protein oligomerization, there are few direct biochemical data bearing on this subject. Here, using mutational analysis, we show that dimerization of the transmembrane alpha-helix of glycophorin A in a detergent environment is spontaneous and highly specific. Very subtle changes in the side-chain structure at certain sensitive positions disrupt the helix-helix association. These sensitive positions occur at approximately every 3.9 residues along the helix, consistent with their comprising the interface of a closely fit transmembranous supercoil of alpha-helices. By contrast with other reported cases of interactions between transmembrane helices, the set of interfacial residues in this case contains no highly polar groups. Amino acids with aliphatic side chains define much of the interface, indicating that precise packing interactions between the helices may provide much of the energy for association. These data highlight the potential general importance of specific interactions between the hydrophobic anchors of integral membrane proteins.     

Correlation of sequence and tertiary structure in globular proteins

The x-ray crystal structures of 19 selected proteins are examined empirically for correlations between the amino acid sequence and long-range, tertiary conformation. There is clear evidence for preferential associations of certain types of amino acids, particularly among the hydrophobic aliphatic, aromatic, and cysteine residues. However, the likelihoods of forming these residue-pair contacts are all less than 12%, so packing and geometric requirements must often take precedent over energetic considerations. The prediction of long-range contacts is not substantially improved by taking into account the sequentially previous residues. The analysis of atom–atom contacts shows a similar lack of predictive ability, but the results show that a good approximation to the interresidue energy function must include different types of interactions at two or three different sites on some amino acids. Backbone–backbone long-range interactions are relatively rare and nonspecific, whereas some “polar” side chains form hydrogen bonds from the polar groups while occasionally forming hydrophobic contacts with the remainder of the chain.     

Protein contacts, inter-residue interactions and side-chain modelling

Three-dimensional structures of proteins are the support of their biological functions. Their folds are stabilized by contacts between residues. Inner protein contacts are generally described through direct atomic contacts, i.e. interactions between side-chain atoms, while contact prediction methods mainly used inter-Calpha distances. In this paper, we have analyzed the protein contacts on a recent high quality non-redundant databank using different criteria. First, we have studied the average number of contacts depending on the distance threshold to define a contact. Preferential contacts between types of amino acids have been highlighted. Detailed analyses have been done concerning the proximity of contacts in the sequence, the size of the proteins and fold classes. The strongest differences have been extracted, highlighting important residues. Then, we studied the influence of five different side-chain conformation prediction methods (SCWRL, IRECS, SCAP, SCATD and SCCOMP) on the distribution of contacts. The prediction rates of these different methods are quite similar. However, using a distance criterion between side chains, the results are quite different, e.g. SCAP predicts 50% more contacts than observed, unlike other methods that predict fewer contacts than observed. Contacts deduced are quite distinct from one method to another with at most 75% contacts in common. Moreover, distributions of amino acid preferential contacts present unexpected behaviours distinct from previously observed in the X-ray structures, especially at the surface of proteins. For instance, the interactions involving Tryptophan greatly decrease.     

Analysis of side-chain orientations in homologous proteins

The side-chain conformations of topologically equivalent residues in seven pairs of proteins ranging in sequence homology from 16% to 60% are compared. Both identical and mutated residues are included. For proteins with greater than 40% homology, it is found that at least 80% of the side-chain orientations of identical residues and 75% or more of the mutated residues in each pair of proteins have matching gamma atom dihedral angles (+/- 40 degrees); the comparison is not based strictly on chi 1 angles. Further, if a match is obtained at the gamma position, there is a high probability of matching for the delta atom(s) of the side-chain. For proteins with less than 25% homology the percentages are somewhat lower. Trends observed for conservative substitutions are essentially the same as those noted for mutated residues in general. Side-chain accessibility does not affect the probability of matches of identical residues; however, less accessible pairs of mutated residues have 10 to 20% higher matching probabilities than do exposed residues. Mismatches can frequently be related to large B-factors, certain types of amino acid substitutions, or the appearance of multiple minima on the side-chain potential energy surfaces and are most likely to occur for certain small residues (Ser, Thr, Val). Analysis of all the results makes possible the formulation of a set of rules for side-chain positioning in the modeling of homologous proteins.     

Hydrophilicity of polar amino acid side-chains is markedly reduced by flanking peptide bonds

Several amino acid side-chain hydropathy scales have been devised on the basis of solubility and water/organic solvent partitioning data obtained with free amino acids or side-chain analogs. In nearly all cases, these scales are based upon the structure-additivity assumption; it has been assumed that the transfer free energies of the amino acid side-chains are the same in these model compounds as they are in a polypeptide. This assumption is probably wrong. In the present study, deviations from additivity for amino acid side-chains are demonstrated by comparing a theoretically derived scale, which N-acetylamino acid amides. The results show that the flanking peptide bonds dramatically reduce the hydrophilicity of the polar side-chains, with deviations up to several kilocalories (1 kcal = 4.184 kJ) for the charged side-chains at pH 7.0. Further calculation shows that these deviations are due to reductions of 40 to 85% in the unfavorable transfer free energy of the polar functional groups. In addition, proximity of the neighboring amide bonds in the parent molecule (N-acetylglycine amide) decreases the hydrophilicity of the -CONH-backbone unit by 36%. This decrease is expected to be twice as large for -CONH- units in the interior of a polypeptide backbone. The significance of these observations is: (1) valid hydropathy scales can be obtained only with model peptides; (2) deviations from additivity are expected in all solvent systems, including non-polar solvents that are thought to mimic the interior of a membrane; (3) the spontaneous insertion of polypeptides into membranes is likely to occur much more readily than has been previously thought. In order to estimate the free energy of transferring the side-chains and the polypeptide backbone from water to the interior of a lipid bilayer, the results of this study are used to construct a hydropathy scale based upon the partitioning of solutes between water and non-polar solvents. The validity of hydropathy scales that are based on criteria other than solubility and water/organic solvent partitioning data is also discussed.     

Non-canonical interactions of porphyrins in porphyrin-containing proteins

In this study we have described the non-canonical interactions between the porphyrin ring and the protein part of porphyrin-containing proteins to better understand their stabilizing role. The analysis reported in this study shows that the predominant type of non-canonical interactions at porphyrins are CH···O and CH···N interactions, with a small percentage of CH···π and non-canonical interactions involving sulfur atoms. The majority of non-canonical interactions are formed from side-chains of charged and polar amino acids, whereas backbone groups are not frequently involved. The main-chain non-canonical interactions might be slightly more linear than the side-chain interactions, and they have somewhat shorter median distances. The analysis, performed in this study, shows that about 44% of the total interactions in the dataset are involved in the formation of multiple (furcated) non-canonical interactions. The high number of porphyrin–water interactions show importance of the inclusion of solvent in protein–ligand interaction studies. Furthermore, in the present study we have observed that stabilization centers are composed predominantly from nonpolar amino acid residues. Amino acids deployed in the environment of porphyrin rings are deposited in helices and coils. The results from this study might be used for structure-based porphyrin protein prediction and as scaffolds for future porphyrin-containing protein design.     

TRANSPORT OF AMINO ACIDS BY THE SOIL ALGA STICHOCOCCUS BACILLARIS1

The green alga Stichococcus bacillaris Naeg. is able to take up at least eleven amino acids. All of these except glutamic and aspartic acids are transported by carrier systems that obey saturation kinetics. The acidic amino acids enter the cell by passive diffusion. Michaelis-Menten parameters (K_s and V_max) were calculated for several amino acids. All obey simple Michaelis-Menten behavior except for 2-methylalanine and leucine which may have double carrier systems of different affinities. Interactions between pairs of amino acids suggest that there is at least one carrier system specific for basic amino acids and probably several systems specific for neutral amino acids. Further analysis of neutral amino acid interactions reveal that the uptake of several amino acids is incompletely inhibited by competitor uptake at infinite concentration. The simplest interpretation of the data is the operation of three carrier systems for neutral amino acids, one of which has higher affinity and broader specificity than the other two. The amino acid carrier systems appear to operate by an active mechanism. The metabolic poison DCCD inhibits uptake up to 99%. The capacities of the neutral amino acid carrier systems are increased when cells are grown in medium containing suboptimal concentrations of nitrogen.     

Derivation and testing of pair potentials for protein folding. When is the quasichemical approximation correct?

Many existing derivations of knowledge-based statistical pair potentials invoke the quasichemical approximation to estimate the expected side-chain contact frequency if there were no amino acid pair-specific interactions. At first glance, the quasichemical approximation that treats the residues in a protein as being disconnected and expresses the side-chain contact probability as being proportional to the product of the mole fractions of the pair of residues would appear to be rather severe. To investigate the validity of this approximation, we introduce two new reference states in which no specific pair interactions between amino acids are allowed, but in which the connectivity of the protein chain is retained. The first estimates the expected number of side-chain contracts by treating the protein as a Gaussian random coil polymer. The second, more realistic reference state includes the effects of chain connectivity, secondary structure, and chain compactness by estimating the expected side-chain contrast probability by placing the sequence of interest in each member of a library of structures of comparable compactness to the native conformation. The side-chain contact maps are not allowed to readjust to the sequence of interest, i.e., the side chains cannot repack. This situation would hold rigorously if all amino acids were the same size. Both reference states effectively permit the factorization of the side-chain contact probability into sequence-dependent and structure-dependent terms. Then, because the sequence distribution of amino acids in proteins is random, the quasichemical approximation to each of these reference states is shown to be excellent. Thus, the range of validity of the quasichemical approximation is determined by the magnitude of the side-chain repacking term, which is, at present, unknown. Finally, the performance of these two sets of pair interaction potentials as well as side-chain contact fraction-based interaction scales is assessed by inverse folding tests both without and with allowing for gaps.     

Analysis of the code relating sequence to conformation in globular proteins. The distribution of residue pairs in turns and kinks in the backbone chain

1. The residue pair is considered as the fundamental unit which differentiates alpha-helix, beta-pleated sheet and the various turns and kink structures of the protein backbone. 2. The HPLG alphabet (Robson & Pain, 1974) is used to group pairs of residues, giving 16 possible conformational pairs, all of which are found with differing frequencies in the nine proteins examined. 3. The frequencies of occurrence of the 16 different types of turn or kink are analysed in relation to the constituent amino acids. Those containing the L or G conformation are of low frequency and are grouped for purposes of this analysis. 4. The distribution of amino acids within all the conformational pairs is non-random, with distinct preferences shown by certain residues. 5. All pairs containing an L or G conformation require the presence of a glycine or a proton-donor side chain. 6. The results are discussed in terms of the determination of these ;random' structures by local interactions.     

Extent and nature of contacts between protein molecules in crystal lattices and between subunits of protein oligomers

A survey was compiled of several characteristics of the intersubunit contacts in 58 oligomeric proteins, and of the intermolecular contacts in the lattice for 223 protein crystal structures. The total number of atoms in contact and the secondary structure elements involved are similar in the two types of interfaces. Crystal contact patches are frequently smaller than patches involved in oligomer interfaces. Crystal contacts result from more numerous interactions by polar residues, compared with a tendency toward nonpolar amino acids at oligomer interfaces. Arginine is the only amino acid prominent in both types of interfaces. Potentials of mean force for residue–residue contacts at both crystal and oligomer interfaces were derived from comparison of the number of observed residue–residue interactions with the number expected by mass action. They show that hydrophobic interactions at oligomer interfaces favor aromatic amino acids and methionine over aliphatic amino acids; and that crystal contacts form in such a way as to avoid inclusion of hydrophobic interactions. They also suggest that complex salt bridges with certain amino acid compositions might be important in oligomer formation. For a protein that is recalcitrant to crystallization, substitution of lysine residues with arginine or glutamine is a recommended strategy. Proteins 28:494–514, 1997. © 1997 Wiley-Liss, Inc.     

Specificity in the genetic code. The role of nucleotide base-amino acid interaction

The mechanism of the unique and specific association of a given amino acid to its t-RNA is investigated by theoretical methods. Several possible schemes are proposed to explain specificity. The physical forces which act within these mechanisms are illustrated by the computer simulation of probable interactions between glycine and nucleotide bases and base pairs. It is demonstrated that glycine has direct and selective affinities for the nucleotide bases and that these interactions are principally determined by the polar groups. Energies have been calculated for the interaction of glycine with several base pairs. From these, the possibility that specificity arises through direct complexing of an amino acid with its anticodon is evaluated.     

Helix propagation and N-cap propensities of the amino acids measured in alanine-based peptides in 40 volume percent trifluoroethanol.

The helix propagation and N-cap propensities of the amino acids have been measured in alanine-based peptides in 40 volume percent trifluoroethanol (40% TFE) to determine if this helix-stabilizing solvent uniformly affects all amino acids. The propensities in 40% TFE are compared with revised values of the helix parameters of alanine-based peptides in water. Revision of the propensities in water is the result of redefining the capping statistical weights and evaluating the helix nucleation constant with N-capping explicitly included in the helix-coil model. The propagation propensities of all amino acids increase in 40% TFE relative to water, but the increases are highly variable. In water, all beta-branched and beta-substituted amino acids are helix breakers. In 40% TFE, the propagation propensities of the nonpolar amino acids increase greatly, leaving charged and neutral polar, beta-substituted amino acids as helix breakers. Glycine and proline are strong helix breakers in both solvents. Free energy differences for helix propagation (delta delta G) between alanine and other nonpolar amino acids are twice as large in water as predicted from side-chain conformational entropies, but delta delta G values in 40% TFE are close to those predicted from side-chain entropies. This dependence of delta delta G on the solvent points to a specific role of water in determining the relative helix propensities of the nonpolar amino acids. The N-cap propensities converge toward a common value in 40% TFE, suggesting that differential solvation by water contributes to the diversity of N-cap values shown by the amino acids.     

Amino acid side-chain partition energies and distribution of residues in soluble proteins

Energies required to transfer amino acid side chains from water to less polar environments were calculated from results of several studies and compared with several statistical analyses of residue distributions in soluble proteins. An analysis that divides proteins into layers parallel with their surfaces is more informative than those that simply classify residues as exposed or buried. Most residues appear to be distributed as a function of the distance from the protein-water interface in a manner consistent with partition energies calculated from partitioning of amino acids between water and octanol phases and from solubilities of amino acids in water, ethanol, and methanol. Lys, Arg, Tyr, and Trp residues tend to concentrate near the water-protein interface where their apolar side-chain components are more buried than their polar side-chain components. Residue distributions calculated in this manner do not correlate well with side-chain solvation energies calculated from vapor pressures of side-chain analogs over a water phase. Results of statistical studies that classify residues as exposed to solvent or buried inside the protein interior appear to depend on the method used to classify residues. Data from some of these studies correlate better with solvation energies, but other data correlate better with partition energies. Most other statistical methods that have been used to evaluate effects of water on residue distributions yield results that correlate better with partition energies than with solvation energies.     

Nearest-neighbour interactions in species-rich shrublands: the roles of abundance, spatial patterns and resources

We explored pairwise nearest‐neighbour interactions in four species‐rich shrubland plant communities, asking the question: how often is an individual of species j the nearest‐neighbour of species i? In the observed data and null models, less than 35% of the maximum possible number of nearest‐neighbour species pairs was present, and at three of the four sites the number of observed nearest‐neighbour pairs were significantly less than those occurring in simulated null communities. Many of the missing pairs included woody shrubs whose absence might be interpreted as evidence of site‐specific competition between larger growth forms for soil resources or space. Less than 5% of the pairs of species that occurred did so at frequencies different from that expected under random mixing, and many of these pairs were conspecific. Of the heterospecific pairs whose frequency differed significantly from random mixing there was a weak bias towards pairs occurring at higher rather than lower frequencies than expected. There was no evidence for asymmetry (interaction of species j with species i but not the reverse) in the frequency of species pairs. Nearest‐neighbour relationships are species‐specific rather than between plant functional types. The four sites form a soil nutrient and water availability gradient, and, according to the stress gradient hypothesis, positive species interactions should be most prevalent at the most stressful sites. However, we found the opposite: the site with the highest availability of resources had both proportionally the most heterospecific pairs, and the most conspecific and heterospecific species pairs with frequencies departing significantly from that expected under random mixing.     

Statistical analysis of predominantly transient protein-protein interfaces

A non-redundant set of 170 protein-protein interfaces of known structure was statistically analyzed for residue and secondary-structure compositions, pairing preferences and side-chain-backbone interaction frequencies. By focussing mainly on transient protein-protein interfaces, the results underline previous findings for protein-protein interfaces but also show some new interesting aspects of transient interfaces. The residue compositions at interfaces found in this study correlate well with the results of other studies. On average, contacts between pairs of hydrophobic and polar residues were unfavorable, and the charged residues tended to pair subject to charge complementarity. Secondary structure composition analysis shows that neither helices nor beta-sheets are dominantly populated at interfaces. Analyzing the pairing preferences of the secondary structure elements revealed a higher affinity within the same elements and alludes to tight packings. In addition, the results for the side-chain and backbone interaction frequencies, which were measured under more stringent conditions, showed a high occurrence of side-chain-backbone interactions. Taking a closer look at the helix and beta-sheet binding frequencies for a given side-chain and backbone interaction underlined the relevance of tight packings. The polarity of interfaces increased with decreasing interface size. These types of information may be useful for scoring complexes in protein-protein docking studies or for prediction of protein-protein interfaces from the sequences alone.     

Excluded volume in protein side-chain packing

The excluded volume occupied by protein side-chains and the requirement of high packing density in the protein interior should severely limit the number of side-chain conformations compatible with a given native backbone. To examine the relationship between side-chain geometry and side-chain packing, we use an all-atom Monte Carlo simulation to sample the large space of side-chain conformations. We study three models of excluded volume and use umbrella sampling to effectively explore the entire space. We find that while excluded volume constraints reduce the size of conformational space by many orders of magnitude, the number of allowed conformations is still large. An average repacked conformation has 20 % of its chi angles in a non-native state, a marked reduction from the expected 67 % in the absence of excluded volume. Interestingly, well-packed conformations with up to 50 % non-native chi angles exist. The repacked conformations have native packing density as measured by a standard Voronoi procedure. Entropy is distributed non-uniformly over positions, and we partially explain the observed distribution using rotamer probabilities derived from the Protein Data Bank database. In several cases, native rotamers that occur infrequently in the database are seen with high probability in our simulation, indicating that sequence-specific excluded volume interactions can stabilize rotamers that are rare for a given backbone. In spite of our finding that 65 % of the native rotamers and 85 % of chi(1) angles can be predicted correctly on the basis of excluded volume only, 95 % of positions can accommodate more than one rotamer in simulation. We estimate that, in order to quench the side-chain entropy observed in the presence of excluded volume interactions, other interactions (hydrophobic, polar, electrostatic) must provide an additional stabilization of at least 0.6 kT per residue in order to single out the native state.     

Orientation of carp, Cyprinm carpio L., to free amino acids from Tubifex extract in an olfactometer

We analysed the effects of L-amino acid combinations, tested at the same relative concentrations as in 500 mgl¹ of a filtrated Tubifex extract, on attraction and exploratory behaviour of juvenile carp, Cyprinus curpio L., in an olfactometer. Tubifex filtrate was significantly active in the range of concentrations tested (5 mgl⁻¹, 50 mgl⁻¹, 500 mgl⁻¹, 5 gl⁻¹). Maximum attraction was obtained in response to a concentration of 500 mg l⁻¹. Soluble TCA fraction and a mixture of the 17 synthetic amino acids found in this fraction also showed significant effects on attraction and exploratory behaviour. Tests of the four chemical groups of amino acids gave the following results: acidic amino acids (aspartic and glutamic) did not produce significant activity. Basic amino acids (histidine, arginine and lysine) and polar, uncharged amino acids (glycine, serine, threonine, tyrosine, asparagine and glutamine) were ineffective as attractants but significantly increased exploration. Non-polar amino acids (alanine, valine, leucine, isoleucine, phenylalanine, and methionine) showed significant effects on both attraction and exploration. In tests where different pairs of these amino acid groups were mixed it was shown that a combination of non-polar amino acids and polar uncharged amino acids was the most effective in inducing both attraction and exploration. The simplest combination to have a significant effect on both attraction and exploration was alanine, valine and glycine.
These results are discussed in the light of recent data on the role of free amino acids in the behaviour of fishes.