Differential helix propensity of small apolar side chains studied by molecular dynamics simulations.

A series of oligoalanine molecules with single amino acid replacements in the middle of the chain has been studied by molecular dynamics simulations. Differences in stability of the alpha-helix (as free energies delta delta G degrees) were estimated for the following series of residues: alpha-aminoisobutyric acid, alanine, alpha-amino-n-butyric acid, valine, glycine, D-alanine, t-leucine (= alpha-amino-beta,beta-dimethyl-n- butyric acid), and proline, arranged here in decreasing order of helix-forming potential. (The results for proline and valine had been reported earlier.) No experimental results were available for alpha-amino-n-butyric acid, D-alanine, and t-leucine at the time these calculations were done. The values of delta delta G degrees, including the three predictions, are in striking agreement with recent experimental results. A combination of free dynamics, dynamics with forced conformational change, and dynamics with forced molecular replacement was used. Conformational distributions were calculated for the peptide backbone of the dipeptides and, where appropriate, for the side chains of the dipeptide and the alpha-helix. The results demonstrate an unexpected level of accuracy for the all-atom model used to represent atomic interactions in the simulations. The simulations permit a detailed analysis of different factors responsible for conformational preferences and differences in stability. These conclusions drawn from this analysis agree with accepted qualitative explanations and allow these explanations to be quantitated to an extent not heretofore possible.     

Proline in alpha-helix: stability and conformation studied by dynamics simulation

Free-energy simulations have been used to estimate the change in the conformational stability of short polyalanine alpha-helices when one of the alanines is replaced by a proline residue. For substituting proline in the middle of the helix the change in free energy of folding (delta delta G degrees) was calculated as 14 kJ/mol (3.4 kcal/mol), in excellent agreement with the one available experimental value. The helix containing proline was found to be strongly kinked; the free energy for reducing the angle of the kink from 40 degrees to 15 degrees was calculated, and found to be small. A tendency to alternate hydrogen bonding schemes was observed in the proline-containing helix. These observations for the oligopeptide agree well with the observation of a range of kink angles (18-35 degrees) and variety of hydrogen bonding schemes, in the rare instances where proline occurs in helices in globular proteins. For substituting proline at the N-terminus of the helix the change in free energy of folding (delta delta G degrees) was calculated as -4 kJ/mol in the first helical position (N1) and +6 kJ/mol in the second helical position (N2). The observed frequent occurrence of proline in position N1 in alpha-helices in proteins therefore has its origin in stability differences of secondary structure. The conclusion reached here that proline may be a better helix former in position N1 than (even) alanine, and thus be a helix initiator may be testable experimentally by measurements of fraction helical conformation of individual residues in oligopeptides of appropriate sequence. The relevance of these results in regards to the frequent occurrence of proline-containing helices in certain membrane proteins is discussed.     

Advantages of fine-grained side chain conformer libraries

We compare the modelling accuracy of two common rotamer libraries, the Dunbrack-Cohen and the 'Penultimate' rotamer libraries, with that of a novel library of discrete side chain conformations extracted from the Protein Data Bank. These side chain conformer libraries are extracted automatically from high-quality protein structures using stringent filters and maintain crystallographic bond lengths and angles. This contrasts with traditional rotamer libraries defined in terms of chi angles under the assumption of idealized covalent geometry. We demonstrate that side chain modelling onto native and near-native main chain conformations is significantly more successful with the conformer libraries than with the rotamer libraries when solely considering excluded-volume interactions. The rotamer libraries are inadequate to model side chains without atomic clashes on over 20% of targets if the backbone is held fixed in the native conformation. An algorithm is described for simultaneously modelling both main chain and side chain atoms during discrete ab initio sampling. The resulting models have equivalent root mean square deviations from the experimentally determined protein loops as models from backbone-only ensembles, indicating that all-atom modelling does not detract from the accuracy of conformational sampling.     

Constructing side chains on near-native main chains for ab initio protein structure prediction

