A hydrogen bonded chain in bacteriorhodopsin by computer modelling approach

The seven alpha-helical segments of Bacteriorhodopsin (BR) passing through the membrane are investigated for a continuous Hydrogen Bonded Chain (HBC). The study is carried out by computer modelling approach. It is assumed that the seven helices are placed as (AGFEDCB), which has been accepted as the best model by several groups. Helices A, D, E and G are considered to be present in right handed alpha-helical conformation. The inter-orientation of these helices are represented by Eulerian angles alpha, beta and gamma. For the helices B, C and F which contain Proline in the middle, several conformational possibilities were considered. In these cases apart from the Eulerian angles alpha, beta and gamma, the dihedral angles phi p-1 and psi p-1 of the residues that are succeeded by Proline residue in the helical regions were also used in fixing the position of the helices with respect to each other. All these parameters were varied to fit with the top, middle and bottom distances reported by electron diffraction studies. Good fit was obtained for all right handed alpha-helical conformations and also for helices B, C and F with a left handed turn at the residue preceeding proline. Hence two structures were analysed for continuous HBC. Structure I which contained all the seven helices in right handed alpha-helical conformation and Structure II, which had the helices A, D, E and G in right handed conformation and the helices B, C and F in right handed alpha-helical conformation with a left handed turn at the residue preceeding proline. All possible staggered conformations were considered for the side chains and the inter atomic distances were analysed for Hydrogen bonds. It was possible to obtain a continuous chain in both the structures which includes most of the residues found to be important by the experiments. However Lys-216 has to be considered in two different conformations to connect the cytoplasmic side with the extra cellular side. The overall height spanned by HBC is about 25A. The chains obtained by both the structures I and II are analysed in terms of the conformational parameters. It has also been possible to place the retinal in the region as predicted by the experiments. The Tryptophan residues which affect the spectral characteristics can be aligned on either side of the retinal.     

Proline-induced constraints in alpha-helices

The disrupting effect of a prolyl residue on an α-helix has been analyzed by means of conformational energy computations. In the preferred, nearly α-helical conformations of Ac-Ala₄-Pro-NHMe and of Ac-Ala₇-Pro-Ala₇-NHMe, only the residue preceding Pro is not α-helical, while all other residues can occur in the α-helical A conformation; i.e., it is sufficient to introduce a conformational change of only one residue in order to accommodate proline in a distorted α-helix. Other low-energy conformations exist in which the conformational state of three residues preceding proline is altered considerably; on the other hand, another conformation in which these three residues retain the near-α-helical A-conformational state (with up to 26° changes of their dihedral angles ? and ψ, and a 48° change in one ω from those of the ideal α-helix) has a considerably higher energy. These conclusions are not altered by the substitution of other residues in the place of the Ala preceding Pro. The conformations of the peptide chain next to prolyl residues in or near an α-helix have been analyzed in 58 proteins of known structure, based on published atomic coordinates. Of 331 α-helices, 61 have a Pro at or next to their N-terminus, 21 have a Pro next to their C-terminus, and 30 contain a Pro inside the helix. Of the latter, 16 correspond to a break in the helix, 9 are located inside distorted first turns of the helix, and 5 are parts of irregular helices. Thus, the reported occurrence of prolyl residues next to or inside observed α-helices in proteins is consistent with the computed steric and energetic requirements of prolyl peptides.     

Directional Information Transfer in Protein Helices

IN the course of studies on the relation between conformation and primary sequence in globular proteins, it has become clear that the choice a residue makes between a right handed α-helical conformation (H?) and any other conformation (H) is determined mainly by that residue and its near neighbours in the primary sequence^1–9. When looking for physical mechanisms to explain these findings it is important to know whether the influence of one residue on the conformational state of a neighbouring residue has any directional characteristics. The possibility of certain residues exerting helix-forming influence in either the COOH-terminal or NH₂-terminal direction preferentially is suggested by the non-random distribution of amino-acid residues between the two ends of helical regions in globular proteins^4,10. Ptitsyn's analysis⁴ suggested that a group of residues containing alanine and leucine tends to occur within helical regions, while positively and negatively charged side chain residues are distributed preferentially at the carboxyl-and amino-terminal ends, respectively, of helices.     

