Analysis and prediction of helix-helix interactions in membrane channels and transporters

Membrane proteins span a large variety of different functions such as cell-surface receptors, redox proteins, ion channels, and transporters. Proteins with functional pores show different characteristics of helix-helix packing as other helical membrane proteins. We found that the helix-helix contacts of 13 nonhomologous high-resolution structures of membrane channels and transporters are mainly accomplished by weakly polar amino acids (G > S > T > F) that preferably create contacts every fourth residue, typical for right-handed helix crossings. There is a strong correlation between the now available biological hydrophobicity scale and the propensities of the weakly polar and hydrophobic residues to be buried at helix-helix interfaces or to be exposed to the lipids in membrane channels and transporters. The polar residues, however, make no major contribution towards the packing of their transmembrane helices, and are therefore subsumed to be primarily exposed to the polar milieu during the folding process. The contact formation of membrane channels and transporters is therefore ruled by the solubility of the residues, which we suppose to be the driving force for the assembly of their transmembrane helices. By contrast, in 14 nonhomologous high-resolution structures of other membrane protein coils, also large and polar amino acids (D > S > M > Q) create characteristic contacts every 3.5th residues, which is a signature for left-handed helix crossings. Accordingly, it seems that dependent on the function, different concepts of folding and stabilization are realized for helical membrane proteins. Using a sequence-based matrix prediction method these differences are exploited to improve the prediction of buried and exposed residues of transmembrane helices significantly. When the sequence motifs typical for membrane channels and transporters were applied for the prediction of helix-helix contacts the quality of prediction rises by 16% to an average value of 76%, compared to the same approach when only single amino acid positions are taken into account.     

Molecular dynamics simulation of the unfolding of individual bacteriorhodopsin helices in sodium dodecyl sulfate micelles

We report molecular dynamics simulations of the trends in the changes in secondary structure of the seven individual helices of bacteriorhodopsin when inserted into sodium dodecyl sulfate (SDS) micelles, and their dependence on the amino acid sequence. The results indicate that the partitioning of the helices in the micelles and their stability are dependent on the hydrophobicity of the transmembrane segments. Helices A, B, and E are stable and retain their initial secondary structure throughout the 100 ns simulation time. In contrast, helices C, D, F, and G show structural perturbations within the first 10 ns. The instabilities are localized near charged residues within the transmembrane segments. The overall structural instability of the helix is correlated with its partitioning to the surface of the micelle and its interaction with polar groups there. The in silico experiments were performed to complement the in vitro experiments that examined the partial denaturation of bacteriorhodopsin in SDS described in the preceding article (DOI 10.1021/bi201769z ). The simulations are consistent with the trends revealed by the experimental results but strongly underestimate the extent of helix to extended coil transformation. The reason may be either that the sampling time was not sufficiently long or, more interestingly, that interhelix residue interactions play a role in the unfolding of the helices.     

The Family of Na+/Cl− Neurotransmitter Transporters

Abstract: The termination of neurotransmission is achieved by rapid uptake of the released neurotransmitter by specific high-affinity neurotransmitter transporters. Most of these transporters are encoded by a family of genes (Na⁺/Cl⁻ transporters) having a similar membrane topography of 12 transmembrane helices. An evolutionary tree revealed five distinct subfamilies: γ-aminobutyric acid transporters, monoamine transporters, amino acid transporters, "orphan" transporters, and the recently discovered bacterial transporters. The bacterial transporters that belong to this family may help to develop heterologous expression systems with the aim of solving the three-dimensional structure of these membrane proteins. Some of the neurotransmitter transporters have been implicated as important sites for drug action. Monoamine transporters, for example, are targeted by major classes of antidepressants, psychostimulants, and antihypertensive drugs. Localization of individual transporters in specific cells and brain areas is pertinent to understanding their contribution to neurotransmission and their potential as targets for drugs. The most important questions in the field include resolving the mechanism of neurotransmitter transport, the structure of the transporters, and the interaction of each transporter in complex neurological activities.     

Amino acid composition analysis of human secondary transport proteins and implications for reliable membrane topology prediction

