α2-Macroglobulin (α2M) is a major carrier of transforming growth factor-β (TGF-β) in vitro and in vivo. By screening glutathione S-transferase (GST) fusion proteins with overlapping sequences, we localized the TGF-β-binding site to aa 700–738 of the mature human α2M subunit. In separate experiments, we screened overlapping synthetic peptides corresponding to aa 696–777 of α2M and identified a single 16-mer (718–733) that binds TGF-β1. Platelet-derived growth factor-BB (PDGF-BB) bound to the same peptide, even though TGF-β and PDGF-BB share almost no sequence identity. The sequence of the growth factor-binding peptide, WDLVVVNSAGVAEVGV, included a high proportion of hydrophobic amino acids. The analogous peptide from murinoglobulin, a human α2M homologue that does not bind growth factors, contained only three nonconservative amino acid substitutions; however, the MUG peptide failed to bind TGF-β1 and PDGF-BB. These results demonstrate that a distinct and highly-restricted site in α2M, positioned near the C-terminal flank of the bait region, mediates growth factor binding. At least part of the growth factor-binding site is encoded by exon 18 of the α2M gene, which is notable for a 5′ splice site polymorphism that has been implicated in Alzheimer's Disease.

MAP30 is a 30 kDa single-stranded, type-I ribosome inactivating protein (RIP) possessing anti-tumor and anti-HIV activities. It binds both ribosomal RNA and the HIV-1 long-terminal repeat DNA. To understand the structural basis for MAP30 activities, we undertook the study of MAP30 by solution NMR spectroscopy. We report nearly complete 1H, 13C, and 15N chemical shift assignments of its 263 amino acids. Based upon an analysis of secondary 13C chemical shifts, 3JHNHA coupling constants, hydrogen exchange data, and nuclear Overhauser effect patterns, we find that the secondary structure and β-sheet topology of MAP30 are very similar to those of the ricin A chain, a subunit of the well-known type-II RIP, even though two proteins display distinct activities. We therefore suggest that MAP30 and ricin A chain share a similar three-dimensional fold, and that the reported functional differences between two proteins arise primarily from differences in local three-dimensional structure and other structural properties such as surface electrostatic potentials.

Sequence analysis of Candida rugosa lipase 1 (LIP1) predicts the presence of three N-linked glycosylation sites at asparagine 291, 314, 351. To investigate the relevance of sugar chains in the activation and stabilization of LIP1, we directed site mutagenesis to replace the above mentioned asparagine with glutamine residues. Comparison of the activity of mutants with that of the wild-type (wt) lipase indicates that both 314 and 351 Asn to Gln substitutions influence, although at a different extent, the enzyme activity both in hydrolysis and esterification reactions, but they do not alter the enzyme water activity profiles in organic solvents or temperature stability. Introduction of Gln to replace Asn351 is likely to disrupt a stabilizing interaction between the sugar chain and residues of the inner side of the lid in the enzyme active conformation. The effect of deglycosylation at position 314 is more difficult to explain and might suggest a more general role of the sugar moiety for the structural stability of lipase 1. Conversely, Asn291Gln substitution does not affect the lipolytic or the esterase activity of the mutant that behaves essentially as the wt enzyme. This observation supports the hypothesis that changes in activity of Asn314Gln and Asn351Gln mutants are specifically due to deglycosylation.

Bile-salt activated lipase (BAL) is a pancreatic enzyme that digests a variety of lipids in the small intestine. A distinct property of BAL is its dependency on bile salts in hydrolyzing substrates of long acyl chains or bulky alcoholic motifs. A crystal structure of the catalytic domain of human BAL (residues 1–538) with two surface mutations (N186D and A298D), which were introduced in attempting to facilitate crystallization, has been determined at 2.3 Å resolution. The crystal form belongs to space group P212121 with one monomer per asymmetric unit, and the protein shows an α/β hydrolase fold. In the absence of bound bile salt molecules, the protein possesses a preformed catalytic triad and a functional oxyanion hole. Several surface loops around the active site are mobile, including two loops potentially involved in substrate binding (residues 115–125 and 270–285).

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.

