The stability of protein is defined not only by the hydrogen bonding, hydrophobic effect, van der Waals interactions, and salt bridges. Additional, much more subtle contributions to protein stability can arise from surface residues that change their properties upon unfolding. The recombinant major cold shock protein of Escherichia coli CspA an all-β protein unfolds reversible in a two-state manner, and behaves in all other respects as typical globular protein. However, the enthalpy of CspA unfolding strongly depends on the pH and buffer composition. Detailed analysis of the unfolding enthalpies as a function of pH and buffers with different heats of ionization shows that CspA unfolding in the pH range 5.5–9.0 is linked to protonation of an amino group. This amino group appears to be the N-terminal α-amino group of the CspA molecule. It undergoes a 1.6 U shift in pKa values between native and unfolded states. Although this shift in pKa is expected to contribute ∼5 kJ/mol to CspA stabilization energy the experimentally observed stabilization is only ∼1 kJ/mol. This discrepancy is related to a strong enthalpy–entropy compensation due, most likely, to the differences in hydration of the protonated and deprotonated forms of the α-amino group.

The Q7K/L33I/L63I HIV-1 protease mutant was expressed in Escherichia coli and the effect of binding a substrate-analog inhibitor, acetyl-pepstatin, was investigated by fluorescence spectroscopy and molecular dynamics. The dimeric enzyme has four intrinsic tryptophans, located at positions 6 and 42 in each monomer. Fluorescence spectra and acrylamide quenching experiments show two differently accessible Trp populations in the apoenzyme with kq1 = 6.85 × 109 M−1 s−1 and kq2 = 1.88 × 109 M−1 s−1, that merge into one in the complex with kq = 1.78 × 109 M−1 s−1.

500 ps trajectory analysis of Trp χ1/χ2 rotameric interconversions suggest a model to account for the observed Trp fluorescence. In the simulations, Trp6/Trp6B rotameric interconversions do not occur on this timescale for both HIV forms. In the apoenzyme simulations, however, both Trp42s and Trp42Bs are flipping between χ1/χ2 states; in the complexed form, no such interconverions occur. A detailed investigation of the local Trp environments sampled during the molecular dynamics simulation suggests that one of the apoenzyme Trp42B rotameric interconversions would allow indole-quencher contact, such as with nearby Tyr59. This could account for the short lifetime component. The model thus interprets the experimental data on the basis of the conformational fluctuations of Trp42s alone. It suggests that the rotameric interconversions of these Trps, located relatively far from the active site and at the very start of the flap region, becomes restrained when the apoenzyme binds the inhibitor. The model is thus consistent with associating components of the fluorescence decay in HIV-1 protease to ground state conformational heterogeneity.

The paucity of detailed X-ray crystallographic structures of integral membrane proteins arises from substantive technical obstacles in the overexpression of multimilligram quantities of protein, and in the crystallization of purified protein-detergent complexes (PDCs). With rare exception, crystal contacts within the lattice are mediated by protein–protein interaction, and the detergent surrounding the protein behaves as a disordered solvent. The addition and use of surfactants that display mesoscopic self-assembly behavior in membrane protein crystallization experiments presents a novel alternative strategy. Well-ordered crystals of the water channel human aquaporin-1 (hAQP1) that diffract to 4 Å resolution have been obtained with this approach.

