Beta-sheet aggregation proposed in sickle cell hemoglobin

A β-sheet conformation is predicted at the N-terminal of β chains in sickle cell hemoglobin (Hb S) as a result of the β6 Glu → Val mutation. Since Glu is the weakest and Val is the strongest β-sheet former in the predictive method of Chou and Fasman [Biochemistry $\text{13}$, 211, 222 (1974)], such a substitution greatly increases the β-sheet potential in the β 1–6 region. The similarity in the concentration and temperature dependence of Hb S gelation to β-sheet formation in polyamino acids suggest that a common aggregation mechanism may be involved. Conditions to cause a β → α trans-formation at the β 1–6 region of Hb S is discussed relative to the treatment of sickle cell disease.     

Infrared study of the secondary structure of lac repressor and its tryptic core

The lac repressor and its tryptic core were studied by ir spectroscopy, and their β-structure content was determined by analysis of the spectra. Using protein-derived reference spectra, we find a β-content for lac repressor of 18% and of 23% for its tryptic core. The higher amount of β-conformation in the tryptic core is confirmed by another type of analysis (decomposition of the spectra in Gaussian curves). These results are discussed with respect to their implications for the structure of the N-terminal “headpiece” of lac repressor and for the mode of interaction of lac repressor with lac operator.     

Conformational structure of bombesin as studied by vibrational and circular dichroism spectroscopy

Raman and Fourier transform infrared (FTIR) spectroscopies and circular dichroism (CD) have been applied to investigate the secondary structure of bombesin in the solid state and in phosphate buffer solution (pH 3.8). At concentrations around 10⁻⁵ M, circular dichroism reveals that bombesin exists as an irregular or disordered conformation. However, the secondary structure of the peptide appears to be a mixture of disordered structure and intermolecular β-sheets in 0.01 M sodium phosphate buffer when the peptide concentrations are higher than around 6.5 mM. The tendency of bombesin to form aggregated β-sheet species seems to be originated mainly in the sequence of the residues 7–14, as supported by the Raman spectra and β-sheet propensities (P_β) of the amino-acid residues. It is the hydrophobic force of this amino-acid sequence, and not a salt bridge effect, that is the factor responsible for the formation of peptide aggregates.     

Analysis and prediction of the packing of α-helices against a β-sheet in the tertiary structure of globular proteins

The packing of α-helices and β-sheets in six $\text{α}\text{β}$ proteins (e.g. flavodoxin) has been analysed. The results provide the basis for a computer algorithm to predict the tertiary structure of an $\text{α}\text{β}$ protein from its amino acid sequence and actual assignment of secondary structure.The packing of an individual α-helix against a β-sheet generally involves two adjacent ± 4 rows of non-polar residues on the α-helix at the positions i, i + 4, i + 8, i + 1, i + 5, i + 9. The pattern of interacting β-sheet residues results from the twisted nature of the sheet surface and the attendant rotation of the side-chains. At a more detailed level, four of the α-helical residues (i + 1, i + 4, i + 5 and i + 8) form a diamond that surrounds one particular β-sheet residue, generally isoleucine, leucine or valine. In general, the α-helix sits 10 Å above the sheet and lies parallel to the strand direction.The prediction follows a combinational approach. First, a list of possible β-sheet structures (10⁶ to 10¹⁴) is constructed by the generation of all β-sheet topologies and β-strand alignments. This list is reduced by constraints on topology and the location of non-polar residues to mediate the sheet/helix packing, and then rank-ordered on the extent of hydrogen bonding. This algorithm was uniformly applied to 16 $\text{α}\text{β}$ domains in 13 proteins. For every structure, one member of the reduced list was close to the crystal structure; the root-mean-square deviation between equivalenced C^α atoms averaged 5.6 Å for 100 residues. For the $\text{α}\text{β}$ proteins with pure parallel β-sheets, the total number of structures comparable to or better than the native in terms of hydrogen bonds was between 1 and 148. For proteins with mixed β-sheets, the worst case is glyceraldehyde-3-phosphate dehydrogenase, where as many as 3800 structures would have to be sampled. The evolutionary significance of these results as well as the potential use of a combinatorial approach to the protein folding problem are discussed.     