We have used X-ray fiber diffraction to probe the structure of fibers of tau and tau fragments. Fibers of fragments from the microtubule binding domain had a cross β-structure that closely resembles that reported both for neurofibrillary tangles found in Alzheimer's disease brain and for fibrous lesions from other protein folding diseases. In contrast, fibers of full-length tau had a different, more complex structure. Despite major differences at the molecular level, all fiber types exhibited very similar morphology by electron microscopy. These results have a number of implications for understanding the etiology of Alzheimer's and other tauopathic diseases. The morphology of the peptide fibers suggests that the region in tau corresponding to the peptides plays a critical role in the nucleation of fiber assembly. The dramatically different structure of the full length tau fibers suggests that some region in tau has enough inherent structure to interfere with the formation of cross β-fibers. Additionally, the similar appearance by electron microscopy of fibrils with varying molecular structure suggests that different molecular arrangements may exist in other samples of fibers formed from tau.

Several recent studies have shown that it is possible to increase protein stability by improving electrostatic interactions among charged groups on the surface of the folded protein. However, the stability increases are considerably smaller than predicted by a simple Coulomb's law calculation, and in some cases, a charge reversal on the surface leads to a decrease in stability when an increase was predicted. These results suggest that favorable charge–charge interactions are important in determining the denatured state ensemble, and that the free energy of the denatured state may be decreased more than that of the native state by reversing the charge of a side chain. We suggest that when the hydrophobic and hydrogen bonding interactions that stabilize the folded state are disrupted, the unfolded polypeptide chain rearranges to compact conformations with favorable long-range electrostatic interactions. These charge–charge interactions in the denatured state will reduce the net contribution of electrostatic interactions to protein stability and will help determine the denatured state ensemble. To support this idea, we show that the denatured state ensemble of ribonuclease Sa is considerably more compact at pH 7 where favorable charge–charge interactions are possible than at pH 3, where unfavorable electrostatic repulsion among the positive charges causes an expansion of the denatured state ensemble. Further support is provided by studies of the ionic strength dependence of the stability of charge–reversal mutants of ribonuclease Sa. These results may have important implications for the mechanism of protein folding.

Presentation of peptides derived from endogenous proteins by class I major histocompatibility complex molecules is essential both for immunological self-tolerance and induction of cytotoxic T-cell responses against intracellular parasites. Despite frequent and diverse post-translational modification of eukaryotic cell proteins, very few class I-bound peptides with post-translationally modified residues are known. Here we describe a natural dodecamer ligand of HLA-B39 (B*3910) derived from an RNA-binding nucleoprotein that carried NG,NG-dimethyl-Arg. Although common among RNA-binding proteins, this modification was not previously known among natural class I ligands. The sequence of this peptide was determined by Edman degradation and electrospray ion trap mass spectrometry. The fragmentation pattern of the dimethyl-Arg side chain observed with this latter technique allowed us to unambiguously assign the isomeric form of the modified residue. The post-translationally modified ligand was a prominent component (1–2%) of the B*3910-bound peptide repertoire. The dimethyl-Arg residue was located in a central position of the peptide, amenable to interacting with T-cell receptors, and most other residues in the middle region of the peptide were Gly. These structural features strongly suggest that the post-translationally modified residue may have a major influence on the antigenic properties of this natural ligand.

We propose a self-consistent approach to analyze knowledge-based atom–atom potentials used to calculate protein–ligand binding energies. Ligands complexed to actual protein structures were first built using the SMoG growth procedure (DeWitte & Shakhnovich, 1996) with a chosen input potential. These model protein–ligand complexes were used to construct databases from which knowledge-based protein–ligand potentials were derived. We then tested several different modifications to such potentials and evaluated their performance on their ability to reconstruct the input potential using the statistical information available from a database composed of model complexes. Our data indicate that the most significant improvement resulted from properly accounting for the following key issues when estimating the reference state: (1) the presence of significant nonenergetic effects that influence the contact frequencies and (2) the presence of correlations in contact patterns due to chemical structure. The most successful procedure was applied to derive an atom–atom potential for real protein–ligand complexes. Despite the simplicity of the model (pairwise contact potential with a single interaction distance), the derived binding free energies showed a statistically significant correlation (∼0.65) with experimental binding scores for a diverse set of complexes.

The role of the switch II region in the conformational transition of activation of Ha-ras-p21 has been investigated by mutating residues predicted to act as hinges for the conformational transition of this loop (Ala59, Gly60, and Gly75) (Diaz JF, Wroblowski B, Schlitter J, Engelborghs Y, 1997, Proteins 28:434–451), as well as mutating the catalytic residue Gln61. The proposed mutations of the hinge residues decrease the rate of the conformational transition of activation as measured by the binding of BeF3− to the GDP-p21 complex. Also, the thermodynamic parameters of the binding reaction are altered by a factor between three and five, depending on the temperature. (Due to changes in activation and reaction enthalpies, partially compensated by entropy changes.) The control mutation Q61H in which only the catalytic residue is changed has only a limited effect on the kinetic rate constants of the conformational transition and on the thermodynamic parameters of the reaction.

