The hypersensitivity resistance response directed by the N′ gene in Nicotiana sylvestris is elicited by the tobacco mosaic virus (TMV) coat protein R46G, but not by the U1 wild-type TMV coat protein. In this study, the structural and hydrodynamic properties of R46G and wild-type coat proteins were compared for variations that may explain N′ gene elicitation. Circular dichroism spectroscopy reveals no significant secondary or tertiary structural differences between the elicitor and nonelicitor coat proteins. Analytical ultracentrifugation studies, however, do show different concentration dependencies of the weight average sedimentation coefficients at 4 °C. Viral reconstitution kinetics at 20 °C were used to determine viral assembly rates and as an initial assay of the rate of 20S formation, the obligate species for viral reconstitution. These kinetic results reveal a decreased lag time for reconstitution performed with R46G that initially lack the 20S aggregate. However, experiments performed with 20S initially present reveal no detectable differences indicating that the mechanism of viral assembly is similar for the two coat protein species. Therefore, an increased rate of 20S formation from R46G subunits may explain the differences in the viral reconstitution lag times. The inferred increase in the rate of 20S formation is verified by direct measurement of the 20S boundary as a function of time at 20 °C using velocity sedimentation analysis. These results are consistent with the interpretation that there may be an altered size distribution and/or lifetime of the small coat protein aggregates in elicitors that allows N. sylvestris to recognize the invading virus.

The usefulness of molecular dynamics to assess the structural integrity of mutants containing several mutations has been investigated. Our goal was to determine whether molecular dynamics would be able to discriminate mutants of a protein having a close–to–wild-type fold, from those that are not folded under the same conditions. We used as a model the B1 domain of protein G in which we replaced the unique central α-helix by the sequence of the second β-hairpin, which has a strong intrinsic propensity to form this secondary structure in solution. In the resulting protein, one-third of the secondary structure has been replaced by a non-native one. Models of the mutants were built based on the three-dimensional structure of the wild-type GB1 domain. During 2 ns of molecular dynamics simulations on these models, mutants containing up to 10 mutations in the helix retained the native fold, while another mutant with an additional mutation unfolded. This result is in agreement with our circular dichroism and NMR experiments, which indicated that the former mutants fold into a structure similar to the wild-type, as opposed to the latter mutant which is partly unfolded. Additionally, a mutant containing six mutations scattered through the surface of the domain, and which is unfolded, was also detected by the simulation. This study suggests that molecular dynamics calculations could be performed on molecular models of mutants of a protein to evaluate their foldability, prior to a mutagenesis experiment.

Cathepsin K is a lysosomal cysteine protease belonging to the papain superfamily. It has been implicated as a major mediator of osteoclastic bone resorption. Wild-type human procathepsin K has been crystallized in a glycosylated and a deglycosylated form. The latter crystals diffract better, to 3.2 Å resolution, and contain four molecules in the asymmetric unit. The structure was solved by molecular replacement and refined to an R-factor of 0.194. The N-terminal fragment of the proregion forms a globular domain while the C-terminal segment is extended and shows substantial flexibility. The proregion interacts with the enzyme along the substrate binding groove and along the proregion binding loop (residues Ser138–Asn156). It binds to the active site in the opposite direction to that of natural substrates. The overall binding mode of the proregion to cathepsin K is similar to that observed in cathepsin L, caricain, and cathepsin B, but there are local differences that likely contribute to the specificity of these proregions for their cognate enzymes. The main observed difference is in the position of the short helix α3p (67p–75p), which occupies the S′ subsites. As in the other proenzymes, the proregion utilizes the S2 subsite for anchoring by placing a leucine side chain there, according to the specificity of cathepsin K toward its substrate.

