Extended retro (reversed) peptide sequences have not previously been accommodated within functional proteins. Here, we show that the entire transmembrane portion of the β-barrel of the pore-forming protein α-hemolysin can be formed by retrosequences comprising a total of 175 amino acid residues, 25 contributed by the central sequence of each subunit of the heptameric pore. The properties of wild-type and retro heptamers in planar bilayers are similar. The single-channel conductance of the retro pore is 15% less than that of the wild-type heptamer and its current-voltage relationship denotes close to ohmic behavior, while the wild-type pore is weakly rectifying. Both wild-type and retro pores are very weakly anion selective. These results and the examination of molecular models suggest that β-barrels may be especially accepting of retro sequences compared to other protein folds. Indeed, the ability to form a retro domain could be diagnostic of a β-barrel, explaining, for example, the activity of the retro forms of many membrane-permeabilizing peptides. By contrast with the wild-type subunits, monomeric retro subunits undergo premature assembly in the absence of membranes, most likely because the altered central sequence fails to interact with the remainder of the subunit, thereby initiating assembly. Despite this difficulty, a technique was devised for obtaining heteromeric pores containing both wild-type and retro subunits. Most probably as a consequence of unfavorable interstrand side-chain interactions, the heteromeric pores are less stable than either the wild-type or retro homoheptamers, as judged by the presence of subconductance states in single-channel recordings. Knowledge about the extraordinary plasticity of the transmembrane β-barrel of α-hemolysin will be very useful in the de novo design of functional membrane proteins based on the β-barrel motif.

Microcalorimetry has been used to measure the stabilities of mutational variants of yeast iso-1 cytochrome c in which F82 and L85 have been replaced by other hydrophobic amino acids. Specifically, F82 has been replaced by Y and L85 by A. The double mutant F82Y, L85A iso-1 has also been studied, and the mutational perturbations are compared to those for the two single mutants, F82Y iso-1 and L85A iso-1. Results are interpreted in terms of known crystallographic structures. The data show that (1) the destabilization of the mutant proteins is similar in magnitude to that which is theoretically predicted by the more obvious mutation-induced structural effects; (2) the free energy of destabilization of the double mutant, F82Y, L85A iso-1, is less than the sum of those of the two single mutants, almost certainly because, in the double mutant, the -OH group of Y82 is able to protrude into the cavity formed by the L85A substitution. The more favorable structural accommodation of the new -OH group in the double mutant leads to additional stability through (1) further decreases in the volumes of internal cavities and (2) formation of an extra protein–protein hydrogen bond.

Production of folded and biologically active protein from Escherichia coli derived inclusion bodies can only be accomplished if a scheme exists for in vitro naturation. Motivated by the need for a rapid and statistically meaningful method of determining and evaluating protein folding conditions, we have designed a new fractional factorial protein folding screen. The screen includes 12 factors shown by previous experiments to enhance protein folding and it incorporates the 12 factors into 16 different folding conditions. By examining a 1/256th fraction of the full factorial, multiple folding conditions were determined for the ligand binding domains from glutamate and kainate receptors, and for lysozyme and carbonic anhydrase B. The impact of each factor on the formation of biologically active material was estimated by calculating factor main effects. Factors and corresponding levels such as pH (8.5) and l-arginine (0.5 M) consistently had a positive effect on protein folding, whereas detergent (0.3 mM lauryl maltoside) and nonpolar additive (0.4 M sucrose) were detrimental to the folding of these four proteins. One of the 16 conditions yielded the most folded material for three out of the four proteins. Our results suggest that this protein folding screen will be generally useful in determining whether other proteins will fold in vitro and, if so, what factors are important. Furthermore, fractional factorial folding screens are well suited to the evaluation of previously untested factors on protein folding.

The relative orientations of adjacent structural elements without many well-defined NOE contacts between them are typically poorly defined in NMR structures. For apo-S100B(ββ) and the structurally homologous protein calcyclin, the solution structures determined by conventional NMR exhibited considerable differences and made it impossible to draw unambiguous conclusions regarding the Ca2+-induced conformational change required for target protein binding. The structure of rat apo-S100B(ββ) was recalculated using a large number of constraints derived from dipolar couplings that were measured in a dilute liquid crystalline phase. The dipolar couplings orient bond vectors relative to a single-axis system, and thereby remove much of the uncertainty in NOE-based structures. The structure of apo-S100B(ββ) indicates a minimal change in the first, pseudo-EF-hand Ca2+ binding site, but a large reorientation of helix 3 in the second, classical EF-hand upon Ca2+ binding.

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.