Is there value in constructing side chains while searching protein conformational space during an ab initio simulation? If so, what is the most computationally efficient method for constructing these side chains? To answer these questions, four published approaches were used to construct side chain conformations on a range of near-native main chains generated by ab initio protein structure prediction methods. The accuracy of these approaches was compared with a naive approach that selects the most frequently observed rotamer for a given amino acid to construct side chains. An all-atom conditional probability discriminatory function is useful at selecting conformations with overall low all-atom root mean square deviation (r.m.s.d.) and the discrimination improves on sets that are closer to the native conformation. In addition, the naive approach performs as well as more sophisticated methods in terms of the percentage of chi(1) angles built accurately and the all-atom r. m.s.d., between the native and near-native conformations. The results suggest that the naive method would be extremely useful for fast and efficient side chain construction on vast numbers of conformations for ab initio prediction of protein structure.     

Free energy surface for rotamers of cis-enol malonaldehyde in aqueous solution studied by molecular dynamics calculations

Two-dimensional free energy surfaces for four rotamers of cis-enol malonaldehyde in water have been investigated by umbrella sampling molecular dynamics (MD) calculations. Biasing potential used in the umbrella sampling calculation was adopted to be the minus of conformational free energy preliminary obtained by the thermodynamic integration MD calculations for the rigid malonaldehyde whose stretching and bending were all fixed. The calculated free energy surface shows that, in water, a rotamer that has an intramolecular hydrogen bond is most stable among the rotamers. This is the same as that in vacuum, while order of relative stability of the other three rotamers is different in water and in vacuum. Inclusion of intramolecular vibrations changed the free energy surface little, i.e. at most 2.6 kJ/mol, which is much smaller than the solvation free energy. Free energy barriers from the most stable intramolecular hydrogen bonded rotamer to the others are lowered by hydration but they are still very high, >50 kJ/mol, such that the malonaldehyde molecule spends most of its time in water taking this conformation. Thus, reaction coordinate for intramolecular proton transfer reaction in water may be constructed assuming this rotamer.     

Designing of peptides with left handed helical structure by incorporating the unusual amino acids

The conformational behaviour of delta Ala has been investigated by quantum mechanical method PCILO in the model dipeptide Ac-delta Ala-NHMe and in the model tripeptides Ac-X-delta Ala-NHMe with X = Gly, Ala, Val, Leu, Abu and Phe and is found to be quite different. The computational results suggest that in the model tripeptides the most stable conformation corresponds to phi 1 = -30 degrees, psi 1 = 120 degrees and phi 2 = psi 2 = 30 degrees in which the > C = 0 of the acetyl group is involved in hydrogen bond formation with N-H of the amide group. Similar results were obtained for the conformational behaviour of D-Ala in Ac-D-Ala-NHMe and Ac-Ala-D-Ala-NHMe. The conformational behaviour of the amino acids delta Ala, D-Ala, Val and Aib in model tripeptides have been utilized in the designing of left handed helical peptides. It is shown that the peptide HCO-(Ala-D-Ala)3-NHMe can adopt both left and right handed helix whereas in the peptide Ac-(Ala-delta Ala)3-NHMe the lowest energy conformer is beta-bend ribbon structure. Left handed helical structure with phi = 30 degrees, psi = 60 degrees for D-Ala residues and phi = psi = 30 degrees for delta Ala is found to be more stable by 4 kcal mole-1 than the corresponding right handed helical structure for the peptide Ac-(D-Ala-delta Ala)3-NHMe. In both the peptides Ac-(Val-delta Ala)3-NHMe and Ac-(D-Val-delta Ala)3-NHMe the most stable conformer is the left handed helix. Comparisons of results for Ac-(Ala-delta Ala)3-NHMe and Ac(Val-delta Ala)3-NHMe and Ac-(D-Ala-delta Ala)3-NHMe and Ac-(D-Val-delta Ala)3-NHMe also reveal that the Val residues facilitate the population of 3(10) left handed helix over the other conformers. It is also shown that the conformational behaviour of Aib residue depends on the chirality of neighbouring amino acids, i.e. Ac-(Aib-Ala)3-NHMe adopts right handed helical structure whereas Ac-(Aib-D-Ala)3-NHMe is found to be in left handed helical structure.     

Energetics of the hairpin to mismatched duplex transition of d(GCCGCAGC) on NaCl solution

