CrmA is an unusual viral serpin that inhibits both cysteine and serine proteinases involved in the regulation of host inflammatory and apoptosis processes. It differs from other members of the serpin superfamily by having a reactive center loop that is one residue shorter, and by its apparent inability to form SDS-stable covalent complexes with cysteine proteinases. To obtain insight into the inhibitory mechanism of crmA, we determined the crystal structure of reactive center loop-cleaved crmA to 2.9 Å resolution. The structure, which is the first of a viral serpin, suggests that crmA can inhibit cysteine proteinases by a mechanism analogous to that used by other serpins against serine proteinases. However, one striking difference from other serpins, which may be significant for in vivo function, is an additional highly charged antiparallel strand for β sheet A, whose sequence and length are unique to crmA.

Decorsin is an antagonist of integrin αIIb β3 and a potent platelet aggregation inhibitor. A synthetic gene encoding decorsin, originally isolated from the leech Macrobdella decora, was designed, constructed, and expressed in Escherichia coli. The synthetic gene was fused to the stII signal sequence and expressed under the transcriptional control of the E. coli alkaline phosphatase promoter. The protein was purified by size-exclusion filtration of the periplasmic contents followed by reversed-phase high-performance liquid chromatography. Purified recombinant decorsin was found to be indistinguishable from leech-derived decorsin based on amino acid composition, mass spectral analysis, and biological activity assays. Complete sequential assignments of 1H and proton bound 13C resonances were established. Stereospecific assignments of 21 of 25 nondegenerate β-methylene groups were determined. The RGD adhesion site recognized by integrin receptors was found at the apex of a most exposed hairpin loop. The dynamic behavior of decorsin was analyzed using several independent NMR parameters. Although the loop containing the RGD sequence is the most flexible one in decorsin, the conformation of the RGD site itself is more restricted than in other proteins with similar activities.

Blue copper proteins are type-I copper-containing redox proteins whose role is to shuttle electrons from an electron donor to an electron acceptor in bacteria and plants. A large amount of experimental data is available on blue copper proteins; however, their functional characterization is hindered by the complexity of redox processes in biological systems. We describe here the application of a semiquantitative method based on a comparative analysis of molecular interaction fields to gain insights into the recognition properties of blue copper proteins. Molecular electrostatic and hydrophobic potentials were computed and compared for a set of 33 experimentally-determined structures of proteins from seven blue copper subfamilies, and the results were quantified by means of similarity indices. The analysis provides a classification of the blue copper proteins and shows that (1) comparison of the molecular electrostatic potentials provides useful information complementary to that highlighted by sequence analysis; (2) similarities in recognition properties can be detected for proteins belonging to different subfamilies, such as amicyanins and pseudoazurins, that may be isofunctional proteins; (3) dissimilarities in interaction properties, consistent with experimentally different binding specificities, may be observed between proteins belonging to the same subfamily, such as cyanobacterial and eukaryotic plastocyanins; (4) proteins with low sequence identity, such as azurins and pseudoazurins, can have sufficient similarity to bind to similar electron donors and acceptors while having different binding specificity profiles.

Human serum albumin (HSA) interacts with a vast array of chemically diverse ligands at specific binding sites. To pinpoint the essential structural elements for the formation of the warfarin binding site on human serum albumin, a defined set of five recombinant proteins comprising combinations of domains and/or subdomains of the N-terminal part were prepared and characterized by biochemical standard procedures, tryptophanyl fluorescence, and circular dichroic measurements, indicating well-preserved secondary and tertiary structures. Affinity constants for binding to warfarin were estimated by fluorescence titration experiments and found to be highest for HSA-DOM I-II and HSA, followed by HSA-DOM IB-II, HSA-DOM II, and HSA-DOM I-IIA. In addition, ultraviolet difference spectroscopy and induced circular dichroism experiments were carried out to get an in depth understanding of the binding mechanism of warfarin to the fragments as stand-alone proteins. This systematic study indicates that the primary warfarin binding site is centered in subdomain IIA with indispensable structural contributions of subdomain IIB and domain I, while domain III is not involved in this binding site, underlining the great potential that lies in the use of combinations of recombinant fragments for the study and accurate localization of ligand binding sites on HSA.

