Thermodynamics of replacing an alpha-helical Pro residue in the P40S mutant of Escherichia coli thioredoxin

Escherichia coli thioredoxin is a 108 amino acid oxidoreductase and contains a single Met residue at position 37. The protein contains a long alpha-helical stretch between residues 32 and 49. The central residue of this helix, Pro40, has been replaced by Ser. The stabilities of the oxidized states of two proteins, the single mutant M37L and the double mutant M37L,P40S, have been characterized by differential scanning calorimetry (DSC) and also by a series of isothermal guanidine hydrochloride (GuHCl) melts in the temperature range of 277 to 333 K. The P40S mutation was found to stabilize the protein at all temperatures upto 340 K though both proteins had similar Tm values of about 356 K. At 298 K, the M37L,P40S mutant was found to be more stable than M37L by 1.5 kcal/mol. A combined analysis of GuHCl and calorimetric data was carried out to determine the enthalpy, entropy, and heat capacity change upon unfolding. At 298 K there was a large, stabilizing enthalpic effect in P40S though significant enthalpy-entropy compensation was observed and the two proteins had similar values of deltaCp. Thus, replacement of a Pro in the interior of an alpha helix can have substantial effects on protein stability.     

The asparagine-stabilized beta-turn of apamin: contribution to structural stability from dynamics simulation and amide hydrogen exchange analysis

Molecular dynamics simulations of bee venom apamin, and an analogue having an Asn to Ala substitution at residue 2 (apamin-N2A), were analyzed to explore the contribution of hydrogen bonds involving Asn2 to local (beta-turn residues N2, C3, K4, A5) and global stability. The wild-type peptide retained a stable conformation during 2.4 ns of simulation at 67 degrees C, with high beta-turn stability characterized by backbone-side chain hydrogen bonds involving beta-turn residues K4 and A5, with the N2 side chain amide carbonyl. The loss of stabilizing interactions involving the N2 side chain resulted in the loss of the beta-turn conformation in the apamin N2A simulations (27 or 67 degrees C). This loss of beta-turn stability propagates throughout the peptide structure, with destabilization of the C-terminal helix connected to the N-terminal region by two disulfide bonds. Backbone stability in a synthetic peptide analogue (apamin-N2A) was characterized by NMR and amide hydrogen exchange measurements. Consistent with the simulations, loss of hydrogen bonds involving the N2 side chain resulted in destabilization of both the N-terminal beta-turn and the C-terminal helix. Amide exchange protection factors in the C-terminal helix were reduced by 9-11-fold in apamin N2A as compared with apamin, corresponding to free energy (deltaDeltaG(uf)) of around 1.5 kcal M(-1) at 20 degrees C. This is equivalent to the contribution of hydrogen bond interactions involving the N2 side chain to the stability of the beta-turn. Together with additional measures of exchange protection factors, the three main contributions to backbone stability in apamin that account for virtually the full thermodynamic stability of the peptide have been quantitated.     

Crystallographic and functional analysis of the ESCRT-I /HIV-1 Gag PTAP interaction

Budding of HIV-1 requires the binding of the PTAP late domain of the Gag p6 protein to the UEV domain of the TSG101 subunit of ESCRT-I. The normal function of this motif in cells is in receptor downregulation. Here, we report the 1.4-1.6?? structures of the human TSG101 UEV domain alone and with wild-type and mutant HIV-1 PTAP and Hrs PSAP nonapeptides. The hydroxyl of the Thr or Ser residue in the P(S/T)AP motif hydrogen bonds with the main chain of Asn69. Mutation of the Asn to Pro, blocking the main-chain amide, abrogates PTAP motif binding in?vitro and blocks budding of HIV-1 from cells. N69P and other PTAP binding-deficient alleles of TSG101 did not rescue HIV-1 budding. However, the mutant alleles did rescue downregulation of endogenous EGF receptor. This demonstrates that the PSAP motif is not rate determining in EGF receptor downregulation under normal conditions.     

Modifications of the active center of T4 thioredoxin by site-directed mutagenesis

