PackHelix: a tool for helix-sheet packing during protein structure prediction

The three‐dimensional structure of a protein is organized around the packing of its secondary structure elements. Predicting the topology and constructing the geometry of structural motifs involving α‐helices and/or β‐strands are therefore key steps for accurate prediction of protein structure. While many efforts have focused on how to pack helices and on how to sample exhaustively the topologies and geometries of multiple strands forming a β‐sheet in a protein, there has been little progress on generating native‐like packings of helices on sheets. We describe a method that can generate the packing of multiple helices on a given β‐sheet for αβα sandwich type protein folds. This method mines the results of a statistical analysis of the conformations of αβ₂ motifs in protein structures to provide input values for the geometric attributes of the packing of a helix on a sheet. It then proceeds with a geometric builder that generates multiple arrangements of the helices on the sheet of interest by sampling through these values and performing consistency checks that guarantee proper loop geometry between the helices and the strands, minimal number of collisions between the helices, and proper formation of a hydrophobic core. The method is implemented as a module of ProteinShop. Our results show that it produces structures that are within 4–6 Å RMSD of the native one, regardless of the number of helices that need to be packed, though this number may increase if the protein has several helices between two consecutive strands in the sequence that pack on the sheet formed by these two strands. Proteins 2011; Published 2011 Wiley‐Liss, Inc.     

Folding of peptide fragments comprising the complete sequence of proteins. Models for initiation of protein folding. I. Myohemerythrin.

In an attempt to delineate potential folding initiation sites for different protein structural motifs, we have synthesized series of peptides that span the entire length of the polypeptide chain of two proteins, and examined their conformational preferences in aqueous solution using proton nuclear magnetic resonance and circular dichroism spectroscopy. We describe here the behavior of peptides derived from a simple four-helix bundle protein, myohemerythrin. The peptides correspond to the sequences of the four long helices (the A, B, C and D helices), the N- and C-terminal loops and the connecting sequences between the helices. The peptides corresponding to the helices of the folded protein all exhibit preferences for helix-like conformations in solution. The conformational ensembles of the A- and D-helix peptides contain ordered helical forms, as shown by extensive series of medium-range nuclear Overhauser effect connectivities, while the B- and C-helix peptides exhibit conformational preferences for nascent helix. All four peptides adopt ordered helical conformations in mixtures of trifluoroethanol and water. The terminal and interconnecting loop peptides also appear to contain appreciable populations of conformers with backbone phi and psi angles in the alpha-region and include highly populated hydrophobic cluster and/or turn conformations in some cases. Trifluoroethanol is unable to drive these peptides towards helical conformations. Overall, the peptide fragments of myohemerythrin have a marked preference towards secondary structure formation in aqueous solution. In contrast, peptide fragments derived from the beta-sandwich protein plastocyanin are relatively devoid of secondary structure in aqueous solution (see accompanying paper). These results suggest that the two different protein structural motifs may require different propensities for formation of local elements of secondary structure to initiate folding, and that there is a prepartitioning of conformational space determined by the local amino acid sequence that is different for the helical and beta-sandwich structural motifs.     

Determination of the secondary structural elements of chicken liver fatty acid binding protein by two‐dimensional homonuclear NMR

A conformational study in solution of the fatty acid binding protein from chicken liver is presented. The nearly complete sequence‐specific ¹H resonance assignment was achieved from homonuclear two‐dimensional nmr experiments using a sample of native protein. The principal elements of secondary structure were identified: 10 antiparallel β‐strands and one helical segment followed by a turn comprising 5 residues. These elements correspond closely with those of the crystal structure of the related protein, and two new secondary structural features obtained from the nmr data are the β‐sheet conformation between the first and the last β‐strand in the protein sequence, as well as a helical loop at the N‐terminus of the polypeptide chain. © 1999 John Wiley & Sons, Inc. Biopoly 50: 1–11, 1999     

Alpha-hairpin stability and folding of transmembrane segments