Secondary transporters in humans are a large group of proteins that transport a wide range of ions, metals, organic and inorganic solutes involved in energy transduction, control of membrane potential and osmotic balance, metabolic processes and in the absorption or efflux of drugs and xenobiotics. They are also emerging as important targets for development of new drugs and as target sites for drug delivery to specific organs or tissues. We have performed amino acid composition (AAC) and phylogenetic analyses and membrane topology predictions for 336 human secondary transport proteins and used the results to confirm protein classification and to look for trends and correlations with structural domains and specific substrates and/or function. Some proteins showed statistically high contents of individual amino acids or of groups of amino acids with similar physicochemical properties. One recurring trend was a correlation between high contents of charged and/or polar residues with misleading results in predictions of membrane topology, which was especially prevalent in Mitochondrial Carrier family proteins. We demonstrate how charged or polar residues located in the middle of transmembrane helices can interfere with their identification by membrane topology tools resulting in missed helices in the prediction. Comparison of AAC in the human proteins with that in 235 secondary transport proteins from Escherichia coli revealed similar overall trends along with differences in average contents for some individual amino acids and groups of similar amino acids that are presumed to result from a greater number of functions and complexity in the higher organism.     

Comparison of helix interactions in membrane and soluble alpha-bundle proteins

Helix-helix interactions are important for the folding, stability, and function of membrane proteins. Here, two independent and complementary methods are used to investigate the nature and distribution of amino acids that mediate helix-helix interactions in membrane and soluble alpha-bundle proteins. The first method characterizes the packing density of individual amino acids in helical proteins based on the van der Waals surface area occluded by surrounding atoms. We have recently used this method to show that transmembrane helices pack more tightly, on average, than helices in soluble proteins. These studies are extended here to characterize the packing of interfacial and noninterfacial amino acids and the packing of amino acids in the interfaces of helices that have either right- or left-handed crossing angles, and either parallel or antiparallel orientations. We show that the most abundant tightly packed interfacial residues in membrane proteins are Gly, Ala, and Ser, and that helices with left-handed crossing angles are more tightly packed on average than helices with right-handed crossing angles. The second method used to characterize helix-helix interactions involves the use of helix contact plots. We find that helices in membrane proteins exhibit a broader distribution of interhelical contacts than helices in soluble proteins. Both helical membrane and soluble proteins make use of a general motif for helix interactions that relies mainly on four residues (Leu, Ala, Ile, Val) to mediate helix interactions in a fashion characteristic of left-handed helical coiled coils. However, a second motif for mediating helix interactions is revealed by the high occurrence and high average packing values of small and polar residues (Ala, Gly, Ser, Thr) in the helix interfaces of membrane proteins. Finally, we show that there is a strong linear correlation between the occurrence of residues in helix-helix interfaces and their packing values, and discuss these results with respect to membrane protein structure prediction and membrane protein stability.     

Genetic mapping of the human C5a receptor. Identification of transmembrane amino acids critical for receptor function

Many hormones and sensory stimuli signal through a superfamily of seven transmembrane-spanning receptors to activate heterotrimeric G proteins. How the seven transmembrane segments of the receptors (a molecular architecture of bundled alpha-helices conserved from yeast to man) work as "on/off" switches remains unknown. Previously, we used random saturation mutagenesis coupled with a genetic selection in yeast to determine the relative importance of amino acids in four of the seven transmembrane segments of the human C5a receptor (Baranski, T. J., Herzmark, P., Lichtarge, O., Gerber, B. O., Trueheart, J., Meng, E. C., Iiri, T., Sheikh, S. P., and Bourne, H. R. (1999) J. Biol. Chem. 274, 15757-15765). In this study, we evaluate helices I, II, and IV, thereby furnishing a complete mutational map of the seven transmembrane helices of the human C5a receptor. Our analysis identified 19 amino acid positions resistant to non-conservative substitutions. When combined with the 25 essential residues previously identified in helices III and V-VII, they delineate two distinct components of the receptor switch: a ligand-binding surface at or near the extracellular surface of the helix bundle and a core cluster in the cytoplasmic half of the bundle. In addition, we found critical amino acids in the first and second helices that are predicted to face the lipid membrane. These residues form an extended surface that might mediate interactions with lipids and other membrane proteins or function as an oligomerization domain with other receptors.     

Helix packing moments reveal diversity and conservation in membrane protein structure

