Well-resolved ion signals of intact large protein assemblies, with molecular masses extending above one million Dalton, have been detected and mass analyzed using electrospray ionization mass spectrometry, with an uncertainty in mass of <0.2%. The mass spectral data seem to reflect known solution-phase behavior of the studied protein assembly and have therefore been directly used to probe the protein assembly topology and stability as a function of ionic strength and pH.

The Ure2 protein from Saccharomyces cerevisiae has been proposed to undergo a prion-like autocatalytic conformational change, which leads to inactivation of the protein, thereby generating the [URE3] phenotype. The first 65 amino acids, which are dispensable for the cellular function of Ure2p in nitrogen metabolism, are necessary and sufficient for [URE3] (Masison & Wickner, 1995), leading to designation of this domain as the Ure2 prion domain (UPD). We expressed both UPD and Ure2 as glutathione-S-transferase (GST) fusion proteins in Escherichia coli and observed both to be initially soluble. Upon cleavage of GST-UPD by thrombin, the released UPD formed ordered fibrils that displayed amyloid-like characteristics, such as Congo red dye binding and green-gold birefringence. The fibrils exhibited high β-sheet content by Fourier transform infrared spectroscopy. Fiber formation proceeded in an autocatalytic manner. In contrast, the released, full-length Ure2p formed mostly amorphous aggregates; a small amount polymerized into fibrils of uniform size and morphology. Aggregation of Ure2p could be seeded by UPD fibrils. Our results provide biochemical support for the proposal that the [URE3] state is caused by a self-propagating inactive form of Ure2p. We also found that the uncleaved GST-UPD fusion protein could polymerize into amyloid fibrils by a strictly autocatalytic mechanism, forcing the GST moiety of the protein to adopt a new, β-sheet-rich conformation. The findings on the GST-UPD fusion protein indicate that the ability of the prion domain to mediate a prion-like conversion process is not specific for or limited to the Ure2p.

This paper presents an analytically tractable model that captures the most elementary aspect of the protein folding problem, namely that both the energy and the entropy decrease as a protein folds. In this model, the system diffuses within a sphere in the presence of an attractive spherically symmetric potential. The native state is represented by a small sphere in the center, and the remaining space is identified with unfolded states. The folding temperature, the time-dependence of the populations, and the relaxation rate are calculated, and the folding dynamics is analyzed for both golf-course and funnel-like energy landscapes. This simple model allows us to illustrate a surprising number of concepts including entropic barriers, transition states, funnels, and the origin of single exponential relaxation kinetics.

In a general approach to the understanding of protein adaptation to high temperature, molecular models of the closely related mesophilic Streptomyces sp. S38 Xyl1 and thermophilic Thermomonospora fusca TfxA family 11 xylanases were built and compared with the three-dimensional (3D) structures of homologous enzymes. Some of the structural features identified as potential contributors to the higher thermostability of TfxA were introduced in Xyl1 by site-directed mutagenesis in an attempt to improve its thermostability and thermophilicity. A new Y11–Y16 aromatic interaction, similar to that present in TfxA and created in Xyl1 by the T11Y mutation, improved both the thermophilicity and thermostability. Indeed, the optimum activity temperature (70 vs. 60 °C) and the apparent Tm were increased by about 9 °C, and the mutant was sixfold more stable at 57 °C. The combined mutations A82R/F168H/N169D/Δ170 potentially creating a R82–D169 salt bridge homologous to that present in TfxA improved the thermostability but not the thermophilicity. Mutations R82/D170 and S33P seemed to be slightly destabilizing and devoid of influence on the optimal activity temperature of Xyl1. Structural analysis revealed that residues Y11 and Y16 were located on β-strands B1 and B2, respectively. This interaction should increase the stability of the N-terminal part of Xyl1. Moreover, Y11 and Y16 seem to form an aromatic continuum with five other residues forming putative subsites involved in the binding of xylan (+3, +2, +1, −1, −2). Y11 and Y16 might represent two additional binding subsites (−3, −4) and the T11Y mutation could thus improve substrate binding to the enzyme at higher temperature and thus the thermophilicity of Xyl1.