Molecular Dynamics (MD) simulations at low dielectric constant have been carried out for peptides matching the double spanning segments of transmembrane proteins. Different folding dynamics have been observed. The peptides folded into the stable helix-turn-helix conformation-alpha-hairpin-with antiparallel-oriented strands or unstable alpha-hairpin conformation that unfolded later into the straight helical structure. The peptide having flexible residues in the TM helices often misfolded into a tangled structure that can be avoided by restricting the flexibility of these residues. General conclusions can be drawn from the observed folding dynamics. The stability and folding of some double spanning transmembrane fragments are self-assembling. The following and/or neighboring peptide chains of the protein may support the stability of the hairpin structure of other fragments. The stability of the TM helices containing flexible residues could be maintained due to contacts with neighboring TM segments.     

Folding of peptide fragments comprising the complete sequence of proteins. Models for initiation of protein folding. II. Plastocyanin.

In an attempt to understand the earliest events in the protein folding pathway, the complete sequence of French bean plastocyanin has been synthesized as a series of short peptide fragments, and the conformational preferences of each peptide examined in aqueous solution using proton n.m.r. methods. Plastocyanin consists largely of beta-sheet, with reverse turns and loops between the strands of the sheet, and one short helix. The n.m.r. experiments indicate that most of the peptides derived from the plastocyanin sequence have remarkably little propensity to adopt folded conformations in aqueous solution, in marked contrast to the peptides derived from the helical protein, myohemerythrin (accompanying paper). For most plastocyanin peptides, the backbone dihedral angles are predominantly in the beta-region of conformational space. Some of the peptides show weak NOE connectivities between adjacent amide protons, indicative of small local populations of backbone conformations in the a region of (phi,psi) space. A conformational preference for a reverse turn is seen in the sequence Ala65-Pro-Gly-Glu68, where a turn structure is found in the folded protein. Significantly, the peptide sequences that populate the alpha-region of (phi,psi) space are mostly derived from turn and loop regions in the protein. The addition of trifluoroethanol does not drive the peptides into helical conformations. In one region of the sequence, the n.m.r. spectra provide evidence of the formation of a hydrophobic cluster involving aromatic and aliphatic side-chains. These results have significance for understanding the initiation of protein folding. From these studies of the fragments of plastocyanin (this paper) and myohemerythrin (accompanying paper), it appears that there is a pre-partitioning of the conformational space sampled by the polypeptide backbone that is related to the secondary structure in the final folded state.     

Three-dimensional solution structure of the reduced form of Escherichia coli thioredoxin determined by nuclear magnetic resonance spectroscopy

The three-dimensional solution structure of reduced (dithiol) thioredoxin from Escherichia coli has been determined with distance and dihedral angle constraints obtained from 1H NMR spectroscopy. Reduced thioredoxin has a well-defined global fold consisting of a central five-strand beta-sheet and three long helices. The beta-strands are packed in the sheet in the order beta 1 beta 3 beta 2 beta 4 beta 5, with beta 1, beta 3, and beta 2 parallel and beta 2, beta 4, and beta 5 arranged in an antiparallel fashion. Two of the helices connect strands of the beta-sheet: alpha 1 between beta 1 and beta 2 and alpha 2 between beta 2 and beta 3. Strands beta 4 and beta 5 are connected by a short loop that contains a beta-bulge. Strands beta 3 and beta 4 are connected by a long loop that contains a series of turn-like or 3(10) helical structures. The active site Cys-Gly-Pro-Cys sequence forms a protruding loop between strand beta 2 and helix alpha 2. The structure is very similar overall to that of oxidized (disulfide) thioredoxin obtained from X-ray crystal structure analysis but differs in the local conformation of the active site loop. The distance between the sulfurs of Cys 32 and Cys 35 increases from 2.05 A in the disulfide bridge to 6.8 +/- 0.6 A in the dithiol of reduced thioredoxin, as a result of a rotation of the side chain of Cys 35 and a significant change in the position of Pro 34. This conformational change has important implications for the mechanism of thioredoxin as a protein disulfide oxidoreductase.     

Conformational analysis and clustering of short and medium size loops connecting regular secondary structures: a database for modeling and prediction.