Helical membrane proteins are more tightly packed and the packing interactions are more diverse than those found in helical soluble proteins. Based on a linear correlation between amino acid packing values and interhelical propensity, we propose the concept of a helix packing moment to predict the orientation of helices in helical membrane proteins and membrane protein complexes. We show that the helix packing moment correlates with the helix interfaces of helix dimers of single pass membrane proteins of known structure. Helix packing moments are also shown to help identify the packing interfaces in membrane proteins with multiple transmembrane helices, where a single helix can have multiple contact surfaces. Analyses are described on class A G protein-coupled receptors (GPCRs) with seven transmembrane helices. We show that the helix packing moments are conserved across the class A family of GPCRs and correspond to key structural contacts in rhodopsin. These contacts are distinct from the highly conserved signature motifs of GPCRs and have not previously been recognized. The specific amino acid types involved in these contacts, however, are not necessarily conserved between subfamilies of GPCRs, indicating that the same protein architecture can be supported by a diverse set of interactions. In GPCRs, as well as membrane channels and transporters, amino acid residues with small side-chains (Gly, Ala, Ser, Cys) allow tight helix packing by mediating strong van der Waals interactions between helices. Closely packed helices, in turn, facilitate interhelical hydrogen bonding of both weakly polar (Ser, Thr, Cys) and strongly polar (Asn, Gln, Glu, Asp, His, Arg, Lys) amino acid residues. We propose the use of the helix packing moment as a complementary tool to the helical hydrophobic moment in the analysis of transmembrane sequences.     

Molecular packing and packing defects in helical membrane proteins

The packing of helices spanning lipid bilayers is crucial for the stability and function of alpha-helical membrane proteins. Using a modified Voronoi procedure, we calculated packing densities for helix-helix contacts in membrane spanning domains. Our results show that the transmembrane helices of protein channels and transporters are significantly more loosely packed compared with helices in globular proteins. The observed packing deficiencies of these membrane proteins are also reflected by a higher amount of cavities at functionally important sites. The cavities positioned along the gated pores of membrane channels and transporters are noticeably lined by polar amino acids that should be exposed to the aqueous medium when the protein is in the open state. In contrast, nonpolar amino acids surround the cavities in those protein regions where large rearrangements are supposed to take place, as near the hinge regions of transporters or at restriction sites of protein channels. We presume that the observed deficiencies of helix-helix packing are essential for the helical mobility that sustains the function of many membrane protein channels and transporters.     

Membrane-integration Characteristics of Two ABC Transporters, CFTR and P-glycoprotein

To what extent do corresponding transmembrane helices in related integral membrane proteins have different membrane-insertion characteristics? Here, we compare, side-by-side, the membrane insertion characteristics of the 12 transmembrane helices in the adenosine triphosphate-binding cassette (ABC) transporters, P-glycoprotein (P-gp) and the cystic fibrosis transmembrane conductance regulator (CFTR). Our results show that 10 of the 12 CFTR transmembrane segments can insert independently into the ER membrane. In contrast, only three of the P-gp transmembrane segments are independently stable in the membrane, while the majority depend on the presence of neighboring loops and/or transmembrane segments for efficient insertion. Membrane-insertion characteristics can thus vary widely between related proteins.     

Comprehensive analysis of the numbers,lengths and amino acid compositions of transmembrane helices in prokaryotic,eukaryotic and viral integral membrane proteins of high-resolution structure

We report a comprehensive analysis of the numbers, lengths and amino acid compositions of transmembrane helices in 235 high-resolution structures of integral membrane proteins. The properties of 1551 transmembrane helices in the structures were compared with those obtained by analysis of the same amino acid sequences using topology prediction tools. Explanations for the 81 (5.2%) missing or additional transmembrane helices in the prediction results were identified. Main reasons for missing transmembrane helices were mis-identification of N-terminal signal peptides, breaks in α-helix conformation or charged residues in the middle of transmembrane helices and transmembrane helices with unusual amino acid composition. The main reason for additional transmembrane helices was mis-identification of amphipathic helices, extramembrane helices or hairpin re-entrant loops. Transmembrane helix length had an overall median of 24 residues and an average of 24.9 ± 7.0 residues and the most common length was 23 residues. The overall content of residues in transmembrane helices as a percentage of the full proteins had a median of 56.8% and an average of 55.7 ± 16.0%. Amino acid composition was analysed for the full proteins, transmembrane helices and extramembrane regions. Individual proteins or types of proteins with transmembrane helices containing extremes in contents of individual amino acids or combinations of amino acids with similar physicochemical properties were identified and linked to structure and/or function. In addition to overall median and average values, all results were analysed for proteins originating from different types of organism (prokaryotic, eukaryotic, viral) and for subgroups of receptors, channels, transporters and others.     