The chaperonin HSP60 (GroEL) proteins are essential in eubacterial genomes and in eukaryotic organelles. Functional regions inferred from mutation studies and the Escherichia coli GroEL 3D crystal complexes are evaluated in a multiple alignment across 43 diverse HSP60 sequences, centering on ATP/ADP and Mg2+ binding sites, on residues interacting with substrate, on GroES contact positions, on interface regions between monomers and domains, and on residues important in allosteric conformational changes. The most evolutionary conserved residues relate to the ATP/ADP and Mg2+ binding sites. Hydrophobic residues that contribute in substrate binding are also significantly conserved. A large number of charged residues line the central cavity of the GroEL–GroES complex in the substrate-releasing conformation. These span statistically significant intra- and inter-monomer three-dimensional (3D) charge clusters that are highly conserved among sequences and presumably play an important role interacting with the substrate. Unaligned short segments between blocks of alignment are generally exposed at the outside wall of the Anfinsen cage complex. The multiple alignment reveals regions of divergence common to specific evolutionary groups. For example, rickettsial sequences diverge in the ATP/ADP binding domain and Gram-positive sequences diverge in the allosteric transition domain. The evolutionary information of the multiple alignment proffers attractive sites for mutational studies.

A combination of epitope excision, epitope extraction, and differential chemical modification followed by mass spectrometric peptide mapping was used for the characterization of a discontinuous epitope that is recognized by the mouse anti-HIV-p24 monoclonal antibody 5E2.A3. In epitope excision, the protein is first bound to an immobilized antibody and then digested with proteolytic enzymes. In epitope extraction, the protein is first digested and subsequently allowed to react with the antibody. After epitope excision of the p24-5E2.A3 complex with endoproteinase Lys-C, a large fragment remained affinity bound corresponding to amino acids 1–158 of HIV-p24 (fragment 1–158). Further digestion, however, resulted in loss of affinity. Moreover, no affinity-bound fragments were observed after an epitope extraction experiment. These data from the epitope excision and extraction experiments suggest that the epitope is discontinuous. For the further characterization of the epitope, amino groups in the epitope-containing fragment were acetylated in both the affinity bound and free states followed by mass spectrometric analysis. Two successive acetylation reactions were performed: (1) the first used a low molar excess of acetic anhydride, and (2) the second, after separation from the antibody, a high molar excess of its hexadeuteroderivative. This isotopic labeling procedure, in combination with high resolution mass spectrometry, allowed the precise determination of relative reactivities of amino groups. In this study, no differences were observed in the ranking of the relative reactivities of five lysine residues. However, the N-terminal amino group was found to be part of the discontinuous epitope.

A single water molecule (w135), buried within the structure of rat intestinal fatty acid binding protein (I-FABP), is investigated by NMR, molecular dynamics simulations, and analysis of known crystal structures. An ordered water molecule was found in structurally analogous position in 24 crystal structures of nine different members of the family of fatty acid binding proteins. There is a remarkable conservation of the local structure near the w135 binding site among different proteins from this family. NMR cross-relaxation measurements imply that w135 is present in the I-FABP:ANS (1-sulfonato-8-(1′)anilinonaphthalene) complex in solution with the residence time of >300 ps. Mean-square positional fluctuations of w135 oxygen observed in MD simulations (0.18 and 0.13 Å2) are comparable in magnitude to fluctuations exhibited by the backbone atoms and result from highly constrained binding pocket as revealed by Voronoi volumes (averages of 27.0 ± 1.8 Å3 and 24.7 ± 2.2 Å3 for the two simulations). Escape of w135 from its binding pocket was observed only in one MD simulation. The escape process was initiated by interactions with external water molecules and was accompanied by large deformations in β-strands D and E. Immediately before the release, w135 assumed three distinct states that differ in hydrogen bonding topology and persisted for about 15 ps each. Computer simulations suggest that escape of w135 from the I-FABP matrix is primarily determined by conformational fluctuations of the protein backbone and interactions with external water molecules.