The active site sequence of T4 thioredoxin, Cys-Val-Tyr-Cys, has been modified in two positions to Cys-Gly-Pro-Cys to mimic that of Escherichia coli thioredoxin. The two point mutants Cys-Gly-Tyr-Cys and Cys-Val-Pro-Cys have also been constructed. The mutant proteins have similar reaction rates with T4 ribonucleotide reductase as has the wild-type T4 thioredoxin. Mutant T4 thioredoxins with Pro instead of Tyr at position 16 in the active site sequence have three to four times lower apparent KM with E. coli ribonucleotide reductase than wild-type T4 thioredoxin. The KM values for these mutant proteins which do not have Tyr in position 16 are thus closer to E. coli thioredoxin than to the wild-type T4 thioredoxin. The bulky tyrosine side chain probably prevents proper interactions to E. coli ribonucleotide reductase. Also the redox potentials of these two mutant thioredoxins are lower than that of the wild-type T4 thioredoxin and are thereby more similar to the redox potential of E. coli thioredoxin. Mutations in position 15 behave more or less like the wild-type protein. The kinetic parameters with E. coli thioredoxin reductase are similar for wild-type and mutant T4 thioredoxins except that the apparent kcat is lower for the mutant protein with Pro instead of Tyr in position 16. The active site sequence of T4 thioredoxin has also been changed to Cys-Pro-Tyr-Cys to mimic that of glutaredoxins. This change does not markedly alter the reaction rate of the mutant protein with T4 ribonucleotide reductase or E. coli thioredoxin reductase, but the redox potential is lower for this mutant protein than for wild-type T4 thioredoxin.     

Hydrogen bonds between short polar side chains and peptide backbone: prevalence in proteins and effects on helix-forming propensities

A survey of 322 proteins showed that the short polar (SP) side chains of four residues, Thr, Ser, Asp, and Asn, have a very strong tendency to form hydrogen bonds with neighboring backbone amides. Specifically, 32% of Thr, 29% of Ser, 26% of Asp, and 19% of Asn engage in such hydrogen bonds. When an SP residue caps the N terminal of a helix, the contribution to helix stability by a hydrogen bond with the amide of the N3 or N2 residue is well established. When an SP residue is in the middle of a helix, the side chain is unlikely to form hydrogen bonds with neighboring backbone amides for steric and geometric reasons. In essence the SP side chain competes with the backbone carbonyl for the same hydrogen-bonding partner (i.e., the backbone amide) and thus SP residues tend to break backbone carbonyl-amide hydrogen bonds. The proposition that this is the origin for the low propensities of SP residues in the middle of alpha helices (relative to those of nonpolar residues) was tested. The combined effects of restricting side-chain rotamer conformations (documented by Creamer and Rose, Proc Acad Sci USA, 1992;89:5937-5941; Proteins, 1994;19:85-97) and excluding side- chain to backbone hydrogen bonds by the helix were quantitatively analyzed. These were found to correlate strongly with four experimentally determined scales of helix-forming propensities. The correlation coefficients ranged from 0.72 to 0.87, which are comparable to those found for nonpolar residues (for which only the loss of side-chain conformational entropy needs to be considered).     

Crystal structure of human thioredoxin revealing an unraveled helix and exposed S-nitrosation site

Thioredoxins reduce disulfide bonds and other thiol modifications in all cells using a CXXC motif. Human thioredoxin 1 is unusual in that it codes for an additional three cysteines in its 105 amino acid sequence, each of which have been implicated in other reductive activities. Cys 62 and Cys 69 are buried in the protein interior and lie at either end of a short helix (helix 3), and yet can disulfide link under oxidizing conditions. Cys 62 is readily S‐nitrosated, giving rise to a SNO modification, which is also buried. Here, we present two crystal structures of the C69S/C73S mutant protein under oxidizing (1.5 Å) and reducing (1.1 Å) conditions. In the oxidized structure, helix 3 is unraveled and displays a new conformation that is stabilized by a series of new hydrogen bonds and a disulfide link with Cys 62 in a neighboring molecule. The new conformation provides an explanation for how a completely buried residue can participate in SNO exchange reactions.     