Loops are regions of nonrepetitive conformation connecting regular secondary structures. We identified 2,024 loops of one to eight residues in length, with acceptable main-chain bond lengths and peptide bond angles, from a database of 223 protein and protein-domain structures. Each loop is characterized by its sequence, main-chain conformation, and relative disposition of its bounding secondary structures as described by the separation between the tips of their axes and the angle between them. Loops, grouped according to their length and type of their bounding secondary structures, were superposed and clustered into 161 conformational classes, corresponding to 63% of all loops. Of these, 109 (51% of the loops) were populated by at least four nonhomologous loops or four loops sharing a low sequence identity. Another 52 classes, including 12% of the loops, were populated by at least three loops of low sequence similarity from three or fewer nonhomologous groups. Loop class suprafamilies resulting from variations in the termini of secondary structures are discussed in this article. Most previously described loop conformations were found among the classes. New classes included a 2:4 type IV hairpin, a helix-capping loop, and a loop that mediates dinucleotide-binding. The relative disposition of bounding secondary structures varies among loop classes, with some classes such as beta-hairpins being very restrictive. For each class, sequence preferences as key residues were identified; those most frequently at these conserved positions than in proteins were Gly, Asp, Pro, Phe, and Cys. Most of these residues are involved in stabilizing loop conformation, often through a positive phi conformation or secondary structure capping. Identification of helix-capping residues and beta-breakers among the highly conserved positions supported our decision to group loops according to their bounding secondary structures. Several of the identified loop classes were associated with specific functions, and all of the member loops had the same function; key residues were conserved for this purpose, as is the case for the parvalbumin-like calcium-binding loops. A significant number, but not all, of the member loops of other loop classes had the same function, as is the case for the helix-turn-helix DNA-binding loops. This article provides a systematic and coherent conformational classification of loops, covering a broad range of lengths and all four combinations of bounding secondary structure types, and supplies a useful basis for modelling of loop conformations where the bounding secondary structures are known or reliably predicted.     

Environmentally induced reversible conformational switching in the yeast cell adhesion protein alpha-agglutinin

The yeast cell adhesion protein alpha-agglutinin is expressed on the surface of a free-living organism and is subjected to a variety of environmental conditions. Circular dichroism (CD) spectroscopy shows that the binding region of alpha-agglutinin has a beta-sheet-rich structure, with only approximately 2% alpha-helix under native conditions (15-40 degrees C at pH 5.5). This region is predicted to fold into three immunoglobulin-like domains, and models are consistent with the CD spectra as well as with peptide mapping and site-specific mutagenesis. However, secondary structure prediction algorithms show that segments comprising approximately 17% of the residues have high alpha-helical and low beta-sheet potential. Two model peptides of such segments had helical tendencies, and one of these peptides showed pH-dependent conformational switching. Similarly, CD spectroscopy of the binding region of alpha-agglutinin showed reversible conversion from beta-rich to mixed alpha/beta structure at elevated temperatures or when the pH was changed. The reversibility of these changes implied that there is a small energy difference between the all-beta and the alpha/beta states. Similar changes followed cleavage of peptide or disulfide bonds. Together, these observations imply that short sequences of high helical propensity are constrained to a beta-rich state by covalent and local charge interactions under native conditions, but form helices under non-native conditions.     

Double helix formation in α‐peptides: a theoretical study

A complete overview on all possible hydrogen bonding patterns of double helices with antiparallel and parallel strand orientation in α‐peptide sequences is provided on the basis of ab initio molecular orbital theory. The most stable representatives belong to the group of antiparallel helices. The study on side chain influence shows that these double helices can only be realized if the strands are composed of L ‐ and D ‐amino acids in alternate order. The stability of the double helices is compared with that of competing single‐stranded helices. The data contribute to an understanding of secondary structure formation in peptides and provide a basis for a rational design of membrane channels. Copyright © 2010 European Peptide Society and John Wiley & Sons, Ltd.     

Cooperative effects in hydrogen-bonding of protein secondary structure elements: a systematic analysis of crystal data using Secbase