The X-ray crystal structure of the human α-thrombin-hirunorm IV complex has been determined at 2.5 Å resolution, and refined to an R-factor of 0.173. The structure reveals an inhibitor binding mode distinctive of a true hirudin mimetic, which justifies the high inhibitory potency and the selectivity of hirunorm IV. This novel inhibitor, composed of 26 amino acids, interacts through the N-terminal end with the α-thrombin active site in a nonsubstrate mode, and binds specifically to the fibrinogen recognition exosite through the C-terminal end. The backbone of the N-terminal tripeptide Chg1"-Arg2"-2Nal3" (Chg, cyclohexyl-glycine; 2Nal, β-(2-naphthyl)-alanine) forms a parallel β-strand to the thrombin main-chain segment Ser214–Gly216. The Chg1" side chain occupies the S2 site, Arg2" penetrates into the S1 specificity site, while the 2Nal3" side chain occupies the aryl binding site. The Arg2" side chain enters the S1 specificity pocket from a position quite apart from the canonical P1 site. This notwithstanding, the Arg2" side chain establishes the typical ion pair with the carboxylate group of Asp189.

Despite decades of intense study, the complementarity of β-lactams for β-lactamases and penicillin binding proteins is poorly understood. For most of these enzymes, β-lactam binding involves rapid formation of a covalent intermediate. This makes measuring the equilibrium between bound and free β-lactam difficult, effectively precluding measurement of the interaction energy between the ligand and the enzyme. Here, we explore the energetic complementarity of β-lactams for the β-lactamase AmpC through reversible denaturation of adducts of the enzyme with β-lactams. AmpC from Escherichia coli was reversibly denatured by temperature in a two-state manner with a temperature of melting (Tm) of 54.6 °C and a van't Hoff enthalpy of unfolding (ΔHVH) of 182 kcal/mol. Solvent denaturation gave a Gibbs free energy of unfolding in the absence of denaturant (ΔGuH2O) of 14.0 kcal/mol. Ligand binding perturbed the stability of the enzyme. The penicillin cloxacillin stabilized AmpC by 3.2 kcal/mol (ΔTm = +5.8 °C); the monobactam aztreonam stabilized the enzyme by 2.7 kcal/mol (ΔTm = +4.9 °C). Both acylating inhibitors complement the active site. Surprisingly, the oxacephem moxalactam and the carbapenem imipenem both destabilized AmpC, by 1.8 kcal/mol (ΔTm = −3.2 °C) and 0.7 kcal/mol (ΔTm = −1.2 °C), respectively. These β-lactams, which share nonhydrogen substituents in the 6(7)α position of the β-lactam ring, make unfavorable noncovalent interactions with the enzyme. Complexes of AmpC with transition state analog inhibitors were also reversibly denatured; both benzo(b)thiophene-2-boronic acid (BZBTH2B) and p-nitrophenyl phenylphosphonate (PNPP) stabilized AmpC. Finally, a catalytically inactive mutant of AmpC, Y150F, was reversibly denatured. It was 0.7 kcal/mol (ΔTm = −1.3 °C) less stable than wild-type (WT) by thermal denaturation. Both the cloxacillin and the moxalactam adducts with Y150F were significantly destabilized relative to their WT counterparts, suggesting that this residue plays a role in recognizing the acylated intermediate of the β-lactamase reaction. Reversible denaturation allows for energetic analyses of the complementarity of AmpC for β-lactams, through ligand binding, and for itself, through residue substitution. Reversible denaturation may be a useful way to study ligand complementarity to other β-lactam binding proteins as well.

Short chain l-3-hydroxyacyl CoA dehydrogenase (SCHAD) is a soluble dimeric enzyme critical for oxidative metabolism of fatty acids. Its primary sequence has been reported to be conserved across numerous tissues and species with the notable exception of the pig heart homologue. Preliminary efforts to solve the crystal structure of the dimeric pig heart SCHAD suggested the unprecedented occurrence of three enzyme subunits within the asymmetric unit, a phenomenon that was thought to have hampered refinement of the initial chain tracing. The recently solved crystal coordinates of human heart SCHAD facilitated a molecular replacement solution to the pig heart SCHAD data. Refinement of the model, in conjunction with the nucleotide sequence for pig heart SCHAD determined in this paper, has demonstrated that the previously published pig heart SCHAD sequence was incorrect. Presented here are the corrected amino acid sequence and the high resolution crystal structure determined for pig heart SCHAD complexed with its NAD+ cofactor (2.8 Å; Rcryst = 22.4%, Rfree = 28.8%). In addition, the peculiar phenomenon of a dimeric enzyme crystallizing with three subunits contained in the asymmetric unit is described.