Both the peptidase activity against small fluorescent peptide substrates and the ATPase activity of Lon (La) proteases are stimulated by unstructured proteins such as α-casein. This stimulation reveals the simultaneous interaction of Lon with two proteolytic substrates—α-casein and the peptide substrate. To understand the cellular function of this stimulation, it is important to determine the physical properties of Lon stimulators. The abilities of compositionally simple random copolymers of amino acids (rcAAs) to stimulate the peptidase and ATPase activities of the Lon protease from Mycobacterium smegmatis (Ms-Lon) and its N-terminal truncation mutant (N-E226) were determined. We report that cationic but not anionic rcAAs stimulated Ms-Lon's peptidase activity but were themselves poor substrates for the enzyme. Peptidase stimulation by rcAAs correlated approximately with the degree of hydrophobicity of these polypeptides and reached levels >10-fold higher than observed previously for Ms-Lon stimulators such as α-casein. In contrast to α-casein, which stimulates Ms-Lon's peptidase activity by 40% and ATPase activity by 150%, rcAAs stimulated peptidase activity without concomitant stimulation of ATPase activity. Active site labeling experiments suggested that both rcAAs and ATP increased peptidase activity by increasing accessibility to the peptidase active site. Peptidase activity assays in the presence of both α-casein and rcAAs revealed that interactions of rcAAs and α-casein with Ms-Lon are extremely complex and not mutually exclusive. Specifically, (1) additions of low concentrations of α-casein (<50 μg/mL) caused a further stimulation of Ms-Lon's rcAA-stimulated peptidase activity; (2) additions of higher concentrations of α-casein inhibited Ms-Lon's rcAA-stimulated peptidase activity; (3) additions of all concentrations of α-casein inhibited N-E226's rcAA-stimulated peptidase activity. We conclude the Ms-Lon can interact with an rcAA, α-casein, and a substrate peptide simultaneously, and that formation of this quaternary complex requires the N-terminal domain of Ms-Lon. These data support models of Ms-Lon that include two allosteric polypeptide binding sites distinct from the catalytic peptidase site.

Viscosities of aqueous solutions of guanidine hydrochloride (GuHCl) were measured in the presence of varying amounts of glucose. At high concentrations of glucose or GuHCl, the measured viscosities showed significant deviation from the values computed using a method proposed by Tanford (1966, J Biol Chem 241:3228–3232). This method was originally derived to allow the calculation of the effects of buffer or low concentrations of salts and other additives on the density and viscosity of aqueous solutions of urea or GuHCl. Recently it has been used to estimate the viscosity of denaturant solutions that contain high concentrations of viscogens. Our results show that the extrapolation of this approach to solutions of highly concentrated viscous co-solutes leads to significant errors. The implications for experimental studies of the viscosity dependence of conformational transitions in proteins is discussed.

We have characterized the thermodynamic stability of the SH3 domain from the Saccharomyces cerevisiae Abp1p protein and found it to be relatively low compared to most other SH3 domains, with a Tm of 60 °C and a ΔGu of 3.08 kcal/mol. Analysis of a large alignment of SH3 domains led to the identification of atypical residues at eight positions in the wild-type Abp1p SH3 domain sequence that were subsequently replaced by the residue seen most frequently at that position in the alignment. Three of the eight mutants constructed in this way displayed increases in Tm ranging from 8 to 15 °C with concomitant increases in ΔGu of up to 1.4 kcal/mol. The effects of these substitutions on folding thermodynamics and kinetics were entirely additive, and a mutant containing all three was dramatically stabilized with a Tm greater than 90 °C and a ΔGu more than double that of the wild-type domain. The folding rate of this hyperstable mutant was 10-fold faster than wild-type, while its unfolding rate was fivefold slower. All of the stabilized mutants were still able to bind a target peptide with wild-type affinity. We have analyzed the stabilizing amino acid substitutions isolated in this study and several other similar sequence alignment based studies. In approximately 25% of cases, increased stability can be explained by enhanced propensity of the substituted residue for the local backbone conformation at the mutagenized site.

High mobility group (HMG) proteins 14 and 17 are nonhistone nuclear proteins that have been implicated in control of transcription and chromatin structure. To examine the posttranslational modifications of HMG-14 and -17 in vivo, HMG proteins were prepared from nuclear vs. cytosolic fractions of human K562 cells treated with 12-O-tetradecanoylphorbol 13-acetate (TPA) or okadaic acid (OA) and examined by electrospray mass spectrometry. Analysis of full-length masses demonstrated mono-, di-, and triphosphorylation of HMG-14 and mono- and diphosphorylation of HMG-17 from OA treated cells, whereas HMG-14 and -17 from TPA treated cells were monophosphorylated. Peptide mass and sequence analysis showed major and minor phosphorylation sites, respectively, at Ser24 and Ser28 in HMG-17, and Ser20 and Ser24 in HMG-14. These sites were found in the consensus sequence RRSARLSAK, within the nucleosomal binding domain of each protein. A third phosphorylation site in HMG-14 was located at either Ser6 or Ser7. Interestingly, the proportion of HMG-14 and -17 found in cytosolic pools increased significantly after 1 h of treatment compared to control cells and showed preferential phosphorylation compared with proteins from nuclear fractions. These results suggest that phosphorylation of HMG-14 and -7 interferes with nuclear localization mechanisms in a manner favoring release from nuclei.

