Who replaced Willie Mays in center field for the Giants or Mickey Mantle in the Yankee outfield? Who succeeded Hornung and Taylor in the backfield of the Green Bay Packers? For those of you who prefer a different type of “football,” who replaced Pele on the Brazilian team? Finally, who took over as Editor of Protein Science after Hans Neurath?

The transmembrane glycoprotein gp130 is the common signal transducing receptor subunit of the interleukin-6-type cytokines. It is a member of the cytokine-receptor superfamily predicted to consist of six domains in its extracellular part. The second and third domain constitute the cytokine-binding module defined by a set of four conserved cysteines and a WSXWS motif, respectively. The three-dimensional structure of the carboxy-terminal domain of this region was determined by multidimensional NMR. The domain consists of seven β-strands constituting a fibronectin type III-like topology. The structure reveals that the WSDWS motif of gp130 is part of an extended tryptophan/arginine zipper which modulates the conformation of the CD loop.

The complete amino acid sequence of the lectin KM+ from Artocarpus integrifolia (jackfruit), which contains 149 residues/mol, is reported and compared to those of other members of the Moraceae family, particularly that of jacalin, also from jackfruit, with which it shares 52% sequence identity. KM+ presents an acetyl-blocked N-terminus and is not posttranslationally modified by proteolytic cleavage as is the case for jacalin. Rather, it possesses a short, glycine-rich linker that unites the regions homologous to the α- and β-chains of jacalin. The results of homology modeling implicate the linker sequence in sterically impeding rotation of the side chain of Asp141 within the binding site pocket. As a consequence, the aspartic acid is locked into a conformation adequate only for the recognition of equatorial hydroxyl groups on the C4 epimeric center (α-d-mannose, α-d-glucose, and their derivatives). In contrast, the internal cleavage of the jacalin chain permits free rotation of the homologous aspartic acid, rendering it capable of accepting hydrogen bonds from both possible hydroxyl configurations on C4. We suggest that, together with direct recognition of epimeric hydroxyls and the steric exclusion of disfavored ligands, conformational restriction of the lectin should be considered to be a new mechanism by which selectivity may be built into carbohydrate binding sites. Jacalin and KM+ adopt the β-prism fold already observed in two unrelated protein families. Despite presenting little or no sequence similarity, an analysis of the β-prism reveals a canonical feature repeatedly present in all such structures, which is based on six largely hydrophobic residues within a β-hairpin containing two classic-type β-bulges. We suggest the term β-prism motif to describe this feature.

NMR spectroscopic analysis of the C-terminal Kunitz domain fragment (α3(VI)) from the human α3-chain of type VI collagen has revealed that the side chain of Trp21 exists in two unequally populated conformations. The major conformation (M) is identical to the conformation observed in the X-ray crystallographic structure, while the minor conformation (m) cannot structurally be resolved in detail by NMR due to insufficient NOE data. In the present study, we have applied: (1) rigid and adiabatic mapping, (2) free energy simulations, and (3) molecular dynamic simulations to elucidate the structure of the m conformer and to provide a possible pathway of the Trp21 side chain between the two conformers. Adiabatic energy mapping of conformations of the Trp21 side chain obtained by energy minimization identified two energy minima: One corresponding to the conformation of Trp21 observed in the X-ray crystallographic structure and solution structure of α3(VI) (the M conformation) and the second corresponding to the m conformation predicted by NMR spectroscopy. A transition pathway between the M and m conformation is suggested. The free-energy difference between the two conformers obtained by the thermodynamic integration method is calculated to 1.77 ± 0.7 kcal/mol in favor of the M form, which is in good agreement with NMR results. Structural and dynamic properties of the major and minor conformers of the α3(VI) molecule were investigated by molecular dynamic. Essential dynamics analysis of the two resulting 800 ps trajectories reveals that when going from the M to the m conformation only small, localized changes in the protein structure are induced. However, notable differences are observed in the mobility of the binding loop (residues Thr13–Ile18), which is more flexible in the m conformation than in the M conformation. This suggests that the reorientation of Trp21 might influence the inhibitory activity against trypsin, despite the relative large distance between the binding loop and Trp21.

