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.

The WW domain is an approximately 38 residue peptide-binding motif that binds a variety of sequences, including the consensus sequence xPPxY. We have displayed hYAP65 WW on the surface of M13 phage and randomized one-third of its three-stranded antiparallel β-sheet. Improved binding to the hydrophobic peptide, GTPPPPYTVG (WW1), was selected in the presence of three different concentrations of proteinase K to simultaneously drive selection for improved stability as well as high-affinity binding. While some of the selected binders show cooperative unfolding transitions, others show noncooperative thermal unfolding curves. Two novel WW consensus sequences have been identified, which bind to the xPPxY motif with higher affinity than the wild-type hYAP65 WW domain. These WW domain sequences are not precedented in any natural WW domain sequence. Thus, there appear to be a large number of motifs capable of recognizing the target peptide sequence, only a subset of which appear to be used in natural proteins.

Human ribonuclease inhibitor (hRI) is a cytosolic protein that protects cells from the adventitious invasion of pancreatic-type ribonucleases. hRI has 32 cysteine residues. The oxidation of these cysteine residues to form disulfide bonds is a rapid, cooperative process that inactivates hRI. The most proximal cysteine residues in native hRI are two pairs that are adjacent in sequence: Cys94 and Cys95, and Cys328 and Cys329. A cystine formed from such adjacent cysteine residues would likely contain a perturbing cis peptide bond within its eight-membered ring, which would disrupt the structure of hRI and could facilitate further oxidation. We find that replacing Cys328 and Cys329 with alanine residues has little effect on the affinity of hRI for bovine pancreatic ribonuclease A (RNase A), but increases its resistance to oxidation by 10- to 15-fold. Similar effects are observed for the single variants, C328A hRI and C329A hRI, suggesting that oxidation resistance arises from the inability to form a Cys328–Cys329 disulfide bond. Replacing Cys94 and Cys95 with alanine residues increases oxidation resistance to a lesser extent, and decreases the affinity of hRI for RNase A. The C328A, C329A, and C328A/C329A variants are likely to be more useful than wild-type hRI for inhibiting pancreatic-type ribonucleases in vitro and in vivo. We conclude that replacing adjacent cysteine residues can confer oxidation resistance in a protein.

The substrate specificity of porcine pepsin has been altered by site-directed mutagenesis in an attempt to selectively cleave bovine hide collagen at only a few sites, similar to cathepsin D, for the production of high quality gelatin. Kinetic parameters were determined using chromogenic peptide substrates based on the sequence Lys-Pro-Xaa-Yaa-Phe*Nph-Arg-Leu (where Xaa is Ile or Pro, Yaa is Glu, Leu, Gln or Lys, Nph is p-nitrophenylalanine, and * is the site of cleavage). Substitution of Thr222 and Glu287 within the S2 subsite of pepsin by Val and Met, respectively, produced a double mutant with a two- to fourfold higher kcat/Km, compared with wild-type pepsin, for the chromogenic peptides with residues Leu, Gln, and Glu at position P2 (Yaa). The results suggest that the functional group of the P2 side chain may be exposed to solvent, while the aliphatic portion interacts with hydrophobic residues comprising S2. Wild-type pepsin cleaved a peptide corresponding to the carboxy-terminal telopeptide region of bovine type I collagen α1 chain, SGGYDLSFLPQPPQE, predominantly at three sites (Asp-Leu, Leu-Ser, and Phe-Leu) and at a significantly lower rate at Ser-Phe. However, Thr222Val/Glu287Met cleaved site Ser-Phe at a rate 20-fold higher than the wild-type. Significantly, enzymes containing the double substitution Phe111Thr/Leu112Phe cleaved this peptide predominantly at one site Leu-Ser (similar to cathepsin D) and at a rate 23-fold higher than the wild-type. These mutants can potentially enhance the rate of solubilization of bovine hide collagen under conditions mild enough to maintain the triple helix structure and hence minimize the rate of subsequent denaturation and proteolytic cleavage.

