Much computational research aimed at understanding ionizable group interactions in proteins has focused on numerical solutions of the Poisson–Boltzmann (PB) equation, incorporating protein exclusion zones for solvent and counterions in a continuum model. Poor agreement with measured pKas and pH-dependent stabilities for a (protein, solvent) relative dielectric boundary of (4,80) has lead to the adoption of an intermediate (20,80) boundary. It is now shown that a simple Debye–Hückel (DH) calculation, removing both the low dielectric and counterion exclusion regions associated with protein, is equally effective in general pKa calculations. However, a broad-based discrepancy to measured pH-dependent stabilities is maintained in the absence of ionizable group interactions in the unfolded state. A simple model is introduced for these interactions, with a significantly improved match to experiment that suggests a potential utility in predicting and analyzing the acid pH-dependence of protein stability. The methods are applied to the relative pH-dependent stabilities of the pore-forming domains of colicins A and N. The results relate generally to the well-known preponderance of surface ionizable groups with solvent-mediated interactions. Although numerical PB solutions do not currently have a significant advantage for overall pKa estimations, development based on consideration of microscopic solvation energetics in tandem with the continuum model could combine the large ΔpKas of a subset of ionizable groups with the overall robustness of the DH model.

A database of functional sites for proteins with known structures, SITE, is constructed and used in conjunction with a simple pattern matching program SiteMatch to evaluate possible function conservation in a recently constructed database of fold predictions for Escherichia coli proteins (Rychlewski L et al., 1999, Protein Sci 8:614–624). In this and other prediction databases, fold predictions are based on algorithms that can recognize weak sequence similarities and putatively assign new proteins into already characterized protein families. It is not clear whether such sequence similarities arise from distant homologies or general similarity of physicochemical features along the sequence. Leaving aside the important question of nature of relations within fold superfamilies, it is possible to assess possible function conservation by looking at the pattern of conservation of crucial functional residues. SITE consists of a multilevel function description based on structure annotations and structure analyses. In particular, active site residues, ligand binding residues, and patterns of hydrophobic residues on the protein surface are used to describe different functional features. SiteMatch, a simple pattern matching program, is designed to check the conservation of residues involved in protein activity in alignments generated by any alignment method. Here, this procedure is used to study conservation of functional features in alignments between protein sequences from the E. coli genome and their optimal structural templates. The optimal templates were identified and alignments taken from the database of genomic structural predictions was described in a previous publication (Rychlewski L et al., 1999, Protein Sci 8:614–624). An automated assessment of function conservation is used to analyze the relation between fold and function similarity for a large number of fold predictions. For instance, it is shown that identifying low significance predictions with a high level of functional residue conservations can be used to extend the prediction sensitivity for fold prediction methods. Over 100 new fold/function predictions in this class were obtained in the E. coli genome. At the same time, about 30% of our previous fold predictions are not confirmed as function predictions, further highlighting the problem of function divergence in fold superfamilies.

The motor domain regions of three novel members of the kinesin superfamily TLKIF1, TLKIFC, and TLBIMC were identified in a thermophilic fungus Thermomyces lanuginosus. Based on sequence similarity, they were classified as members of the known kinesin families Unc104/KIF1, KAR3, and BIMC. TLKIF1 was subsequently expressed in Escherichia coli. The expression level was high, and the protein was mostly soluble, easy to purify, and enzymatically active. TLKIF1 is a monomeric kinesin motor, which in a gliding motility assay displays a robust plus-directed microtubule movement up to 2 μm/s. The discovery of TLKIF1 also demonstrates that a family of kinesin motors not previously found in fungi may in fact be used in this group of organisms.

Outer membrane protein A (OmpA) of Escherichia coli is a β-barrel membrane protein that unfolds in 8 M urea to a random coil. OmpA refolds upon urea dilution in the presence of certain detergents or lipids. To examine the minimal requirements for secondary and tertiary structure formation in β-barrel membrane proteins, folding of OmpA was studied as a function of the hydrophobic chain length, the chemical structure of the polar headgroup, and the concentration of a large array of amphiphiles. OmpA folded in the presence of detergents only above a critical minimal chain length of the apolar chain as determined by circular dichroism spectroscopy and a SDS-PAGE assay that measures tertiary structure formation. Details of the chemical structure of the polar headgroup were unimportant for folding. The minimal chain length required for folding correlated with the critical micelle concentration in each detergent series. Therefore, OmpA requires preformed detergent micelles for folding and does not adsorb monomeric detergent to its perimeter after folding. Formation of secondary and tertiary structure is thermodynamically coupled and strictly dependent on the interaction with aggregated amphiphiles.