Circular dichroism of the parallel β helical proteins pectate lyase C and E

The pectate lyases, PelC and PelE, have an unusual folding motif, known as a parallel β-helix, in which the polypeptide chain is coiled into a larger helix composed of three parallel β-sheets connected by loops having variable lengths and conformations. Since the regular secondary structure consists almost entirely of parallel β-sheets these proteins provide a unique opportunity to study the effect of parallel β-helical structure on circular dichroism (CD). We report here the CD spectra of PelC and PelE in the presence and absence of Ca²⁺, derive the parallel β-helical components of the spectra, and compare these results with previous CD studies of parallel β-sheet structure. The shape and intensity of the parallel β-sheet spectrum is distinctive and may be useful in identifying other proteins that contain the parallel β-helical folding motif. © 1995 Wiley-Liss, Inc.     

Studied on gene control regions IX. The effect of hypoxanthine-substituted lac operators on the lac operator--lac repressor interaction.

Hypoxanthine was substituted for guanine at specific sites in the lac operator DNA by a combination of chemical and enzymatic procedures. The stability of these modified lac operators with wild-type (SQ) and tight binding (QX86) lac repressors was measured. Effects were variable. At some sites insertion of hypoxanthine significantly reduced the stability of the complex whereas at other sites substitution with hypoxanthine did not alter the repressor—operator interaction. In addition, insertion of this analog at two sites increased the stability of the complexes. These changes were used to partially map regions of the lac operator that are in contact with lac represser. The results suggest that lac repressor recognizes the guanine 2-amino group at specific sites in the minor groove of lac operator.     

Simulation Studies on the Stabilities of Aggregates Formed by Fibril-Forming Segments of α-Synuclein

Abstract

We performed molecular dynamics simulations for various oligomers with different β-sheet conformations consisting of α-Synuclein 71–82 residues using an all atom force field and explicit water model. Tetramers of antiparallel β-sheet are shown to be stable, whereas parallel sheets are highly unstable due to the repulsive interactions between bulky and polar side chains as well as the weaker backbone hydrogen bonds. We also investigated the stabilities of double antiparallel β-sheets stacked with asymmetric and symmetric geometries. Our results show that this 12 amino acid residue peptide can form stable β-sheet conformers at 320K and higher temperatures. The backbone hydrogen bonds in β-sheet and the steric packing between hydrophobic side chains between β-sheets are shown to give conformational stabilities.     

The possible roles of residues 79 and 80 of the Trp repressor from Escherichia coli K-12 in trp operator recognition

We constructed mutants of the Trp repressor from Escherichia coli K-12 with all possible single amino acid exchanges at positions 79 and 80 (residues 1 and 2 of the recognition helix). We tested these mutants in vivo by measuring the repression of synthesis of β-galactosidase with symmetric variants of α- and β-centered trp operators, which replace the lac operator in a synthetic lac system. The Trp repressor carrying a substitution of isoleucine 79 by lysine, showed a marked specificity change with respect to base pair 7 of the α-centered trp operator. Gel retardation experiments confirmed this result. Trp repressor mutant IR79 specifically recognizes a trp operator variant with substitutions in positions 7 and 8. Another mutant, with glycine in position 79, exhibited loss of contact at base pair 7. We speculate that the side chain of Ile79 interacts with the AT base pairs 7 and 8 of the α-centered trp operator, possibly with the methyl groups of thymines. Replacement of thymine in position 7 or 8 by uracil confirms the involvement of the methyl group of thymine 8 in repressor binding. Several Trp repressor mutants in position 80 (i.e. AI80, AL80, AM80 and AP80) broaden the specificity of the Trp repressor for α-centered trp operator variants with exchanges in positions 3, 4 and 5.     

Papain Does Not Cleave Operator-Bound Lambda Repressor: Structural Characterization of the Carboxy Terminal Domain and the Hinge

Abstract

