The change in heat capacity ΔCp for the folding of ribonuclease A was determined using differential scanning calorimetry and thermal denaturation curves. The methods gave equivalent results, ΔCp = 1.15 ± 0.08 kcal mol−1 K−1. Estimates of the conformational stability of ribonuclease A based on these results from thermal unfolding are in good agreement with estimates from urea unfolding analyzed using the linear extrapolation method.

The monomer–dimer equilibrium for the human immunodeficiency virus type 1 (HIV-1) protease has been investigated under physiological conditions. Dimer dissociation at pH 7.0 was correlated with a loss in β-sheet structure and a lower degree of ANS binding. An autolysis-resistant mutant, Q7K/L33I/L63I, was used to facilitate sedimentation equilibrium studies at neutral pH where the wild-type enzyme is typically unstable in the absence of bound inhibitor. The dimer dissociation constant (KD) of the triple mutant was 5.8 μM at pH 7.0 and was below the limit of measurement (∼100 nM) at pH 4.5. Similar studies using the catalytically inactive D25N mutant yielded a KD value of 1.0 μM at pH 7.0. These values differ significantly from a previously reported value of 23 nM obtained indirectly from inhibitor binding measurements (Darke et al., 1994). We show that the discrepancy may result from the thermodynamic linkage between the monomer–dimer and inhibitor binding equilibria. Under conditions where a significant degree of monomer is present, both substrates and competitive inhibitors will shift the equilibrium toward the dimer, resulting in apparent increases in dimer stability and decreases in ligand binding affinity. Sedimentation equilibrium studies were also carried out on several drug-resistant HIV-1 protease mutants: V82F, V82F/I84V, V82T/I84V, and L90M. All four mutants exhibited reduced dimer stability relative to the autolysis-resistant mutant at pH 7.0. Our results indicate that reductions in drug affinity may be due to the combined effects of mutations on both dimer stability and inhibitor binding.

The effect of methylurea, N,N′-dimethylurea, ethylurea, and butylurea as well as guanidine hydrochloride (GuHCl), urea and pH on the thermal stability, structural properties, and preferential solvation changes accompanying the thermal unfolding of ribonuclease A (RNase A) has been investigated by differential scanning calorimetry (DSC), UV, and circular dichroism (CD) spectroscopy. The results show that the thermal stability of RNase A decreases with increasing concentration of denaturants and the size of the hydrophobic group substituted on the urea molecule. From CD measurements in the near- and far-UV range, it has been observed that the tertiary structure of RNase A melts at about 3 °C lower temperature than its secondary structure, which means that the hierarchy in structural building blocks exists for RNase A even at conditions at which according to DSC and UV measurements the RNase A unfolding can be interpreted in terms of a two-state approximation. The far-UV CD spectra also show that the final denatured states of RNase A at high temperatures in the presence of different denaturants including 4.5 M GuHCl are similar to each other but different from the one obtained in 4.5 M GuHCl at 25 °C. The concentration dependence of the preferential solvation change ΔΓ23, expressed as the number of cosolvent molecules entering or leaving the solvation shell of the protein upon denaturation and calculated from DSC data, shows the same relative denaturation efficiency of alkylureas as other methods.

We have made a comparative structure based analysis of the thermodynamics of lectin-carbohydrate (L-C) binding and protein folding. Examination of the total change in accessible surface area in those processes revealed a much larger decrease in free energy per unit of area buried in the case of L-C associations. According to our analysis, this larger stabilization of L-C interactions arises from a more favorable enthalpy of burying a unit of polar surface area, and from higher proportions of polar areas. Hydrogen bonds present at 14 L-C interfaces were identified, and their overall characteristics were compared to those reported before for hydrogen bonds in protein structures. Three major factors might explain why polar–polar interactions are stronger in L-C binding than in protein folding: (1) higher surface density of hydrogen bonds; (2) better hydrogen-bonding geometry; (3) larger proportion of hydrogen bonds involving charged groups. Theoretically, the binding entropy can be partitioned into three main contributions: entropy changes due to surface desolvation, entropy losses arising from freezing rotatable bonds, and entropic effects that result from restricting translation and overall rotation motions. These contributions were estimated from structural information and added up to give calculated binding entropies. Good correlation between experimental and calculated values was observed when solvation effects were treated according to a parametrization developed by other authors from protein folding studies. Finally, our structural parametrization gave calculated free energies that deviate from experimental values by 1.1 kcal/mol on the average; this amounts to an uncertainty of one order of magnitude in the binding constant.

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.