Structures of wild-type and mutant signal sequences of Escherichia coli ribose binding protein.

The structure of a chemically synthesized 25-residue-long functional signal peptide of Escherichia coli ribose binding protein was compared with that of a nonfunctional mutant-signal peptide using circular dichroism and two-dimensional 1H NMR in solvents mimicking the amphiphilic environments. The functional peptide forms an 18-residue-long alpha-helix starting from the NH2-terminal region and reaching to the hydrophobic stretch in a solvent consisting of 10% dimethylsulfoxide, 40% water, and 50% trifluoroethanol (v/v). The nonfunctional mutant peptide, which contains a Pro at position 9 instead of a Leu in the wild-type peptide, does not have any secondary structure in that solvent but forms a 12-residue-long alpha-helix within the hydrophobic stretch in water/trifluoroethanol (50:50, v/v) solvent. It seems that the Pro-9 residue in the nonfunctional peptide disturbs the helix propagation from the hydrophobic stretch to the NH2-terminal region. Because both of these peptides have stable helices within the hydrophobic stretch, it may be concluded that the additional 2 turns of the alpha-helix in the NH2-terminal region of the wild-type signal peptide is important for its function.     

High-resolution structure of the temperature-sensitive mutant of phage lysozyme, Arg 96----His

The structure of the temperature-sensitive mutant lysozyme of bacteriophage T4 in which arginine 96 is replaced by histidine has been determined crystallographically and refined to a residual of 17.6% at 1.9-A resolution. Overall, the three-dimensional structure of the mutant protein is extremely similar to that of wild type. There are local distortions in the mutant structure suggesting that the substituted His 96 residue is under strain. This appears to be one of the major reasons for the decreased thermostability. In wild-type lysozyme the guanidinium of Arg 96 is located at the carboxy terminus of alpha-helix 82-90 and makes a pair of hydrogen bonds to two of the carbonyl groups in the last turn of the helix. The loss of this "helix dipole" interaction also appears to contribute to the destabilization. The pKa* of His 96 in the mutant lysozyme has been determined by nuclear magnetic resonance and found to be 6.8 at 10 degrees C. This relatively normal value of the histidine pKa* suggests that the protonated and unprotonated forms of the imidazole ring are perturbed equally by the protein environment or, what is equivalent, the mutant lysozyme is equally stable with either histidine species.     

Crystal structure and molecular dynamics simulation of ubiquitin-like domain of murine parkin

Parkin is the gene product identified as the major cause of autosomal recessive juvenile Parkinsonism (AR-JP). Parkin, a ubiquitin ligase E3, contains a unique ubiquitin-like domain in its N-terminus designated Uld which is assumed to be a interaction domain with the Rpn 10 subunit of 26S proteasome. To elucidate the structural and functional role of Uld in parkin at the atomic level, the X-ray crystal structure of murine Uld was determined and a molecular dynamics simulation of wild Uld and its five mutants (K27N, R33Q, R42P, K48A and V56E) identified from AR-JP patients was performed. Murine Uld consists of two alpha helices [Ile23-Arg33 (alpha1) and Val56-Gln57 (alpha2)] and five beta strands [Met1-Phe7 (beta1), Tyr11-Asp18 (beta2), Leu41-Phe45 (beta3), Lys48-Pro51 (beta4) and Ser65-Arg72 (beta5)] and its overall structure is essentially the same as that of human ubiquitin with a 1.22 A rmsd for the backbone atoms of residues 1-76; however, the sequential identity and similarity between both molecules are 32% and 63%, respectively. This close resemblance is due to the core structure built by same hydrogen bond formations between and within the backbone chains of alpha1 and beta1/2/5 secondary structure elements and by nearly the same hydrophobic interactions formed between the nonpolar amino acids of their secondary structures. The side chain NetaH of Lys27 on the alpha1 helix was crucial to the stabilization of the spatial orientations of beta3 and beta4 strands, possible binding region with Rpn 10 subunit, through three hydrogen bonds. The MD simulations showed the K27N and R33Q mutations increase the structural fluctuation of these beta strands including the alpha1 helix. Reversely, the V56E mutant restricted the spatial flexibility at the periphery of the short alpha2 helix by the interactions between the polar atoms of Glu56 and Ser19 residues. However, a large fluctuation of beta4 strand with respect to beta5 strand was induced in the R42P mutant, because of the impossibility of forming paired hydrogen bonds of Pro for Arg42 in wild Uld. The X-ray structure showed that the side chains of Asp39, Gln40 and Arg42 at the N-terminal periphery of beta3 strand protrude from the molecular surface of Uld and participate in hydrogen bonds with the polar residues of neighboring Ulds. Thus, the MD simulation suggests that the mutation substitution of Pro for Arg42 not only causes the large fluctuation of beta3 strand in the Uld but also leads to the loss of the ability of Uld to trap the Rpn 10 subunit. In contrast, the MD simulation of K48A mutant showed little influence on the beta3-beta4 loop structure, but a large fluctuation of Lys48 side chain, suggesting the importance of flexibility of this side chain for the interaction with the Rpn 10 subunit. The present results would be important in elucidating the impaired proteasomal binding mechanism of parkin in AR-JP.     