The circular dichroism spectra of three different purified carboxy terminal fragments 93–236, 112–236 and 132–236 of the bacteriophage γ cI repressor have been measured and compared with those of the intact repressor and the amino terminal fragment 1–92. All three carboxy terminal fragments contain mostly β-strands and loops, a minor helix content increasing with the size of the fragment, showing that the 93–131 region previously called a hinge is structured. Fourier transformed infrared spectra also showed that fragment 93–236 contains α-helices, β-sheets and turns but fragment 132–236 contains no detectable α-helix, only β-sheets and turns. Papain is known to cleave the γ repressor, but it is shown here that it cannot cleave the operator-bound repressor dimer. For the 132–236 fragment, both the wt and the SN228 mutant previously shown to be dimerization defective in the intact, gave similar dimerization properties as investigated by HPLC at 2 to 100 µM protein concentration, with a KD of 13.2 µM and 19.1 µM respectively. The papain cleavage for wt and SN228 proceed at equal rates for the first cleavage at 92–93; however, the subsequent cleavages are faster for SN228. The three Cys residues in the 132–236 fragment were found to be unreactive upon incubation with DTNB, indicating the thiol sulfur atoms are buried in the repressor carboxy terminal domain. Denaturation of the 132–236 fragment studied by tryptophan fluorescence shows two transitions centered at 1.5 M and 4.5 M of urea.     

Local hydrophobicity stabilizes secondary structures in proteins

The probability of occurrence of helix and β-sheet residues in 47 globular proteins was determined as a function of local hydrophobicity, which was defined by the sum of the Nozaki-Tanford transfer free energies at two nearest-neighbors on both sides of the amino acid sequence. In general, hydrophilic amino acids favor neither helix nor β-sheet formations when neighbor residues are also hydrophilic but favor helix formation at higher local hydrophobicity. On the other hand, some hydrophobic amino acids such as Met, Leu, and Ile favor helix formation when neighbor residues are hydrophilic. None of the hydrophobic amino acids favor β-sheet formation with hydrophilic neighbors, but most of them strongly favor β-sheet formation at high local hydrophobicity. When the average of 20 amino acids is taken, both helix and β-sheet residue probabilities are higher at higher local hydrophobicity, although the increase is steeper for β-sheets. Therefore, β-sheet formation is more influenced by local hydrophobicity than helix formation. Generally, helices are nearer the surface and tend to have hydrophilic and hydrophobic faces at opposite sides. The tendency of alternating regions of hydrophilic and hydrophobic residues in a helical sequence was revealed by calculating the correlation of the Nozaki-Tanford values. Such amphipathic helices may be important in protein–protein and protein–lipid interactions and in forming hydrophilic channels in the membrane. The choice of 30 nonhomologous proteins as the data set did not alter the above results.     

Logical analysis of the mechanism of protein folding III. Prediction of the strong long-range interactions.

It has been found that strong long-range interactions occur in regions having large β-structural potentials. As has been described previously (Nagano, 1974), interactions among regions having both helical and β-structural potentials (αβ-gaβ interactions) are also very important. Accordingly, an idea is presented in this paper that the relative stability of a protein conformation could be estimated by a relatively simple mathematical function of sequence and conformation. The function P(p,q) is called the non-energy part of pseudo-free energy, because minimization of the sum of P(p,q) and energy functions (cf. Levitt, 1974; Warme &; Scheraga, 1974) can be expected to lead to a plausible model of a protein. A merit of the function is that it can help us decide which way to go in manipulating a temporarily built model, e.g. towards a helix-rich protein or towards a β-structure-rich protein. The estimation of P(p,q) as an artificial potential does not use much computer time because only the co-ordinates of the β-carbon atoms (α-carbon atoms if the residue is Gly) are used. It is composed of terms of the long-range interactions P_L and short-range interactions P_S. The term P_L represents the relative strength of helix-helix interactions, helix-β-candidate interactions and β-candidate-β-candidate interactions. It is assumed that both helical and β-structural potentials can be measured as the differences between the predicted function for helix and β-structure, respectively, as defined previously (Nagano, 1973), and the corresponding largest values ever found. A hypothesis that two residues distantly separated in the primary sequence contribute less to the stability of the whole molecule is finally discarded because the true conformation of concanavalin A becomes very unstable compared with its false conformation folded like the main part of subtilisin. The parameters thus determined indicate that the helix-β-candidate interactions are almost as important as the β-candidate-β-candidate interactions. Both helix and loop prediction functions are combined to give the short-range interactions term, P_S, according to whether the region is really helical or not, and to whether it is really looped or not. The function P(p,q) can be used as a criterion for judging whether the predicted conformation is realistic or false, because the parameters can be adjusted to give, within limits, reasonable values of −10 kcal/residue for true conformations and higher than −5 kcal/residue for false conformations.As an application of the present theory of protein folding, the tertiary structure of bacteriophage T4 lysozyme is predicted and presented in Figure 1, prior to the X-ray structure becoming available.     