β-Lactamases are the major resistance mechanism to β-lactam antibiotics and pose a growing threat to public health. Recently, bacteria have become resistant to β-lactamase inhibitors, making this problem pressing. In an effort to overcome this resistance, non-β-lactam inhibitors of β-lactamases were investigated for complementarity to the structure of AmpC β-lactamase from Escherichia coli. This led to the discovery of an inhibitor, benzo(b)thiophene-2-boronic acid (BZBTH2B), which inhibited AmpC with a Ki of 27 nM. This inhibitor is chemically dissimilar to β-lactams, raising the question of what specific interactions are responsible for its activity. To answer this question, the X-ray crystallographic structure of BZBTH2B in complex with AmpC was determined to 2.25 Å resolution. The structure reveals several unexpected interactions. The inhibitor appears to complement the conserved, R1-amide binding region of AmpC, despite lacking an amide group. Interactions between one of the boronic acid oxygen atoms, Tyr150, and an ordered water molecule suggest a mechanism for acid/base catalysis and a direction for hydrolytic attack in the enzyme catalyzed reaction. To investigate how a non-β-lactam inhibitor would perform against resistant bacteria, BZBTH2B was tested in antimicrobial assays. BZBTH2B significantly potentiated the activity of a third-generation cephalosporin against AmpC-producing resistant bacteria. This inhibitor was unaffected by two common resistance mechanisms that often arise against β-lactams in conjunction with β-lactamases. Porin channel mutations did not decrease the efficacy of BZBTH2B against cells expressing AmpC. Also, this inhibitor did not induce expression of AmpC, a problem with many β-lactams. The structure of the BZBTH2B/AmpC complex provides a starting point for the structure-based elaboration of this class of non-β-lactam inhibitors.

To explore the structure of the pore-forming fragment of colicin E1 in membranes, a series of 23 consecutive single cysteine substitution mutants was prepared in the sequence 402–424. Each mutant was reacted with a sulfhydryl-specific reagent to generate a nitroxide labeled side chain, and the mobility of the side chain and its accessibility to collision with paramagnetic reagents was determined from the electron paramagnetic resonance spectrum. Individual values of these quantities were used to identify tertiary contact sites and the nature of the surrounding solvent, while their periodic dependence on sequence position was used to identify secondary structure. In solution, the data revealed a regular helix of 11 residues in the region 406–416, consistent with helix IV of the crystal structure. Upon binding to negatively charged membranes at pH 4.0, helix IV apparently grows to a length of 19 residues, extending from 402–420. One face of the helix is solvated by the lipid bilayer, and the other by an environment of a polar nature. Surprisingly, a conserved charged pair, D408–R409, is located on the lipid-exposed face. Evidence is presented to suggest a transmembrane orientation of this new helix, although other topographies may exist in equilibrium.

The phage 434 Cro protein, the N-terminal domain of its repressor (R1–69) and that of phage λ (λ6–85) constitute a group of small, monomeric, single-domain folding units consisting of five helices with striking structural similarity. The intrinsic helix stabilities in λ6–85 have been correlated to its rapid folding behavior, and a residual hydrophobic cluster found in R1–69 in 7 M urea has been proposed as a folding initiation site. To understand the early events in the folding of 434 Cro, and for comparison with R1–69 and λ6–85, we examined the conformational behavior of five peptides covering the entire 434 Cro sequence in water, 40% (by volume) TFE/water, and 7 M urea solutions using CD and NMR. Each peptide corresponds to a helix and adjacent residues as identified in the native 434 Cro NMR and crystal structures. All are soluble and monomeric in the solution conditions examined except for the peptide corresponding to the 434 Cro helix 4, which has low water solubility. Helix formation is observed for the 434 Cro helix 1 and helix 2 peptides in water, for all the peptides in 40% TFE and for none in 7 M urea. NMR data indicate that the helix limits in the peptides are similar to those in the native protein helices. The number of side-chain NOEs in water and TFE correlates with the helix content, and essentially none are observed in 7 M urea for any peptide, except that for helix 5, where a hydrophobic cluster may be present. The low intrinsic folding propensities of the five helices could account for the observed stability and folding behavior of 434 Cro and is, at least qualitatively, in accord with the results of the recently described diffusion-collision model incorporating intrinsic helix propensities.