The refolding of four disulfide lysozyme (at pH 5.2, 20 °C) involves parallel pathways, which have been proposed to merge at a near-native state. This species contains stable structure in the α- and β-domains but lacks a functional active site. Although previous experiments have demonstrated that the near-native state is populated on the fast refolding pathway, its relevance to slow refolding molecules could not be directly determined from previous experiments. In this paper, we describe experiments that investigate the effect of added salts on the refolding pathway of lysozyme at pH 5.2, 20 °C. We show, using stopped flow tryptophan fluorescence, inhibitor binding, and circular dichroism (CD), that the rate of formation of native lysozyme on the slow refolding track is significantly reduced in solutions of high ionic strength in a manner dependent on the position of the anion in the Hofmeister series. By contrast, the rate of evolution of hydrogen exchange (HX) protection monitored by electrospray ionization mass spectrometry (ESI MS) is unchanged under the refolding conditions studied. The data show, therefore, that at high ionic strengths β-domain stabilization and native state formation on the slow refolding pathway become kinetically decoupled such that the near-native state becomes significantly populated. Thus, by changing the energy landscape with the addition of salts new insights into the relevance of intermediate states in lysozyme refolding are revealed.

Folding of apomyoglobin is characterized by formation of a compact intermediate that contains substantial helicity. To determine whether this intermediate is obligatory or whether the protein can fold directly into the native state via an alternate parallel pathway, we have combined quench-flow hydrogen-exchange pulse labeling techniques with electrospray ionization mass spectrometry. The mass spectra of apomyoglobin obtained at various refolding times suggest that apomyoglobin indeed folds through a single pathway containing an obligatory intermediate with a significant hydrogen-bonded secondary structure content.

Binding of the protein Raf to the active form of Ras promotes activation of the MAP kinase signaling pathway, triggering cell growth and differentiation. Raf/Arg89 in the center of the binding interface plays an important role determining Ras–Raf binding affinity. We have investigated experimentally and computationally the Raf-R89K mutation, which abolishes signaling in vivo. The binding to [γ-35S]GTP-Ras of a fusion protein between the Raf-binding domain (RBD) of Raf and GST was reduced at least 175-fold by the mutation, corresponding to a standard binding free energy decrease of at least 3.0 kcal/mol. To compute this free energy and obtain insights into the microscopic interactions favoring binding, we performed alchemical simulations of the RBD, both complexed to Ras and free in solution, in which residue 89 is gradually mutated from Arg into Lys. The simulations give a standard binding free energy decrease of 2.9 ± 1.9 kcal/mol, in agreement with experiment. The use of numerous runs with three different force fields allows insights into the sources of uncertainty in the free energy and its components. The binding decreases partly because of a 7 kcal/mol higher cost to desolvate Lys upon binding, compared to Arg, due to better solvent interactions with the more concentrated Lys charge in the unbound state. This effect is expected to be general, contributing to the lower propensity of Lys to participate in protein–protein interfaces. Large contributions to the free energy change also arise from electrostatic interactions with groups up to 8 Å away, namely residues 37–41 in the conserved effector domain of Ras (including 4 kcal/mol from Ser39 which loses a bifurcated hydrogen bond to Arg89), the conserved Lys84 and Lys87 of Raf, and 2–3 specific water molecules. This analysis will provide insights into the large experimental database of Ras–Raf mutations.