The use of Escherichia coli asparaginase II as a drug for the treatment of acute lymphoblastic leukemia is complicated by the significant glutaminase side activity of the enzyme. To develop enzyme forms with reduced glutaminase activity, a number of variants with amino acid replacements in the vicinity of the substrate binding site were constructed and assayed for their kinetic and stability properties. We found that replacements of Asp248 affected glutamine turnover much more strongly than asparagine hydrolysis. In the wild-type enzyme, N248 modulates substrate binding to a neighboring subunit by hydrogen bonding to side chains that directly interact with the substrate. In variant N248A, the loss of transition state stabilization caused by the mutation was 15 kJ mol−1 for l-glutamine compared to 4 kJ mol−1 for l-aspartic β-hydroxamate and 7 kJ mol−1 for l-asparagine. Smaller differences were seen with other N248 variants. Modeling studies suggested that the selective reduction of glutaminase activity is the result of small conformational changes that affect active-site residues and catalytically relevant water molecules.

Electrospray ionization (ESI) mass spectra of both well-characterized and novel metallothioneins (MTs) from various species were recorded to explore their metal-ion-binding modes and stoichiometries. The ESI mass spectra of the zinc- and cadmium-binding MTs showed a single main peak corresponding to metal-to-protein ratios of 4, 6, or 7. These findings combined with data obtained by other methods suggest that these MTs bind zinc or cadmium in a single predominant form and are consistent with the presence of three- and four-metal clusters. An unstable copper-specific MT isoform from Roman snails (Helix pomatia) could be isolated intact and was shown to preferentially bind 12 copper ions. To obtain additional information on the formation and relative stability of metal-thiolate clusters in MTs, a mass spectrometric titration study was conducted. One to seven molar equivalents of zinc or of cadmium were added to metal-free human MT-2 at neutral pH, and the resulting complexes were measured by ESI mass spectrometry. These experiments revealed that the formation of the four-metal cluster and of the thermodynamically less stable three-metal cluster is sequential and largely cooperative for both zinc and cadmium. Minor intermediate forms between metal-free MT, Me4MT, and fully reconstituted Me7MT were also observed. The addition of increasing amounts of cadmium to metal-free blue crab MT-I resulted in prominent peaks whose masses were consistent with apoMT, Cd3MT, and Cd6MT, reflecting the known structure of this MT with two Me3Cys9 centers. In a similar reconstitution experiment performed with Caenorhabditits elegans MT-II, a series of signals corresponding to apoMT and Cd3MT to Cd6MT species were observed.

It is now possible to identify over 30 functional subfamilies among the WD-repeat-containing proteins found in the completed genomes. The majority of these subfamilies have at least one member for which experimental data allow assignment to a cellular pathway or process. Half of the 63 WD-repeat-containing proteins in Saccharomyces cerevisiae, half of the 70 in Caenorhabditis elegans, and a third of the 100 plus predicted in Drosophila can be assigned to 23 of these functional subfamilies. Perhaps indicative of the future, 33 WD-repeat-containing proteins from the partial genome of Arabidopsis thaliana can now be assigned to 18 of these subfamilies. These assignments have been made possible by combining traditional sequence similarity with an implied common beta propeller structural context to obtain measures of protein–protein surface similarity. The beta propeller structural context is represented in the form of a Hidden Markov Model. The procedure is completely automated.

Comparisons of protein sequence via cyclic training of Hidden Markov Models (HMMs) in conjunction with alignments of three-dimensional structure, using the Combinatorial Extension (CE) algorithm, reveal two putative EF-hand metal binding domains in acetylcholinesterase. Based on sequence similarity, putative EF-hands are also predicted for the neuroligin family of cell surface proteins. These predictions are supported by experimental evidence. In the acetylcholinesterase crystal structure from Torpedo californica, the first putative EF-hand region binds the Zn2+ found in the heavy metal replacement structure. Further, the interaction of neuroligin 1 with its cognate receptor neurexin depends on Ca2+. Thus, members of the α,β hydrolase fold family of proteins contain potential Ca2+ binding sites, which in some family members may be critical for heterologous cell associations.