We have performed a detailed study of methanol-induced conformational transitions of horse heart apomyoglobin (apoMb) to investigate the existence of the compact and expanded denatured states. A combination of far- and near-ultraviolet circular dichroism, NMR spectroscopy, and small-angle X-ray scattering (SAXS) was used, allowing a phase diagram to be constructed as a function of pH and the methanol concentration. The phase diagram contains four conformational states, the native (N), acid-denatured (UA), compact denatured (IM), and expanded helical denatured (H) states, and indicates that the compact denatured state (IM) is stable under relatively mild denaturing conditions, whereas the expanded denatured states (UA and H) are realized under extreme conditions of pH (strong electric repulsion) or alcohol concentration (weak hydrophobic interaction). The results of this study, together with many previous studies in the literature, indicate the general existence of the compact denatured states not only in the salt-pH plane but also in the alcohol-pH plane. Furthermore, to determine the general feature of the H conformation we used several proteins including ubiquitin, ribonuclease A, α-lactalbumin, β-lactoglobulin, and Streptomyces subtilisin inhibitor (SSI) in addition to apoMb. SAXS studies of these proteins in 60% methanol showed that the H states of these all proteins have expanded and nonglobular conformations. The qualitative agreement of the experimental data with computer-simulated Kratky profiles also supports this structural feature of the H state.

Human pancreatic α-amylase (HPA) was expressed in the methylotrophic yeast Pichia pastoris and two mutants (D197A and D197N) of a completely conserved active site carboxylic acid were generated. All recombinant proteins were shown by electrospray ionization mass spectrometry (ESI-MS) to be glycosylated and the site of attachment was shown to be Asn461 by peptide mapping in conjunction with ESI-MS. Treatment of these proteins with endoglycosidase F demonstrated that they contained a single N-linked oligosaccharide and yielded a protein product with a single N-acetyl glucosamine (GlcNAc), which could be crystallized. Solution of the crystal structure to a resolution of 2.0 Å confirmed the location of the glycosyl group as Asn461 and showed that the recombinant protein had essentially the same conformation as the native enzyme. The kinetic parameters of the glycosylated and deglycosylated wild-type proteins were the same while the kcat/Km values for D197A and D197N were 106–107 times lower than the wild-type enzyme. The decreased kcat/Km values for the mutants confirm that D197 plays a crucial role in the hydrolytic activity of HPA, presumably as the catalytic nucleophile.

Staphylococcal nuclease forms three different partially-folded intermediates at low pH in the presence of low to moderate concentration of anions, differing in the amount of secondary structure, globularity, stability, and compactness. Although these intermediates are monomeric at low protein concentration (≤0.25 mg/mL), increasing concentrations of protein result in the formation of dimers and soluble oligomers, ultimately leading to larger insoluble aggregates. Unexpectedly, increasing protein concentration not only led to association, but also to increased structure of the intermediates. The secondary structure, stability, and globularity of the two less-ordered partially-folded intermediates (A1 and A2) were substantially increased upon association, suggesting that aggregation induces structure. An excellent correlation was found between degree of association and amount of structure measured by different techniques, including circular dichroism, fluorescence, Fourier transform infrared spectroscopy (FTIR), and small-angle X-ray scattering. The associated states were also substantially more stable toward urea denaturation than the monomeric forms. A mechanism is proposed, in which the observed association of monomeric intermediates involves intermolecular interactions which correspond to those found intramolecularly in normal folding to the native state.

In the recently determined crystal structure of human thioredoxin (hTRX), the protein is a dimer, covalently linked via an intermolecular disulfide bond involving Cys73 from each monomer (Weichsel et al., 1996). The active site loop, comprising Trp31-Cys-Gly-Pro-Cys35, forms part of the dimer interface and is blocked. Surprisingly, even a mutant protein in which Cys73 is replaced by Ser, is dimeric in the crystal, which led to the suggestion that the dimer of hTRX may be of physiological importance (Weichsel et al., 1996). In addition, the role of Asp60 in dimer formation was probed both crystallographically and with a dimerization assay using diamide as the oxidant (Andersen et al., 1997). Since no evidence for the presence of a dimer was found in the solution NMR structure of hTRX (Qin et al., 1994), it was suggested that the threonine for methionine substitution at position 74 in the hTRX used for the NMR studies might be responsible for this difference between crystal and solution states (Weichsel et al., 1996). To resolve this issue, we have examined the molecular weight and rotational correlation time of hTRX by analytical ultracentrifugation and NMR spectroscopy, respectively. Two variants were investigated, namely hTRX identical in amino acid sequence to the one whose NMR structure we previously determined (C62A, C69A, C73A, M74T) and hTRX (C62A, C69A, C73A, M74) containing the wild-type amino acid methionine at position 74.