Studies of the mechanism of phenol hydroxylase: effect of mutation of proline 364 to serine

An active site residue in phenol hydroxylase (PHHY), Pro364, was mutated to serine to investigate its role in enzymatic catalysis. In the presence of phenol, the reaction between the reduced flavin of P364S and oxygen is very fast, but only 13% of the flavin is utilized to hydroxylate the substrate, compared to nearly 100% for the wild-type enzyme. The oxidative half-reaction of PHHY using m-cresol as a substrate is similarly affected by the mutation. Pro364 was suggested to be important in stabilizing the transition state of the oxygen transfer step by forming a hydrogen bond between its carbonyl oxygen and the C4a-hydroperoxyflavin [Ridder, L., Mullholland, A. J., Rietjens, I. M. C. M., and Vervoort, J. (2000) J. Am. Chem. Soc. 122, 8728-8738]. The P364S mutation may weaken this interaction by increasing the flexibility of the peptide chain; hence, the transition state would be destabilized to result in a decreased level of hydroxylation of phenol. However, when the oxidative half-reaction was studied using resorcinol as a substrate, the P364S mutant form was not significantly different from the wild-type enzyme. The rate constants for all the reaction steps as well as the hydroxylation efficiency (coupling between NADPH oxidation and resorcinol consumption) are comparable to those of the wild-type enzyme. It is suggested that the function of Pro364 in catalysis, stabilization of the transition state, is not as important in the reaction with resorcinol, possibly because the position of hydroxylation is different with resorcinol than with phenol and m-cresol.     

Role of proline residues in the structure and function of a membrane transport protein

By use of site-directed mutagenesis, each prolyl residue in the lac permease of Escherichia coli at positions 28 (putative helix I), 31 (helix I), 61 (helix II), 89 (helix III), 97 (helix III), 123 (helix IV), 192 (putative hydrophilic region 7), 220 (helix VII), 280 (helix VIII), and 327 [helix X; Lolkema, J. S., et al. (1988) Biochemistry 27, 8307] was systematically replaced with Gly, Ala, or Leu or deleted by truncation of the C-terminus [i.e., Pro403 and Pro405; Roepe, P.D., et al. (1989) Proc. Natl. Acad. Sci. U.S.A. 86, 3992]. Replacements were chosen on the basis of side-chain helical propensity: Gly, like Pro, is thought to be a "helix breaker", while Ala and Leu are "helix makers". With the exception of Pro28, each prolyl residue can be replaced with Gly or Ala, and Pro403 and -405 can be deleted with the C-terminal tail, and significant lac permease activity is retained. In contrast, when Pro28 is replaced with Gly, Ala, or Ser, lactose transport is abolished, but permease with Ser28 binds p-nitrophenyl alpha-D-galactopyranoside and catalyzes active transport of beta-galactopyranosyl-1-thio-beta-D- galactopyranoside. Replacement of Pro28, -31, -123, -280, or -327 with Leu abolishes lactose transport, while replacement of Pro61, -89, -97, or -220 with Leu has relatively minor effects. None of the alterations in permease activity is due to inability of the mutant proteins to insert into the membrane or to diminished lifetimes after insertion, since the concentration of each mutant permease in the membrane is comparable to that of wild-type permease as judged by immunological analyses.(ABSTRACT TRUNCATED AT 250 WORDS)     