Circular dichroism and laser Raman spectroscopic analysis of the secondary structure ofCerebratulus lacteus toxin B-IV

The secondary structure ofCerebratulus lacteus toxin B-IV, a neurotoxic polypeptide containing 55 amino acid residues and four disulfide bonds, was experimentally estimated by computer analyses of toxin circular dichroism (CD) and laser Raman spectra. The CD spectrum of the toxin displayed typical α-helical peaks at 191, 208, and 222 nm. At neutralpH, the α-helix estimates from CD varied between 49 and 55%, when nonrepresentative spectrum analytical methods were used. Analysis of the laser Raman spectrum obtained at a much higher toxin concentration yielded a 78% α-helix estimate. Both CD and Raman spectroscopic methods failed to detect any β-sheet structure. The spectroscopic analyses revealed significantly more α-helix and less β-sheet for toxin B-IV than was predicted from its sequence. To account for the difference between the 49–55% helix estimate from CD spectra and the 78% helix estimate from the Raman spectrum, we postulate that some terminal residues are unfolded at the low toxin concentrations used for CD measurements but form helix at the high toxin concentration used for Raman measurements. Our CD observations showing thatCerebatulus toxin B-IV helix content increases about 15% in trifluoroethanol or at highpH are consistent with this interpretation.     

DNA recognition by β-sheets

The modes of DNA recognition by β-sheets are analyzed by using the known crystal and solution three-dimensional structures of DNA-protein complexes. Close fitting of the protein surface and the DNA surface determines the binding geometry. Interaction takes place so that essentially the N-to-C direction of the β-strands either follows or crosses the DNA groove. Upon following the major groove a two-stranded antiparallel β-sheet dives into the groove and contacts DNA bases with its convex side facing the DNA, while upon following the minor groove, it binds around the sugar-phosphate backbones, with its opposite concave side shielding the DNA. In order for the β-strands crossing the minor groove to interact with the DNA, the dinucleotide steps need to almost totally helically untwist and roll around major groove. The β-sheet, on the other hand, needs to adopt a concave curvature on the binding surface in the direction that follows the DNA minor groove, and a convex surface in the direction that bridges the sugar-phosphate backbones across the groove. The result is to produce a hyperbolic paraboloidal DNA-binding surface. © 1998 John Wiley & Sons, Inc. Biopoly 44: 335–359, 1997     

The influence of nearest-neighboring amino acid residues on aspects of secondary structure of proteins. Attempts to locate -helices and -sheets

The influence of nearest-neighbor pairs of amino acids (n ? 1) and (n + 1) on the conformation of amino acid (n) in proteins has been studied. From experimental data on eleven proteins of known three dimensional structures, our definition of an α-helical domain in the Φ,Ψ plot has been reexamined and found to be satisfactory. On the same principle, a regular β-sheet domain has been delineated. We then revised our 20 × 20 table of frequencies of occurrences of various conformations tabulating three values: α-helical, β-sheet, and neither. These frequencies were then used to locate the helixbreaking positions in cytochrome b₅, papain, thermolysin, and calcium-binding protein. In conjuction with the helical wheel method, they were useful for predicting the locations of most α-helical segments. Similarly the β-sheet breaking positions in papain were located and most of the β-sheets found by X-ray diffraction were close to or between these positions. Data on β-sheets are extremely sparse so that extensive tests were not possible. The application of this method to abnormal hemoglobins suggested possible distortions of helices and in several instances correlated with abnormal properties of the hemoglobins and association with disease. The variable region of human immunoglobin heavy chains was found to have a very low α-helical content though β-sheet structures might exist.     

