The in vitro refolding of hen egg-white lysozyme is studied in the presence of various osmolytes. Proline is found to prevent aggregation during protein refolding. However, other osmolytes used in this study fail to exhibit a similar property. Experimental evidence suggests that proline inhibits protein aggregation by binding to folding intermediate(s) and trapping the folding intermediate(s) into enzymatically inactive, “aggregation-insensitive” state(s). However, elimination of proline from the refolded protein mixture results in significant recovery of the bacteriolytic activity. At higher concentrations (>1.5 M), proline is shown to form loose, higher-order molecular aggregate(s). The supramolecular assembly of proline is found to possess an amphipathic character. Formation of higher-order aggregates is believed to be crucial for proline to function as a protein folding aid. In addition to its role in osmoregulation under water stress conditions, the results of this study hint at the possibility of proline behaving as a protein folding chaperone.

We have modeled the structure of human lymphotactin (hLpnt), by homology modeling and molecular dynamics simulations. This chemokine is unique in having a single disulfide bond and a long C-terminal tail. Because other structural classes of chemokines have two pairs of Cys residues, compared to one in Lpnt, and because it has been shown that both disulfide bonds are required for stability and function, the question arises how the Lpnt maintains its structural integrity. The initial structure of hLpnt was constructed by homology modeling. The first 63 residues in the monomer of hLpnt were modeled using the structure of the human CC chemokine, RANTES, whose sequence appeared most similar. The structure of the long C-terminal tail, missing in RANTES, was taken from the human muscle fatty-acid binding protein. In a Protein Data Bank search, this protein was found to contain a sequence that was most homologous to the long tail. Consequently, the modeled hLpnt C-terminal tail consisted of both α-helical and β-motifs. The complete model of the hLpnt monomer consisted of two α-helices located above the five-stranded β-sheet. Molecular dynamics simulations of the solvated initial model have indicated that the stability of the predicted fold is related to the geometry of Pro78. The five-stranded β-sheet appeared to be preserved only when Pro78 was modeled in the cis conformation. Simulations were also performed both for the C-terminal truncated forms of the hLpnt that contained one or two (CC chemokine-like) disulfide bonds, and for the chicken Lpnt (cLpnt). Our MD simulations indicated that the turn region (T30–G34) in hLpnt is important for the interactions with the receptor, and that the long C-terminal region stabilizes both the turn (T30–G34) and the five-stranded β-sheet. The major conclusion from our theoretical studies is that the lack of one disulfide bond and the extension of the C-terminus in hLptn are mutually complementary. It is very likely that removal of two Cys residues sufficiently destabilizes the structure of a chemokine molecule, particularly the core β-sheet, to abolish its biological function. However, this situation is rectified by the long C-terminal segment. The role of this long region is most likely to stabilize the first β-turn region and α-helix H1, explaining how this chemokine can function with a single disulfide bond.

Comparative protein structure prediction is limited mostly by the errors in alignment and loop modeling. We describe here a new automated modeling technique that significantly improves the accuracy of loop predictions in protein structures. The positions of all nonhydrogen atoms of the loop are optimized in a fixed environment with respect to a pseudo energy function. The energy is a sum of many spatial restraints that include the bond length, bond angle, and improper dihedral angle terms from the CHARMM-22 force field, statistical preferences for the main-chain and side-chain dihedral angles, and statistical preferences for nonbonded atomic contacts that depend on the two atom types, their distance through space, and separation in sequence. The energy function is optimized with the method of conjugate gradients combined with molecular dynamics and simulated annealing. Typically, the predicted loop conformation corresponds to the lowest energy conformation among 500 independent optimizations. Predictions were made for 40 loops of known structure at each length from 1 to 14 residues. The accuracy of loop predictions is evaluated as a function of thoroughness of conformational sampling, loop length, and structural properties of native loops. When accuracy is measured by local superposition of the model on the native loop, 100, 90, and 30% of 4-, 8-, and 12-residue loop predictions, respectively, had <2 Å RMSD error for the mainchain N, Cα, C, and O atoms; the average accuracies were 0.59 ± 0.05, 1.16 ± 0.10, and 2.61 ± 0.16 Å, respectively. To simulate real comparative modeling problems, the method was also evaluated by predicting loops of known structure in only approximately correct environments with errors typical of comparative modeling without misalignment. When the RMSD distortion of the main-chain stem atoms is 2.5 Å, the average loop prediction error increased by 180, 25, and 3% for 4-, 8-, and 12-residue loops, respectively. The accuracy of the lowest energy prediction for a given loop can be estimated from the structural variability among a number of low energy predictions. The relative value of the present method is gauged by (1) comparing it with one of the most successful previously described methods, and (2) describing its accuracy in recent blind predictions of protein structure. Finally, it is shown that the average accuracy of prediction is limited primarily by the accuracy of the energy function rather than by the extent of conformational sampling.