We report Potential of Mean Force studies to describe the relative thermodynamic stabilities of d(GCCGCAGC) in a mismatched duplex and a hairpin monomer conformation in NaCl solution. The PMF calculations are combined with previous molecular mechanics and normal mode analysis in order to estimate the role of different components of the free energy in determining the relative stability of the duplex and hairpin structures. The high entropy associated with the loop region and the lack of minor groove phosphate-phosphate interactions in the hairpin compete against the gain in enthalpic contribution to the free energy due to base pairing in the mismatched duplex. The combined free energy calculations show that the hairpin is the most stable conformation at low salt and that a hairpin to duplex transition takes place at approximately 0.47 M NaCl. In addition, we studied the hairpin to partially stacked single helical conformation equilibrium at low salt. We found a small variation in transition temperature in salt concentration, delta Tm/delta log10(cs) approximately 2-3 degrees K/decade, in contrast to the duplex to hairpin or duplex to partially stacked single helix transition where the transition temperature exhibited marked dependence on salt concentration. This is in qualitative agreement with experimental data. Based on the Potential of Mean Force free energy calculation, the order of relative stability of the three-conformations studied varies with salt concentration. We observed the following orders of stability: stacked single helix greater than hairpin greater than duplex for cs less than 0.77 M NaCl; single helix greater than duplex greater than hairpin for 0.77 less than Cs less than 2.1 M; and duplex greater than hairpin greater than single strand for cs greater than 2.1 M. From the calculated PMF free energy curves in the NaCl concentration range, 0.012 less than cs less than 5.0 M, we can assign upper and lower bounds for the non-ionic differences in free energy between the duplex, hairpin, and stacked single helical states (at standard conditions: cs = 1.0 M, T = 25 degrees C, and 1 M oligomer concentration). We found that for delta G duplex single helix = G duplex - 2 x G single helix less than -7.38 Kcal/mol, the single helix is the least stable state. For the duplex-to-hairpin free energy difference in the range, -1.87 less than delta G duplex-hairpin less than 0.03 Kcal/mol, there will always be a salt-induced hairpin-to-duplex transition for 0.01 less than cs less than 1.6 M NaCl. If delta G duplex-hairpin less than -1.87, the duplex is always more stable than the hairpin; and for delta G duplex-hairpin greater than Kcal/mol, the hairpin state is always more stable than the duplex, for all salt concentrations.     

IRECS: a new algorithm for the selection of most probable ensembles of side-chain conformations in protein models

We introduce a new algorithm, IRECS (Iterative REduction of Conformational Space), for identifying ensembles of most probable side-chain conformations for homology modeling. On the basis of a given rotamer library, IRECS ranks all side-chain rotamers of a protein according to the probability with which each side chain adopts the respective rotamer conformation. This ranking enables the user to select small rotamer sets that are most likely to contain a near-native rotamer for each side chain. IRECS can therefore act as a fast heuristic alternative to the Dead-End-Elimination algorithm (DEE). In contrast to DEE, IRECS allows for the selection of rotamer subsets of arbitrary size, thus being able to define structure ensembles for a protein. We show that the selection of more than one rotamer per side chain is generally meaningful, since the selected rotamers represent the conformational space of flexible side chains. A knowledge-based statistical potential ROTA was constructed for the IRECS algorithm. The potential was optimized to discriminate between side-chain conformations of native and rotameric decoys of protein structures. By restricting the number of rotamers per side chain to one, IRECS can optimize side chains for a single conformation model. The average accuracy of IRECS for the chi1 and chi1+2 dihedral angles amounts to 84.7% and 71.6%, respectively, using a 40 degrees cutoff. When we compared IRECS with SCWRL and SCAP, the performance of IRECS was comparable to that of both methods. IRECS and the ROTA potential are available for download from the URL http://irecs.bioinf.mpi-inf.mpg.de.     

Relation of the progestagenic activity and conformation of the 17beta-acetyl side chain in the D'6-pentarane series