The hydrophilic subunit of the mannose transporter (IIABMan) of Escherichia coli is a homodimer that contains four tryptophans per monomer, three in the N-terminal domain (Trp12, Trp33, and Trp69) and one in the C-terminal domain (Trp182). Single and double Trp-Phe mutants of IIABMan and of the IIA domain were produced. Fluorescence emission studies revealed that Trp33 and Trp12 are the major fluorescence emitters, Trp69 is strongly quenched in the native protein and Trp182 strongly blue shifted, indicative of a hydrophobic environment. Stabilities of the Trp mutants of dimeric IIAMan and IIABMan were estimated from midpoints of the GdmHCl-induced unfolding transitions and from the amount of dimers that resisted dissociation by SDS (sodium dodecyl sulfate), respectively. W12F exhibited increased stability, but only 6% of the wild-type phosphotransferase activity, whereas W33F was marginally and W69F significantly destabilized, but fully active. Second site mutations W33F and W69F in the background of the W12F mutation reduced protein stability and suppressed the functional defect of W12F. These results suggest that flexibility is required for the adjustments of protein–protein contacts necessary for the phosphoryltransfer between the phosphorylcarrier protein HPr, IIAMan, IIBMan, and the incoming mannose bound to the transmembrane IICMan–IIDMan complex.

The effects of histidine residue placement in a de novo-designed four-α-helix bundle are investigated by placement of histidine residues at coiled coil heptad a positions in two distinct heptads and at each position within a single heptad repeat of our prototype heme protein maquette, [H10H24]2 [{Ac-CGGGELWKL·HEELLKK·FEELLKL·HEERLKK·L-CONH2}2]2 composed of a generic (α-SS-α)2 peptide architecture. The heme to peptide stoichiometry of variants of [H10H24]2 with either or both histidines on each helix replaced with noncoordinating alanine residues ([H10A24]2, [A10H24]2, and [A10A24]2) demonstrates the obligate requirement of histidine for biologically significant heme affinity. Variants of [A10A24]2, [{Ac-CGGGELWKL·AEELLKK·FEELLKL·AEERLKK·L-CONH2}2]2, containing a single histidine per helix in positions 9 to 15 were evaluated to verify the design based on molecular modeling. The bis-histidine site formed between heptad positions a at 10 and 10′ bound ferric hemes with the highest affinity, Kd1 and Kd2 values of 15 and 800 nM, respectively. Placement of histidine at position 11 (heptad position b) resulted in a protein that bound a single heme with moderate affinity, Kd1 of 9.5 μM, whereas the other peptides had no measurable apparent affinity for ferric heme with Kd1 values >200 μM. The bis-histidine ligation of heme to [H10A24]2 and [H11A24]2 was confirmed by electron paramagnetic resonance spectroscopy. The protein design rules derived from this study, together with the narrow tolerances revealed, are applicable for improving future heme protein designs, for analyzing the results of randomized heme protein combinatorial libraries, as well as for implementation in automated protein design.

The N-terminal 17 residues of ubiquitin have been shown by 1H NMR to fold autonomously into a β-hairpin structure in aqueous solution. This structure has a specific, native-like register, though side-chain contacts differ in detail from those observed in the intact protein. An autonomously folding hairpin has previously been identified in the case of streptococcal protein G, which is structurally homologous with ubiquitin, but remarkably, the two are not in topologically equivalent positions in the fold. This suggests that the organization of folding may be quite different for proteins sharing similar tertiary structures. Two smaller peptides have also been studied, corresponding to the isolated arms of the N-terminal hairpin of ubiquitin, and significant differences from simple random coil predictions observed in the spectra of these subfragments, suggestive of significant limitation of the backbone conformational space sampled, presumably as a consequence of the strongly β-structure favoring composition of the sequences. This illustrates the ability of local sequence elements to express a propensity for β-structure even in the absence of actual sheet formation. Attempts were made to estimate the population of the folded state of the hairpin, in terms of a simple two-state folding model. Using published “random coil” values to model the unfolded state, and values derived from native ubiquitin for the putative unique, folded state, it was found that the apparent population varied widely for different residues and with different NMR parameters. Use of the spectra of the subfragment peptides to provide a more realistic model of the unfolded state led to better agreement in the estimates that could be obtained from chemical shift and coupling constant measurements, while making it clear that some other approaches to population estimation could not give meaningful results, because of the tendency to populate the β-region of conformational space even in the absence of the hairpin structure.