Calmodulin (CaM), the ubiquitous, eukaryotic, bilobal calcium-binding regulatory protein, has been cleaved by thrombin to create two fragments, TM1 (1–106) and TM2 (107–148). NMR and CD results indicate that TM1 and TM2 can associate in the presence of Ca2+ to form a complex similar to native CaM, even though the cleavage site is not in the linker region between two helix-loop-helix domains, but rather within an α-helix. Cadmium-113 NMR results show that this complex has enhanced metal-ion binding properties when compared to either TM1 or TM2 alone. This complex can bind several CaM-binding target peptides, as shown by gel bandshift assays, circular dichroism spectra, and 13C NMR spectra of biosynthetically methyl-13C-Met-labeled TM1 and TM2; moreover, gel bandshift assays show that the addition of a target peptide strengthens the interactions between TM1 and TM2 and increases the stability of the complex. Cadmium-113 NMR spectra indicate that the TM1:TM2 complex can also bind the antipsychotic drug trifluoperazine. However, in contrast to CaM:peptide complexes, the TM1:TM2:peptide complexes are disrupted by 4 M urea; moreover, TM1 and TM2 in combination are unable to activate CaM-dependent enzymes. This suggests that TM1:TM2 mixtures cannot bind target molecules as tightly as intact CaM, or perhaps that binding occurs but additional interactions with the target enzymes that are necessary for proper activation are perturbed by the proteolytic cleavage. The results presented here reflect the importance of the existence of helix-loop-helix Ca2+-binding domains in pairs in proteins such as CaM, and extend the understanding of the association of such domains in this class of proteins in general.

Domains belonging to the immunoglobulin-like fold are responsible for a wide variety of molecular recognition processes. Here we describe a new family of domains, the HYR family, which is predicted to belong to this fold, and which appears to be involved in cellular adhesion. HYR domains were identified in several eukaryotic proteins, often associated with Complement Control Protein (CCP) modules or arranged in multiple copies. Our analysis provides a sequence and structural basis for understanding the role of these domains in interaction mechanisms and leads to further characterization of heretofore undescribed repeated domains with similar folds found in several bacterial proteins involved in enzymatic activities (some chitinases) or in cell surface adhesion (streptococcal C-alpha antigen).

A simple method for growing protein crystals in the metastable zone using the vapor diffusion technique is described. The coverslips holding the hanging drops are transferred, after being incubated for some time at conditions normally giving many small crystals, over reservoirs at concentrations that normally yield clear drops. Fewer, much larger and better diffracting crystals are obtained, compared with conventional crystallization at similar conditions. To our knowledge, this is the first report of a significant crystal improvement due to “backing off” from nucleation conditions, using the hanging drop method.

A correlation of the transfer time with published results for vapor diffusion equilibration of poly(ethylene glycol) solutions is also presented.

The fluorescence time decay parameters of the β-lactoglobulin-1-anilinonaphthalene-8-sulfonate complex have been investigated under physical and chemical perturbations (2 < pH < 8 and added electrolyte 0 < NaCl < 0.5 M) to obtain new insight on the nature of the protein binding interactions. A double exponential decay of the bound probe lifetime has been confirmed by the presence of a longer component, 11 to 14.5 ns, and a shorter component, 2.5 to 3.5 ns. The two lifetimes are ascribed to different binding modes associated also with different exposure to the solvent; in particular, the longer component is attributed to binding inside the hydrophobic beta barrel, while a “surface” site is suggested for the shorter component. A detailed analysis of the lifetime fractional intensities correlates the binding constants with ionic strength and supports the presence of electrostatic effects at both sites. A Debye–Hückel approach, applied to extrapolate the electrostatic free energy contribution vs. pH at vanishing ionic strength, gives interesting clues on the effective charge felt by the ANS ligands in the proximity of each site. In particular, binding is found to parallel the aspartate and glutamate titrations between pH 3 and pH 4.5; the “surface” site mainly responds to the presence of these local titrating charges while the “internal” site more closely follows the overall protein net charge.