The folding kinetics of human common-type acylphosphatase (cAcP) from its urea- and TFE-denatured states have been determined by stopped-flow fluorescence techniques. The refolding reaction from the highly unfolded state formed in urea is characterized by double exponential behavior that includes a slow phase associated with isomerism of the Gly53–Pro54 peptide bond. However, this slow phase is absent when refolding is initiated by dilution of the highly α-helical denatured state formed in the presence of 40% trifluoroethanol (TFE). NMR studies of a peptide fragment corresponding to residues Gly53–Gly69 of cAcP indicate that only the native-like trans isomer of the Gly–Pro peptide bond is significantly populated in the presence of TFE, whereas both the cis and trans isomers are found in an ∼1:9 ratio for the peptide bond in aqueous solution. Molecular modeling studies in conjunction with NMR experiments suggest that the trans isomer of the Gly53–Pro54 peptide bond is stabilized in TFE by the formation of a nonnative-like hydrogen bond between the CO group of Gly53 and the NH group of Lys57. These results therefore reveal that a specific nonnative interaction in the denatured state can increase significantly the overall efficiency of refolding.

The rubredoxin from the cryptomonad Guillardia theta is one of the first examples of a rubredoxin encoded in a eukaryotic organism. The structure of a soluble zinc-substituted 70-residue G. theta rubredoxin lacking the membrane anchor and the thylakoid targeting sequence was determined by multidimensional heteronuclear NMR, representing the first three-dimensional (3D) structure of a eukaryotic rubredoxin. For the structure calculation a strategy was applied in which information about hydrogen bonds was directly inferred from a long-range HNCO experiment, and the dynamics of the protein was deduced from heteronuclear nuclear Overhauser effect data and exchange rates of the amide protons. The structure is well defined, exhibiting average root-mean-square deviations of 0.21 Å for the backbone heavy atoms and 0.67 Å for all heavy atoms of residues 7–56, and an increased flexibility toward the termini. The structure of this core fold is almost identical to that of prokaryotic rubredoxins. There are, however, significant differences with respect to the charge distribution at the protein surface, suggesting that G. theta rubredoxin exerts a different physiological function compared to the structurally characterized prokaryotic rubredoxins. The amino-terminal residues containing the putative signal peptidase recognition/cleavage site show an increased flexibility compared to the core fold, but still adopt a defined 3D orientation, which is mainly stabilized by nonlocal interactions to residues of the carboxy-terminal region. This orientation might reflect the structural elements and charge pattern necessary for correct signal peptidase recognition of the G. theta rubredoxin precursor.

Several recent publications illustrated advantages of using sequence profiles in recognizing distant homologies between proteins. At the same time, the practical usefulness of distant homology recognition depends not only on the sensitivity of the algorithm, but also on the quality of the alignment between a prediction target and the template from the database of known proteins. Here, we study this question for several supersensitive protein algorithms that were previously compared in their recognition sensitivity (Rychlewski et al., 2000). A database of protein pairs with similar structures, but low sequence similarity is used to rate the alignments obtained with several different methods, which included sequence–sequence, sequence–profile, and profile–profile alignment methods. We show that incorporation of evolutionary information encoded in sequence profiles into alignment calculation methods significantly increases the alignment accuracy, bringing them closer to the alignments obtained from structure comparison.

In general, alignment quality is correlated with recognition and alignment score significance. For every alignment method, alignments with statistically significant scores correlate with both correct structural templates and good quality alignments. At the same time, average alignment lengths differ in various methods, making the comparison between them difficult. For instance, the alignments obtained by FFAS, the profile–profile alignment algorithm developed in our group are always longer that the alignments obtained with the PSI-BLAST algorithms. To address this problem, we develop methods to truncate or extend alignments to cover a specified percentage of protein lengths. In most cases, the elongation of the alignment by profile–profile methods is reasonable, adding fragments of similar structure. The examples of erroneous alignment are examined and it is shown that they can be identified based on the model quality.

The native form of some proteins such as strained plasma serpins (serine protease inhibitors) and the spring-loaded viral membrane fusion proteins are in a metastable state. The metastable native form is thought to be a folding intermediate in which conversion into the most stable state is blocked by a very high kinetic barrier. In an effort to understand how the spontaneous conversion of the metastable native form into the most stable state is prevented, we designed mutations of α1-antitrypsin, a prototype serpin, which can bypass the folding barrier. Extending the reactive center loop of α1-antitrypsin converts the molecule into a more stable state. Remarkably, a 30-residue loop extension allows conversion into an extremely stable state, which is comparable to the relaxed cleaved form. Biochemical data strongly suggest that the strain release is due to the insertion of the reactive center loop into the major β-sheet, A sheet, as in the known stable conformations of serpins. Our results clearly show that extending the reactive center loop is sufficient to bypass the folding barrier of α1-antitrypsin and suggest that the constrain held by polypeptide connection prevents the conversion of the native form into the lowest energy state.