The fact that mutations of the hinge residues of the switch II region affect both the binding of the phosphate analog and the conformational transition of activation indicates that the switch II is implicated both in the early and the late states of the transition.

The backbone dynamics and overall tumbling of protein G have been investigated using 15N relaxation. Comparison of measured R2/R1 relaxation rate ratios with known three-dimensional coordinates of the protein show that the rotational diffusion tensor is significantly asymmetric, exhibiting a prolate axial symmetry. Extensive Monte Carlo simulations have been used to estimate the uncertainty due to experimental error in the relaxation rates to be D∥/D⊥ = 1.68 ± 0.08, while the dispersion in the NMR ensemble leads to a variation of D∥/D⊥ = 1.65 ± 0.03. Incorporation of this tensorial description into a Lipari–Szabo type analysis of internal motion has allowed us to accurately describe the local dynamics of the molecule. This analysis differs from an earlier study where the overall rotational diffusion was described by a spherical top. In this previous analysis, exchange parameters were fitted to many of the residues in the alpha helix. This was interpreted as reflecting a small motion of the alpha helix with respect to the beta sheet. We propose that the differential relaxation properties of this helix compared to the beta sheet are due to the near-orthogonality of the NH vectors in the two structural motifs with respect to the unique axis of the diffusion tensor. Our analysis shows that when anisotropic rotational diffusion is taken into account NH vectors in these structural motifs appear to be equally rigid. This study underlines the importance of a correct description of the rotational diffusion tensor if internal motion is to be accurately investigated.

Two new sets of scoring matrices are introduced: H2 for the protein sequence comparison and T2 for the protein sequence–structure correlation. Each element of H2 or T2 measures the frequency with which a pair of amino acid types in one protein, k-residues apart in the sequence, is aligned with another pair of residues, of given amino acid types (for H2) or in given structural states (for T2), in other structurally homologous proteins. There are four types, corresponding to the k-values of 1 to 4, for both H2 and T2. These matrices were set up using a large number of structurally homologous protein pairs, with little sequence homology between the pair, that were recently generated using the structure comparison program SHEBA.

The two scoring matrices were incorporated into the main body of the sequence alignment program SSEARCH in the FASTA package and tested in a fold recognition setting in which a set of 107 test sequences were aligned to each of a panel of 3,539 domains that represent all known protein structures. Six procedures were tested; the straight Smith-Waterman (SW) and FASTA procedures, which used the Blosum62 single residue type substitution matrix; BLAST and PSI-BLAST procedures, which also used the Blosum62 matrix; PASH, which used Blosum62 and H2 matrices; and PASSC, which used Blosum62, H2, and T2 matrices. All procedures gave similar results when the probe and target sequences had greater than 30% sequence identity. However, when the sequence identity was below 30%, a similar structure could be found for more sequences using PASSC than using any other procedure. PASH and PSI-BLAST gave the next best results.

Elongation factor Tu (EF-Tu) promotes the binding of aminoacyl-tRNA (aa-tRNA) to the acceptor site of the ribosome. During the elongation cycle, EF-Tu interacts with guanine nucleotides, aa-tRNA and its nucleotide exchange factor (EF-Ts). Quantitative determination of the equilibrium dissociation constants that govern the interactions of mammalian mitochondrial EF-Tu (EF-Tumt) with guanine nucleotides was the focus of the work reported here. Equilibrium dialysis with [3H]GDP was used to measure the equilibrium dissociation constant of the EF-Tumt·GDP complex (KGDP = 1.0 ± 0.1 μM). Competition of GTP with a fluorescent derivative of GDP (mantGDP) for binding to EF-Tumt was used to measure the dissociation constant of the EF-Tumt·GTP complex (KGTP = 18 ± 9 μM). The analysis of these data required information on the dissociation constant of the EF-Tumt·mantGDP complex (KmGDP = 2.0 ± 0.5 μM), which was measured by equilibrium dialysis. Both KGDP and KGTP for EF-Tumt are quite different (about two orders of magnitude higher) than the dissociation constants of the corresponding complexes formed by Escherichia coli EF-Tu. The forward and reverse rate constants for the association and dissociation of the EF-Tumt·GDP complex were determined using the change in the fluorescence of mantGDP upon interaction with EF-Tumt. These values are in agreement with a simple equilibrium binding interaction between EF-Tumt and GDP. The results obtained are discussed in terms of the recently described crystal structure of the EF-Tumt·GDP complex.