We report the solution structure of the chemotactic cytokine (chemokine) vMIP-II. This protein has unique biological activities in that it blocks infection by several different human immunodeficiency virus type 1 (HIV-1) strains. This occurs because vMIP-II binds to a wide range of chemokine receptors, some of which are used by HIV to gain cell entry. vMIP-II is a monomeric protein, unlike most members of the chemokine family, and its structure consists of a disordered N-terminus, followed by a helical turn (Gln25–Leu27), which leads into the first strand of a three-stranded antiparallel β-sheet (Ser29–Thr34; Gly42–Thr47; Gln52–Asp56). Following the sheet is a C-terminal α-helix, which extends from residue Asp60 until Gln68. The final five residues beyond the C-terminal helix (Pro70–Arg74) are in an extended conformation, but several of these C-terminal residues contact the first β-strand. The structure of vMIP-II is compared to other chemokines that also block infection by HIV-1, and the structural basis of its lack of ability to form a dimer is discussed.

Residue Asn47 at position L1 of a type II′ β-turn of the α-spectrin SH3 domain is located in a disallowed region of the Ramachandran plot (φ = 56 ± 12, ψ = −118 ± 17). Therefore, it is expected that replacement of Asn47 by Gly should result in a considerable stabilization of the protein. Thermodynamic analysis of the N47G and N47A mutants shows that the change in free energy is small (∼0.7 kcal/mol; ∼3 kJ/mol) and comparable to that found when mutating a Gly to Ala in a α-helix or β-sheet. X-ray structural analysis of these mutants shows that the conformation of the β-turn does not change upon mutation and, therefore, that there is no relaxation of the structure, nor is there any gain or loss of interactions that could explain the small energy change. Our results indicate that the energetic definition of II′ region of the Ramachandran plot (φ = 60 ± 30, ψ = −115 ± 15) should be revised for at least Ala and Asn in structure validation and protein design.

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.

The X-ray crystallographic structures of two mutants (K206Q and H207E) of the N-lobe of human transferrin (hTF/2N) have been determined to high resolution (1.8 and 2.0 Å, respectively). Both mutant proteins bind iron with greater affinity than native hTF/2N. The structures of the K206Q and H207E mutants show interactions (both H-bonding and electrostatic) that stabilize the interaction of Lys296 in the closed conformation, thereby stabilizing the iron bound forms.

The small heat shock proteins (sHsps), which counteract heat and oxidative stress in an unknown way, belong to a protein family of sHsps and α-crystallins whose members form large oligomeric complexes. The chloroplast-localized sHsp, Hsp21, contains a conserved methionine-rich sequence, predicted to form an amphipatic helix with the methionines situated along one of its sides. Here, we report how methionine sulfoxidation was detected by mass spectrometry in proteolytically cleaved peptides that were produced from recombinant Arabidopsis thaliana Hsp21, which had been treated with varying concentrations of hydrogen peroxide. Sulfoxidation of the methionine residues in the conserved amphipatic helix coincided with a significant conformational change in the Hsp21 protein oligomer.

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.

To investigate the role of some tertiary interactions, the disulfide bonds, in the early stages of refolding of hen lysozyme, we report the kinetics of reoxidation of denatured and reduced lysozyme under the same refolding conditions as those previously used to investigate the kinetics of regain of its circular dichroism (CD), fluorescence, and activity. At different stages of the refolding, the oxidation of the protein was blocked by alkylation of the free cysteines with iodoacetamide and the various oxidation states present in the samples were identified by electrospray-mass spectrometry. Thus, it was possible to monitor the appearance and/or disappearance of the species with 0 to 4 disulfide bonds. Using a simulation program, these kinetics were compared with those of regain of far-UV CD, fluorescence, and enzymatic activity and were discussed in terms of a refined model for the refolding of reduced hen egg white lysozyme.