Helix kinks are equally prevalent in soluble and membrane proteins

Helix kinks are a common feature of α‐helical membrane proteins, but are thought to be rare in soluble proteins. In this study we find that kinks are a feature of long α‐helices in both soluble and membrane proteins, rather than just transmembrane α‐helices. The apparent rarity of kinks in soluble proteins is due to the relative infrequency of long helices (≥20 residues) in these proteins. We compare length‐matched sets of soluble and membrane helices, and find that the frequency of kinks, the role of Proline, the patterns of other amino acid around kinks (allowing for the expected differences in amino acid distributions between the two types of protein), and the effects of hydrogen bonds are the same for the two types of helices. In both types of protein, helices that contain Proline in the second and subsequent turns are very frequently kinked. However, there are a sizeable proportion of kinked helices that do not contain a Proline in either their sequence or sequence homolog. Moreover, we observe that in soluble proteins, kinked helices have a structural preference in that they typically point into the solvent. Proteins 2014; 82:1960–1970. © 2014 The Authors. Proteins: Structure, Function, and Bioinformatics Published by Wiley Periodicals, Inc.     

A simple method for predicting transmembrane alpha helices with better accuracy.

The prediction of a protein's structure from its amino acid sequence has been a long-standing goal of molecular biology. In this work, a new set of conformational parameters for membrane spanning alpha helices was developed using the information from the topology of 70 membrane proteins. Based on these conformational parameters, a simple algorithm has been formulated to predict the transmembrane alpha helices in membrane proteins. A FORTRAN program has been developed which takes the amino acid sequence as input and gives the predicted transmembrane alpha-helices as output. The present method correctly identifies 295 transmembrane helical segments in 70 membrane proteins with only two overpredictions. Furthermore, this method predicts all 45 transmembrane helices in the photosynthetic reaction center, bacteriorhodopsin and cytochrome c oxidase to an 86% level of accuracy and so is better than all other methods published to date.     

A conserved amino acid motif (R-X-G-R-R) in the Glut1 glucose transporter is an important determinant of membrane topology.

The Glut1 glucose transporter is one of over 300 members of the major facilitator superfamily of membrane transporters. These proteins are extremely diverse in substrate specificity and differ in their transport mechanisms. The two most common features shared by many members of this superfamily are the presence of 12 predicted transmembrane segments and an amino acid motif, R-X-G-R-R, present at equivalent positions within the cytoplasmic loops joining transmembrane segments 2-3 and 8-9. The structural and functional roles of the arginine residues within these motifs in Glut1 were investigated by expression of site-directed mutant transporters in Xenopus oocytes followed by analyses of intrinsic transport activity and the membrane topology of mutant glycosylation-scanning reporter Glut1 molecules. Substitution of lysine residues for the cluster of 3 arginine residues in each of the 2 cytoplasmic pentameric motifs of Glut1 revealed no absolute requirement for arginine side chains at any of the 6 positions for transport of 2-deoxyglucose. However, removal of the 3 positive charges at either site by substitution of glycines for the arginines completely abolished transport activity as the result of a local perturbation in the membrane topology in which the cytoplasmic loop was aberrantly translocated into the exoplasm along with the two flanking transmembrane segments. Substitution of lysines for the arginines had no affect on membrane topology. We conclude that the positive charges in the R-X-G-R-R motif form critical local cytoplasmic anchor points involved in determining the membrane topology of Glut1. These data provide a simple explanation for the presence of this conserved amino acid motif in hundreds of functionally diverse membrane transporters that share a common predicted membrane topology.     

蛋白质折叠速率与其氨基酸序列极性的关系