X-ray diffraction studies of a silkmoth chorion

High- and low-angle diffraction studies have been performed on mature chorion (eggshell) of the silkmoth, Antheraea polyphemus. The results confirm the prevalence of β-sheet structure, previously suggested by predictions based on known primary structure and by results of laser Raman spectroscopy. The patterns obtained with different irradiation geometries suggest that a significant proportion of β-sheets are stacked and oriented with respect to the chorion surface and the ultrastructurally evident fibrillar components. Strong similarities are evident with the organization of β-sheets in chicken scale keratin.     

Solvent-induced tuning of internal structure in a protein amyloid protofibril

An important goal in studies of protein aggregation is to obtain an understanding of the structural diversity that is characteristic of amyloid fibril and protofibril structures at the molecular level. In this study, what to our knowledge are novel assays based on time-resolved fluorescence anisotropy decay and dynamic quenching measurements of a fluorophore placed at different specific locations in the primary structure of a small protein, barstar, have been used to determine the extent to which the protein sequence participates in the structural core of protofibrils. The fluorescence measurements reveal the structural basis of how modulating solvent polarity results in the tuning of the protofibril conformation from a pair of parallel β-sheets in heat-induced protofibrils to a single parallel β-sheet in trifluorethanol-induced protofibrils. In trifluorethanol-induced protofibrils, the single β-sheet is shown to be built up from in-register β-strands formed by nearly the entire protein sequence, while in heat-induced protofibrils, the pair of β-sheets motif is built up from β-strands formed by only the last two-third of the protein sequence.     

Interplay of Protein and DNA Structure Revealed in Simulations of the lac Operon

The E. coli Lac repressor is the classic textbook example of a protein that attaches to widely spaced sites along a genome and forces the intervening DNA into a loop. The short loops implicated in the regulation of the lac operon suggest the involvement of factors other than DNA and repressor in gene control. The molecular simulations presented here examine two likely structural contributions to the in-vivo looping of bacterial DNA: the distortions of the double helix introduced upon association of the highly abundant, nonspecific nucleoid protein HU and the large-scale deformations of the repressor detected in low-resolution experiments. The computations take account of the three-dimensional arrangements of nucleotides and amino acids found in crystal structures of DNA with the two proteins, the natural rest state and deformational properties of protein-free DNA, and the constraints on looping imposed by the conformation of the repressor and the orientation of bound DNA. The predicted looping propensities capture the complex, chain-length-dependent variation in repression efficacy extracted from gene expression studies and in vitro experiments and reveal unexpected chain-length-dependent variations in the uptake of HU, the deformation of repressor, and the folding of DNA. Both the opening of repressor and the presence of HU, at levels approximating those found in vivo, enhance the probability of loop formation. HU affects the global organization of the repressor and the opening of repressor influences the levels of HU binding to DNA. The length of the loop determines whether the DNA adopts antiparallel or parallel orientations on the repressor, whether the repressor is opened or closed, and how many HU molecules bind to the loop. The collective behavior of proteins and DNA is greater than the sum of the parts and hints of ways in which multiple proteins may coordinate the packaging and processing of genetic information.     

Complexes of Escherichia coli lac-repressor with non-operator DNA revealed by electron microscopy: two repressor molecules can share the same segment of DNA.

Complexes of Escherichia coli lac-repressor with non-operator DNA have been visualized in the electron microscope using high-resolution metal shadowing and negative staining. Under conditions of a high ratio of repressor to DNA, all the DNA molecules are covered by repressor molecules and the resulting complexes appear as flattened ribbons with a width of approximately 200 Å. The overall dimensions of these complexes and their substructure indicate that it is very likely that repressor molecules are tightly packed on both “sides” of the DNA helix. Thus two repressor molecules can share the same segment of non-operator DNA by binding to opposite sides of the DNA helix.