Calmodulin (CaM) is a 148-residue regulatory calcium-binding protein that activates a wide range of target proteins and enzymes. Calcium-saturated CaM has a bilobal structure, and each domain has an exposed hydrophobic surface region where target proteins are bound. These two “active sites” of calmodulin are remarkably rich in Met residues. Here we have biosynthetically substituted (up to 90% incorporation) the unnatural amino acids ethionine (Eth) and norleucine (Nle) for the nine Met residues of CaM. The substituted proteins bind in a calcium-dependent manner to hydrophobic matrices and a synthetic peptide, encompassing the CaM-binding domain of myosin light-chain kinase (MLCK). Infrared and circular dichroism spectroscopy show that there are essentially no changes in the secondary structure of these proteins compared to wild-type CaM (WT-CaM). One- and two-dimensional NMR studies of the Eth-CaM and Nle-CaM proteins reveal that, while the core of the proteins is relatively unaffected by the substitutions, the two hydrophobic interaction surfaces adjust to accommodate the Eth and Nle residues. Enzyme activation studies with MLCK show that Eth-CaM and Nle-CaM activate the enzyme to 90% of its maximal activity, with little changes in dissociation constant. For calcineurin only 50% activation was obtained, and the KD for Nle-CaM also increased 3.5-fold compared with WT-CaM. These data show that the “active site” Met residues of CaM play a distinct role in the activation of different target enzymes, in agreement with site-directed mutagenesis studies of the Met residues of CaM.

Lumazine synthase catalyzes the penultimate step in the synthesis of riboflavin in plants, fungi, and microorganisms. The enzyme displays two quaternary structures, the pentameric forms in yeast and fungi and the 60-meric icosahedral capsids in plants and bacteria. To elucidate the structural features that might be responsible for differences in assembly, we have determined the crystal structures of lumazine synthase, complexed with the inhibitor 5-nitroso-6-ribitylamino-2,4-pyrimidinedione, from spinach and the fungus Magnaporthe grisea to 3.3 and 3.1 Å resolution, respectively. The overall structure of the subunit and the mode of inhibitor binding are very similar in these enzyme species. The core of the subunit consists of a four-stranded parallel β-sheet sandwiched between two helices on one side and three helices on the other. The packing of the five subunits in the pentameric M. grisea lumazine synthase is very similar to the packing in the pentameric substructures in the icosahedral capsid of the plant enzyme. Two structural features can be correlated to the differences in assembly. In the plant enzyme, the N-terminal β-strand interacts with the β-sheet of the adjacent subunit, thus extending the sheet from four to five strands. In fungal lumazine synthase, an insertion of two residues after strand β1 results in a completely different orientation of this part of the polypeptide chain and this conformational difference prevents proper packing of the subunits at the trimer interface in the icosahedron. In the spinach enzyme, the β-hairpin connecting helices α4 and α5 participates in the packing at the trimer interface of the icosahedron. Another insertion of two residues at this position of the polypeptide chain in the fungal enzyme disrupts the hydrogen bonding in the hairpin, and the resulting change in conformation of this loop also interferes with proper intrasubunit contacts at the trimer interface.

Guanylin is a guanylyl cyclase (GC)-activating peptide that is mainly secreted as the corresponding prohormone of 94 amino acid residues. In this study, we show that the originally isolated 15-residue guanylin, representing the COOH-terminal part of the prohormone, is released from the prohormone by cleavage of an Asp-Pro amide bond under conditions applied during the isolation procedures. Thus, the 15-residue guanylin is probably a non-native, chemically induced GC-activating peptide. This guanylin molecule contains two disulfide bonds that are absolutely necessary for receptor activation. We demonstrate that the folding of the reduced 15-residue guanylin results almost completely in the formation of the two inactive disulfide isomers. In contrast, the reduced form of proguanylin containing the entire prosequence folds to a product with the native cysteine connectivity. Because proguanylin lacking the 31 NH2-terminal residues of the prosequence folds only to a minor extent to guanylin with the native disulfide bonds, it is evident that this NH2-terminal region contributes significantly to the correct disulfide-coupled folding. Structural studies using CD and NMR spectroscopy show that native proguanylin contains a considerable amount of α-helical and, to a lesser extent, β-sheet structural elements. In addition, a close proximity of the NH2- and the COOH-terminal regions was found by NOESY. It appears that this interaction is important for the constitution of the correct conformation and provides an explanation of the minor guanylyl cyclase activity of proguanylin by shielding the bioactive COOH-terminal domain from the receptor.