The synthesis of 2' beta-methyl-16 alpha,17 alpha-cyclohexanoprogesterone and its MM2 conformational analysis have been performed. The acetyl side chain was shown to have an unusual conformation with the torsion angle C13-C17-C20-O20 being -32.1 degrees. This conformation is by 5.4 kJ.mol-1 more stable than the usual one with the torsion angle 130.3 degrees. 2' beta-Methyl-16 alpha,17 alpha-cyclohexanoprogesterone proved to be inactive as a progestogen (pregnancy maintenance and McPhail tests). The lack of the activity may be due to the additional methyl group in D'-ring causing a change of the conformation of the 17 beta-acetyl side chain, thus hindering the formation of the conformation necessary for binding to the progesterone receptor.     

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.     

Crystal structure and ligand-binding studies of a screened peptide complexed with streptavidin.

The thermodynamic binding parameters and crystal structure for streptavidin-peptide complexes where the peptide sequences were obtained by random screening methods are reported. The affinities between streptavidin and two heptapeptides were determined by titrating calorimetric methods [Phe-Ser-His-Pro-Gln-Asn-Thr, Ka = 7944 (+/- 224) M-1, delta G degrees = -5.32 (+/- 0.01) kcal/mol, and delta H degrees = -19.34 (+/- 0.48) kcal/mol; His-Asp-His-Pro-Gln-Asn-Leu, Ka = 3542 (+/- 146) M-1, delta G degrees = -4.84 (+/- 0.03) kcal/mol, and delta H degrees = -19.00 (+/- 0.64) kcal/mol]. The crystal structure of streptavidin complexed with one of these peptides has been determined at 2.0-A resolution. The peptide (Phe-Ser-His-Pro-Gln-Asn-Thr) binds in a turn conformation with the histidine, proline, and glutamine side chains oriented inward at the biotin-binding site. A water molecule is immobilized between the histidine and glutamine side chains of the peptide and an aspartic acid side chain of the protein. Although some of the residues that participate in binding biotin also interact with the screened peptide, the peptide adopts an alternate method of utilizing binding determinants in the biotin-binding site of streptavidin.     

Addition of side chain interactions to modified Lifson-Roig helix-coil theory: application to energetics of phenylalanine-methionine interactions.

We introduce here i, i + 3 and i, i + 4 side chain interactions into the modified Lifson-Roig helix-coil theory of Doig et al. (1994, Biochemistry 33:3396-3403). The helix/coil equilibrium is a function of initiation, propagation, capping, and side chain interaction parameters. If each of these parameters is known, the helix content of any isolated peptide can be predicted. The model considers every possible conformation of a peptide, is not limited to peptides with only a single helical segment, and has physically meaningful parameters. We apply the theory to measure the i, i + 4 interaction energies between Phe and Met side chains. Peptides with these residues spaced i, i + 4 are significantly more helical than controls where they are spaced i, i + 5. Application of the model yields delta G for the Phe-Met orientation to be -0.75 kcal.mol-1, whereas that for the Met-Phe orientation is -0.54 kcal.mol-1. These orientational preferences can be explained, in part, by rotamer preferences for the interacting side chains. We place Phe-Met i, i + 4 at the N-terminus, the C-terminus, and in the center of the host peptide. The model quantitatively predicts the observed helix contents using a single parameter for the side chain-side chain interaction energy. This result indicates that the model works well even when the interaction is at different locations in the helix.     

New tubular single-stranded helix of poly-L-amino acids suggested by molecular mechanics calculations: I. Homopolypeptides in isolated environments.

A search was made in terms of molecular mechanics calculation for tubular, or pore-forming, single-stranded helices of poly-L-amino acids. Such a helix was found in the vicinity of (phi, psi, omega) = 81 degrees, 98 degrees, 170 degrees) in the conformational space. It was the 6.2(20) helix of right-handedness. The 6.2(20) helix, here named the "mu helix," had a cylindrical pore along its helical axis. The diameter of the pore was 6.6 A on the basis of the atom centers of carbonyl carbons and amino nitrogens. The left-handed mu helix was less stable than the right-handed counterpart. The conformation energy of the mu helix, expressed relative to that of the alpha helix of the same polypeptide, depended to a great extent on amino acid composition. It varied over a range of a few kilocalories per mol per residue above and below the conformation energy of the alpha helix of the same polypeptide. This marked diversity in the relative conformation energy of the mu helix can be ascribed primarily to the difference in the relative position of alpha-carbons between the mu and the alpha helices. The conformational entropy made only a small contribution, if any, to the relative stability of the mu helix. There was a hydrogen-bonded network of side chains in the mu helices of poly-L-glutamine and poly-L-asparagine.     