Transmembrane helices are the most readily predictable secondary structure components of proteins. They can be predicted to a high degree of accuracy in a variety of ways. Many of these methods compare new sequence data with the sequence characteristics of known transmembrane domains. However, the known transmembrane sequences are not necessarily representative of a particular organism. We attempt to demonstrate that parameters optimized for the known transmembrane domains are far from optimal when predicting transmembrane regions in a given genome. In particular, we have tested the effect of nucleotide bias upon the composition and hence the prediction characteristics of transmembrane helices. Our analysis shows that nucleotide bias of a genome has a strong and predictable influence upon the occurrences of several of the most important hydrophobic amino acids found within transmembrane helices. Thus, we show that nucleotide bias should be taken into account when determining putative transmembrane domains from sequence data.

NMR spectroscopy was used to search for mechanistically significant differences in the local mobility of the main-chain amides of Bacillus circulans xylanase (BCX) in its native and catalytically competent covalent glycosyl-enzyme intermediate states. 15N T1, T2, and 15N{1H} NOE values were measured for ~120 out of 178 peptide groups in both the apo form of the protein and in BCX covalently modified at position Glu78 with a mechanism-based 2-deoxy-2-fluoro-β-xylobioside inactivator. Employing the model-free formalism of Lipari and Szabo, the measured relaxation parameters were used to calculate a global correlation time (τm) for the protein in each form (9.2 ± 0.2 ns for apo-BCX; 9.8 ± 0.3 ns for the modified protein), as well as individual order parameters for the main-chain NH bond vectors. Average values of the order parameters for the protein in the apo and complexed forms were S2 = 0.86 ± 0.04 and S2 = 0.91 ± 0.04, respectively. No correlation is observed between these order parameters and the secondary structure, solvent accessibility, or hydrogen bonding patterns of amides in either form of the protein. These results demonstrate that the backbone of BCX is well ordered in both states and that formation of the glycosyl-enzyme intermediate leads to little change, in any, in the dynamic properties of BCX on the time scales sampled by 15N-NMR relaxation measurements.

The oxidative refolding of ribonuclease A has been investigated in several experimental conditions using a variety of redox systems. All these studies agree that the formation of disulfide bonds during the process occurs through a nonrandom mechanism with a preferential coupling of certain cysteine residues. We have previously demonstrated that in the presence of glutathione the refolding process occurs through the reiteration of two sequential reactions: a mixed disulfide with glutathione is produced first which evolves to form an intramolecular S-S bond. In the same experimental conditions, protein disulfide isomerase (PDI) was shown to catalyze formation and reduction of mixed disulfides with glutathione as well as formation of intramolecular S-S bonds.

This paper reports the structural characterization of the one-disulfide intermediate population during the oxidative refolding of Ribonuclease A under the presence of PDI and glutathione with the aim of defining the role of the enzyme at the early stages of the reaction. The one-disulfide intermediate population occurring at the early stages of both the uncatalyzed and the PDI-catalyzed refolding was purified and structurally characterized by proteolytic digestion followed by MALDI-MS and LC/ESIMS analyses. In the uncatalyzed refolding, a total of 12 disulfide bonds out of the 28 theoretical possible cysteine couplings was observed, confirming a nonrandom distribution of native and nonnative disulfide bonds. Under the presence of PDI, only two additional nonnative disulfides were detected. Semiquantitative LC/ESIMS analysis of the distribution of the S-S bridged peptides showed that the most abundant species were equally populated in both the uncatalyzed and the catalyzed process.

This paper shows the first structural characterization of the one-disulfide intermediate population formed transiently during the refolding of ribonuclease A in quasi-physiological conditions that mimic those present in the ER lumen. At the early stages of the process, three of the four native disulfides are detected, whereas the Cys26–Cys84 pairing is absent. Most of the nonnative disulfide bonds identified are formed by nearest-neighboring cysteines. The presence of PDI does not significantly alter the distribution of S-S bonds, suggesting that the ensemble of single-disulfide species is formed under thermodynamic control.