Fructose 1,6-bisphosphate aldolase catalyzes the reversible cleavage of fructose 1,6-bisphosphate and fructose 1-phosphate to dihydroxyacetone phosphate and either glyceraldehyde 3-phosphate or glyceraldehyde, respectively. Catalysis involves the formation of a Schiff's base intermediate formed at the ε-amino group of Lys229. The existing apo-enzyme structure was refined using the crystallographic free-R-factor and maximum likelihood methods that have been shown to give improved structural results that are less subject to model bias. Crystals were also soaked with the natural substrate (fructose 1,6-bisphosphate), and the crystal structure of this complex has been determined to 2.8 Å. The apo structure differs from the previous Brookhaven-deposited structure (1ald) in the flexible C-terminal region. This is also the region where the native and complex structures exhibit differences. The conformational changes between native and complex structure are not large, but the observed complex does not involve the full formation of the Schiff's base intermediate, and suggests a preliminary hydrogen-bonded Michaelis complex before the formation of the covalent complex.

The crystal structure of the FMN-binding domain of human NADPH-cytochrome P450 reductase (P450R-FMN), a key component in the cytochrome P450 monooxygenase system, has been determined to 1.93 Å resolution and shown to be very similar both to the global fold in solution (Barsukov I et al., 1997, J Biomol NMR 10:63–75) and to the corresponding domain in the 2.6 Å crystal structure of intact rat P450R (Wang M et al., 1997, Proc Nat Acad Sci USA 94:8411–8416). The crystal structure of P450R-FMN reported here confirms the overall similarity of its α-β-α architecture to that of the bacterial flavodoxins, but reveals differences in the position, number, and length of the helices relative to the central β-sheet. The marked similarity between P450R-FMN and flavodoxins in the interactions between the FMN and the protein, indicate a striking evolutionary conservation of the FMN binding site. The P450R-FMN molecule has an unusual surface charge distribution, leading to a very strong dipole, which may be involved in docking cytochrome P450 into place for electron transfer near the FMN. Several acidic residues near the FMN are identified by mutagenesis experiments to be important for electron transfer to P4502D6 and to cytochrome c, a clear indication of the part of the molecular surface that is likely to be involved in substrate binding. Somewhat different parts are found to be involved in binding cytochrome P450 and cytochrome c.

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.

Hydrophobic substitutions at solvent-exposed positions in two α-helical regions of the bacteriophage P22 Arc repressor were introduced by combinatorial mutagenesis. In helix A, hydrophobic residues were tolerated individually at each of the five positions examined, but multiple substitutions were poorly tolerated as shown by the finding that mutants with more than two additional hydrophobic residues were biologically inactive. Several inactive helix A variants were purified and found to have reduced thermal stability relative to wild-type Arc, with a rough correlation between the number of polar-to-hydrophobic substitutions and the magnitude of the stability defect. Quite different results were obtained in helix B, where variants with as many as five polar-to-hydrophobic substitutions were found to be biologically active and one variant with three hydrophobic substitutions had a tm 6 °C higher than wild-type. By contrast, a helix A mutant with three similar polar-to-hydrophobic substitutions was 23 °C less stable than wild-type. Also, one set of three polar-to-hydrophobic substitutions in helix B was tolerated when introduced into the wild-type background but not when introduced into an equally active mutant having a nearly identical structure. Context effects occur both when comparing different regions of the same protein and when comparing the same region in two different homologues.

A statistical analysis of a representative data set of 169 known protein structures was used to analyze the specificity of residue interactions between spatial neighboring strands in β-sheets. Pairwise potentials were derived from the frequency of residue pairs in nearest contact, second nearest and third nearest contacts across neighboring β-strands compared to the expected frequency of residue pairs in a random model. A pseudo-energy function based on these statistical pairwise potentials recognized native β-sheets among possible alternative pairings. The native pairing was found within the three lowest energies in 73% of the cases in the training data set and in 63% of β-sheets in a test data set of 67 proteins, which were not part of the training set. The energy function was also used to detect tripeptides, which occur frequently in β-sheets of native proteins. The majority of native partners of tripeptides were distributed in a low energy range. Self-correcting distance geometry (SECODG) calculations using distance constraints sets derived from possible low energy pairing of β-strands uniquely identified the native pairing of the β-sheet in pancreatic trypsin inhibitor (BPTI). These results will be useful for predicting the structure of proteins from their amino acid sequence as well as for the design of proteins containing β-sheets.