Rotamer libraries and probabilities of transition between rotamers for the side chains in protein-protein binding

Conformational changes in the side chains are essential for protein-protein binding. Rotameric states and unbound- to-bound conformational changes in the surface residues were systematically studied on a representative set of protein complexes. The side-chain conformations were mapped onto dihedral angles space. The variable threshold algorithm was developed to cluster the dihedral angle distributions and to derive rotamers, defined as the most probable conformation in a cluster. Six rotamer libraries were generated: full surface, surface noninterface, and surface interface-each for bound and unbound states. The libraries were used to calculate the probabilities of the rotamer transitions upon binding. The stability of amino acids was quantified based on the transition maps. The noninterface residues' stability was higher than that of the interface. Long side chains with three or four dihedral angles were less stable than the shorter ones. The transitions between the rotamers at the interface occurred more frequently than on the noninterface surface. Most side chains changed conformation within the same rotamer or moved to an adjacent rotamer. The highest percentage of the transitions was observed primarily between the two most occupied rotamers. The probability of the transition between rotamers increased with the decrease of the rotamer stability. The analysis revealed characteristics of the surface side-chain conformational transitions that can be utilized in flexible docking protocols.     

Conformational stability of bovine holo and apo adrenodoxin--a scanning calorimetric study.

Holo and apo adrenodoxin were studied by differential scanning calorimetry, absorption spectroscopy, limited proteolysis, and size-exclusion chromatography. To determine the conformational stability of adrenodoxin, a method was found that prevents the irreversible destruction of the iron-sulfur center. The approach makes use of a buffer solution that contains sodium sulfide and mercaptoethanol. The thermal transition of adrenodoxin takes place at Ttrs = 46-57 degrees C, depending on the Na2S concentration with a denaturation enthalpy of delta H = 300-380 kJ/mol. From delta H versus Ttrs a heat capacity change was determined as delta Cp = 7.5 +/- 1.2 kJ/mol/K. The apo protein is less stable than the holo protein as judged by the lower denaturation enthalpy (delta H = 93 +/- 14 kJ/mol at Ttrs = 37.4 +/- 3.3 degrees C) and the higher proteolytic susceptibility. The importance of the iron-sulfur cluster for the conformational stability of adrenodoxin and some conditions for refolding of the thermally denatured protein are discussed.     

Prediction of conformation of rat galanin in the presence and absence of water with the use of monte Carlo methods and the ECEPP/3 force field

The conformation of the 29-residue rat galanin neuropeptide was studied using the Monte Carlo with energy minimization (MCM) and electrostatically driven Monte Carlo (EDMC) methods. According to a previously elaborated procedure, the polypeptide chain was first treated in a united-residue approximation, in order to enable extensive exploration of the conformational space to be carried out (with the use of MCM), Then the low-energy united-residue conformations were converted to the all-atom representations, and EDMC simulations were carried out for the all-atom polypeptide chains, using the ECEPP/3 force field with hydration included. In order to estimate the effect of environment on galanin conformation, the low-energy conformations obtained as a result of these simulations were taken as starting structures for further EDMC runs that did not include hydration. The lowest-energy conformation obtained in aqueous solution calculations had a nonhelical N-terminal part packed against the nonpolar face of a residual helix that extended from Pro¹³ toward the C-terminus. One next lowest-energy structure was a nearly-all-helical conformation, but with a markedly higher energy. In contrast, all of the low-energy conformations in the absence of water were all-helical differing only by the extent to which the helix was kinked around Pro¹³. These results are in qualitative agreement with the available NMR and CD data of galanin in aqueous and nonaqueous solvents.     

Molecular dynamics study of structure and stability of a model coiled coil