A dominant feature of folding of cytochrome c is the presence of nonnative His-heme kinetic traps, which either pre-exist in the unfolded protein or are formed soon after initiation of folding. The kinetically trapped species can constitute the majority of folding species, and their breakdown limits the rate of folding to the native state. A temperature jump (T-jump) relaxation technique has been used to compare the unfolding/folding kinetics of yeast iso-2 cytochrome c and a genetically engineered double mutant that lacks His-heme kinetic traps, H33N,H39K iso-2. The results show that the thermodynamic properties of the transition states are very similar. A single relaxation time τobs is observed for both proteins by absorbance changes at 287 nm, a measure of solvent exclusion from aromatic residues. At temperatures near Tm, the midpoint of the thermal unfolding transitions, τobs is four to eight times faster for H33N,H39K iso-2 (τobs ∼ 4–10 ms) than for iso-2 (τobs ∼ 20–30 ms). T-jumps show that there are no kinetically unresolved (τ < 1–3 μs T-jump dead time) “burst” phases for either protein. Using a two-state model, the folding (kf) and unfolding (ku) rate constants and the thermodynamic activation parameters ΔGf[Dagger], ΔGu[Dagger], ΔHf[Dagger], ΔHu[Dagger], ΔSf[Dagger], ΔSu[Dagger] are evaluated by fitting the data to a function describing the temperature dependence of the apparent rate constant kobs (= τobs−1) = kf + ku. The results show that there is a small activation enthalpy for folding, suggesting that the barrier to folding is largely entropic. In the “new view,” a purely entropic kinetic barrier to folding is consistent with a smooth funnel folding landscape.

The X-ray structure of chitinase from the fungal pathogen Coccidioides immitis has been solved to 2.2 Å resolution. Like other members of the class 18 hydrolase family, this 427 residue protein is an eight-stranded β/α-barrel. Although lacking an N-terminal chitin anchoring domain, the enzyme closely resembles the chitinase from Serratia marcescens. Among the conserved features are three cis peptide bonds, all involving conserved active site residues. The active site is formed from conserved residues such as tryptophans 47, 131, 315, 378, tyrosines 239 and 293, and arginines 52 and 295. Glu171 is the catalytic acid in the hydrolytic mechanism; it was mutated to a Gln, and activity was abolished. Allosamidin is a substrate analog that strongly inhibits the class 18 enzymes. Its binding to the chitinase hevamine has been observed, and we used conserved structural features of the two enzymes to predict the inhibitors binding to the fungal enzyme.

The refolding kinetics of a single-chain Fv (scFv) fragment, derived from a stabilized mutant of the phosphorylcholine binding antibody McPC603, was investigated by H/D exchange and ESI-MS and compared with the folding kinetics of its constituting domains VH and VL. Both VH and VL adopt essentially native-like exchange protection within the dead time of the manual-mixing H/D exchange experiment (10 s) and in the case of VL, which contains two cis-prolines in the native conformation, this fast protection is independent of proline cis/trans isomerization. At the earliest time point resolvable by manual mixing, fewer deuterons are protected in the scFv fragment than in the two isolated domains together, despite the fact that the scFv fragment is significantly more stable than VL and VH. Full H/D exchange protection in the scFv fragment is gained on a time scale of minutes. This means that the domains in the scFv fragment do not refold independently. Rather, they associate prematurely and in nonnative form, a kinetic trap. Unproductive domain association is observed both after equilibrium- and short-term denaturation. For the equilibrium-denatured scFv fragment, whose native structure formation is dependent on a cis conformation of an interface proline in VL, this cis/trans isomerization reaction proceeds about one order in magnitude more slowly than the escape from the trap to a conformation where full H/D exchange protection is already achieved. We interpret these data in terms of a general kinetic scheme involving intermediates with and without domain association.