The effects of 1,1,1,3,3,3-hexafluoro-isopropanol (HFIP) on the conformation of cytochrome c (cyt c) at pH 1.9 were studied using a combination of spectroscopic and physical methods. Analysis varying the HFIP concentration showed that a compact denatured conformation (MHF) accumulates in a low concentration range of HFIP in the middle of structural transition from the highly unstructured acid-denatured state to the highly helical alcohol-denatured state of cyt c. This contrasts clearly with the effect of isopropanol (IP), in which no compact conformation accompanied with the transition. Analysis varying concentrations of HFIP and NaCl concurrently showed that the MHF state of cyt c is essentially identical to the salt-induced molten-globule (MG) state, and the MG state in the presence of salt was also stabilized by a low concentration of HFIP. Furthermore, 2,2,2-trifluoroethanol stabilized MHF similarly to HFIP, supporting the proposition that the specific effect observed for HFIP is caused by fluorination of alcohol. The mechanism stabilizing compact conformation by HFIP remains unclear, but is probably distinct from that of salts and polyols, which are also known to stabilize the MG-like state.

The issue of specificity in tyrosine kinase intracellular signaling mediated by src homology 2 (SH2) domains has great importance in the understanding how individual signals maintain their mutual exclusivity and affect downstream responses. Several proteins contain tandem SH2 domains that, on interacting with their ligand, provide a higher level of specificity than can be afforded by the interaction of a single SH2 domain. In this study, we focus on the comparison of two proteins ZAP70 and the p85 subunit of PI 3-kinase, which although distinctly different in function and general structure, possess tandem SH2 domains separated by a linker region and which bind to phosphorylated receptor molecules localized to the cell membrane. Binding studies using isothermal titration calorimetry show that these two proteins interact with peptides mimicking their physiological ligands in very different ways. In the case of the SH2 domains from ZAP70, they interact with a stoichiometry of unity, while p85 is able to make two distinct interactions, one with a stoichiometry of 1:1 and the other with two p85 molecules interacting with one receptor. The observation of two different modes of binding of p85 might be important in providing different cellular responses based on fluctuating intracellular concentration regimes of this protein. Thermodynamic data on both proteins suggest that a conformational change occurs on binding. On investigation of this structural change using a truncated form of p85 (including just the two SH2 domains and the inter-SH2 region), both NMR and circular dichroism spectroscopic studies failed to show significant changes in secondary structure. This suggests that any conformational change associated with binding is small and potentially limited to loop regions of the protein.

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

Six helix surface positions of protein G (Gβ1) were redesigned using a computational protein design algorithm, resulting in the five fold mutant Gβ1m2. Gβ1m2 is well folded with a circular dichroism spectrum nearly identical to that of Gβ1, and a melting temperature of 91 °C, ∼6 °C higher than that of Gβ1. The crystal structure of Gβ1m2 was solved to 2.0 Å resolution by molecular replacement. The absence of hydrogen bond or salt bridge interactions between the designed residues in Gβ1m2 suggests that the increased stability of Gβ1m2 is due to increased helix propensity and more favorable helix dipole interactions.

Metabotropic glutamate receptors (mGluRs) belong to the family 3 of G-protein-coupled receptors. On these proteins, agonist binding on the extracellular domain leads to conformational changes in the 7-transmembrane domains required for G-protein activation. To elucidate the structural features that might be responsible for such an activation mechanism, we have generated models of the amino terminal domain (ATD) of type 4 mGluR (mGlu4R). The fold recognition search allowed the identification of three hits with a low sequence identity, but with high secondary structure conservation: leucine isoleucine valine-binding protein (LIVBP) and leucine-binding protein (LBP) as already known, and acetamide-binding protein (AmiC). These proteins are characterized by a bilobate structure in an open state for LIVBP/LBP and a closed state for AmiC, with ligand binding in the cleft. Models for both open and closed forms of mGlu4R ATD have been generated. ACPT-I (1-aminocyclopentane 1,3,4-tricarboxylic acid), a selective agonist, has been docked in the two models. In the open form, ACPT-I is only bound to lobe I through interactions with Lys74, Arg78, Ser159, and Thr182. In the closed form, ACPT-I is trapped between both lobes with additional binding to Tyr230, Asp312, Ser313, and Lys317 from lobe II. These results support the hypothesis that mGluR agonists bind a closed form of the ATDs, suggesting that such a conformation of the binding domain corresponds to the active conformation.