The binding free energies of four inhibitors to bovine β-trypsin are calculated. The inhibitors use either ornithine, lysine, or arginine to bind to the S1 specificity site. The electrostatic contribution to binding free energy is calculated by solving the finite difference Poisson–Boltzmann equation, the contribution of nonpolar interactions is calculated using a free energy-surface area relationship and the loss of conformational entropy is estimated both for trypsin and ligand side chains. Binding free energy values are of a reasonable magnitude and the relative affinity of the four inhibitors for trypsin is correctly predicted. Electrostatic interactions are found to oppose binding in all cases. However, in the case of ornithine- and lysine-based inhibitors, the salt bridge formed between their charged group and the partially buried carboxylate of Asp189 is found to stabilize the complex. Our analysis reveals how the molecular architecture of the trypsin binding site results in highly specific recognition of substrates and inhibitors. Specifically, partially burying Asp189 in the inhibitor-free enzyme decreases the penalty for desolvation of this group upon complexation. Water molecules trapped in the binding interface further stabilize the buried ion pair, resulting in a favorable electrostatic contribution of the ion pair formed with ornithine and lysine side chains. Moreover, all side chains that form the trypsin specificity site are partially buried, and hence, relatively immobile in the inhibitor-free state, thus reducing the entropic cost of complexation. The implications of the results for the general problem of recognition and binding are considered. A novel finding in this regard is that like charged molecules can have electrostatic contributions to binding that are more favorable than oppositely charged molecules due to enhanced interactions with the solvent in the highly charged complex that is formed.

The disulfide bond pattern of catrocollastatin-C was determined by N-terminal sequencing and mass spectrometry. The N-terminal disintegrin-like domain is a compact structure including eight disulfide bonds, seven of them in the same pattern as the disintegrin bitistatin. The protein has two extra cysteine residues (XIII and XVI) that form an additional disulfide bond that is characteristically found in the disintegrin-like domains of cellular metalloproteinases (ADAMs) and PIII snake venom Zn-metalloproteinases (SVMPs). The C-terminal cysteine-rich domain of catrocollastatin-C contains five disulfide bonds between nearest-neighbor cysteines and a long range disulfide bridge between CysV and CysX. These results provide structural evidence for a redefinition of the disintegrin-like and cysteine-rich domain boundaries. An evolutionary pathway for ADAMs, PIII, and PII SVMPs based on disulfide bond engineering is also proposed.

A 12-residue peptide designed to form an α-helix and self-associate into an antiparallel 4-α-helical bundle yields a 0.9 Å crystal structure revealing unanticipated features. The structure was determined by direct phasing with the “Shake-and-Bake” program, and contains four crystallographically distinct 12-mer peptide molecules plus solvent for a total of 479 atoms. The crystal is formed from nearly ideal α-helices hydrogen bonded head-to-tail into columns, which in turn pack side-by-side into sheets spanning the width of the crystal. Within each sheet, the α-helices run antiparallel and are closely spaced (9–10 Å center-to-center). The sheets are more loosely packed against each other (13–14 Å between helix centers). Each sheet is amphiphilic: apolar leucine side chains project from one face, charged lysine and glutamate side chains from the other face. The sheets are stacked with two polar faces opposing and two apolar faces opposing. The result is a periodic biomaterial composed of packed protein bilayers, with alternating polar and apolar interfaces. All of the 30 water molecules in the unit cell lie in the polar interface or between the stacked termini of helices. A section through the sheet reveals that the helices packed at the apolar interface resemble the four-α-helical bundle of the design, but the helices overhang parts of the adjacent bundles, and the helix crossing angles are less steep than intended (7–11° rather than 18°).