A simple method is presented for the determination of individual rate constants for substrate hydrolysis by serine proteases and other enzymes with similar catalytic mechanism. The method does not require solvent perturbation like viscosity changes, or solvent isotope effects, that often compromise nonspecifically the activity of substrate and enzyme. The rates of substrate diffusion into the active site (k1), substrate dissociation (k−1), acylation (k2), and deacylation (k3) in the accepted mechanism of substrate hydrolysis by serine proteases are derived from the temperature dependence of the Michaelis–Menten parameters kcat/Km and kcat. The method also yields the activation energies for these molecular events. Application to wild-type and mutant thrombins reveals how the various steps of the catalytic mechanism are affected by Na+ binding and site-directed mutations of the important residues Y225 in the Na+ binding environment and L99 in the S2 specificity site. Extension of this method to other proteases should enable the derivation of detailed information on the kinetic and energetic determinants of protease function.

The development of an energy or scoring function for protein structure prediction is greatly enhanced by testing the function on a set of computer-generated conformations (decoys) to determine whether it can readily distinguish native-like conformations from nonnative ones. We have created “Decoys ‘R’ Us,” a database containing many such sets of conformations, to provide a resource that allows scoring functions to be improved.

Escherichia coli CspA is a member of the cold shock protein family. All cold shock proteins studied to date fold rapidly by an apparent two-state mechanism. CspA contains an unusual cluster of aromatic amino acids on its surface that is necessary for nucleic acid binding and also provides stability to CspA (Hillier et al., 1998). To elucidate the role this aromatic cluster plays in the determining the folding rate and pathway of CspA, we have studied the folding kinetics of mutants containing either leucine or serine substituted for Phe18, Phe20, and/or Phe31. The leucine substitutions are found to accelerate folding and the serine substitutions to decelerate folding. Because these residues exert effects on the free energy of the folding transition state, they may be necessary for nucleating folding. They are not responsible, however, for the very compact, native-like transition state ensemble seen in the cold shock proteins, as the refolding rates of the mutants all show a similar, weak dependence of unfolding rate on denaturant concentration. Using mutant cycle analysis, we show that there is energetic coupling among the three residues between the unfolded and transition states, suggesting that the cooperative nature of these interactions helps to determine the unfolding rate. Overall, our results suggest that separate evolutionary pressures can act simultaneously on the same group of residues to maintain function, stability, and folding rate.

Guanine-containing mono- and dinucleotides bind to the active site of ribonuclease A in a nonproductive mode (retro-binding) (Aguilar CF, Thomas PJ, Mills A, Moss DS, Palmer RA, 1992, J Mol Biol 224:265–267). Guanine binds to the highly specific pyrimidine site by forming hydrogen bonds with Thr45 and with the sulfate anion located in the P1 site. To investigate the influence of the anion present in the P1 site on retro-binding, we determined the structure of two new complexes of RNase A with uridylyl(2′,5′)guanosine obtained by soaking two different forms of pre-grown RNase A crystals. In one case, RNase A was crystallized without removing the sulfate anion strongly bound to the active site; in the other, the protein was first equilibrated with a basic solution to displace the anion from the P1 site. The X-ray structures of the complexes with and without sulfate in P1 were refined using diffraction data up to 1.8 Å (R-factor 0.192) and 2.0 Å (R-factor 0.178), respectively. The binding mode of the substrate analogue to the protein differs markedly in the two complexes. When the sulfate is located in P1, we observe retro-binding; whereas when the anion is removed from the active site, the uridine is productively bound at the B1 site. In the productive complex, the electron density is very well defined for the uridine moiety, whereas the downstream guanine is disordered. This finding indicates that the interactions of guanine in the B2 site are rather weak and that this site is essentially adenine preferring. In this crystal form, there are two molecules per asymmetric unit, and due to crystal packing, only the active site of one molecule is accessible to the ligand. Thus, in the same crystal we have a ligand-bound and a ligand-free RNase A molecule. The comparison of these two structures furnishes a detailed and reliable picture of the structural alterations induced by the binding of the substrate. These results provide structural information to support the hypotheses on the role of RNase A active site residues that have recently emerged from site-directed mutagenesis studies.