The thioredoxin action upon the 2-oxoacid dehydrogenase complexes is investigated by using different thioredoxins, both wild-type and mutated. The attacking cysteine residue of thioredoxin is established to be essential for the thioredoxin-dependent activation of the complexes. Mutation of the buried cysteine residue to serine is not crucial for the activation, but prevents inhibition of the complexes, exhibited by the Clamydomonas reinhardtii thioredoxin m disulfide. Site-directed mutagenesis of D26, W31, F/W12, and Y/A70 (the Escherichia coli thioredoxin numbering is employed for all the thioredoxins studied) indicates that both the active site and remote residues of thioredoxin are involved in its interplay with the 2-oxoacid dehydrogenase complexes. Sequences of 11 thioredoxin species tested biochemically are aligned. The thioredoxin residues at the contact between the α3/310 and α1 helices, the length of the α1 helix and the charges in the α2-β3 and β4-β5 linkers are found to correlate with the protein influence on the 2-oxoacid dehydrogenase complexes (the secondary structural elements of thioredoxin are defined according to Eklund H et al., 1991, Proteins 11:13–28). The distribution of the charges on the surface of the thioredoxin molecules is analyzed. The analysis reveals the species specific polarization of the thioredoxin active site surroundings, which corresponds to the efficiency of the thioredoxin interplay with the 2-oxoacid dehydrogenase systems. The most effective mitochondrial thioredoxin is characterized by the strongest polarization of this area and the highest value of the electrostatic dipole vector of the molecule. Not only the magnitude, but also the orientation of the dipole vector show correlation with the thioredoxin action. The dipole direction is found to be significantly influenced by the charges of the residues 13/14, 51, and 83/85, which distinguish the activating and inhibiting thioredoxin disulfides.

The structure of the trigonal crystal form of bovine β-lactoglobulin variant B at pH 7.1 has been determined by X-ray diffraction methods at a resolution of 2.22 Å and refined to values for R and Rfree of 0.239 and 0.286, respectively. By comparison with the structure of the trigonal crystal form of bovine β-lactoglobulin variant A at pH 7.1, which was determined previously [Qin BY et al., 1998, Biochemistry 37:14014–14023], the structural consequences of the sequence differences D64G and V118A of variants A and B, respectively, have been investigated. Only minor differences in the core calyx structure occur. In the vicinity of the mutation site D64G on loop CD (residues 61–67), there are small changes in main-chain conformation, whereas the substitution V118A on β-strand H is unaccompanied by changes in the surrounding structure, thereby creating a void volume and weakened hydrophobic interactions with a consequent loss of thermal stability relative to variant A. A conformational difference is found for the loop EF, implicated in the pH-dependent conformational change known as the Tanford transition, but it is not clear whether this reflects differences intrinsic to the variants in solution or differences in crystallization.

Electrostatic interactions are often critical for determining the specificity of protein-protein complexes. To study the role of electrostatic interactions for assembly of helical bundles, we previously designed a thermostable, heterotrimeric coiled coil, ABC, in which charged residues were employed to drive preferential association of three distinct, 34-residue helices. To investigate the basis for heterotrimer specificity, we have used multiwavelength anomalous diffraction (MAD) analysis to determine the 1.8 Å resolution crystal structure of ABC. The structure shows that ABC forms a heterotrimeric coiled coil with the intended arrangement of parallel chains. Over half of the ion pairs engineered to restrict helix associations were apparent in the experimental electron density map. As seen in other trimeric coiled coils, ABC displays acute knobs-into-holes packing and a buried anion coordinated by core polar amino acids. These interactions validate the design strategy and illustrate how packing and polar contacts determine structural uniqueness.

The X-ray crystal structure of the human α-thrombin-hirunorm IV complex has been determined at 2.5 Å resolution, and refined to an R-factor of 0.173. The structure reveals an inhibitor binding mode distinctive of a true hirudin mimetic, which justifies the high inhibitory potency and the selectivity of hirunorm IV. This novel inhibitor, composed of 26 amino acids, interacts through the N-terminal end with the α-thrombin active site in a nonsubstrate mode, and binds specifically to the fibrinogen recognition exosite through the C-terminal end. The backbone of the N-terminal tripeptide Chg1"-Arg2"-2Nal3" (Chg, cyclohexyl-glycine; 2Nal, β-(2-naphthyl)-alanine) forms a parallel β-strand to the thrombin main-chain segment Ser214–Gly216. The Chg1" side chain occupies the S2 site, Arg2" penetrates into the S1 specificity site, while the 2Nal3" side chain occupies the aryl binding site. The Arg2" side chain enters the S1 specificity pocket from a position quite apart from the canonical P1 site. This notwithstanding, the Arg2" side chain establishes the typical ion pair with the carboxylate group of Asp189.