The effects of 1,1,1,3,3,3-hexafluoro-isopropanol (HFIP) on the conformation of cytochrome c (cyt c) at pH 1.9 were studied using a combination of spectroscopic and physical methods. Analysis varying the HFIP concentration showed that a compact denatured conformation (MHF) accumulates in a low concentration range of HFIP in the middle of structural transition from the highly unstructured acid-denatured state to the highly helical alcohol-denatured state of cyt c. This contrasts clearly with the effect of isopropanol (IP), in which no compact conformation accompanied with the transition. Analysis varying concentrations of HFIP and NaCl concurrently showed that the MHF state of cyt c is essentially identical to the salt-induced molten-globule (MG) state, and the MG state in the presence of salt was also stabilized by a low concentration of HFIP. Furthermore, 2,2,2-trifluoroethanol stabilized MHF similarly to HFIP, supporting the proposition that the specific effect observed for HFIP is caused by fluorination of alcohol. The mechanism stabilizing compact conformation by HFIP remains unclear, but is probably distinct from that of salts and polyols, which are also known to stabilize the MG-like state.

The issue of specificity in tyrosine kinase intracellular signaling mediated by src homology 2 (SH2) domains has great importance in the understanding how individual signals maintain their mutual exclusivity and affect downstream responses. Several proteins contain tandem SH2 domains that, on interacting with their ligand, provide a higher level of specificity than can be afforded by the interaction of a single SH2 domain. In this study, we focus on the comparison of two proteins ZAP70 and the p85 subunit of PI 3-kinase, which although distinctly different in function and general structure, possess tandem SH2 domains separated by a linker region and which bind to phosphorylated receptor molecules localized to the cell membrane. Binding studies using isothermal titration calorimetry show that these two proteins interact with peptides mimicking their physiological ligands in very different ways. In the case of the SH2 domains from ZAP70, they interact with a stoichiometry of unity, while p85 is able to make two distinct interactions, one with a stoichiometry of 1:1 and the other with two p85 molecules interacting with one receptor. The observation of two different modes of binding of p85 might be important in providing different cellular responses based on fluctuating intracellular concentration regimes of this protein. Thermodynamic data on both proteins suggest that a conformational change occurs on binding. On investigation of this structural change using a truncated form of p85 (including just the two SH2 domains and the inter-SH2 region), both NMR and circular dichroism spectroscopic studies failed to show significant changes in secondary structure. This suggests that any conformational change associated with binding is small and potentially limited to loop regions of the protein.

Search and study the general principles that govern kinetics and thermodynamics of protein folding generates new insight into the factors that control this process. Here, we demonstrate based on the known experimental data and using theoretical modeling of protein folding that side-chain entropy is one of the general determinants of protein folding. We show for proteins belonging to the same structural family that there exists an optimal relationship between the average side-chain entropy and the average number of contacts per residue for fast folding kinetics. Analysis of side-chain entropy for proteins that fold without additional agents demonstrates that there exists an optimal region of average side-chain entropy for fast folding. Deviation of the average side-chain entropy from the optimal region results in an anomalous protein folding process (prions, α-lytic protease, subtilisin, some DNA-binding proteins). Proteins with high or low side-chain entropy would have extended unfolded regions and would require some additional agents for complete folding. Such proteins are common in nature, and their structure properties have biological importance.

The genome of the unicellular cyanobacterium Synechocystis sp. PCC 6803 contains a gene (slr2097, glbN) encoding a 123 amino-acid product with sequence similarity to globins. Related proteins from cyanobacteria, ciliates, and green algae bind oxygen and have a pronounced tendency to coordinate the heme iron with two protein ligands. To study the structural and functional properties of Synechocystis sp. PCC 6803 hemoglobin, slr2097 was cloned and overexpressed in Escherichia coli. Purification of the hemoglobin was performed after addition of hemin to the clarified cell lysate. Recombinant, heme-reconstituted ferric Synechocystis sp. PCC 6803 hemoglobin was found to be a stable helical protein, soluble to concentrations higher than 500 μM. At neutral pH, it yielded an electronic absorption spectrum typical of a low-spin ferric species, with maxima at 410 and 546 nm. The proton NMR spectrum revealed sharp lines spread over a chemical shift window narrower than 40 ppm, in support of low-spin hexacoordination of the heme iron. Nuclear Overhauser effects demonstrated that the heme is inserted in the protein matrix to produce one major equilibrium form. Addition of dithionite resulted in an absorption spectrum with maxima at 426, 528, and 560 nm. This reduced form appeared capable of carbon monoxide binding. Optical data also suggested that cyanide ions could bind to the heme in the ferric state. The spectral properties of the putative Synechocystis sp. PCC 6803 hemoglobin confirmed that it can be used for further studies of an ancient hemoprotein structure.