以大肠杆菌60个蛋白酶以及几种常见病毒（SARS病毒、艾滋病病毒、丙型肝炎病毒及乙型肝炎病毒）各蛋白质序列中的所有α-螺旋和β-折叠片段为研究对象,计算了各片段的折叠速率和平均极性,分别在各物种的α-螺旋和β-折叠两类二级结构片段中分析了二者的相关性。得到结论：不论是大肠杆菌中的蛋白酶还是病毒蛋白,其中的两类氨基酸片段的平均极性与折叠速率都是极显著相关的：对于所有的α片段,二者呈线性正相关,而对于所有的β片段,二者成线性负相关。结果证实了在蛋白质折叠中,氨基酸的极性起着重要的作用。     

Gene sequence and predicted amino acid sequence of the motA protein, a membrane-associated protein required for flagellar rotation in Escherichia coli.

The motA and motB gene products of Escherichia coli are integral membrane proteins necessary for flagellar rotation. We determined the DNA sequence of the region containing the motA gene and its promoter. Within this sequence, there is an open reading frame of 885 nucleotides, which with high probability (98% confidence level) meets criteria for a coding sequence. The 295-residue amino acid translation product had a molecular weight of 31,974, in good agreement with the value determined experimentally by gel electrophoresis. The amino acid sequence, which was quite hydrophobic, was subjected to a theoretical analysis designed to predict membrane-spanning alpha-helical segments of integral membrane proteins; four such hydrophobic helices were predicted by this treatment. Additional amphipathic helices may also be present. A remarkable feature of the sequence is the existence of two segments of high uncompensated charge density, one positive and the other negative. Possible organization of the protein in the membrane is discussed. Asymmetry in the amino acid composition of translated DNA sequences was used to distinguish between two possible initiation codons. The use of this method as a criterion for authentication of coding regions is described briefly in an Appendix.     

Structural Features of the Glutamate Transporter Family

Neuronal and glial glutamate transporters remove the excitatory neurotransmitter glutamate from the synaptic cleft and thus prevent neurotoxicity. The proteins belong to a large and widespread family of secondary transporters, including bacterial glutamate, serine, and C₄-dicarboxylate transporters; mammalian neutral-amino-acid transporters; and an increasing number of bacterial, archaeal, and eukaryotic proteins that have not yet been functionally characterized. Sixty members of the glutamate transporter family were found in the databases on the basis of sequence homology. The amino acid sequences of the carriers have diverged enormously. Homology between the members of the family is most apparent in a stretch of approximately 150 residues in the C-terminal part of the proteins. This region contains four reasonably well-conserved sequence motifs, all of which have been suggested to be part of the translocation pore or substrate binding site. Phylogenetic analysis of the C-terminal stretch revealed the presence of five subfamilies with characterized members: (i) the eukaryotic glutamate transporters, (ii) the bacterial glutamate transporters, (iii) the eukaryotic neutral-amino-acid transporters, (iv) the bacterial C₄-dicarboxylate transporters, and (v) the bacterial serine transporters. A number of other subfamilies that do not contain characterized members have been defined. In contrast to their amino acid sequences, the hydropathy profiles of the members of the family are extremely well conserved. Analysis of the hydropathy profiles has suggested that the glutamate transporters have a global structure that is unique among secondary transporters. Experimentally, the unique structure of the transporters was recently confirmed by membrane topology studies. Although there is still controversy about part of the topology, the most likely model predicts the presence of eight membrane-spanning α-helices and a loop-pore structure which is unique among secondary transporters but may resemble loop-pores found in ion channels. A second distinctive structural feature is the presence of a highly amphipathic membrane-spanning helix that provides a hydrophilic path through the membrane. Recent data from analysis of site-directed mutants and studies on the mechanism and pharmacology of the transporters are discussed in relation to the structural model.     

Conformational preference functions for predicting helices in membrane proteins