In addition to the Cys-Xaa-Xaa-Cys motif at position 30–33, DsbA, the essential catalyst for disulfide bond formation in the bacterial periplasm shares with other oxidoreductases of the thioredoxin family a cis-proline in proximity of the active site residues. In the variant DsbAP151A, this residue has been changed to an alanine, an almost isosteric residue which is not disposed to adopt the cis conformation. The substitution strongly destabilized the structure of DsbA, as determined by the decrease in the free energy of folding. The pKa of the thiol of Cys30 was only marginally decreased. Although in vivo the variant appeared to be correctly oxidized, it exhibited an activity less than half that of the wild-type enzyme with respect to the folding of alkaline phosphatase, used as a reporter of the disulfide bond formation in the periplasm. DsbAP151A crystallized in a different crystal form from the wild-type protein, in space group P21 with six molecules in the asymmetric unit. Its X-ray structure was determined to 2.8 Å resolution. The most significant conformational changes occurred at the active site. The loop 149–152 adopted a new backbone conformation with Ala151 in a trans conformation. This rearrangement resulted in the loss of van der Waals interactions between this loop and the disulfide bond. His32 from the Cys-Xaa-Xaa-Cys sequence presented in four out of six molecules in the asymmetric unit a gauche− conformation not observed in the wild-type protein. The X-ray structure and folding studies on DsbAP151A were consistent with the cis-proline playing a major role in the stabilization of the protein. A role for the positioning of the substrate is discussed. These important properties for the enzyme function might explain the conservation of this residue in DsbA and related proteins possessing the thioredoxin fold.

The mechanism of the disulfide-bond forming enzyme DsbA depends on the very low pKa of a cysteine residue in its active-site and on the relative instability of the oxidized enzyme compared to the reduced one. A thermodynamic cycle has been used to correlate its redox properties to the difference in the free energies of folding (ΔΔGred/ox) of the oxidized and reduced forms. However, the relation was proved unsatisfied for a number of DsbA variants. In this study, we investigate the thermodynamic and redox properties of a highly destabilized variant DsbAP151A (substitution of cis-Pro151 by an alanine) by the means of intrinsic tryptophan fluorescence and by high-sensitivity differential scanning calorimetry (HS-DSC). When the value of ΔΔGred/ox obtained fluorimetrically for DsbAP151A does not correlate with the value expected from its redox potential, the value of ΔΔGred/ox provided by HS-DSC are in perfect agreement with the predicted thermodynamic cycle for both wild-type and variant. HS-DSC data indicate that oxidized wild-type enzyme and the reduced forms of both wild-type and variant unfold according to a two-state mechanism. Oxidized DsbAP151A shows a deviation from two-state behavior that implies the loss of interdomain cooperativity in DsbA caused by Pro151 substitution. The presence of chaotrope in fluorimetric measurements could facilitate domain uncoupling so that the fluorescence probe (Trp76) does not reflect the whole unfolding process of DsbAP151A anymore. Thus, theoretical thermodynamic cycle is respected when an appropriate method is applied to DsbA unfolding under conditions in which protein domains still conserve their cooperativity.

Calmodulin (CaM) is a 148-residue regulatory calcium-binding protein that activates a wide range of target proteins and enzymes. Calcium-saturated CaM has a bilobal structure, and each domain has an exposed hydrophobic surface region where target proteins are bound. These two “active sites” of calmodulin are remarkably rich in Met residues. Here we have biosynthetically substituted (up to 90% incorporation) the unnatural amino acids ethionine (Eth) and norleucine (Nle) for the nine Met residues of CaM. The substituted proteins bind in a calcium-dependent manner to hydrophobic matrices and a synthetic peptide, encompassing the CaM-binding domain of myosin light-chain kinase (MLCK). Infrared and circular dichroism spectroscopy show that there are essentially no changes in the secondary structure of these proteins compared to wild-type CaM (WT-CaM). One- and two-dimensional NMR studies of the Eth-CaM and Nle-CaM proteins reveal that, while the core of the proteins is relatively unaffected by the substitutions, the two hydrophobic interaction surfaces adjust to accommodate the Eth and Nle residues. Enzyme activation studies with MLCK show that Eth-CaM and Nle-CaM activate the enzyme to 90% of its maximal activity, with little changes in dissociation constant. For calcineurin only 50% activation was obtained, and the KD for Nle-CaM also increased 3.5-fold compared with WT-CaM. These data show that the “active site” Met residues of CaM play a distinct role in the activation of different target enzymes, in agreement with site-directed mutagenesis studies of the Met residues of CaM.