Hsp70 molecular chaperones contain three distinct structural domains, a 44 kDa N-terminal ATPase domain, a 17 kDa peptide-binding domain, and a 10 kDa C-terminal domain. The ATPase and peptide binding domains are conserved in sequence and are functionally well characterized. The function of the 10 kDa variable C-terminal domain is less well understood. We have characterized the secondary structure and dynamics of the C-terminal domain from the Escherichia coli Hsp70, DnaK, in solution by high-resolution NMR. The domain was shown to be comprised of a rigid structure consisting of four helices and a flexible C-terminal subdomain of approximately 33 amino acids. The mobility of the flexible region is maintained in the context of the full-length protein and does not appear to be modulated by the nucleotide state. The flexibility of this region appears to be a conserved feature of Hsp70 architecture and may have important functional implications. We also developed a method to analyze 15N nuclear spin relaxation data, which allows us to extract amide bond vector directions relative to a unique diffusion axis. The extracted angles and rotational correlation times indicate that the helices form an elongated, bundle-like structure in solution.

Haloalkane dehalogenase (DhlA) hydrolyzes short-chain haloalkanes to produce the corresponding alcohols and halide ions. Release of the halide ion from the active-site cavity can proceed via a two-step and a three-step route, which both contain slow enzyme isomerization steps. Thermodynamic analysis of bromide binding and release showed that the slow unimolecular isomerization steps in the three-step bromide export route have considerably larger transition state enthalpies and entropies than those in the other route. This suggests that the three-step route involves different and perhaps larger conformational changes than the two-step export route. We propose that the three-step halide export route starts with conformational changes that result in a more open configuration of the active site from which the halide ion can readily escape. In addition, we suggest that the two-step route for halide release involves the transfer of the halide ion from the halide-binding site in the cavity to a binding site somewhere at the protein surface, where a so-called collision complex is formed in which the halide ion is only weakly bound. No large structural rearrangements are necessary for this latter process.

We examine the similarities and differences between two widely used knowledge-based potentials, which are expressed as contact matrices (consisting of 210 elements) that gives a scale for interaction energies between the naturally occurring amino acid residues. These are the Miyazawa–Jernigan contact interaction matrix M and the potential matrix S derived by Skolnick J et al., 1997, Protein Sci 6:676–688. Although the correlation between the two matrices is good, there is a relatively large dispersion between the elements. We show that when Thr is chosen as a reference solvent within the Miyazawa and Jernigan scheme, the dispersion between the M and S matrices is reduced. The resulting interaction matrix B gives hydrophobicities that are in very good agreement with experiment. The small dispersion between the S and B matrices, which arises due to differing reference states, is shown to have dramatic effect on the predicted native states of lattice models of proteins. These findings and other arguments are used to suggest that for reliable predictions of protein structures, pairwise additive potentials are not sufficient. We also establish that optimized protein sequences can tolerate relatively large random errors in the pair potentials. We conjecture that three body interaction may be needed to predict the folds of proteins in a reliable manner.

A simple approach to estimate the number of α-helical and β-strand segments from protein circular dichroism spectra is described. The α-helix and β-sheet conformations in globular protein structures, assigned by DSSP and STRIDE algorithms, were divided into regular and distorted fractions by considering a certain number of terminal residues in a given α-helix or β-strand segment to be distorted. The resulting secondary structure fractions for 29 reference proteins were used in the analyses of circular dichroism spectra by the SELCON method. From the performance indices of the analyses, we determined that, on an average, four residues per α-helix and two residues per β-strand may be considered distorted in proteins. The number of α-helical and β-strand segments and their average length in a given protein were estimated from the fraction of distorted α-helix and β-strand conformations determined from the analysis of circular dichroism spectra. The statistical test for the reference protein set shows the high reliability of such a classification of protein secondary structure. The method was used to analyze the circular dichroism spectra of four additional proteins and the predicted structural characteristics agree with the crystal structure data.

Unfolded apocytochrome c acquires an α-helical conformation upon interaction with lipid. Folding kinetic results below and above the lipid's CMC, together with energy transfer measurements of lipid bound states, and salt-induced compact states in solution, show that the folding transition of apocytochrome c from the unfolded state in solution to a lipid-inserted helical conformation proceeds via a collapsed intermediate state (IC). This initial compact state is driven by a hydrophobic collapse of the polypeptide chain in the absence of the heme group and may represent a heme-free analogue of an early compact intermediate detected on the folding pathway of cytochrome c in solution. Insertion into the lipid phase occurs via an unfolding step of IC through a more extended state associated with the membrane surface (IS). While IC appears to be as compact as salt-induced compact states in solution with substantial α-helix content, the final lipid-inserted state (Hmic) is as compact as the unfolded state in solution at pH 5 and has an α-helix content which resembles that of native cytochrome c.