The intermolecular contact regions between monomers of the homodimeric DNA binding protein ParR and the interaction between the glycoproteins CD28 and CD80 were investigated using a strategy that combined chemical cross-linking with differential MALDI-MS analyses. ParR dimers were modified in vitro with the thiol-cleavable cross-linker 3,3′-dithio-bis(succinimidylproprionate) (DTSSP), proteolytically digested with trypsin and analyzed by MALDI-MS peptide mapping. Comparison of the peptide maps obtained from digested cross-linked ParR dimers in the presence and absence of a thiol reagent strongly supported a “head-to-tail” arrangement of the monomers in the dimeric complex. Glycoprotein fusion constructs CD28-IgG and CD80-Fab were cross-linked in vitro by DTSSP, characterized by nonreducing SDS-PAGE, digested in situ with trypsin and analyzed by MALDI-MS peptide mapping (± thiol reagent). The data revealed the presence of an intermolecular cross-link between the receptor regions of the glycoprotein constructs, as well as a number of unexpected but nonetheless specific interactions between the fusion domains of CD28-IgG and the receptor domain of CD80-Fab. The strategy of chemical cross-linking combined with differential MALDI-MS peptide mapping (± thiol reagent) enabled localization of the interface region(s) of the complexes studied and clearly demonstrates the utility of such an approach to obtain structural information on interacting noncovalent complexes.

Chemical and thermal denaturation of calmodulin has been monitored spectroscopically to determine the stability for the intact protein and its two isolated domains as a function of binding of Ca2+ or Mg2+. The reversible urea unfolding of either isolated apo-domain follows a two-state mechanism with relatively low ΔG°20 values of ∼2.7 (N-domain) and ∼1.9 kcal/mol (C-domain). The apo-C-domain is significantly unfolded at normal temperatures (20–25 °C). The greater affinity of the C-domain for Ca2+ causes it to be more stable than the N-domain at [Ca2+] ≥0.3 mM. By contrast, Mg2+ causes a greater stabilization of the N- rather than the C-domain, consistent with measured Mg2+ affinities.

For the intact protein (±Ca2+), the bimodal denaturation profiles can be analyzed to give two ΔG°20 values, which differ significantly from those of the isolated domains, with one domain being less stable and one domain more stable. The observed stability of the domains is strongly dependent on solution conditions such as ionic strength, as well as specific effects due to metal ion binding. In the intact protein, different folding intermediates are observed, depending on the ionic composition. The results illustrate that a protein of low intrinsic stability is liable to major perturbation of its unfolding properties by environmental conditions and liganding processes and, by extension, mutation. Hence, the observed stability of an isolated domain may differ significantly from the stability of the same structure in a multidomain protein. These results address questions involved in manipulating the stability of a protein or its domains by site directed mutagenesis and protein engineering.

The biotin holoenzyme synthetases (BHS) are essential enzymes in all organisms that catalyze post-translational linkage of biotin to biotin-dependent carboxylases. The primary sequences of a large number of these enzymes are now available and homologies are found among all. The glycine-rich sequence, GRGRXG, constitutes one of the homologous regions in these enzymes and, based on its similarity to sequences found in a number of mononucleotide binding enzymes, has been proposed to function in ATP binding in the BHSs. In the Escherichia coli enzyme, the only member of the family for which a three-dimensional structure has been determined, the conserved sequence is found in a partially disordered surface loop. Mutations in the sequence have previously been isolated and characterized in vivo. In this work these single-site mutants, G115S, R118G, and R119W, of the E. coli BHS have been purified and biochemically characterized with respect to binding of small molecule substrates and the intermediate in the biotinylation reaction. Results of this characterization indicate that, rather than functioning in ATP binding, this glycine-rich sequence is required for binding the substrate biotin and the intermediate in the biotinylation reaction, biotinyl-5′-AMP. These results are of general significance for understanding structure-function relationships in biotin holoenzyme synthetases.

We have used the homonuclear Overhauser effect (NOE) to characterize a model protein: carbonic anhydrase B. We have obtained NOE difference spectra for this protein, centering the on-resonance signals either at the methyl-proton or at the water-proton signals. The spin-diffusion spectra obtained as a function of protein concentration and temperature provide direct evidence of much greater protein–water interaction in the molten-globule state than in the native and denatured states. Furthermore, although the protein loses its gross tertiary structure in both the molten-globule and denatured states, it remains almost as compact in its molten-globule state as it is in the native state. The spin-diffusion spectra, obtained as a function of a variable delay time after the saturation pulse, allowed us to measure the relaxation times of several types of proton in the solution. These spectra contain enough information to distinguish between those water molecules solvating the protein and the free ones present as bulk water.