A systematic analysis of the hydrogen-bonding geometry in helices and beta sheets has been performed. The distances and angles between the backbone carbonyl O and amide N atoms were correlated considering more than 1500 protein chains in crystal structures determined to a resolution better than 1.5 A. They reveal statistically significant trends in the H-bond geometry across the different secondary structural elements. The analysis has been performed using Secbase, a modular extension of Relibase (Receptor Ligand Database) which integrates information about secondary structural elements assigned to individual protein structures with the various search facilities implemented into Relibase. A comparison of the mean hydrogen-bond distances in alpha helices and 3(10) helices of increasing length shows opposing trends. Whereas in alpha helices the mean H-bond distance shrinks with increasing helix length and turn number, the corresponding mean dimension in 3(10) helices expands in a comparable series. Comparing similarly the hydrogen-bond lengths in beta sheets there is no difference to be found between the mean H-bond length in antiparallel and parallel beta sheets along the strand direction. In contrast, an interesting systematic trend appears to be given for the hydrogen bonds perpendicular to the strands bridging across an extended sheet. With increasing number of accumulated strands, which results in a growing number of back-to-back piling hydrogen bonds across the strands, a slight decrease of the mean H-bond distance is apparent in parallel beta sheets whereas such trends are obviously not given in antiparallel beta sheets. This observation suggests that cooperative effects mutually polarizing spatially well-aligned hydrogen bonds are present either in alpha helices and parallel beta sheets whereas such influences seem to be lacking in 3(10) helices and antiparallel beta sheets.     

Conformational analysis corresponding to intra-chain disulfide bridged peptides in proteins of known three-dimensional structure

We have carried out a systematic analysis in order to evaluate whether Intra-Chain Disulfide Bridged Peptides (ICDBPs) observed in proteins of known three-dimensional structure adopt structurally similar conformations as they may correspond to structural/functional motifs. 406 representative ICDBPs comprising between 3 to 17 amino acid residues could be classified according to peptide sequence length and main-chain secondary structure conformation into 146 classes. ICDBPs comprising 6 amino acid residues are maximally represented in the Protein Data Bank. They also represent the maximum number of main-chain secondary structure conformational classes. Individual ICDBPs in each class represent different protein superfamilies and correspond to different amino acid sequences. We identified 145 ICDBP pairs that had not less-than 0.5 A root mean square deviation value corresponding to their equivalent peptide backbone atoms. We believe these ICDBPs represent structural motifs and possible candidates in order to further explore their structure/function role in the corresponding proteins. The common conformational classes observed for ICDBPs defined according to the main-chain secondary structure conformations; H (helix), B (residue in a isolated beta bridge), C (coil), E (extended beta strand), G (3(10) helix), I (pi helix), S (bend), T (hydrogen-bonded turn) were; "CHHH", "CTTC", "CSSS" and "CSSC" (for ICDBP length 4), "CSSCC" (length 5), "EETTEE", "CCSSCC", "CCSSSC" (length 6), "EETTTEE" (length 7), "EETTTTEE" (length 8), "EEEETTEEEE" (length 10), "EEEETTTEEEE" (length 11) and "EEEETTTTEEEE" (length 12).     

The dual role of a loop with low loop contact distance in folding and domain swapping

Alpha helices, beta strands, and loops are the basic building blocks of protein structure. The folding kinetics of alpha helices and beta strands have been investigated extensively. However, little is known about the formation of loops. Experimental studies show that for some proteins, the formation of a single loop is the rate-determining step for folding, whereas for others, a loop (or turn) can misfold to serve as the hinge loop region for domain-swapped species. Computer simulations of an all-atom model of fragment B of Staphylococcal protein A found that the formation of a single loop initiates the dominant folding pathway. On the other hand, the stability analysis of intermediates suggests that the same loop is a likely candidate to serve as a hinge loop for domain swapping. To interpret the simulation result, we developed a simple structural parameter: the loop contact distance (LCD), or the sequence distance of contacting residues between a loop and the rest of the protein. The parameter is applied to a number of other proteins, including SH3 domains and prion protein. The results suggest that a locally interacting loop (low LCD) can either promote folding or serve as the hinge region for domain swapping. Thus, there is an intimate connection between folding and domain swapping, a possible cause of misfolding and aggregation.     