All the protein sequences from SWISS-PROT database were analyzed for occurrence of single amino acid repeats, tandem oligo-peptide repeats, and periodically conserved amino acids. Single amino acid repeats of glutamine, serine, glutamic acid, glycine, and alanine seem to be tolerated to a considerable extent in many proteins. Tandem oligo-peptide repeats of different types with varying levels of conservation were detected in several proteins and found to be conspicuous, particularly in structural and cell surface proteins. It appears that repeated sequence patterns may be a mechanism that provides regular arrays of spatial and functional groups, useful for structural packing or for one to one interactions with target molecules. To facilitate further explorations, a database of Tandem Repeats in Protein Sequences (TRIPS) has been developed and is available at URL: http://www.ncl-india.org/trips.

The rate of macromolecular surface formation in yeast iso-2 cytochrome c and its site-specific mutant, N52I iso-2, has been studied using a monoclonal antibody that recognizes a tertiary epitope including K58 and H39. The results indicate that epitope refolding occurs after fast folding but prior to slow folding, in contrast to horse cytochrome c where surface formation occurs early. The antibody-detected (ad) kinetic phase accompanying epitope formation has kad = 0.2 s−1 and is ∼40-fold slower than the fastest detectable event in the folding of yeast iso-2 cytochrome c (k2f ∼ 8 s−1), but occurs prior to the absorbance- and fluorescence-detected slow folding steps (k1a ∼ 0.06 s−1; k1b ∼ 0.09 s−1). N52I iso-2 cytochrome c exhibits similar kinetic behavior with respect to epitope formation. A detailed dissection of the mechanistic differences between the folding pathways of horse and yeast cytochromes c identifies possible reasons for the slow surface formation in the latter. Our results suggest that non-native ligation involving H33 or H39 during refolding may slow down the formation of the tertiary epitope in iso-2 cytochrome c. This study illustrates that surface formation can be coupled to early events in protein folding. Thus, the rate of macromolecular surface formation is fine tuned by the residues that make up the surface and the interactions they entertain during refolding.

Prolyl oligopeptidase, an enzyme implicated in memory disorders, is a member of a new serine peptidase family. Crystallographic studies (Fülöp et al., 1998) revealed a novel oxyanion binding site containing a tyrosine residue, Tyr473. To study the importance of Tyr473 OH, we have produced prolyl oligopeptidase and its Tyr473Phe variant in Escherichia coli. The specificity rate constant, kcat/Km, for the modified enzyme decreased by a factor of 8–40 with highly specific substrates, Z-Gly-Pro-Nap, and a fluorogenic octapeptide. With these compounds, the decline in kcat was partly compensated for by reduction in Km, a difference from the extensively studied subtilisin. With the less specific suc-Gly-Pro-Nap, the Km value, which approximates Ks, was not significantly changed, resulting in greater diminution (∼500-fold) in kcat/Km. The second-order rate constant for the reaction with Z-Pro-prolinal, a slow tight-binding transition-state analogue inhibitor, and the Ki values for a slow substrate and two product-like inhibitors were not significantly affected by the Tyr473 OH group. The mechanism of transition-state stabilization was markedly dependent upon the nature of substrate and varied with pH as the enzyme interconverted between its two catalytically competent forms.

Using a cross-linking approach, we recently demonstrated that radiolabeled peptides or misfolded proteins specifically interact in vitro with two luminal proteins in crude extracts from pancreas microsomes. The proteins were the folding catalysts protein disulfide isomerase (PDI) and PDIp, a glycosylated, PDI-related protein, expressed exclusively in the pancreas. In this study, we explore the specificity of these proteins in binding peptides and related ligands and show that tyrosine and tryptophan residues in peptides are the recognition motifs for their binding by PDIp. This peptide-binding specificity may reflect the selectivity of PDIp in binding regions of unfolded polypeptide during catalysis of protein folding.

We report here that in vitro exposure of monomeric actin to hydrogen peroxide leads to a conversion of 6 of the 16 methionine residues to methionine sulfoxide residues. Although the initial effect of H2O2 on actin is the oxidation of Cys374, we have found that Met44, Met47, Met176, Met190, Met269, and Met355 are the other sites of the oxidative modification. Met44 and Met47 are the methionyl sites first oxidized. The methionine residues that are oxidized are not simply related to their accessibility to the external medium and are found in all four subdomains of actin. The conformations of subdomain 1, a region critical for the functional binding of different actin-binding proteins, and subdomain 2, which plays important roles in the polymerization process and stabilization of the actin filament, are changed upon oxidation. The conformational changes are deduced from the increased exposure of hydrophobic residues, which correlates with methionine sulfoxide formation, from the perturbations in tryptophan fluorescence, and from the decreased susceptibility to limited proteolysis of oxidized actin.