Homologs of the tumor suppressor p53, called p63 and p73, have been identified. The p63 and p73 family members possess a domain structure similar to p53, but contain variable C-terminal extensions. We find that some of the C-terminal extensions contain Sterile Alpha Motif (SAM) domains. SAM domains are protein modules that are involved in protein–protein interactions. Consistent with this role, the C-terminal SAM domains of the p63 and p73 may regulate function by recruiting other protein effectors.

Catalysis by thioredoxin reductase (TrxR) from Escherichia coli requires alternation between two domain arrangements. One of these conformations has been observed by X-ray crystallography (Waksman G, Krishna TSR, Williams CH Jr, Kuriyan J, 1994, J Mol Biol 236:800–816). This form of TrxR, denoted FO, permits the reaction of enzyme-bound reduced FAD with a redox-active disulfide on TrxR. As part of an investigation of conformational changes and intermediates in catalysis by TrxR, an X-ray structure of the FO form of TrxR with both the FAD and active site disulfide reduced has been determined. Reduction after crystallization resulted in significant local conformation changes. The isoalloxazine ring of the FAD cofactor, which is essentially planar in the oxidized enzyme, assumes a 34° “butterfly” bend about the N(5)–N(10) axis in reduced TrxR. Theoretical calculations reported by others predict ring bending of 15–28° for reduced isoalloxazines protonated at N(1). The large bending in reduced TrxR is attributed in part to steric interactions between the isoalloxazine ring and the sulfur of Cys138, formed by reduction of the active site disulfide, and is accompanied by changes in the positions and interactions of several of the ribityl side-chain atoms of FAD. The bending angle in reduced TrxR is larger than that for any flavoprotein in the Protein Data Bank. Distributions of bending angles in published oxidized and reduced flavoenzyme structures are different from those found in studies of free flavins, indicating that the protein environment has a significant effect on bending.

Evidence already available is used to demonstrate that although prostaglandin G/H synthase hydroxylates arachidonic acid through radical intermediates, it effects cyclizations through a carbocation center at C-10. This is produced following migration of H to the initial radical at C-13 and a 1ε oxidation. Under orbital symmetry control, the cyclizations can give only the ring size and trans stereochemistry actually observed. After cyclization, the H-shift reverses to take the sequence back into current radical theory for hydroxylation at C-15. Thus 10,10-difluoroarachidonic acid cannot be cyclized, although it can be hydroxylated. Acetylation of Ser516 in the isoform synthase-2 is considered to oppose carbocation formation and/or H-migration and so prevent cyclizations while permitting hydroxylations; the associated inversion of chirality at C-15 can then readily be accommodated without the change in conformation required by other schemes. Suicide inhibition occurs when carbocations form stable bonds upon (thermal) contact with adjacent heteroatoms, etc. Because the cyclooxygenase and peroxidase functions operate simultaneously through the same heme, phenol acts as reducing cosubstrate for the cyclooxygenase, thus enabling it to promote PGG2 production and protect the enzyme from oxidative destruction.