The monomer–dimer equilibrium of the glycophorin A (GpA) transmembrane (TM) fragment has been used as a model system to investigate the amino acid sequence requirements that permit an appropriate helix–helix packing in a membrane-mimetic environment. In particular, we have focused on a region of the helix where no crucial residues for packing have been yet reported. Various deletion and replacement mutants in the C-terminal region of the TM fragment showed that the distance between the dimerization motif and the flanking charged residues from the cytoplasmic side of the protein is important for helix packing. Furthermore, selected GpA mutants have been used to illustrate the rearrangement of TM fragments that takes place when leucine repeats are introduced in such protein segments. We also show that secondary structure of GpA derivatives was independent from dimerization, in agreement with the two-stage model for membrane protein folding and oligomerization.

Nitrile hydratase from Rhodococcus sp. N-771 is an αβ heterodimer with a nonheme ferric iron in the catalytic center. In the catalytic center, αCys112 and αCys114 are modified to a cysteine sulfinic acid (Cys-SO2H) and a cysteine sulfenic acid (Cys-SOH), respectively. To understand the function and the biogenic mechanism of these modified residues, we reconstituted the nitrile hydratase from recombinant unmodified subunits. The αβ complex reconstituted under argon exhibited no activity. However, it gradually gained the enzymatic activity through aerobic incubation. ESI-LC/MS analysis showed that the anaerobically reconstituted αβ complex did not have the modification of αCys112-SO2H and aerobic incubation induced the modification. The activity of the reconstituted αβ complex correlated with the amount of αCys112-SO2H. Furthermore, ESI-LC/MS analyses of the tryptic digest of the reconstituted complex, removed of ferric iron at low pH and carboxamidomethylated without reduction, suggested that αCys114 is modified to Cys-SOH together with the sulfinic acid modification of αCys112. These results suggest that αCys112 and αCys114 are spontaneously oxidized to Cys-SO2H and Cys-SOH, respectively, and αCys112-SO2H is responsible for the catalytic activity solely or in combination with αCys114-SOH.

The binding of the methylanthraniloyl derivatives of ATP (mant-ATP), ADP (mant-ADP), 2′deoxyATP (mant-2′deoxyATP), and 3′deoxyATP (mant-3′deoxyATP) to the catalytic subunit of protein kinase A was studied to gain insights into the mechanism of nucleotide binding. The binding of the mant nucleotides leads to a large increase in fluorescence energy transfer at 440 nm, allowing direct measurements of nucleotide affinity. The dissociation constant of mant-ADP is identical to that for ADP, while that for mant-ATP is approximately threefold higher than that for ATP. The dissociation constant for mant-3′deoxyATP is approximately fivefold higher than that for 3′deoxyATP while derivatization of 2′deoxyATP does not affect affinity. The time-dependent binding of mant-ATP, mant-2′deoxyATP, and mant-ADP, measured using stopped-flow fluorescence spectroscopy, is best fit to three exponentials. The fast phase is ligand dependent, while the two slower phases are ligand independent. The slower phases are similar but not identical in rate, and have opposite fluorescence amplitudes. Both isomers of mant-ATP are equivalent substrates, as judged by reversed-phase chromatography, although the rate of phosphorylation is approximately 20-fold lower than the natural nucleotide. The kinetic data are consistent with a three-step binding mechanism in which initial association of the nucleotide derivatives produces a highly fluorescent complex. Either one or two conformational changes can occur after the formation of this binary species, but one of the isomerized forms must have low fluorescence compared to the initial binary complex. These data soundly attest to the structural plasticity within the kinase core that may be essential for catalysis. Overall, the mant nucleotides present a useful reporter system for gauging these conformational changes in light of the prevailing three-dimensional models for the enzyme.