Effect of side chains on the conformational energy and rotational strength of the n-pi transition for some alpha-helical poly-alpha-amino acids

Calculations of the dependence of the conformational energy and the rotational strength of the amide n–π* electronic transition (in a series of α-helical polyhel-α- amino acids with different side chains) on conformation have been carried out. The conformational energies were computed by procedures developed in this laboratory; the computation of rotational strengths was carried out by the method of Schellman and Oriel, with a slight modification. Polyamino acids with both nonpolar and polar side chains were considered; in the latter case, it was assumed that the only influence of the polar side chain was on the backbone conformation and on the electrostatic field which perturbs the amide chromophore of the backbone. Only conformations in the range of backbone dihedral angles of the right- and left-handed a-helices were considered, and the assumption of regularity (i.e., uniformity of dihedral angles in every residue) was made. The rotational strength per residue was found to vary markedly with chain length (in oligomers of up to 40 residues long); both the conformational energy per residue and the rotational strength per residue were found to vary significantly with the backbone conformation, which in turn depends on the nature of the side chain. The geometry of the hydrogen bond in the α-helical backbone is the most important factor which influences the dependence of the rotational strength on conformation. The implications of these results, for the interpretation of experimental circular dichroism and optical rotatory dispersion data, are discussed.     

Crystal Structure and Molecular Conformation of the Cyclic Hexapeptide cyclo-(Gly-Aib-Gly)2

We have synthesized and crystallized the cyclic peptide (Gly-Aib-Gly) ₂. Its structure has been determined by conventional X-ray diffracti on methods. In the crystal it adopts a conformation with one β-turn (type I) and its mirror image at the other side of the ring. All conformation al angles are similar to those reported for these amino acid residues. In particular the Aib residue has a conformation intermediate between α- and 3₁₀-helical conformations. The ring is an adequate model for the β-turn conformation. A molecule of formic acid is found in the crystal which shows a very short hydrogen bond with one of the glycine carbonyl groups.     

Studies on the conformation of amino acids. IX. Conformations of butyl, seryl, threonyl, cystennyl, and valyl residues in a dipeptide unit

Potential energies of conformation of a dipeptide unit with butyl, seryl, threonyl, eysteinyl, and valyl side groups have been computed by using classical energy expressions. The presence of a γ-atom introduces characteristic restrictions on the backbone rotational angles ? and ψ the γ-atom itself is restricted to three staggered positions about the C^α—C^β bond. The important results are that a γ-carbon in position I (χ¹ ? 60°) cannot be accommodated in the standard right-and left-handed α-helices, whereas a γ-oxygen or sulfur could easily be accommodated in the right-handed α-helix. Further, a γ-carbon or a heteroatom in position II (χ¹ ? 180°) does not favor a conformation ψ ? 180°, compared to two other positions. The valyl side group significantly reduces the allowed ? and ψ values and energetically prefers a β-conformation compared to right-or left-handed α-helical conformations. The less favorable α-helical conformation is possible only for γ (III, II) combination of the valyl residue. The observed ?, ψ, and χ¹ values of all the amino acid residues in the three protein molecules, lysozyme, myoglobin, and chymotrypsin are compared with the theoretical predictions and the agreement is excellent. The results bring out the important fact that even in large molecules, the conformation of local segments are predominantly governed by the short-range intramolecular interactions.     

Folding and stability of helical proteins: Carp myogen

In this work we use our very simple general representation of protein structures to study the mainly helical protein carp myogen. The representation, which treats the amino acid side-chains as simple spheres, is further simplified by rigidly fixing residues in α-helices. With this model we are able to reproduce the geometry and energetic stability of the native myogen conformation. Studies of the formation of α-helical sub-assemblies showed that the simulated folding of two and four-helix systems worked well, reaching compact native-like conformations with a good rate of success. Greater problems were encountered with the whole molecule (six helices), possibly due to the omission of entropic effects or to simulating the folding too rapidly. Finally, studies of the conformation of a pair of helices when isolated and when part of the whole molecule native conformation showed that long-range interactions have an unexpectedly strong influence on the conformation of the pair of helices.     