β-Galactosidase (lacZ) from Escherichia coli is a 464 kDa homotetramer. Each subunit consists of five domains, the third being an α/β barrel that contains most of the active site residues. A comparison is made between each of the domains and a large set of proteins representative of all structures from the protein data bank. Many structures include an α/β barrel. Those that are most similar to the α/β barrel of E. coli β-galactosidase have similar catalytic residues and belong to the so-called “4/7 superfamily” of glycosyl hydrolases. The structure comparison suggests that β-amylase should also be included in this family. Of three structure comparison methods tested, the “ProSup” procedure of Zu-Kang and Sippl and the “Superimpose” procedure of Diederichs were slightly superior in discriminating the members of this superfamily, although all procedures were very powerful in identifying related protein structures. Domains 1, 2, and 4 of E. coli β-galactosidase have topologies related to “jelly-roll barrels” and “immunoglobulin constant” domains. This fold also occurs in the cellulose binding domains (CBDs) of a number of glycosyl hydrolases. The fold of domain 1 of E. coli β-galactosidase is closely related to some CBDs, and the domain contributes to substrate binding, but in a manner unrelated to cellulose binding by the CBDs. This is typical of domains 1, 2, 4, and 5, which appear to have been recruited to play roles in β-galactosidase that are unrelated to the functions that such domains provide in other contexts. It is proposed that β-galactosidase arose from a prototypical single domain α/β barrel with an extended active site cleft. The subsequent incorporation of elements from other domains could then have reduced the size of the active site from a cleft to a pocket to better hydrolyze the disaccharide lactose and, at the same time, to facilitate the production of inducer, allolactose.

One hundred-forty-five full-length aldehyde dehydrogenase-related sequences were aligned to determine relationships within the aldehyde dehydrogenase (ALDH) extended family. The alignment reveals only four invariant residues: two glycines, a phenylalanine involved in NAD binding, and a glutamic acid that coordinates the nicotinamide ribose in certain E-NAD binary complex crystal structures, but which may also serve as a general base for the catalytic reaction. The cysteine that provides the catalytic thiol and its closest neighbor in space, an asparagine residue, are conserved in all ALDHs with demonstrated dehydrogenase activity. Sixteen residues are conserved in at least 95% of the sequences; 12 of these cluster into seven sequence motifs conserved in almost all ALDHs. These motifs cluster around the active site of the enzyme. Phylogenetic analysis of these ALDHs indicates at least 13 ALDH families, most of which have previously been identified but not grouped separately by alignment. ALDHs cluster into two main trunks of the phylogenetic tree. The largest, the “Class 3” trunk, contains mostly substrate-specific ALDH families, as well as the class 3 ALDH family itself. The other trunk, the “Class 1/2” trunk, contains mostly variable substrate ALDH families, including the class 1 and 2 ALDH families. Divergence of the substrate-specific ALDHs occurred earlier than the division between ALDHs with broad substrate specificities. A site on the World Wide Web has also been devoted to this alignment project.

The structural organization of the B1 domain of streptococcal protein G (PGA) has been probed using molecular dynamics simulations, with a particular emphasis on the role of the solvent exposed Ile6 residue. In addition to the native protein (WT-PGA), three single-mutants (I6G-PGA, I6F-PGA, and I6T-PGA), one double-mutant (I6T,T53G-PGA), and three isolated peptide fragments (corresponding to the helix and the two β-hairpins) were studied in the presence of explicit water molecules. Comparative analysis of the various systems showed that the level of perturbation was directly related to the hydrophobicity and the size of the side chain of residue 6, the internal rigidity of the proteins decreasing in the order I6T-PGA > I6G-PGA > WT-PGA > I6F-PGA. The results emphasized the importance of residue 6 in controlling both the integrity of the sheet's surface and the orientation of the helix in relation to the sheet by modulation of surface/core interactions. The effects of mutations were delocalized across the structure, and glycine residues, in particular, absorbed most of the introduced strain. A qualitative structural decomposition of the native fold into elementary building-blocks was achieved using principal component analysis and mechanical response matrices. Within this framework, internal motions of the protein were described as coordinated articulations of these structural units, mutations affecting mostly the amplitude of the motions rather than the structure/location of the building-blocks. Analysis of the isolated peptidic fragments suggested that packing did not play a determinant role in defining the elementary building-blocks, but that chain topology was mostly responsible.