Eotaxin is a member of the chemokine family of about 40 proteins that induce cell migration. Eotaxin binds the CC chemokine receptor CCR3 that is highly expressed by eosinophils, and it is considered important in the pathology of chronic respiratory disorders such as asthma. The high resolution structure of eotaxin is known. The 74 amino acid protein has two disulfide bridges and shows a typical chemokine fold comprised of a core of three antiparallel β-strands and an overlying α-helix. In this paper, we report the backbone dynamics of eotaxin determined through 15N-T1, T2, and {1H}-15N nuclear Overhauser effect heteronuclear multidimensional NMR experiments. This is the first extensive study of the dynamics of a chemokine derived from 600, 500, and 300 MHz NMR field strengths. From the T1, T2, and NOE relaxation data, parameters that describe the internal motions of eotaxin were derived using the Lipari–Szabo model free analysis. The most ordered regions of the protein correspond to the known secondary structure elements. However, surrounding the core, the regions known to be functionally important in chemokines show a range of motions on varying timescales. These include extensive subnanosecond to picosecond motions in the N-terminus, C-terminus, and the N-loop succeeding the disulfides. Analysis of rotational diffusion anisotropy of eotaxin and chemical exchange terms at multiple fields also allowed the confident identification of slow conformational exchange through the “30s” loop, disulfides, and adjacent residues. In addition, we show that these motions may be attenuated in the dimeric form of a synthetic eotaxin. The structure and dynamical basis for eotaxin receptor binding is discussed in light of the dynamics data.

The objective of this study was to evaluate the suitability of the WW domain as a desirable model system to understand the folding and stability of an isolated three-stranded antiparallel β-sheet structure. The WW domain was subjected to thermal and chaotropic denaturation/reconstitution utilizing a variety of biophysical methods. This three-stranded sheet folds reversibly and cooperatively utilizing both urea and GdnHCl as denaturants; however, the denatured state retains structure in the form of a hydrophobic cluster involving at least one aromatic side chain. In contrast to chaotropic denaturation, thermal denaturation appears to be more complete and may be a two state process. The suitability of the WW domain for future studies aimed at understanding the kinetics and thermodynamics of antiparallel β-sheet folding clearly emerges from this initial study. The most exciting and significant result in this manuscript is the finding that the chaotropic denatured state of WW has a hydrophobic cluster as discerned by near-UV CD evidence. The role that the denatured state plays in the folding and stability of a three-stranded β-sheets, and its capacity for preventing aggregation may be particularly important and is the subject of ongoing studies.

β-Galactosidase (lacZ) from Escherichia coli is a 464 kDa homotetramer. Each subunit consists of five domains, the third being an α/β barrel that contains most of the active site residues. A comparison is made between each of the domains and a large set of proteins representative of all structures from the protein data bank. Many structures include an α/β barrel. Those that are most similar to the α/β barrel of E. coli β-galactosidase have similar catalytic residues and belong to the so-called “4/7 superfamily” of glycosyl hydrolases. The structure comparison suggests that β-amylase should also be included in this family. Of three structure comparison methods tested, the “ProSup” procedure of Zu-Kang and Sippl and the “Superimpose” procedure of Diederichs were slightly superior in discriminating the members of this superfamily, although all procedures were very powerful in identifying related protein structures. Domains 1, 2, and 4 of E. coli β-galactosidase have topologies related to “jelly-roll barrels” and “immunoglobulin constant” domains. This fold also occurs in the cellulose binding domains (CBDs) of a number of glycosyl hydrolases. The fold of domain 1 of E. coli β-galactosidase is closely related to some CBDs, and the domain contributes to substrate binding, but in a manner unrelated to cellulose binding by the CBDs. This is typical of domains 1, 2, 4, and 5, which appear to have been recruited to play roles in β-galactosidase that are unrelated to the functions that such domains provide in other contexts. It is proposed that β-galactosidase arose from a prototypical single domain α/β barrel with an extended active site cleft. The subsequent incorporation of elements from other domains could then have reduced the size of the active site from a cleft to a pocket to better hydrolyze the disaccharide lactose and, at the same time, to facilitate the production of inducer, allolactose.

Refolding of b*C40A/C82A/P27A is comprised of several kinetically detectable folding phases. The slowest phase in refolding originates from trans → cis isomerization of the Tyr47–Pro48 peptide bond being in cis conformation in the native state. This refolding phase can be accelerated by the peptidyl-prolyl cis/trans isomerase human cytosolic cyclophilin (Cyp18) with a kcat/KM of 254,000 M−1 s−1. The fast refolding phase is not influenced by the enzyme.