Molecular structure for microbial polysaccharides: conformation of the Klebsiella serotype K9 capsular polysaccharide

Fibre type X-ray diffraction patterns have been obtained from oriented, semicrystalline films prepared from the sodium salt form of the bacterial capsular polysaccharide of Klebsiella serotype K9. The molecule has a pentasaccharide repeating sequence, with four neutral residues in the backbone and a glucoronic acid side chain. A novel feature of the molecule is the incorporation of α-l-rhamnose residues, one 1,2 linked and two 1,3 linked in the backbone. Analysis of the X-ray diffraction results indicate an extended three-fold helical conformation with an axially projected chemical repeat of 1.377 nm. Both left and right handed helices have been examined using linked atom least squares techniques to optimize the stereochemistry while simultaneously meeting the observed helical parameters.     

Analysis of the code relating sequence to conformation in globular proteins. An informational analysis of the role of the residue in determining the conformation of its neighbours in the primary sequence

1. The effect exerted by a residue on the conformation of neighbouring residues was analysed by using data from nine globular proteins of known sequence and conformation. 2. An information measure was used which estimated the role of a residue in influencing neighbouring conformations and also its tendency to influence the lengths of runs of residues in that conformation. This measure was estimated for each residue in all conformations defined by domains on the varphi, psi diagram. 3. Plots of the information measure yielded an intercept, which was a measure of intra-residue information for a residue. The slope was a measure of the statistical co-operativity or tendency of the residue to influence the occurrence of its neighbours in runs of a particular conformation. Both parameters are a function of the residue type. Statistical co-operativity is found in the alpha(1)-helical (H(1)) and beta-pleated-sheet (P(2)) conformations and, to a lesser extent, in their distorted variants H(2) and P(1). 4. The directional nature of these influences for H(1) and P(2) conformations is illustrated by plots of the information measure against the distance m from the residue, for m=-10 to +10. 5. The results for statistical co-operativity are discussed in relation to theories of helix-coil and pleated-sheet-coil transitions. The value of the information-theory-derived parameters in obtaining s parameters for the Zimm & Bragg (1959) equations is illustrated. 6. Directional effects are discussed with particular relation to mechanisms of the termination of helices and the involvement of the alpha(II) conformation and also to discontinuities in pleated-sheet conformations.     

Molecular orbital calculations on the conformation of polypeptides and proteins. V. Conformational energy maps and stereochemical rotational states of aliphatic residues

The quantum-mechanical calculations by the PCILO method on the conformation of amino acid residues of proteins have been extended to the valyl, leucyl, and isoleucyl residues. In distinction to the earlier “empirical” computations, the quantum-mechanical results indicate very similar energy contours for the stable conformations of the three residues. Their general outline is also similar to that of the alanyl residue, although reduced by about 25%. Contrary to the “empirical” computations, the present results predict that the region corresponding to the α-helix should be one of great stability for the three residues and in particular for the valyl residue. The quantum-mechanical results are in excellent agreement with the experimental conformations of the aliphatic residues in lysozyme and myoglobin. Their prediction as to the ready availability of the valyl residue in the α-helical conformation agrees moreover with Ptitsyn's statistical evaluation of the participation of this residue in the inner turns of the helical regions in six globular proteins. The maximum conformational space allowed for the aliphatic residues is somewhat smaller than that allowed for the aromatic ones, while the minimum conformational space (region of stability common to all the residues) is similar in both groups.     

Left handed �� helix models for mammalian prion fibrils