Thermotoga maritima (Tm) expresses a 7 kDa monomeric protein whose 18 N-terminal amino acids show 81% identity to N-terminal sequences of cold shock proteins (Csps) from Bacillus caldolyticus and Bacillus stearothermophilus. There were only trace amounts of the protein in Thermotoga cells grown at 80 °C. Therefore, to perform physicochemical experiments, the gene was cloned in Escherichia coli. A DNA probe was produced by PCR from genomic Tm DNA with degenerated primers developed from the known N-terminus of TmCsp and the known C-terminus of CspB from Bacillus subtilis. Southern blot analysis of genomic Tm DNA allowed to produce a partial gene library, which was used as a template for PCRs with gene- and vector-specific primers to identify the complete DNA sequence. As reported for other csp genes, the 5′ untranslated region of the mRNA was anomalously long; it contained the putative Shine–Dalgarno sequence. The coding part of the gene contained 198 bp, i.e., 66 amino acids. The sequence showed 61% identity to CspB from B. caldolyticus and high similarity to all other known Csps. Computer-based homology modeling allowed the conclusion that TmCsp represents a β-barrel similar to CspB from B. subtilis and CspA from E. coli.

As indicated by spectroscopic analysis, analytical gel permeation chromatography, and mass spectrometry, overexpression of the recombinant protein yielded authentic TmCsp with a molecular weight of 7,474 Da. This was in agreement with the results of analytical ultracentrifugation confirming the monomeric state of the protein. The temperature-induced equilibrium transition at 87 °C exceeds the maximum growth temperature of Tm and represents the maximal Tm-value reported for Csps so far.

β-Lactamases are responsible for bacterial resistance to β-lactams and are thus of major clinical importance. However, the identity of the general base involved in their mechanism of action is still unclear. Two candidate residues, Glu166 and Lys73, have been proposed to fulfill this role. Previous studies support the proposal that Glu166 acts during the deacylation, but there is no consensus on the possible role of this residue in the acylation step. Recent experimental data and theoretical considerations indicate that Lys73 is protonated in the free β-lactamases, showing that this residue is unlikely to act as a proton abstractor. On the other hand, it has been proposed that the pKa of Lys73 would be dramatically reduced upon substrate binding and would thus be able to act as a base. To check this hypothesis, we performed continuum electrostatic calculations for five wild-type and three β-lactamase mutants to estimate the pKa of Lys73 in the presence of substrates, both in the Henri–Michaelis complex and in the tetrahedral intermediate. In all cases, the pKa of Lys73 was computed to be above 10, showing that it is unlikely to act as a proton abstractor, even when a β-lactam substrate is bound in the enzyme active site. The pKa of Lys234 is also raised in the tetrahedral intermediate, thus confirming a probable role of this residue in the stabilization of the tetrahedral intermediate. The influence of the β-lactam carboxylate on the pKa values of the active-site lysines is also discussed.

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.

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.

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.

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.

We have identified a ∼140 amino acid domain that is shared by a variety of proteins in budding and fission yeast, nematode, rat, mouse, frog, oat, and man. Typically, this domain is located within 20 residues of the N-terminus of the various proteins. The percent identity among the domains in the 12 proteins ranges from 42 to 93%, with 16 absolutely conserved residues: N-x11–13-V-x2-A-T-x34–36-R-x7–8-W-R-x3-K-x12-G-x-E-x15-L-x11–12-D-x-G-R-x11-D-x7-R. Even though these proteins share little beyond their segment of homology, data are emerging that several of the proteins are involved in endocytosis and or regulation of cytoskeletal organization. We have named this protein segment the ENTH domain, for Epsin N-terminal Homology domain, and hypothesize that it is a candidate for binding specific ligands and/or enzymatic activity in the cell.