This work experimentally confirms the pathway of activation of Ha-ras-p21, which was calculated by the method of Targeted Molecular Dynamics (TMD) (Díaz JF, Wroblowski B, Schlitter J, Engelborghs Y, 1997a, Proteins Struct Funct Genet 28:434–451). The process can be studied experimentally by analyzing the binding of BeF3− to the GDP complex of the active fluorescent mutant Y32W (Díaz JF, Sillen A, Engelborghs Y, 1997b, J Biol Chem 227:23138–23143). Two mutants, V29G and I36G, have been constructed at both sides of the effector loop of the active fluorescent mutant. This was done to check the proposed reaction pathway and to provide further insight into the mechanism of the activation of ras proteins. Both mutations accelerate the conformational isomerization with two orders of magnitude, demonstrating convincingly the role of these residues as hinges of the effector loop in one or more of the transitions of the conformational change. These results provide experimental support to the pathway calculated by TMD analysis.

The monoclonal antibody 1696, directed against the HIV-1 protease, displays strong inhibitory effects toward the catalytic activity of the enzyme of both the HIV-1 and HIV-2 isolates. This antibody cross-reacts with peptides that include the N-terminus of the enzyme, a region that is well conserved in sequence among different viral strains and which, furthermore, is crucial for homodimerization to the active enzymatic form. This observation, as well as antigen-binding studies in the presence of an active site inhibitor, suggest that 1696 inhibits the HIV protease by destabilizing its active homodimeric form. To characterize further how the antibody 1696 inhibits the HIV-1 and HIV-2 proteases, we have solved the crystal structure of its Fab fragment by molecular replacement and refined it at 3.0 Å resolution. The antigen binding site has a deep cavity at its center, which is lined mainly by acidic and hydrophobic residues, and is large enough to accommodate several antigen residues. The structure of the Fab 1696 could form a starting basis for the design of alternative HIV protease-inhibiting molecules of broad specificity.

The protein–protein interaction energy of 12 nonhomologous serine protease-inhibitor and 15 antibody-antigen complexes is calculated using a molecular mechanics formalism and dissected in terms of the main-chain vs. side-chain contribution, nonrotameric side-chain contributions, and amino acid residue type involvement in the interface interaction.

There are major differences in the interactions of the two types of protein–protein complex. Protease-inhibitor complexes interact predominantly through a main-chain–main-chain mechanism while antibody-antigen complexes interact predominantly through a side-chain–side-chain or a side-chain–main-chain mechanism. However, there is no simple correlation between the main-chain–main-chain interaction energy and the percentage of main-chain surface area buried on binding. The interaction energy is equally effected by the presence of nonrotameric side-chain conformations, which constitute ∼20% of the interaction energy. The ability to reproduce the interface interaction energy of the crystal structure if original side-chain conformations are removed from the calculation is much greater in the protease-inhibitor complexes than the antibody-antigen complexes. The success of a rotameric model for protein–protein docking appears dependent on the extent of the main-chain–main-chain contribution to binding.

Analysis of (1) residue type and (2) residue pair interactions at the interface show that antibody-antigen interactions are very restricted with over 70% of the antibody energy attributable to just six residue types (Tyr > Asp > Asn > Ser > Glu > Trp) in agreement with previous studies on residue propensity. However, it is found here that 50% of the antigen energy is attributable to just four residue types (Arg = Lys > Asn > Asp). On average just 12 residue pair interactions (6%) contribute over 40% of the favorable interaction energy in the antibody-antigen complexes, with charge-charge and charge/polar-tyrosine interactions being prominent. In contrast protease inhibitors use a diverse set of residue types and residue pair interactions.

β-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.

BBA1 is a designed protein that has only 23 residues. It is the smallest protein without disulfide bridges that has a well-defined tertiary structure in solution. We have performed unfolding molecular dynamics simulations on BBA1 and some of its mutants at 300, 330, 360, and 400 K to study their kinetic stability as well as the unfolding mechanism of BBA1. It was shown that the unfolding simulations can provide insights into the forces that stabilize the protein. Packing, hydrophobic interactions, and a salt bridge between Asp12 and Lys16 were found to be important to the protein's stability. The unfolding of BBA1 goes through two major steps: (1) disruption of the hydrophobic core and (2) unfolding of the helix. The β-hairpin remains stable in the unfolding because of the high stability of the type II′ turn connecting the two β-strands.