We propose models for in vitro grown mammalian prion protein fibrils based upon left handed beta helices formed both from the N-terminal and C-terminal regions of the proteinase resistant infectious prion core. The C-terminal threading onto a β-helical structure is almost uniquely determined by fixing the cysteine disulfide bond on a helix corner. In comparison to known left handed helical peptides, the resulting model structures have similar stability attributes including relatively low root mean square deviations in all atom molecular dynamics, substantial side-chain-to-side-chain hydrogen bonding, good volume packing fraction, and low hydrophilic/hydrophobic frustration. For the N-terminus, we propose a new threading of slightly more than two turns, which improves upon the above characteristics relative to existing three turn β-helical models. The N-terminal and C-terminal beta helices can be assembled into eight candidate models for the fibril repeat units, held together by large hinge (order 30 residues) domain swapping, with three amenable to fibril promoting domain swapping via a small (five residue) hinge on the N-terminal side. Small concentrations of the metastable C-terminal β helix in vivo might play a significant role in templating the infectious conformation and in enhancing conversion kinetics for inherited forms of the disease and explain resistance (for canines) involving hypothesized coupling to the methionine 129 sulfur known to play a role in human disease.Key words: prion, amyloid fibril, domain swap, beta helix, computational biology     

The effects of changes in ring geometry on computer models of amylose

A review of single-crystal studies shows that α-D-glucopyranose, residues of which constitute the monomeric units of amylose, is flexible within the constraints of the Cl conformation, and that the internal differences among the rings are most clearly indicated by the variety in ring-torsion angles (or conformation angles). An index of the cumulative effect of changes in these angles is provided by the length of the virtual bond, O-1—O-4, and classification of residue geometries by virtual bond-length permits a systematic selection of suitable residues for the construction of models of amylose. By the use of D-glucopyranose residues having different geometries, it is possible to build models of (a) V-amylose helices having 6, 7, and 8 residues per turn, (b) single and double helical B-amyloses, and (c) KBr-amylose, all of which satisfy reasonable stereochemical criteria. Because no single residue can satisfactorily model all of the well known polymorphs of amylose, it is suggested that structural determinations that utilize a rigid residue approximation should make use of the full range of known, residue geometries.     

Structure and localization of an essential transmembrane segment of the proton translocation channel of yeast H-V-ATPase

Vacuolar (H⁺)-ATPase (V-ATPase) is a proton pump present in several compartments of eukaryotic cells to regulate physiological processes. From biochemical studies it is known that the interaction between arginine 735 present in the seventh transmembrane (TM7) segment from subunit a and specific glutamic acid residues in the subunit c assembly plays an essential role in proton translocation. To provide more detailed structural information about this protein domain, a peptide resembling TM7 (denoted peptide MTM7) from Saccharomyces cerevisiae (yeast) V-ATPase was synthesized and dissolved in two membrane-mimicking solvents: DMSO and SDS. For the first time the secondary structure of the putative TM7 segment from subunit a is obtained by the combined use of CD and NMR spectroscopy. SDS micelles reveal an α-helical conformation for peptide MTM7 and in DMSO three α-helical regions are identified by 2D ¹H-NMR. Based on these conformational findings a new structural model is proposed for the putative TM7 in its natural environment. It is composed of 32 amino acid residues that span the membrane in an α-helical conformation. It starts at the cytoplasmic side at residue T719 and ends at the luminal side at residue W751. Both the luminal and cytoplasmatic regions of TM7 are stabilized by the neighboring hydrophobic transmembrane segments of subunit a and the subunit c assembly from V-ATPase.     

A survey of left-handed polyproline II helices