Folding propensities of peptide fragments of myoglobin.

Myoglobin has been studied extensively as a paradigm for protein folding. As part of an ongoing study of potential folding initiation sites in myoglobin, we have synthetized a series of peptides covering the entire sequence of sperm whale myoglobin. We report here on the conformation preferences of a series of peptides that cover the region from the A helix to the FG turn. Structural propensities were determined using circular dichroism and nuclear magnetic resonance spectroscopy in aqueous solution, trifluoroethanol, and methanol. Peptides corresponding to helical regions in the native protein, namely the B, C, D, and E helices, populate the alpha region of (phi, psi) space in water solution but show no measurable helix formation except in the presence of trifluoroethanol. The F-helix sequence has a much lower propensity to populate helical conformations even in TFE. Despite several attempts, we were not successful in synthesizing a peptide corresponding to the A-helix region that was soluble in water. A peptide termed the AB domain was constructed spanning the A- and B-helix sequences. The AB domain is not soluble in water, but shows extensive helix formation throughout the peptide when dissolved in methanol, with a break in the helix at a site close to the A-B helix junction in the intact folded myoglobin protein. With the exception of one local preference for a turn conformation stabilized by hydrophobic interactions, the peptides corresponding to turns in the folded protein do not measurably populate beta-turn conformations in water, and the addition of trifluoroethanol does not enhance the formation of either helical or turn structure. In contrast to the series of peptides described here, either studies of peptides from the GH region of myoglobin show a marked tendency to populate helical structures (H), nascent helical structures (G), or turn conformations (GH peptide) in water solution. This region, together with the A-helix and part of the B-helix, has been shown to participate in an early folding intermediate. The complete analysis of conformational properties of isolated myoglobin peptides supports the hypothesis that spontaneous secondary structure formation in local regions of the polypeptide may play an important role in the initiation of protein folding.     

Effect of the N-terminal glycine on the secondary structure, orientation, and interaction of the influenza hemagglutinin fusion peptide with lipid bilayers.

The amino-terminal segment of the membrane-anchored subunit of influenza hemagglutinin (HA) plays a crucial role in membrane fusion and, hence, has been termed the fusion peptide. We have studied the secondary structure, orientation, and effects on the bilayer structure of synthetic peptides corresponding to the wild-type and several fusogenic and nonfusogenic mutants with altered N-termini of the influenza HA fusion peptide by fluorescence, circular dichroism, and Fourier transform infrared spectroscopy. All peptides contained segments of alpha-helical and beta-strand conformation. In the wild-type fusion peptide, 40% of all residues were in alpha-secondary and 30% in beta-secondary structures. By comparison, the nonfusogenic peptides exhibited larger beta/alpha secondary structure ratios. The order parameters of the helices and the amide carbonyl groups of the beta-strands of the wild-type fusion peptide were measured separately, based on the infrared dichroism of the respective absorption bands. Order parameters in the range 0.1-0.7 were found for both segments of the wild-type peptide, which indicates that they are most likely aligned at oblique angles to the membrane normal. The nonfusogenic but not the fusogenic peptides induced splitting of the infrared absorption band at 1735 cm(-1), which is assigned to stretching vibrations of the lipid ester carbonyl bond. This splitting, which reports on an alteration of the hydrogen bonds formed between the lipid ester carbonyls and water and/or hydrogen-donating groups of the fusion peptides, correlated with the beta/alpha ratio of the peptides, suggesting that unpaired beta-strands may replace water molecules and hydrogen-bond to the lipid ester carbonyl groups. The profound structural changes induced by single amino acid replacements at the extreme N-terminus of the fusion peptide further suggest that tertiary or quaternary structural interactions may be important when fusion peptides bind to lipid bilayers.     

Ligand-assisted inhibition in cytochrome P450 158A2 from Streptomyces coelicolor A3(2)