Intramolecular catalysis of a proline isomerization reaction in the folding of dihydrofolate reductase.

The cis/trans isomerization of the peptide bond preceding proline residues in proteins can limit the rate at which a protein folds to its native conformation. Mutagenic analyses of dihydrofolate reductase (DHFR) from Escherichia coli show that this isomerization reaction can be intramolecularly catalyzed by a side chain from an amino acid which is distant in sequence but adjacent in the native conformation. The guanidinium NH2 nitrogen of Arg 44 forms one hydrogen bond to the imide nitrogen and a second to the carbonyl oxygen of Pro 66 in wild-type DHFR. Replacement of Arg 44 with Leu results in a change of the nature of the two slow steps in refolding from being limited by the acquisition of secondary and/or tertiary structure to being limited by isomerization. The simultaneous replacement of Pro 66 with Ala (i.e., the Leu 44/Ala 66 double mutant) eliminates this isomerization reaction and once again makes protein folding the limiting process. Apparently, one or both of the hydrogen bonds between Arg 44 and Pro 66 accelerate the isomerization of the Gln 65-Pro 66 peptide bond. The replacement of Arg 44 with Leu affects the kinetics of the slow folding reactions in a fashion which indicates that the crucial hydrogen bonds form in the transition states for the rate-limiting steps in folding.     

The P23T cataract mutation causes loss of solubility of folded gammaD-crystallin

Mutations in the human gammaD-crystallin gene have been linked to several types of congenital cataracts. In particular, the Pro23 to Thr (P23T) mutation of human gammaD crystallin has been linked to cerulean, lamellar, coralliform, and fasciculiform congenital cataracts. We have expressed and purified wild-type human gammaD, P23T, and the Pro23 to Ser23 (P23S) mutant. Our measurements show that P23T is significantly less soluble than wild-type human gammaD, with P23S having an intermediate solubility. Using synchrotron radiation circular dichroism spectroscopy, we have determined that the P23T mutant has a slightly increased content of beta-sheet, which may be attributed to the extension of an edge beta-strand due to the substitution of Pro23 with a residue able to form hydrogen bonds. Neither of the point mutations appears to have reduced the thermal stability of the protein significantly, nor its resistance to guanidine hydrochloride-induced unfolding. These results suggest that insolubility, rather than loss of stability, is the primary basis for P23T congenital cataracts.     

Structure of oxidized bacteriophage T4 glutaredoxin (thioredoxin). Refinement of native and mutant proteins.

The structure of wild-type bacteriophage T4 glutaredoxin (earlier called thioredoxin) in its oxidized form has been refined in a monoclinic crystal form at 2.0 A resolution to a crystallographic R-factor of 0.209. A mutant T4 glutaredoxin gives orthorhombic crystals of better quality. The structure of this mutant has been solved by molecular replacement methods and refined at 1.45 A to an R-value of 0.175. In this mutant glutaredoxin, the active site residues Val15 and Tyr16 have been substituted by Gly and Pro, respectively, to mimic that of Escherichia coli thioredoxin. The main-chain conformation of the wild-type protein is similar in the two independently determined molecules in the asymmetric unit of the monoclinic crystals. On the other hand, side-chain conformations differ considerably between the two molecules due to heterologous packing interactions in the crystals. The structure of the mutant protein is very similar to the wild-type protein, except at mutated positions and at parts involved in crystal contacts. The active site disulfide bridge between Cys14 and Cys17 is located at the first turn of helix alpha 1. The torsion angles of these residues are similar to those of Escherichia coli thioredoxin. The torsion angle around the S-S bond is smaller than that normally observed for disulfides: 58 degrees, 67 degrees and 67 degrees for wild-type glutaredoxin molecule A and B and mutant glutaredoxin, respectively. Each sulfur atom of the disulfide cysteines in T4 glutaredoxin forms a hydrogen bond to one main-chain nitrogen atom. The active site is shielded from solvent on one side by the beta-carbon atoms of the cysteine residues plus side-chains of residues 7, 9, 21 and 33. From the opposite side, there is a cleft where the sulfur atom of Cys14 is accessible and can be attacked by a nucleophilic thiolate ion in the initial step of the reduction reaction.     