A neural network has been used to predict both the location and the type of β-turns in a set of 300 nonhomologous protein domains. A substantial improvement in prediction accuracy compared with previous methods has been achieved by incorporating secondary structure information in the input data. The total percentage of residues correctly classified as β-turn or not-β-turn is around 75% with predicted secondary structure information. More significantly, the method gives a Matthews correlation coefficient (MCC) of around 0.35, compared with a typical MCC of around 0.20 using other β-turn prediction methods. Our method also distinguishes the two most numerous and well-defined types of β-turn, types I and II, with a significant level of accuracy (MCCs 0.22 and 0.26, respectively).

In addition to the Cys-Xaa-Xaa-Cys motif at position 30–33, DsbA, the essential catalyst for disulfide bond formation in the bacterial periplasm shares with other oxidoreductases of the thioredoxin family a cis-proline in proximity of the active site residues. In the variant DsbAP151A, this residue has been changed to an alanine, an almost isosteric residue which is not disposed to adopt the cis conformation. The substitution strongly destabilized the structure of DsbA, as determined by the decrease in the free energy of folding. The pKa of the thiol of Cys30 was only marginally decreased. Although in vivo the variant appeared to be correctly oxidized, it exhibited an activity less than half that of the wild-type enzyme with respect to the folding of alkaline phosphatase, used as a reporter of the disulfide bond formation in the periplasm. DsbAP151A crystallized in a different crystal form from the wild-type protein, in space group P21 with six molecules in the asymmetric unit. Its X-ray structure was determined to 2.8 Å resolution. The most significant conformational changes occurred at the active site. The loop 149–152 adopted a new backbone conformation with Ala151 in a trans conformation. This rearrangement resulted in the loss of van der Waals interactions between this loop and the disulfide bond. His32 from the Cys-Xaa-Xaa-Cys sequence presented in four out of six molecules in the asymmetric unit a gauche− conformation not observed in the wild-type protein. The X-ray structure and folding studies on DsbAP151A were consistent with the cis-proline playing a major role in the stabilization of the protein. A role for the positioning of the substrate is discussed. These important properties for the enzyme function might explain the conservation of this residue in DsbA and related proteins possessing the thioredoxin fold.

In this paper, elevated pressures up to 750 atm (1 atm = 101 kPa) were found to have a strong stabilizing effect on two extremely thermophilic glutamate dehydrogenases (GDHs): the native enzyme from the hyperthermophile Pyrococcus furiosus (Pf), and a recombinant GDH mutant containing an extra tetrapeptide at the C-terminus (rGDHt). The presence of the tetrapeptide greatly destabilized the recombinant mutant at ambient pressure; however, the destabilizing effect was largely reversed by the application of pressure. Electron spin resonance (ESR) spectroscopy of a spin-label attached to the terminal cysteine of rGDHt revealed a high degree of mobility, suggesting that destabilization is due to weakened intersubunit ion-pair interactions induced by thermal fluctuations of the tetrapeptide. For both enzymes, the stabilizing effect of pressure increased with temperature as well as pressure, reaching 36-fold for rGDHt at 105 °C and 750 atm, the largest pressure-induced thermostabilization of an enzyme reported to date. Stabilization of both native GDH and rGDHt was also achieved by adding glycerol. Based on the kinetics of thermal inactivation and the known effects of glycerol on protein structure, a mechanism of pressure-induced thermostabilization is proposed.

Maltose binding protein (MBP) is a large, monomeric two domain protein containing 370 amino acids. In the absence of denaturant at neutral pH, the protein is in the native state, while at pH 3.0 it forms a molten globule. The molten globule lacks a tertiary circular dichroism signal but has secondary structure similar to that of the native state. The molten globule binds 8-anilino-1-naphthalene sulfonate (ANS). The unfolding thermodynamics of MBP at both pHs were measured by carrying out a series of isothermal urea melts at temperatures ranging from 274–329 K. At 298 K, values of ΔG°, ΔCp, and Cm were 3.1 ± 0.2 kcal mol−1, 5.9 ± 0.8 kcal mol−1 K−1 (15.9 cal (mol-residue)−1 K−1), and 0.8 M, respectively, at pH 3.0 and 14.5 ± 0.4 kcal mol−1, 8.3 ± 0.7 kcal mol−1 K−1 (22.4 kcal (mol-residue)−1 K−1), and 3.3 M, respectively, at pH 7.1. Guanidine hydrochloride denaturation at pH 7.1 gave values of ΔG° and ΔCp similar to those obtained with urea. The m values for denaturation are strongly temperature dependent, in contrast to what has been previously observed for small globular proteins. The value of ΔCp per mol-residue for the molten globule is comparable to corresponding values of ΔCp for the unfolding of typical globular proteins and suggests that it is a highly ordered structure, unlike molten globules of many small proteins. The value of ΔCp per mol-residue for the unfolding of the native state is among the highest currently known for any protein.

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.