The prediction of binding energies from the three-dimensional (3D) structure of a protein–ligand complex is an important goal of biophysics and structural biology. Here, we critically assess the use of empirical, solvent-accessible surface area-based calculations for the prediction of the binding of Src-SH2 domain with a series of tyrosyl phosphopeptides based on the high-affinity ligand from the hamster middle T antigen (hmT), where the residue in the pY+3 position has been changed. Two other peptides based on the C-terminal regulatory site of the Src protein and the platelet-derived growth factor receptor (PDGFR) are also investigated. Here, we take into account the effects of proton linkage on binding, and test five different surface area-based models that include different treatments for the contributions to conformational change and protein solvation. These differences relate to the treatment of conformational flexibility in the peptide ligand and the inclusion of proximal ordered solvent molecules in the surface area calculations. This allowed the calculation of a range of thermodynamic state functions (ΔCp, ΔS, ΔH, and ΔG) directly from structure. Comparison with the experimentally derived data shows little agreement for the interaction of SrcSH2 domain and the range of tyrosyl phosphopeptides. Furthermore, the adoption of the different models to treat conformational change and solvation has a dramatic effect on the calculated thermodynamic functions, making the predicted binding energies highly model dependent. While empirical, solvent-accessible surface area based calculations are becoming widely adopted to interpret thermodynamic data, this study highlights potential problems with application and interpretation of this type of approach. There is undoubtedly some agreement between predicted and experimentally determined thermodynamic parameters; however, the tolerance of this approach is not sufficient to make it ubiquitously applicable.

The three-dimensional structure of the 56 residue polypeptide Apis mellifera chymotrypsin/cathepsin G inhibitor 1 (AMCI-1) isolated from honey bee hemolymph was calculated based on 730 experimental NMR restraints. It consists of two approximately perpendicular β-sheets, several turns, and a long exposed loop that includes the protease binding site. The lack of extensive secondary structure features or hydrophobic core is compensated by the presence of five disulfide bridges that stabilize both the protein scaffold and the binding loop segment. A detailed analysis of the protease binding loop conformation reveals that it is similar to those found in other canonical serine protease inhibitors. The AMCI-1 structure exhibits a common fold with a novel family of inhibitors from the intestinal parasitic worm Ascaris suum. The pH-induced conformational changes in the binding loop region observed in the Ascaris inhibitor ATI are absent in AMCI-1. Similar binding site sequences and structures strongly suggest that the lack of the conformational change can be attributed to a Glu → Gln substitution at the P1′ position in AMCI-1, compared to ATI. Analysis of amide proton temperature coefficients shows very good correlation with the presence of hydrogen bond donors in the calculated AMCI-1 structure.

To provide a framework for understanding the hyperthermostability of some rubredoxins, a comprehensive analysis of the thermally induced denaturation of rubredoxin (Rd) from the mesophile, Clostridium pasteurianum was undertaken. Rds with three different metals in its M(SCys)4 site (M = Fe3+/2+, Zn2+, or Cd2+) were examined. Kinetics of metal ion release were monitored anaerobically at several fixed temperatures between 40 and 100 °C, and during progressive heating of the iron-containing protein. Both methods gave a thermal stability of metal binding in the order Fe2+ [Lt ] Fe3+ < Zn2+ < Cd2+. The temperature at which half of the iron was released from the protein in temperature ramp experiments was 69 °C for Fe2+Rd and 83 °C for Fe3+Rd. Temperature-dependent changes in the protein structure were monitored by differential scanning calorimetry, tryptophan fluorescence, binding of a fluorescent hydrophobic probe, and 1H NMR. Major but reversible structural changes, consisting of swelling of the hydrophobic core and opening of a loop region, were found to occur at temperatures (50–70 °C) much lower than those required for loss of the metal ion. For the three divalent metal ions, the results suggest that the onset of the reversible, lower-temperature structural changes is dependent on the size of the MS4 site, whereas the final, irreversible loss of metal ion is dependent on the inherent M-SCys bond strength. In the case of Fe3+Rd, stoichiometric Fe3+/cysteine–ligand redox chemistry also occurs during metal ion loss. The results indicate that thermally induced unfolding of the native Cp Rd must surmount a significant kinetic barrier caused by stabilizing interactions both within the protein and within the M(SCys)4 site.