In the past decade or so, the universal dominance of l- over d-amino acids in living organisms has been questioned. This interesting and wide-ranging book summarizes the current state of knowledge on the occurrence, biosynthesis, and roles of d-amino acid containing peptides, with a strong emphasis on those in multicellular organisms. The opening chapter by Scaloni et al. provides a practical guide to the detection and analysis of d-amino acids. The subsequent four chapters describe studies related to the presence of d-amino acids in antimicrobial and opioid peptides in frogs, molluscan neuropeptides, and crustacean hyperglycemic hormones. Although the work could have been summarized more succinctly and there is, as is common with this type of text, some overlap between the chapters, the inclusion of experimental details will be of interest to the specialist.

Automated methods for identifying and characterizing regular β-barrels from coordinate data have been developed to analyze and classify various kinds of barrel structures based on geometric parameters such as the barrel strand number (n) and shear number (S). In total, we find 1,316 barrels in the January 1998 release of Protein Data Bank. Of 1,316 barrels, 1,277 barrels had an even shear number, corresponding to 50 nonhomologous families. The (βα)8 triose phosphate isomerase (TIM) barrel (n = 8, S = 8) fold has the largest number of apparently nonhomologous entries, 16, although the trypsin like antiparallel (n = 6, S = 8) barrels (representing only three families) are the most common with 527 barrels. Of all the protein families that exhibit barrel structures, 68% are found to be various kinds of enzymes, the remainder being binding proteins and transport membrane proteins. In addition, the layers of side chains, which form the cores of barrels with S = n and S = 2n, are also analyzed. More sophisticated methods were developed for detecting TIM barrels specifically, including consideration of the amino acid propensities for the side chains that form the layers. We found that the residues on the outside of the eight stranded parallel β-barrel, buried by the α-helices, are much more hydrophobic than the residues inside the barrel.

In the fold recognition approach to structure prediction, a sequence is tested for compatibility with an already known fold. For membrane proteins, however, few folds have been determined experimentally. Here the feasibility of computing the vast majority of likely membrane protein folds is tested. The results indicate that conformation space can be effectively sampled for small numbers of helices. The vast majority of potential monomeric membrane protein structures can be represented by about 30-folds for three helices, but increases exponentially to about 1,500,000 folds for seven helices. The generated folds could serve as templates for fold recognition or as starting points for conformational searches that are well distributed throughout conformation space.

Urea-induced unfolding of Escherichia coli citrate synthase occurs in two phases, as monitored by circular dichroism at 222 nm (measuring secondary structure) or by tryptophan fluorescence. In this paper we characterize the intermediate state, which retains about 40% of the ellipticity of the native state, and is stable between 2.5 M and 5.5 M urea, approximately. This intermediate binds significant amounts of the probe for hydrophobic surfaces, anilinonaphthalene sulfonate, but forms aggregates at least as high as an octamer, as shown by transverse urea gradient polyacrylamide electrophoresis. Thermal denaturation of E. coli citrate synthase also produces an intermediate at temperatures near 60 °C, which also retains about 40% of the native ellipticity and forms aggregates, as measured by electrospray-ionization/time-of-flight mass spectrometry. We have used a collection of “cavity-forming” mutant proteins, in which bulky buried hydrophobic residues are replaced by alanines, to explore the nature of the intermediate state further. A certain amount of these mutant proteins shows a destabilized intermediate, as measured by the urea concentration range in which the intermediate is observed. These mutants are found in parts of the citrate synthase sequence that, in a native state, form helices G, M, N, Q, R, and S. From this and other evidence, it is argued that the intermediate state is an aggregated state in which these six helices, or parts of them, remain folded, and that formation of this intermediate is also likely to be a key step in the folding of E. coli citrate synthase.