Endo-β-N-acetylglucosaminidase H hydrolyzes the β-(1–4)-glycosidic link of the N,N′-diacetylchitobiose core of high-mannose and hybrid asparagine-linked oligosaccharides. Seven mutants of the active site residues, Asp130 and Glu132, have been prepared, assayed, and crystallized. They include single site mutants of each residue to the corresponding amide, to Ala and to the alternate acidic residue, and to the double amide mutant. The mutants of Asp130 are more active than the corresponding Glu132 mutants, consistent with the assignment of the latter residue as the primary catalytic residue. The amide mutants are more active than the alternate acidic residue mutants, which in turn are more active than the Ala mutants. The structures of the Asn mutant of Asp130 and the double mutant are very similar to that of the wild-type enzyme. Several residues surrounding the mutated residues, including some that form part of the core of the β-barrel and especially Tyr168 and Tyr244, adopt a very different conformation in the structures of the other two mutants of Asp130 and in the Asp mutant of Glu132. The results show that the residues in the upper layers of the β-barrel can organize into two very distinct packing arrangements that depend on subtle electrostatic and steric differences and that greatly affect the geometry of the substrate-binding cleft. Consequently, the relative activities of several of the mutants are defined by structural changes, leading to impaired substrate binding, in addition to changes in functionality.

Apomyoglobin from sperm whale is often used for studies of ligand binding, protein folding, and protein stability. In an effort to describe its conformational properties in solution, homonuclear and heteronuclear (13C and 15N) NMR methods were applied to the protein in its native state. Assignments were confirmed for nuclear Overhauser effects (NOEs) involving side chain and backbone protons in the folded regions of the structure. These NOEs were used to derive distance restraints. The shifts induced by the hydrophobic dye 8-anilino-1-naphthalenesulfonic acid (ANS) were inspected in the regions remote from its binding site and served as an indicator of conformational flexibility. 3JαH-NH values were obtained to assess dihedral angle averaging and to provide additional restraints. A family of structures was calculated with X-PLOR and an ab initio simulated annealing protocol using holomyoglobin as a template. Where the structure appeared well defined by chemical shift, line width, ANS perturbation, and density of NOEs, the low resolution model of apomyoglobin provides a valid approximation for the structure. The new model offers an improved representation of the folded regions of the protein, which encompass the A, B, E, helices as well as parts of the G and H helices. Regions that are less well defined at this stage of calculations include the CD corner and the end of the H-helix. The EF-F-FG segment remains uncharacterized.

An α-helix and a β-strand are said to be interactively packed if at least one residue in each of the secondary structural elements loses 10% of its solvent accessible contact area on association with the other secondary structural element. An analysis of all such 5,975 nonidentical α/β units in protein structures, defined at ≤2.5 Å resolution, shows that the interaxial distance between the α-helix and the β-strand is linearly correlated with the residue-dependent function, log[(V/nda)/n-int], where V is the volume of amino acid residues in the packing interface, nda is the normalized difference in solvent accessible contact area of the residues in packed and unpacked secondary structural elements, and n-int is the number of residues in the packing interface. The β-sheet unit (βu), defined as a pair of adjacent parallel or antiparallel hydrogen-bonded β-strands, packing with an α-helix shows a better correlation between the interaxial distance and log(V/nda) for the residues in the packing interface. This packing relationship is shown to be useful in the prediction of interaxial distances in α/β units using the interacting residue information of equivalent α/β units of homologous proteins. It is, therefore, of value in comparative modeling of protein structures.

The crystal structure of the S642A mutant of mitochondrial aconitase (mAc) with citrate bound has been determined at 1.8 Å resolution and 100 K to capture this binding mode of substrates to the native enzyme. The 2.0 Å resolution, 100 K crystal structure of the S642A mutant with isocitrate binding provides a control, showing that the Ser → Ala replacement does not alter the binding of substrates in the active site. The aconitase mechanism requires that the intermediate product, cis-aconitate, flip over by 180° about the Cα–Cβ double bond. Only one of these two alternative modes of binding, that of the isocitrate mode, has been previously visualized. Now, however, the structure revealing the citrate mode of binding provides direct support for the proposed enzyme mechanism.