Cytochrome P450 158A2 (CYP158A2) can polymerize flaviolin to red-brown pigments, which may afford physical protection to the organism, possibly against the deleterious effects of UV radiation. We have found that the small molecule malonic acid enables cocrystallization of this mixed function oxidase with the azole inhibitor 4-phenylimidazole. The presence of malonate molecules affects the behavior of the binding of 4-phenylimidazole to CYP158A2 and increases inhibition potency up to 2-fold compared to 4-phenylimidazole alone. We report here the crystal structure of the 4-phenylimidazole/malonate complex of CYP158A2 at 1.5 A. Two molecules of malonate used in crystallization are found above the single inhibitor molecule in the active site. Those two molecules are linked between the BC loop and beta 1-4/beta 6-1 strands via hydrogen bond interactions to stabilize the conformational changes of the BC loop and beta strands that take place upon inhibitor binding compared to the ligand-free structure we have reported previously. 4-Phenylimidazole can launch an extensive hydrogen-bonding network in the region of the F/G helices which may stabilize the conformational changes. Our findings clearly show that two molecules of malonate assist the inhibitor 4-phenylimidazole to assume a specific location producing more inhibition in the enzyme catalytic activity.     

Understanding peptide competitive inhibition of botulinum neurotoxin a binding to SV2 protein via molecular dynamics simulations

Botulinum neurotoxins (BoNTs) are known as the most toxic natural substances. Synaptic vesicle protein 2 (SV2) has been proposed to be a protein receptor for BoNT/A. Recently, two short peptides (BoNT/A‐A2 and SV2C‐A3) were designed to inhibit complex formation between the BoNT/A receptor‐binding domain (BoNT/A‐RBD) and the synaptic vesicle protein 2C luminal domain (SV2C‐LD). In this article, the two peptide complex systems are studied by molecular dynamics (MD) simulations. The structural stability analysis indicates that BoNT/A‐A2 system is more stable than SV2C‐A3 system. The conformational analysis implies that the β‐sheet in BoNT/A‐A2 system maintains its secondary structure but the two β‐strands in SV2C‐A3 system have remarkable conformational changes. Based on the calculation of hydrogen bonds, hydrophobic interactions and cation‐π interactions, it is found that the internal hydrogen bonds play crucial roles in the structural stability of the peptides. Because of the stable secondary structure, the β‐sheet in BoNT/A‐A2 system establishes effective interactions at the interface and inhibits BoNT/A‐RBD binding to SV2C‐LD. In contrast, without other β‐strands forming internal hydrogen bonds, the two isolated β‐strands in SV2C‐A3 system become the random coil. This conformational change breaks important hydrogen bonds and weakens cation‐π interaction in the interface, so the complex formation is only partially inhibited by the two β‐strands. These results are consistent with experimental studies and may be helpful in understanding the inhibition mechanisms of peptide inhibitors. © 2015 Wiley Periodicals, Inc. Biopolymers 103: 597–608, 2015.     

alpha beta Spectrin coiled coil association at the tetramerization site

On the basis of sequence homology studies, it has been suggested that the association of human erythrocytes alpha and beta spectrin at the tetramerization site involves interactions between helices. However, no empirical details are available, presumably due to the experimental difficulties in studying spectrin molecules because of its size and/or its structural flexibility. It has been speculated that erythrocyte tetramerization involves helical bundling rather than coiled coil association. We have used recombinant spectrin peptides to model alpha and beta spectrin to study their association at the tetramerization site. Two alpha peptides, Sp alpha 1-156 and Sp alpha 1-368, and one beta peptide, Sp beta 1898-2083, were used as model peptides to demonstrate the formation of the alpha beta complex. We also found that the replacement of R28 in Sp alpha 1-368 to give Sp alpha 1-368R28C abolished complex formation with the beta peptide. Circular dichroism techniques were used to monitor the secondary structures of the individual peptides and of the complex, and the results showed that both Sp alpha 1-156 and Sp beta 1898-2083 peptides in solution, separately, included helices that were not paired with other helices in the absence of their binding partners. However, in a mixture of Sp alpha 1-156 and Sp beta 1898-2083 and formation of the alpha beta complex, the unpaired helices associated to form coiled coils. Since the sequences of these two peptides that are involved in the coiled coil association are derived from a native protein, the information obtained from this study also provides insight toward a better understanding of naturally occurring coiled coil subunit-subunit association.     