This paper employs methods used earlier to study helix propensity in a model α-helix. The methods are extended to simulations of a motif structure of the α-helical coiled coil, i.e., a structure with a simple amino acid sequence, containing only alanine, leucine, and valine, with leucine and valine forming hydrophobic contacts in the helix interface (positions “d” and “a”). Dynamic simulations of the model coiled-coil structure reproduce characteristic features of the coiled-coil motif seen in experimental studies. Free energy simulations were used to assess the change in stability of the model when a leucine pair or a valine pair in the helix interface was replaced with an alanine pair. A leucine pair at position d was found to contribute 3.4 kcal/mol to the stability of the coiled coil relative to an alanine pair, and a valine pair at postion a was found to contribute 0.8 kcal/mol relative to an alanine pair. The value for the leucine pair agrees with reports in two experimental studies with molecules having different amino sequence. The value for the valine pair is reasonable given the smaller size of the valine side chain and the intrinsic low helix propensity of valine. No experimental value was available for comparison. © 1993 Wiley-Liss, Inc.     

Effect of the side group on the helix-forming tendency of α-alkyl-β-L-aspartamyl residues

The conformational preferences of the monomeric units of a series of poly(α-alkyl-β-L-aspartate)s were examined by quantum mechanical calculations. α-Alkyl-β-aspartamyl m-oligopeptides with a number of residues m ranging from 1 to 7 and arranged in the conformations experimentally observed for their corresponding polymers were computed. Both their total relative energies and their cooperative energy differences were compared as a function of the length of the oligopeptide and the nature of the alkyl side group. Results revealed that the 13/4 helical arrangement is the most stable structure for the isolated polymer chain for any side group, although a 17/4 helix becomes favored in the case of methyl and ethyl groups due to the packing effects. On the other hand, the stability of the 4/1 helix appears to be the preferred conformation for side groups with a branched constitution. © 1997 John Wiley & Sons, Inc. Biopoly 41: 721–729, 1997     

Molecular dynamics simulations of synthetic peptide folding

Because the time scale of protein folding is much greater than that of the widely used simulations of native structures, a detailed report of molecular dynamics simulations of folding has not been available. In this study, we Included the average solvent effect in the potential functions to simplify the calculation of the solvent effect and carried out long molecular dynamics simulations of the alanine-based synthetic peptides at 274 K. From either an extended or a randomly generated conformation, the simulations approached a helix-coil equilibrium in about 3 ns. The multiple minima problem did not prevent helix folding. The calculated helical ratio of Ac-AAQAAAAQAAAAQAAY-NH₂ was 47%, in good agreement with the circular dichroism measurement (about 50%). A helical segment with frayed ends was the most stable conformation, but the hydrophobic interaction favored the compact, distorted helix-turn-helix conformations. The transition between the two types of conformations occurred in a much larger time scale than helix propagation. The transient hydrogen bonds between the glutamine side chain and the backbone carbonyl group could reduce the free energy barrier of helix folding and unfolding. The substitution of a single alanine residue in the middle of the peptide with valine or glycine decreased the average helical ratio significantly, in agreement with experimental observations. © 1996 Wiley-Liss, Inc.     

Serine and threonine residues bend alpha-helices in the chi(1) = g(-) conformation

The relationship between the Ser, Thr, and Cys side-chain conformation (chi(1) = g(-), t, g(+)) and the main-chain conformation (phi and psi angles) has been studied in a selection of protein structures that contain alpha-helices. The statistical results show that the g(-) conformation of both Ser and Thr residues decreases their phi angles and increases their psi angles relative to Ala, used as a control. The additional hydrogen bond formed between the O(gamma) atom of Ser and Thr and the i-3 or i-4 peptide carbonyl oxygen induces or stabilizes a bending angle in the helix 3-4 degrees larger than for Ala. This is of particular significance for membrane proteins. Incorporation of this small bending angle in the transmembrane alpha-helix at one side of the cell membrane results in a significant displacement of the residues located at the other side of the membrane. We hypothesize that local alterations of the rotamer configurations of these Ser and Thr residues may result in significant conformational changes across transmembrane helices, and thus participate in the molecular mechanisms underlying transmembrane signaling. This finding has provided the structural basis to understand the experimentally observed influence of Ser residues on the conformational equilibrium between inactive and active states of the receptor, in the neurotransmitter subfamily of G protein-coupled receptors.