Mutant human lysozymes (Ile56Thr & Asp67His) have been reported to form amyloid deposits in the viscera. From the standpoint of understanding the mechanism of amyloid formation, we searched for conditions of amyloid formation in vitro using hen egg lysozyme, which has been extensively studied from a physicochemical standpoint. It was found that the circular dichroism spectra in the far-ultraviolet region of the hen egg lysozyme changed to those characteristic of a β-structure from the native α-helix rich spectrum in 90% ethanol solution. When the concentration of protein was increased to 10 mg/mL, the protein solution formed a gel in the presence of 90% ethanol, and precipitated on further addition of 10 mM NaCl. The precipitates were examined by electron microscopy, their ability to bind Congo red, and X-ray diffraction to determine whether amyloid fibrils were formed in the precipitates. Electron micrographs displayed unbranched protofilament with a diameter of ∼70 Å. The peak point of the difference spectrum for the Congo red binding assay was 541 nm, which is characteristic of amyloid fibrils. The X-ray diffraction pattern showed a sharp and intense diffraction ring at 4.7 Å, a reflection that arises from the interstrand spacing in β-sheets. These results indicate that the precipitates of hen egg lysozyme are amyloid protofilament, and that the amyloid protofilament formation of hen egg lysozyme closely follows upon the destruction of the helical and tertiary structures.

The denaturant-induced (un)folding of apoflavodoxin from Azotobacter vinelandii has been followed at the residue level by NMR spectroscopy. NH groups of 21 residues of the protein could be followed in a series of 1H–15N heteronuclear single-quantum coherence spectra recorded at increasing concentrations of guanidinium hydrochloride despite the formation of protein aggregate. These NH groups are distributed throughout the whole apoflavodoxin structure. The midpoints of unfolding determined by NMR coincide with the one obtained by fluorescence emission spectroscopy. Both techniques give rise to unfolding curves with transition zones at significantly lower denaturant concentrations than the one obtained by circular dichroism spectroscopy. The NMR (un)folding data support a mechanism for apoflavodoxin folding in which a relatively stable intermediate is involved. Native apoflavodoxin is shown to cooperatively unfold to a molten globule-like state with extremely broadened NMR resonances. This initial unfolding step is slow on the NMR chemical shift timescale. The subsequent unfolding of the molten globule is faster on the NMR chemical shift timescale and the limited appearance of 1H–15N HSQC cross peaks of unfolded apoflavodoxin in the denaturant range studied indicates that it is noncooperative.

The nucleotide sequence of the gene ( pepA) of a zymogen of an aspartic proteinase from Penicillium janthinellum with a 71% identity in the deduced amino acid sequence to penicillopepsin (which we propose to call penicillopepsin-JT1) has been determined. The gene consists of 60 codons for a putative leader sequence of 20 amino acid residues, a sequence of about 150 nucleotides that probably codes for an activation peptide and a sequence with two introns that codes for the active aspartic proteinase. This gene, inserted into the expression vector pGPT-pyrG1, was expressed in an aspartic proteinase-free strain of Aspergillus niger var. awamori in high yield as a glycosylated form of the active enzyme that we call penicillopepsin-JT2. After removal of the carbohydrate component with endoglycosidase H, its relative molecular mass is between 33,700 and 34,000. Its kinetic properties, especially the rate-enhancing effects of the presence of alanine residues in positions P3 and P′2 of substrates, are similar to those of penicillopepsin-JT1, endothiapepsin, rhizopuspepsin, and pig pepsin. Earlier findings suggested that this rate-enhancing effect was due to a hydrogen bond between the -NH- of P3 and the hydrogen bond accepting oxygen of the side chain of the fourth amino acid residue C-terminal to Asp215. Thr219 of penicillopepsin-JT2 was mutated to Ser, Val, Gly, and Ala. Thr219Ser showed an increase in kcat when a P3 residue was present in the substrate, which was similar to that of the wild-type, whereas the mutants Thr219Val, Thr219Gly, and Thr219Ala showed no significant increase when a P3 residue was added. The results show that the putative hydrogen bond alone is responsible for the increase. We propose that by locking the -NH- of P3 to the enzyme, the scissile peptide bond between P1 and P′1 becomes distorted toward a tetrahedral conformation and becomes more susceptible to nucleophilic attack by the catalytic apparatus without the need of a conformational change in the enzyme.