Thermodynamic effects of replacements of Pro residues in helix interiors of maltose-binding protein

Introduction of Pro residues into helix interiors results in protein destabilization. It is currently unclear if the converse substitution (i.e., replacement of Pro residues that naturally occur in helix interiors would be stabilizing). Maltose-binding protein is a large 370-amino acid protein that contains 21 Pro residues. Of these, three nonconserved residues (P48, P133, and P159) occur at helix interiors. Each of the residues was replaced with Ala and Ser. Stabilities were characterized by differential scanning calorimetry (DSC) as a function of pH and by isothermal urea denaturation studies as a function of temperature. The P48S and P48A mutants were found to be marginally more stable than the wild-type protein. In the pH range of 5-9, there is an average increase in T(m) values of P48A and P48S of 0.4 degrees C and 0.2 degrees C, respectively, relative to the wild-type protein. The other mutants are less stable than the wild type. Analysis of the effects of such Pro substitutions in MBP and in three other proteins studied to date suggests that substitutions are more likely to be stabilizing if the carbonyl group i-3 or i-4 to the mutation site is not hydrogen bonded in the wild-type protein.     

Effect of the N2 residue on the stability of the α-helix for all 20 amino acids

N2 is the second position in the alpha-helix. All 20 amino acids were placed in the N2 position of a synthetic helical peptide (CH(3)CO-[AXAAAAKAAAAKAAGY]-NH(2)) and the helix content was measured by circular dichroism spectroscopy at 273K. The dependence of peptide helicity on N2 residue identity has been used to determine a free-energy scale by analysis with a modified Lifson-Roig helix-coil theory that includes a parameter for the N2 energy (n2). The rank order of DeltaDeltaG((relative to Ala)) is Glu(-), Asp(-) > Ala > Glu(0), Leu, Val, Gln, Thr, Ile, Ser, Met, Asp(0), His(0), Arg, Cys, Lys, Phe > Asn, > Gly, His(+), Pro, Tyr. The results correlate very well with N2 propensities in proteins, moderately well with N1 and helix interior preferences, and not at all with N-cap preferences. The strongest energetic effects result from interactions with the helix dipole, which favors negative charges at the helix N terminus. Hydrogen bonds to side chains at N2, such as Gln, Ser, and Thr, are weak, despite occurring frequently in protein crystal structures, in contrast to the N-cap position. This is because N-cap hydrogen bonds are close to linear, whereas N2 hydrogen bonds have poor geometry. These results can be used to modify protein stability rationally, help design helices, and improve prediction of helix location and stability.     

Multiple binding modes for palmitate to barley lipid transfer protein facilitated by the presence of proline 12

Molecular dynamics simulations have been used to characterise the binding of the fatty acid ligand palmitate in the barley lipid transfer protein 1 (LTP) internal cavity. Two different palmitate binding modes (1 and 2), with similar protein–ligand interaction energies, have been identified using a variety of simulation strategies. These strategies include applying experimental protein–ligand atom–atom distance restraints during the simulation, or protonating the palmitate ligand, or using the vacuum GROMOS 54B7 force‐field parameter set for the ligand during the initial stages of the simulations. In both the binding modes identified the palmitate carboxylate head group hydrogen bonds with main chain amide groups in helix A, residues 4 to 19, of the protein. In binding mode 1 the hydrogen bonds are to Lys 11, Cys 13, and Leu 14 and in binding mode 2 to Thr 15, Tyr 16, Val 17, Ser 24 and also to the OH of Thr 15. In both cases palmitate binding exploits irregularity of the intrahelical hydrogen‐bonding pattern in helix A of barley LTP due to the presence of Pro 12. Simulations of two variants of barley LTP, namely the single mutant Pro12Val and the double mutant Pro12Val Pro70Val, show that Pro 12 is required for persistent palmitate binding in the LTP cavity. Overall, the work identifies key MD simulation approaches for characterizing the details of protein–ligand interactions in complexes where NMR data provide insufficient restraints.     