We describe the results of a procedure for maximizing the number of sequences that can be reliably linked to a protein of known three-dimensional structure. Unlike other methods, which try to increase sensitivity through the use of fold recognition software, we only use conventional sequence alignment tools, but apply them in a manner that significantly increases the number of relationships detected. We analyzed 11 genomes and found that, depending on the genome, between 23 and 32% of the ORFs had significant matches to proteins of known structure. In all cases, the aligned region consisted of either >100 residues or >50% of the smaller sequence. Slightly higher percentages could be attained if smaller motifs were also included. This is significantly higher than most previously reported methods, even those that have a fold-recognition component. We survey the biochemical and structural characteristics of the most frequently occurring proteins, and discuss the extent to which alignment methods can realistically assign function to gene products.

The lipases from Rhizopus and Rhizomucor are members of the family of Mucorales lipases. Although they display high sequence homology, their stereoselectivity toward triradylglycerols (sn-2 substituted triacylglycerols) varies. Four different triradylglycerols were investigated, which were classified into two groups: flexible substrates with rotatable O′-C1′ ether or ester bonds adjacent to C2 of glycerol and rigid substrates with a rigid N′-C1′ amide bond or a phenyl ring in sn-2. Although Rhizopus lipase shows opposite stereopreference for flexible and rigid substrates (hydrolysis in sn-1 and sn-3, respectively), Rhizomucor lipase hydrolyzes both groups of triradylglycerols preferably in sn-1. To explain these experimental observations, computer-aided molecular modeling was applied to study the molecular basis of stereoselectivity. A generalized model for both lipases of the Mucorales family highlights the residues mediating stereoselectivity: (1) L258, the C-terminal neighbor of the catalytic histidine, and (2) G266, which is located in a loop contacting the glycerol backbone of a bound substrate. Interactions with triradylglycerol substrates are dominated by van der Waals contacts. Stereoselectivity can be predicted by analyzing the value of a single substrate torsion angle that discriminates between sn-1 and sn-3 stereopreference for all substrates and lipases investigated here. This simple model can be easily applied in enzyme and substrate engineering to predict Mucorales lipase variants and synthetic substrates with desired stereoselectivity.

The reaction mechanism of the catalytic phosphoryl transfer of cAMP-dependent protein kinase (cAPK) was investigated by semi-empirical AM1 molecular orbital computations of an active site model system derived from the crystal structure of the catalytic subunit of the enzyme. The activation barrier is calculated as 20.7 kcal mol−1 and the reaction itself to be exothermic by 12.2 kcal mol−1. The active site residue Asp166, which was often proposed to act as a catalytic base, does not accept a proton in any of the reaction steps. Instead, the hydroxyl hydrogen of serine is shifted to the simultaneously transferred phosphate group of ATP. Although the calculated transition state geometry indicates an associative phosphoryl transfer, no concentration of negative charge is found. To study the influence of protein mutations on the reaction mechanism, we compared two-dimensional energy hypersurfaces of the protein kinase wild-type model and a corresponding mutant in which Asp166 was replaced by alanine. Surprisingly, they show similar energy profiles despite the experimentally known decrease of catalytic activity for corresponding mutants. Furthermore, a model structure was examined, where the charged NH3 group of Lys168 was replaced by a neutral methyl group. The energetic hypersurface of this hypothetical mutant shows two possible pathways for phosphoryl transfer, which both require significantly higher activation energies than the other systems investigated, while the energetic stabilization of the reaction product is similar in all systems. As the position of the amino acid side chains and the substrate peptide is virtually unchanged in all model systems, our results suggest that the exchange of Asp166 by other amino acid is less important to the phosphoryl transfer itself, but crucial to maintain the configuration of the active site in vivo. The positively charged side chain of Lys168, however, is necessary to stabilize the intermediate reaction states, particularly the side chain of the substrate peptide.