A suite of FORTRAN programs, PREF, is described for calculating preference functions from the data base of known protein structures and for comparing smoothed profiles of sequence-dependent preferences in proteins of unknown structure. Amino acid preferences for a secondary structure are considered as functions of a sequence environment. Sequence environment of amino acid residue in a protein is defined as an average over some physical, chemical, or statistical property of its primary structure neighbors. The frequency distribution of sequence environments in the data base of soluble protein structures is approximately normal for each amino acid type of known secondary conformation. An analytical expression for the dependence of preferences on sequence environment is obtained after each frequency distribution is replaced by corresponding Gaussian function. The preference for the α-helical conformation increases for each amino acid type with the increase of sequence environment of buried solvent-accessible surface areas. We show that a set of preference functions based on buried surface area is useful for predicting folding motifs in α-class proteins and in integral membrane proteins. The prediction accuracy for helical residues is 79% for 5 integral membrane proteins and 74% for 11 α-class soluble proteins. Most residues found in transmembrane segments of membrane proteins with known α-helical structure are predicted to be indeed in the helical conformation because of very high middle helix preferences. Both extramembrane and transmembrane helices in the photosynthetic reaction center M and L subunits are correctly predicted. We point out in the discussion that our method of conformational preference functions can identify what physical properties of the amino acids are important in the formation of particular secondary structure elements. © 1993 John Wiley & Sons, Inc.     

Structural features of the glutamate transporter family.

Neuronal and glial glutamate transporters remove the excitatory neurotransmitter glutamate from the synaptic cleft and thus prevent neurotoxicity. The proteins belong to a large and widespread family of secondary transporters, including bacterial glutamate, serine, and C4-dicarboxylate transporters; mammalian neutral-amino-acid transporters; and an increasing number of bacterial, archaeal, and eukaryotic proteins that have not yet been functionally characterized. Sixty members of the glutamate transporter family were found in the databases on the basis of sequence homology. The amino acid sequences of the carriers have diverged enormously. Homology between the members of the family is most apparent in a stretch of approximately 150 residues in the C-terminal part of the proteins. This region contains four reasonably well-conserved sequence motifs, all of which have been suggested to be part of the translocation pore or substrate binding site. Phylogenetic analysis of the C-terminal stretch revealed the presence of five subfamilies with characterized members: (i) the eukaryotic glutamate transporters, (ii) the bacterial glutamate transporters, (iii) the eukaryotic neutral-amino-acid transporters, (iv) the bacterial C4-dicarboxylate transporters, and (v) the bacterial serine transporters. A number of other subfamilies that do not contain characterized members have been defined. In contrast to their amino acid sequences, the hydropathy profiles of the members of the family are extremely well conserved. Analysis of the hydropathy profiles has suggested that the glutamate transporters have a global structure that is unique among secondary transporters. Experimentally, the unique structure of the transporters was recently confirmed by membrane topology studies. Although there is still controversy about part of the topology, the most likely model predicts the presence of eight membrane-spanning alpha-helices and a loop-pore structure which is unique among secondary transporters but may resemble loop-pores found in ion channels. A second distinctive structural feature is the presence of a highly amphipathic membrane-spanning helix that provides a hydrophilic path through the membrane. Recent data from analysis of site-directed mutants and studies on the mechanism and pharmacology of the transporters are discussed in relation to the structural model.     

Characteristics affecting expression and solubilization of yeast membrane proteins

Biochemical and structural analysis of membrane proteins often critically depends on the ability to overexpress and solubilize them. To identify properties of eukaryotic membrane proteins that may be predictive of successful overexpression, we analyzed expression levels of the genomic complement of over 1000 predicted membrane proteins in a recently completed Saccharomyces cerevisiae protein expression library. We detected statistically significant positive and negative correlations between high membrane protein expression and protein properties such as size, overall hydrophobicity, number of transmembrane helices, and amino acid composition of transmembrane segments. Although expression levels of membrane and soluble proteins exhibited similar negative correlations with overall hydrophobicity, high-level membrane protein expression was positively correlated with the hydrophobicity of predicted transmembrane segments. To further characterize yeast membrane proteins as potential targets for structure determination, we tested the solubility of 122 of the highest expressed yeast membrane proteins in six commonly used detergents. Almost all the proteins tested could be solubilized using a small number of detergents. Solubility in some detergents depended on protein size, number of transmembrane segments, and hydrophobicity of predicted transmembrane segments. These results suggest that bioinformatic approaches may be capable of identifying membrane proteins that are most amenable to overexpression and detergent solubilization for structural and biochemical analyses. Bioinformatic approaches could also be used in the redesign of proteins that are not intrinsically well-adapted to such studies.     

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.