Pokeweed antiviral protein (PAP) is a ribosome-inactivating protein (RIP), which enzymatically removes a single adenine base from a conserved, surface exposed loop sequence of ribosomal rRNA. We now present unprecedented experimental evidence that PAP can release not only adenine but guanine as well from Escherichia coli rRNA, albeit at a rate 20 times slower than for adenine. We also report X-ray structure analysis and supporting modeling studies for the interactions of PAP with guanine. Our modeling studies indicated that PAP can accommodate a guanine base in the active site pocket without large conformational changes. This prediction was experimentally confirmed, since a guanine base was visible in the active site pocket of the crystal structure of the PAP-guanine complex.

NMR offers the possibility of accurate secondary structure for proteins that would be too large for structure determination. In the absence of an X-ray crystal structure, this information should be useful as an adjunct to protein fold recognition methods based on low resolution force fields. The value of this information has been tested by adding varying amounts of artificial secondary structure data and threading a sequence through a library of candidate folds. Using a literature test set, the threading method alone has only a one-third chance of producing a correct answer among the top ten guesses. With realistic secondary structure information, one can expect a 60–80% chance of finding a homologous structure. The method has then been applied to examples with published estimates of secondary structure. This implementation is completely independent of sequence homology, and sequences are optimally aligned to candidate structures with gaps and insertions allowed. Unlike work using predicted secondary structure, we test the effect of differing amounts of relatively reliable data.

Hepatitis C virus (HCV) nonstructural protein 3 (NS3) has been shown to possess protease and helicase activities and has also been demonstrated to spontaneously associate with nonstructural protein NS4A (NS4A) to form a stable complex. Previous attempts to produce the NS3/NS4A complex in recombinant baculovirus resulted in a protein complex that aggregated and precipitated in the absence of nonionic detergent and high salt. A single-chain form of the NS3/NS4A complex (His-NS4A21–32–GSGS-NS33–631) was constructed in which the NS4A core peptide is fused to the N-terminus of the NS3 protease domain as previously described (Taremi et al., 1998). This protein contains a histidine tagged NS4A peptide (a.a. 21–32) fused to the full-length NS3 (a.a. 3–631) through a flexible tetra amino acid linker. The recombinant protein was expressed to high levels in Escherichia coli, purified to homogeneity, and examined for NTPase, nucleic acid unwinding, and proteolytic activities. The single-chain recombinant NS3-NS4A protein possesses physiological properties equivalent to those of the NS3/NS4A complex except that this novel construct is stable, soluble and sixfold to sevenfold more active in unwinding duplex RNA. Comparison of the helicase activity of the single-chain recombinant NS3-NS4A with that of the full-length NS3 (without NS4A) and that of the helicase domain alone suggested that the presence of the protease domain and at least the NS4A core peptide are required for optimal unwinding activity.

We propose a new approach that permits correlation of specific domains defined by their primary sequence with their location in the structure of complex macromolecular aggregates. It is based on the combination of well-established structural analysis methods that incorporate the use of overlapping peptides on cellulose membranes for the isolation and purification of specific antibodies from a polyclonal antiserum. Monospecific antibodies to the connector protein of bacteriophage φ29 were isolated from polyclonal antisera using a new development of the spotscan method. These antibodies can be purified in quantities that allow antigenicity testing in enzyme-linked immunosorbent assays, Western blotting and immunoprecipitations, demonstrating the specificity of this isolation procedure. This approach has allowed us to generate direct antibody probes for immunoelectron microscopy mapping of different connector protein domains in a low resolution three-dimensional epitope map.

A novel strategy to obtain high-level production of mature proteins exported to the periplasm of Escherichia coli is described. It is based on a modified signal sequence generated by insertion of a coding sequence of the polypeptide precursor of interest at the BamHI site of the commercial vector pQE-30 resulting in an addition of a dodeca-peptide (MRGSH6GS) at the N-terminus of the precursor. The modification does not affect correct processing of the modified signal nor proper folding of the target protein, resulting in an untagged native product. The method is simple for avoiding onerous optimization of translation initiation and screening of host stains. The usefulness of this method is illustrated by overexpression of DsbC and DsbA. Induced by 0.01 mM IPTG at 37 °C, proteins were overproduced to comprise 20–30% of the total cellular proteins, and more than 95% of the expressed proteins were correctly processed and exported into the periplasm with yields of more than 100 mg per liter culture.