Controls exerted by the Aib residue: helix formation and helix reversal

The helix forming properties of the achiral alpha-amino isobutyric residue (Aib) have been demonstrated by numerous crystal structure analyses of designed and naturally occurring peptides containing one or more Aib residues in the sequence. Experimental and computational results concerning the type of helix obtained, whether the 3(10)-helix with 4 --> 1 type hydrogen bonds or the alpha-helix with 5 --> 1 hydrogen bonds or mixtures of the two, have been published. This paper deals with residues that, if inserted into a sequence, could perturb the helix-forming propensity afforded by the presence of Aib residues. Examples of structures will be presented in which Pro, Hyp, Gly-Gly, d-Ala-Gly, and Lac have been centrally placed in the sequence. In addition to the formation of helices, detailed experimentally obtained conformation information is presented for the role of the Aib residue in reversing the sense of the helix (the Schellman motif) with the consequent formation of the 6 --> 1 type hydrogen bond or a solvated 6 --> 1 hydrogen bond. Data are presented for 13 molecules with helix reversals at the C-terminus or near the center of the sequence.     

Sequence determinants of the capping box, a stabilizing motif at the N-termini of alpha-helices.

The capping box, a recurrent hydrogen bonded motif at the N-termini of alpha-helices, caps 2 of the initial 4 backbone amide hydrogen donors of the helix (Harper ET, Rose GD, 1993, Biochemistry 32:7605-7609). In detail, the side chain of the first helical residue forms a hydrogen bond with the backbone of the fourth helical residue and, reciprocally, the side chain of the fourth residue forms a hydrogen bond with the backbone of the first residue. We now enlarge the earlier definition of this motif to include an accompanying hydrophobic interaction between residues that bracket the capping box sequence on either side. The expanded box motif--in which 2 hydrogen bonds and a hydrophobic interaction are localized within 6 consecutive residues--resembles a glycine-based capping motif found at helix C-termini (Aurora R, Srinivasan R, Rose GD, 1994, Science 264:1126-1130).     

NMR studies of structure, hydrogen exchange, and main-chain dynamics in a disrupted-core mutant of thioredoxin.

Core-packing mutants of proteins often approach molten globule states, and hence may have attributes of folding intermediates. We have studied a core-packing mutant of thioredoxin, L78K, in which a leucine residue is substituted by lysine, using 15N heteronuclear two- and three-dimensional NMR. Chemical shift differences between the mutant and wild-type main-chain resonances reveal that structural changes caused by the mutation are localized within 12 A of the altered side chain. The majority of resonances are unchanged, as are many 1H-1H NOEs indicative of the main-chain fold, suggesting that the structure of L78K is largely similar to wild type. Hydrogen exchange studies reveal that residues comprising the central beta-sheet of both mutant and wild-type proteins constitute a local unfolding unit, but with the unfolding/folding equilibrium approximately 12 times larger in L78K. The dynamics of main-chain NH bonds in L78K were studied by 15N spin relaxation and compared with a previous study of wild type. Order parameters for angular motion of NH bonds in the mutant are on average lower than in wild type, suggesting greater spatial freedom on a rapid time scale, but may also be related to different rotational correlation times in the two proteins. There is also evidence of greater conformational exchange in the mutant. Differences between mutant and wild type in hydrogen exchange and main-chain dynamics are not confined to the vicinity of the mutation. We infer that mispacking of the protein core in one location affects local dynamics and stability throughout.