Effects of guanidine hydrochloride (GdnHCl) on the structure and dynamics of wild-type Humicola lanuginosa lipase (HLL) and its two mutants were studied. The latter were S146A (with the active site Ser replaced by Ala) and the single Trp mutant W89m, with substitutions W117F, W221H, and W260H. Steady-state, stopped-flow, and time-resolved laser-induced fluorescence spectroscopy were carried out as a function of [GdnHCl]. The maximum emission wavelength and fluorescence lifetimes revealed the microenvironment of the tryptophan(s) in these lipases to become more polar upon increasing [GdnHCl]. However, significant extent of tertiary structure in GdnHCl is suggested by the observation that both wild-type HLL and W89m remain catalytically active at rather high GdnHCl concentrations of >6 and 4.0 M, respectively. Changes in steady-state emission anisotropy, as well as variation in rotational correlation times and residual anisotropy values, demonstrate that upon increasing [GdnHCl] the structure of the lipases became more loose, with increasing amplitude of structural fluctuations. Finally, intermediate states in the course of exposure of the proteins to GdnHCl were revealed by stopped-flow fluorescence measurements.

Staphylococcal aureus epidermolytic toxins (ETs) A and B are responsible for the induction of staphylococcal scalded skin syndrome, a disease of neonates and young children. The clinical features of this syndrome vary from localized blisters to severe exfoliation affecting most of the body surface. Comparison of the crystal structures of two subtypes of ETs-rETA (at 2.0 Å resolution), rETB (at 2.8 Å resolution), and an active site variant of rETA, Ser195Ala at 2.0 Å resolution has demonstrated that their overall topology resembles that of a “trypsin-like” serine protease, but with significant differences at the N- and C-termini and loop regions. The details of the catalytic site in both ET structures are very similar to those in glutamate-specific serine proteases, suggesting a common catalytic mechanism. However, the “oxyanion hole,” which is part of the catalytic sites of glutamate specific serine proteases, is in the closed or inactive conformation for rETA, yet in the open or active conformation for rETB. The ETs contain a unique amphipathic helix at the N-terminus, and it appears to be involved in optimizing the conformation of the catalytic site residues. Determination of the structure of the rETA catalytic site variant, Ser195Ala, showed no significant perturbation at the active site, establishing that the loss of biological and esterolytic activity can be attributed solely to disruption of the catalytic serine residue. Finally, the crystal structure of ETs, together with biochemical data and mutagenesis studies, strongly confirms the classification of these molecules as “serine proteases” rather than “superantigens.”

The crystal structure of the fibrinolytic enzyme tissue plasminogen activator (tPA) shows that the bulky side chain of Y99 hinders access to the active site by partially occluding the S2 site and may be responsible for the low catalytic activity of tPA toward plasminogen. We have tested the role of Y99 by replacing it with Leu, the residue found in more proficient proteases like trypsin and thrombin. The Y99L replacement results in an increase in the kcat/Km for chromogenic substrates due to enhanced diffusion into the active site. The increase is modest (threefold) for substrates specific for tPA that carry Pro or Gly at P2, but reaches 80-fold for less specific substrates carrying Arg at P2. On the other hand, the Y99L mutation has no effect on the activity of tPA toward the natural substrate plasminogen, that carries Gly at P2, and reduces more than 10-fold the inhibition of tPA by plasminogen activator inhibitor-1 (PAI-1), that carries Ala at P2. We conclude that the steric hindrance provided by Y99 in the crystal structure affects mostly nonphysiological substrates with bulky residues at P2. In addition, residue Y99 plays an active role in the recognition of PAI-1, but not plasminogen. Mutations of Y99 could therefore afford a resistance to inhibition by PAI-1 without compromising the fibrinolytic potency of tPA, a result of potential therapeutic relevance.