Insights into the domains required for dimerization and assembly of human alphaB crystallin

Protein pin array technology was used to identify subunit-subunit interaction sites in the small heat shock protein (sHSP) alphaB crystallin. Subunit-subunit interaction sites were defined as consensus sequences that interacted with both human alphaA crystallin and alphaB crystallin. The human alphaB crystallin protein pin array consisted of contiguous and overlapping peptides, eight amino acids in length, immobilized on pins that were in a 96-well ELISA plate format. The interaction of alphaB crystallin peptides with physiological partner proteins, alphaA crystallin and alphaB crystallin, was detected using antibodies and recorded using spectrophotometric absorbance. Five peptide sequences including 37LFPTSTSLSPFYLRPPSF54 in the N terminus, 75FSVNLDVK82)(beta3), 131LTITSSLS138 (beta8) and 141GVLTVNGP148 (beta9) that form beta strands in the conserved alpha crystallin core domain, and 155PERTIPITREEK166 in the C-terminal extension were identified as subunit-subunit interaction sites in human alphaB crystallin using the novel protein pin array assay. The subunit-subunit interaction sites were mapped to a three-dimensional (3D) homology model of wild-type human alphaB crystallin that was based on the crystal structure of wheat sHSP16.9 and Methanococcus jannaschi sHSP16.5 (Mj sHSP16.5). The subunit-subunit interaction sites identified and mapped onto the homology model were solvent-exposed and had variable secondary structures ranging from beta strands to random coils and short alpha helices. The subunit-subunit interaction sites formed a pattern of hydrophobic patches on the 3D surface of human alphaB crystallin.     

Conformational analysis of peptide surrogates. Reduced and retro-amide links in blocked alanine and in secondary structures.

We have investigated the conformational effects of modifying the amide link of peptides. We studied a reverse amide bond psi [NHCO], a reduced amide bond psi [CH2NH] and a retro-reduced bond psi [NHCH2] as surrogates for the amide link [CONH] in native peptides. A complete search of the conformational space available to residues with these modified links was carried out. The local minima and the rotational barriers were described and compared to the minima of the native residue. The results are compatible with the available observed structural data. These modified links have been incorporated in secondary structure units such as beta turns, alpha helices, and parallel and anti-parallel beta sheets. It was found that a reduced amide link can lead to stabilised beta turns, while the retro modification can be incorporated in stable beta sheets. A significant reduction in the stability of alpha helices is caused by the retro links, while a reduced amide link results in only a small destabilisation.     

Constraints on the conformation of the cytoplasmic face of dark-adapted and light-excited rhodopsin inferred from antirhodopsin antibody imprints

Rhodopsin is the best-understood member of the large G protein-coupled receptor (GPCR) superfamily. The G-protein amplification cascade is triggered by poorly understood light-induced conformational changes in rhodopsin that are homologous to changes caused by agonists in other GPCRs. We have applied the "antibody imprint" method to light-activated rhodopsin in native membranes by using nine monoclonal antibodies (mAbs) against aqueous faces of rhodopsin. Epitopes recognized by these mAbs were found by selection from random peptide libraries displayed on phage. A new computer algorithm, FINDMAP, was used to map the epitopes to discontinuous segments of rhodopsin that are distant in the primary sequence but are in close spatial proximity in the structure. The proximity of a segment of the N-terminal and the loop between helices VI and VIII found by FINDMAP is consistent with the X-ray structure of the dark-adapted rhodopsin. Epitopes to the cytoplasmic face segregated into two classes with different predicted spatial proximities of protein segments that correlate with different preferences of the antibodies for stabilizing the metarhodopsin I or metarhodopsin II conformations of light-excited rhodopsin. Epitopes of antibodies that stabilize metarhodopsin II indicate conformational changes from dark-adapted rhodopsin, including rearrangements of the C-terminal tail and altered exposure of the cytoplasmic end of helix VI, a portion of the C-3 loop, and helix VIII. As additional antibodies are subjected to antibody imprinting, this approach should provide increasingly detailed information on the conformation of light-excited rhodopsin and be applicable to structural studies of other challenging protein targets.