Left-handed polyproline II helices (PPII) are contiguous elements of protein secondary structure in which the phi and psi angles of constituent residues are restricted to around -75 degrees and 145 degrees, respectively. They are important in structural proteins, in unfolded states and as ligands for signaling proteins. Here, we present a survey of 274 nonhomologous polypeptide chains from proteins of known structure for regions that form these structures. Such regions are rare, but the majority of proteins contain at least one PPII helix. Most PPII helices are shorter than five residues, although the longest found contained 12 amino acids. Proline predominates in PPII, but Gln and positively charged residues are also favored. The basis of Gln's prevalence is its ability to form an i, i + 1 side-chain to main-chain hydrogen bond with the backbone carbonyl oxygen of the proceeding residue; this helps to fix the psi angle of the Gln and the phi and psi of the proceeding residue in PPII conformations and explains why Gln is favored at the first position in a PPII helix. PPII helices are highly solvent exposed, which explains why apolar amino acids are disfavored despite preferring this region of phi/psi space when in isolation. PPII helices have perfect threefold rotational symmetry and within these structures we find significant correlation between the hydrophobicity of residues at i and i + 3; thus, PPII helices in globular proteins can be considered to be amphipathic.     

Far-infrared spectra of sequential polypeptides with -helical and polyglycine II structures

The sequential copolymers of glycine and L -alanine, L-valine and L -alanine, L-leucine and L -alanine, and L-phenylalanine and L -alanine and those containing the L-proline residues were synthesized. The infrared spectra in the region from 700 to 200 cm^-1 were measured for these polypeptides with the α-helical conformation or the polyglycine II structure and compared with the spectra of the β-form structures. The results showed that several infrared bands observed in the region from 600 to 200 cm^-1 clearly reflect not only the backbone conformations but also the local conformations of component amino acid residues of polypeptides with the α-helical, β-form and polyglycine II structures.     

Raman spectra of copoly(D,L-alanines)

Raman spectra in the region 1000–150 cm^?1 were measured for copoly(D ,L -alanines) with the D -residue contents, 3, 7, 10, and 20%, and compared with the spectrum of the α-helical poly-L -alanine. The 532- and 378-cm^?1 peaks were assigned to the L -residues with a right-handed α-helix-like local conformation or to the D -residues with a left-handed α-helix-like local conformation. From the intensity of the latter peak the contents of these local conformations were estimated as a function of the D -residue contents for the copolymers. The 264-cm^?1 peak, which has been assigned to the breathing vibration of the α-helical poly-L -alanine, shows a marked decrease in its intensity upon the introduction of the D residues. This result suggests that the overall deformation vibration of the α-helix arises from rather long sequences of the L - and D -alanine residues with the α-helical conformation and that the intensity of this vibration depends on the content of these sequences in the copolymers.     

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

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

Conformational 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.     

An experimental Ramachandran plot for retropeptide derivatives: Conformational features of derivatives of GEM-diamino and malonyl amino acids

Malonic and diaminomethane residues, equivalent to the two possible retro modifications of a glycine unit, with an inverted peptide group, present particular conformations that differ from those found in glycine and, in general, in α-amino acids. In both cases the φ_i and ψ_i torsional angles have restricted values as deduced from inspection of the Cambridge Structural Data Bank and from compounds studied by us. Thus, both ψ_i angles tend to be equal to 115° (or −115°) in the malonyl residues, whereas the φ_i angles tend to be equal to 88° (or −88°) in the diaminomethane residues. These results are in agreement with previous experimental data on polymers, but in the case of malonyl residues they differ from theoretical calculations on isolated molecules. The experimental data for both residues can be represented in a way similar to the usual Ramachandran plot, which will be useful in analyzing the incorporation of these residues into proteins. When side chains are present in either type of residue, they are similar to conventional α-amino acids, although the orientation of the peptide groups is different. In such cases they acquire conformations similar to those found in peptide residues in the α-helix and β-sheet conformations, although other conformations are also possible. © 1998 John Wiley & Sons, Inc. Biopoly 45: 149–155, 1998     

Conformation of amino acid residue contiguous to helical polypeptide

The calculation of the conformational energy of the terminal D - or L -alanine residue contiguous to an α-helical polypeptide, polyalanine, was made. Both L -and D -residues contiguous to the carboxyl terminal of α-helical poly(L -alanine) are considered to prefer the α-helical conformation due to the effect of the α-helical structure of the polymer. The residue at the amino terminal is found to be less affected by the α-helical structure of the polymer.