Staphylococcal nuclease forms three different partially-folded intermediates at low pH in the presence of low to moderate concentration of anions, differing in the amount of secondary structure, globularity, stability, and compactness. Although these intermediates are monomeric at low protein concentration (≤0.25 mg/mL), increasing concentrations of protein result in the formation of dimers and soluble oligomers, ultimately leading to larger insoluble aggregates. Unexpectedly, increasing protein concentration not only led to association, but also to increased structure of the intermediates. The secondary structure, stability, and globularity of the two less-ordered partially-folded intermediates (A1 and A2) were substantially increased upon association, suggesting that aggregation induces structure. An excellent correlation was found between degree of association and amount of structure measured by different techniques, including circular dichroism, fluorescence, Fourier transform infrared spectroscopy (FTIR), and small-angle X-ray scattering. The associated states were also substantially more stable toward urea denaturation than the monomeric forms. A mechanism is proposed, in which the observed association of monomeric intermediates involves intermolecular interactions which correspond to those found intramolecularly in normal folding to the native state.

Hsp16.3, the α-crystallin-related small heat shock protein of Mycobacterium tuberculosis that is maximally expressed during the stationary phase and is a major membrane protein, has been reported to form specific trimer-of-trimers structure and to act as an effective molecular chaperone (Chang Z et al., 1996, J Biol Chem 271:7218–7223). However, little is known about its action mechanism. In this study, Hsp16.3 conformational intermediates with dramatically increased chaperone activities were detected after treatment with very low concentrations of guanidine hydrochloride (0.05 M), urea (0.3 M), or mild heating (30°C). The intermediates showed a significant increase in their capacity to bind the hydrophobic probe 1-anilino-8-naphthalene sulfonate (ANS), indicating an increased exposure of hydrophobic surfaces. Interestingly, the greatest chaperone activities of Hsp16.3 were observed in the presence of 0.3 M guanidine HCl or when heated to 35°C. CD spectroscopy studies revealed no significant changes in protein secondary and tertiary structures at these mild treatments. Our in vitro studies also indicate that long-time-heated Hsp16.3, heated even to temperatures as high as 85°C, has almost the same, if not a slightly greater, chaperone activities as the native protein when cooled to room temperature and its secondary structures also almost recovered. Together, these results suggest that Hsp16.3 modulates its chaperone activity by exposing hydrophobic surfaces and that the protein structure is highly stable and flexible, thus highly adapted for its function.

The aspartyl dyad of free HIV-1 protease has apparent pKas of ∼3 and ∼6, but recent NMR studies indicate that the aspartyl dyad is fixed in the doubly protonated form over a wide pH range when cyclic urea inhibitors are bound, and in the monoprotonated form when the inhibitor KNI-272 is bound. We present computations and measurements related to these changes in protonation and to the thermodynamic linkage between protonation and inhibition. The Poisson–Boltzmann model of electrostatics is used to compute the apparent pKas of the aspartyl dyad in the free enzyme and in complexes with four different inhibitors. The calculations are done with two parameter sets. One assigns ε = 4 to the solute interior and uses a detailed model of ionization; the other uses ε = 20 for the solute interior and a simplified representation of ionization. For the free enzyme, both parameter sets agree well with previously measured apparent pKas of ∼3 and ∼6. However, the calculations with an internal dielectric constant of 4 reproduce the large pKa shifts upon binding of inhibitors, but the calculations with an internal dielectric constant of 20 do not. This observation has implications for the accurate calculation of pKas in complex protein environments.

Because binding of a cyclic urea inhibitor shifts the pKas of the aspartyl dyad, changing the pH is expected to change its apparent binding affinity. However, we find experimentally that the affinity is independent of pH from 5.5 to 7.0. Possible explanations for this discrepancy are discussed.