Vascular endothelial growth factor (VEGF)-A interacts with the receptor tyrosine kinases VEGF-R1 and R2, and the importance of this interaction in endothelial cell (EC) function and blood vessel development has been well documented. Other ligands that interact differentially with these receptors and that are structurally related to VEGF-A include VEGF-B, VEGF-C, VEGF-D, and placenta growth factor (PLGF). Compared with VEGF-A, relatively little is known about the biological role of the VEGF-R1 specific ligand, VEGF-B. Two splice variant isoforms that differ at the COOH-terminus and which retain unique solubility characteristics are widely expressed throughout embryonic and postnatal development. Recent analysis of mice with a targeted deletion of the VEGF-B gene has revealed a defect in heart development and function consistent with an important role in vascularization of the myocardium (Bellomo D et al., 2000, Circ Res 86:E29–E35). To facilitate further characterization of VEGF-B, we have developed a protocol for expression and purification of refolded recombinant protein from Escherichia coli inclusion bodies (IBs). The approach developed resolves a number of significant issues associated with VEGF-B, including the ability to heterodimerize with endogenous VEGF-A when co-expressed in mammalian cells, a complex secondary structure incorporating inter- and intrachain disulfide bonds and hydrophobic characteristics that preclude the use of standard chromatographic resins. The resulting purified disulfide-linked homodimer was demonstrated to bind to VEGF-R1 and to compete with VEGF-A for binding to this receptor.

The Figure 2 caption that was published in the above referenced article (Volume 9, Number 2, 2000, pp. 344–352) is incorrect. We are reproducing Figure 2 below with its correct caption. As a consequence of this error, we are reprinting the following section from Results and discussion.

Coiled coils consist of two or more amphipathic α-helices wrapped around each other to form a superhelical structure stabilized at the interhelical interface by hydrophobic residues spaced in a repeating 3–4 sequence pattern. Dimeric coiled coils have been shown to often form in a single step reaction in which association and folding of peptide chains are tightly coupled. Here, we ask whether such a simple folding mechanism may also apply to the formation of a three-stranded coiled coil. The designed 29-residue peptide LZ16A was shown previously to be in a concentration-dependent equilibrium between unfolded monomer (M), folded dimer (D), and folded trimer (T ). We show by time-resolved fluorescence change experiments that folding of LZ16A to D and T can be described by [equation] and [equation]. The following rate constants were determined (25 °C, pH 7): k1 = 7.8 × 104 M−1 s−1, k−1 = 0.015 s−1, k2 = 6.5 × 105 M−1 s−1, and k−2 = 1.1 s−1. In a separate experiment, equilibrium binding constants were determined from the change with concentration of the far-ultraviolet circular dichroism spectrum of LZ16A and were in good agreement with the kinetic rate constants according to KD = k1/2k−1 and KT = k2/k−2. Furthermore, pulsed hydrogen-exchange experiments indicated that only unfolded M and folded D and T were significantly populated during folding. The results are compatible with a two-step reaction in which a subpopulation of association competent (e.g., partly helical) monomers associate to dimeric and trimeric coiled coils.

We present an analysis of the assumptions behind some of the most commonly used methods for evaluating the goodness of the fit between a sequence and a structure. Our studies on a lattice model show that methods based on statistical considerations are easy to use and can capture some of the features of protein-like sequences and their corresponding native states, but unfortunately are incapable of recognizing, with certainty, the native-like conformation of a sequence among a set of decoys. Meanwhile, an optimization method, entailing the determination of the parameters of an effective free energy of interaction, is much more reliable in recognizing the native state of a sequence. However, the statistical method is shown to perform quite well in tests of protein design.

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.

Conformational changes are essential for the activity of many proteins. If, or how fast, internal fluctuations are related to slow conformational changes that mediate protein function is not understood. In this study, we measure internal fluctuations of the transport protein lactose permease in the presence and absence of substrate by tryptophan fluorescence spectroscopy. We demonstrate that nanosecond fluctuations of α-helices are enhanced when the enzyme transports substrate. This correlates with previously published kinetic data from transport measurements showing that millisecond conformational transitions of the substrate-loaded carrier are faster than those in the absence of substrate. These findings corroborate the hypothesis of the hierarchical model of protein dynamics that predicts that slow conformational transitions are based on fast, thermally activated internal motions.