A 38-residue protein associated with cholesteryl ester transfer inhibition has been identified in baboons (Papio sp.). The cholesteryl ester transfer inhibitor protein (CETIP) corresponds to the N-terminus of baboon apoC-I. Relative to CETIP, baboon apoC-I is a weak inhibitor of baboon cholesteryl ester transferase (CET). To study the structural features responsible for CET inhibition, CETIP was synthesized by solid-phase methods. Using sodium dodecyl sulfate (SDS) to model the lipoprotein environment, the solution structure of CETIP was probed by optical and 1H NMR spectroscopy. Circular dichroism data show that the protein lacks a well-defined structure in water but, upon the addition of SDS, becomes helical (56%). A small blue shift of 8 nm was observed in the intrinsic tryptophan fluorescence of CETIP in the presence of saturating amounts of SDS, suggesting that tryptophan-23 is not buried deeply in the lipid environment. The helical nature of CETIP in the presence of SDS was confirmed by upfield 1Hα secondary shifts and an average solution structure determined by distance geometry/simulated annealing calculations using 476 NOE-based distance restraints. The backbone (N-Cα-C=O) root-mean-square deviation of an ensemble of 17 out of 25 calculated structures superimposed on the average structure was 1.06 ± 0.30 Å using residues V4–P35 and 0.51 ± 0.17 Å using residues A7–S32. Although the side-chain orientations fit the basic description of a class A amphipathic helix, both intramolecular salt bridge formation and “snorkeling” of basic side chains toward the polar face play minor, if any, roles in stabilizing the lipid-bound amphipathic structure. Conformational features of the calculated structures for CETIP are discussed relative to models of CETIP inhibition of cholesteryl ester transferase.

Integral membrane proteins carry out some of the most important functions of living cells, yet relatively few details are known about their structures. This is due, in large part, to the difficulties associated with preparing membrane protein crystals suitable for X-ray diffraction analysis. Mechanistic studies of membrane protein crystallization may provide insights that will aid in determining future membrane protein structures. Accordingly, the solution behavior of the bacterial outer membrane protein OmpF porin was studied by static light scattering under conditions favorable for crystal growth. The second osmotic virial coefficient (B22) was found to be a predictor of the crystallization behavior of porin, as has previously been found for soluble proteins. Both tetragonal and trigonal porin crystals were found to form only within a narrow window of B22 values located at approximately −0.5 to −2 × 10−4 mol mL g−2, which is similar to the “crystallization slot” observed for soluble proteins. The B22 behavior of protein-free detergent micelles proved very similar to that of porin-detergent complexes, suggesting that the detergent's contribution dominates the behavior of protein-detergent complexes under crystallizing conditions. This observation implies that, for any given detergent, it may be possible to construct membrane protein crystallization screens of general utility by manipulating the solution properties so as to drive detergent B22 values into the crystallization slot. Such screens would limit the screening effort to the detergent systems most likely to yield crystals, thereby minimizing protein requirements and improving productivity.

The interdomain disulfide bond present in the C-lobe of all the transferrins was postulated to restrict the domain movement resulting in the slow rate of iron uptake and release. In the present study, the conformational stability and iron binding properties of a derivative of the isolated C-lobe of ovotransferrin in which the interdomain disulfide bond, Cys478-Cys671 was selectively reduced and alkylated with iodoacetamide were compared with the disulfide intact form at the endosomal pH of 5.6. Pyrophosphate and chloride mediated iron release kinetics showed no difference between the disulfide-intact and disulfide-reduced/alkylated forms; the two protein forms yielded similar observed rate constants showing an apparent hyperbolic dependency for anion concentrations. The conformational stability evaluated by unfolding and refolding experiments was greater for the disulfide-intact form than for the disulfide-reduced/alkylated form: the ΔGDH2O values at 30 °C obtained by using urea were 9.0 ± 0.8 and 6.0 ± 0.4 kJ/mol for the former and latter protein forms, respectively, and the corresponding values obtained by using guanidine hydrochloride were 6.2 ± 0.9 and 4.3 ± 0.5 kJ/mol. The dissociation constant of iron (kd) was almost the same for the two protein forms, and it varied only subtly with urea concentrations but increased markedly with GdnHCl concentrations. The nonidentical values of ΔGDH2O and kd for urea and GdnHCl can be attributed to the ionic nature of the later denaturant, in which chloride anion may influence the structure and iron uptake-release properties of the ovotransferrin C-lobe. Taken together, we conclude that the interdomain disulfide bond has no effect on the iron uptake and release function but significantly decreases the conformational stability in the C-lobe.

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.

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.

L6, IL-TMP, and TM4SF5 are cell surface proteins predicted to have four transmembrane domains. Previous sequence analysis led to their assignment as members of the tetraspanin superfamily. In this paper, we identify a new sequence (L6D) that is strikingly similar to L6, IL-TMP, and TM4SF5. Analyses of these four sequences indicate that they are not significantly related to genuine tetraspanins, but instead constitute their own L6 superfamily.

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.

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.