The potential for engineering stable proteins with multiple amino acid substitutions was explored. Eleven lysine, five methionine, two tryptophan, one glycine, and three threonine substitutions were simultaneously made in barley chymotrypsin inhibitor-2 (CI-2) to substantially improve the essential amino acid content of the protein. These substitutions were chosen based on the three-dimensional structure of CI-2 and an alignment of homologous sequences. The initial engineered protein folded into a wild-type-like structure, but had a free energy of unfolding of only 2.2 kcal/mol, considerably less than the wild-type value of 7.5 kcal/mol. Restoration of the lysine mutation at position 67 to the wild-type arginine increased the free energy of unfolding to 3.1 kcal/mol. Subsequent cysteine substitutions at positions 22 and 82 resulted in disulfide bond formation and a protein with nearly wild-type thermodynamic stability (7.0 kcal/mol). None of the engineered proteins retained inhibitory activity against chymotrypsin or elastase, and all had substantially reduced inhibitory activity against subtilisin. The proteolytic stabilities of the proteins correlated with their thermodynamic stabilities. Reduction of the disulfide bond resulted in substantial loss of both thermodynamic and proteolytic stabilities, confirming that the disulfide bond, and not merely the cysteine substitutions, was responsible for the increased stability. We conclude that it is possible to replace over a third of the residues in CI-2 with minimal disruption of stability and structural integrity.

The biotin carboxyl carrier protein (BCCP) is a subunit of acetyl-CoA carboxylase, a biotin-dependent enzyme that catalyzes the first committed step of fatty acid biosynthesis. In its functional cycle, this protein engages in heterologous protein–protein interactions with three distinct partners, depending on its state of post-translational modification. Apo-BCCP interacts specifically with the biotin holoenzyme synthetase, BirA, which results in the post-translational attachment of biotin to a single lysine residue on BCCP. Holo-BCCP then interacts with the biotin carboxylase subunit of acetyl-CoA carboxylase, which leads to the addition of the carboxylate group of bicarbonate to biotin. Finally, the carboxy-biotinylated form of BCCP interacts with transcarboxylase in the transfer of the carboxylate to acetyl-CoA to form malonyl-CoA. The determinants of protein–protein interaction specificity in this system are unknown. The NMR solution structure of the unbiotinylated form of an 87 residue C-terminal domain fragment (residue 70–156) of BCCP (holoBCCP87) and the crystal structure of the biotinylated form of a C-terminal fragment (residue 77–156) of BCCP from Escherichia coli acetyl-CoA carboxylase have previously been determined. Comparative analysis of these structures provided evidence for small, localized conformational changes in the biotin-binding region upon biotinylation of the protein. These structural changes may be important for regulating specific protein–protein interactions. Since the dynamic properties of proteins are correlated with local structural environments, we have determined the relaxation parameters of the backbone 15N nuclear spins of holoBCCP87, and compared these with the data obtained for the apo protein. The results indicate that upon biotinylation, the inherent mobility of the biotin-binding region and the protruding thumb, with which the biotin group interacts in the holo protein, are significantly reduced.

The serine protease subtilisin BPN′ is a useful catalyst for peptide synthesis when dissolved in high concentrations of a water-miscible organic co-solvent such as N,N-dimethylformamide (DMF). However, in 50% DMF, the kcat for amide hydrolysis is two orders of magnitude lower than in aqueous solution. Surprisingly, the kcat for ester hydrolysis is unchanged in 50% DMF. To explain this alteration in activity, the structure of subtilisin 8397+1 was determined in 20, 35, and 50% (v/v) DMF to 1.8 Å resolution. In 50% DMF, the imidazole ring of His64, the central residue of the catalytic triad, has rotated approximately 180° around the Cβ-Cγ bond. Two new water molecules in the active site stabilize the rotated conformation. This rotation places His64 in an unfavorable geometry to interact with the other members of the catalytic triad, Ser221 and Asp32. NMR experiments confirm that the characteristic resonance due to the low barrier hydrogen bond between the His64 and Asp32 is absent in 50% DMF. These experiments provide a clear structural basis for the change in activity of serine proteases in organic co-solvents.