The amino acid sequence of the small copper protein auracyanin A isolated from the thermophilic photosynthetic green bacterium Chloroflexus aurantiacus has been determined to be a polypeptide of 139 residues. His58, Cys123, His128, and Met132 are spaced in a way to be expected if they are the evolutionary conserved metal ligands as in the known small copper proteins plastocyanin and azurin. Secondary structure prediction also indicates that auracyanin has a general β-barrel structure similar to that of azurin from Pseudomonas aeruginosa and plastocyanin from poplar leaves. However, auracyanin appears to have sequence characteristics of both small copper protein sequence classes. The overall similarity with a consensus sequence of azurin is roughly the same as that with a consensus sequence of plastocyanin, namely 30.5%. We suggest that auracyanin A, together with the B forms, is the first example of a new class of small copper proteins that may be descendants of an ancestral sequence to both the azurin proteins occurring in prokaryotic nonphotosynthetic bacteria and the plastocyanin proteins occurring in both prokaryotic cyanobacteria and eukaryotic algae and plants. The N-terminal sequence region 1–18 of auracyanin is remarkably rich in glycine and hydroxy amino acids, and required mass spectrometric analysis to be determined. The nature of the blocking group X is not yet known, although its mass has been determined to be 220 Da. The auracyanins are the first small blue copper proteins found and studied in anoxygenic photosynthetic bacteria and are likely to mediate electron transfer between the cytochrome bc1 complex and the photosynthetic reaction center.

Alanine scanning mutagenesis of the HyHEL-10 paratope of the HyHEL-10/HEWL complex demonstrates that the energetically important side chains (hot spots) of both partners are in contact. A plot of ΔΔGHyHEL-10_mutant vs. ΔΔGHEWL_mutant for the five of six interacting side-chain hydrogen bonds is linear (Slope = 1). Only 3 of the 13 residues in the HEWL epitope contribute >4 kcal/mol to the free energy of formation of the complex when replaced by alanine, but 6 of the 12 HyHEL-10 paratope amino acids do. Double mutant cycle analysis of the single crystallographically identified salt bridge, D32H/K97, shows that there is a significant energetic penalty when either partner is replaced with a neutral side-chain amino acid, but the D32HN/K97M complex is as stable as the WT. The role of the disproportionately high number of Tyr residues in the CDR was evaluated by comparing the ΔΔG values of the Tyr → Phe vs. the corresponding Tyr → Ala mutations. The nonpolar contacts in the light chain contribute only about one-half of the total ΔΔG observed for the Tyr → Ala mutation, while they are significantly more important in the heavy chain. Replacement of the N31L/K96 hydrogen bond with a salt bridge, N31DL/K96, destabilizes the complex by 1.4 kcal/mol. The free energy of interaction, ΔΔGint, obtained from double mutant cycle analysis showed that ΔΔGint for any complex for which the HEWL residue probed is a major immunodeterminant is very close to the loss of free energy observed for the HyHEL-10 single mutant. Error propagation analysis of double mutant cycles shows that data of atypically high precision are required to use this method meaningfully, except where large ΔΔG values are analyzed.

The accumulation of hundreds of olfactory receptor (OR) sequences, along with the recent availability of detailed models of other G-protein-coupled receptors, allows us to analyze the OR amino acid variability patterns in a structural context. A Fourier analysis of 197 multiply aligned olfactory receptor sequences showed an α-helical periodicity in the variability profile. This was particularly pronounced in the more variable transmembranal segments 3, 4, and 5. Rhodopsin-based homology modeling demonstrated that the inferred variable helical faces largely point to the interior of the receptor barrel. We propose that a set of 17 hypervariable residues, which point to the barrel interior and are more extracellularly disposed, constitute the odorant complementarity determining regions. While 12 of these residues coincide with established ligand-binding contact postions in other G-protein-coupled receptors, the rest are suggested to form an olfactory-unique aspect of the binding pocket. Highly conserved olfactory receptor-specific sequence motifs, found in the second and third intracellular loops, may comprise the G-protein recognition epitope. The prediction of olfactory receptor functional sites provides concrete suggestions of site-directed mutagenesis experiments for altering ligand and G-protein specificity.

We present a neural network based method (ChloroP) for identifying chloroplast transit peptides and their cleavage sites. Using cross-validation, 88% of the sequences in our homology reduced training set were correctly classified as transit peptides or nontransit peptides. This performance level is well above that of the publicly available chloroplast localization predictor PSORT. Cleavage sites are predicted using a scoring matrix derived by an automatic motif-finding algorithm. Approximately 60% of the known cleavage sites in our sequence collection were predicted to within ±2 residues from the cleavage sites given in SWISS-PROT. An analysis of 715 Arabidopsis thaliana sequences from SWISS-PROT suggests that the ChloroP method should be useful for the identification of putative transit peptides in genome-wide sequence data. The ChloroP predictor is available as a web-server at http://www.cbs.dtu.dk/services/ChloroP/.

Following hints from X-ray data (Ostermeier C et al., 1997, Proc Natl Acad Sci USA 94:10547–10553; Yoshikawa S et al., 1998, Science 280:1723–1729), chemical evidence is presented from four distantly related cytochrome-c oxidases for the existence of a copperB-coordinated His240–Tyr244) cross-link at the O2-activating Heme Fea3–CuB center in the catalytic subunit I of the enzyme. The early evolutionary invention of this unusual structure may have prevented demaging [bull ]OH-radical release at e−-transfer to dioxygen and thus have enabled O2 respiration.

An expanded, highly dynamic denatured state of staphylococcal nuclease exhibits a native-like topology in the apparent absence of tight packing and fixed hydrogen bonds (Gillespie JR, Shortle D, 1997, J Mol Biol 268:158–169, 170–184). To address the physical basis of the long-range spatial ordering of this molecule, we probe the effects of perturbations of the sequence and solution conditions on the local chain dynamics of a denatured 101-residue fragment that is missing the first three beta strands. Structural interactions between chain segments are inferred from correlated changes in the motional behavior of residues monitored by 15N NMR relaxation measurements. Restoration of the sequence corresponding to the first three beta strands significantly increases the average order of all chain segments that form the five strand beta barrel including loops but has no effect on the carboxy terminal 30 residues. Addition of the denaturing salt sodium perchlorate enhances ordering over the entire sequence of this fragment. Analysis of seven different substitution mutants points to a complex set of interactions between the hydrophobic segment corresponding to beta strand 5 and the remainder of the chain. General patterns in the data suggest there is a hierarchy of native-like interactions that occur transiently in the denatured state and are consistent with the overall topology of the denatured state ensemble being determined by many coupled local interactions rather than a few highly specific long-range interactions.

Methods for protein structure (3D)–sequence (1D) compatibility evaluation (threading) have been developed during the past decade. The protocol in which a sequence can recognize its compatible structure in the structural library (i.e., the fold recognition or the forward-folding search) is available for the structure prediction of new proteins. However, the reverse protocol, in which a structure recognizes its homologous sequences among a sequence database, named the inverse-folding search, is a more difficult application. In this study, we have investigated the feasibility of the latter approach. A structural library, composed of about 400 well-resolved structures with mutually dissimilar sequences, was prepared, and 163 of them had remote homologs in the library. We examined whether they could correctly seek their homologs by both forward- and inverse-folding searches. The results showed that the inverse-folding protocol is more effective than the forward-folding protocol, once the reference states of the compatibility functions are appropriately adjusted. This adjustment only slightly affects the ability of the forward-folding search. We noticed that the scoring, in which a given sequence is re-mounted onto a structure according to the 3D–1D alignment determined by the dynamic programming method, is only effective in the forward-folding protocol and not in the inverse-folding protocol. Namely, the inverse-folding search works significantly better with the score given by the 3D–1D alignment per se, rather than that obtained by the re-mounting. The implications of these results are discussed.

Human growth hormone (hGH) binds to its receptor (hGHr) in a three-body interaction: one molecule of the hormone and two identical monomers of the receptor form a trimer. Curiously, the hormone-receptor interactions in the trimer are not equivalent and the formation of the complex occurs in a specific kinetic order (Cunningham BC, Ultsch M, De Vos AM, Mulkerrin MG, Clauser KR, Wells JA, 1991, Science 254:821–825). In this paper, we model the recognition of hGH to the hGHr using shape complementarity of the three-dimensional structures and macromolecular docking to explore possible binding modes between the receptor and hormone. The method, reported previously (Hendrix DK, Kuntz ID, 1998, Pacific symposium on biocomputing 1998, pp 1234–1244), is based upon matching complementary-shaped strategic sites on the molecular surface. We modify the procedure to examine three-body systems. We find that the order of binding seen experimentally is also essential to our model. We explore the use of mutational data available for hGH to guide our model. In addition to docking hGH to the hGHr, we further test our methodology by successfully reproducing 16 macromolecular complexes from X-ray crystal structures, including enzyme-inhibitor, antibody-antigen, protein dimer, and protein-DNA complexes.

Site-directed mutagenesis was used to replace Lys68 of the human hypoxanthine phosphoribosyltransferase (HGPRTase) with alanine to exploit this less reactive form of the enzyme to gain additional insights into the structure activity relationship of HGPRTase. Although this substitution resulted in only a minimal (one- to threefold) increase in the Km values for binding pyrophosphate or phosphoribosylpyrophosphate, the catalytic efficiencies (kcat/Km) of the forward and reverse reactions were more severely reduced (6- to 30-fold), and the mutant enzyme showed positive cooperativity in binding of α-d-5-phosphoribosyl-1-pyrophosphate (PRPP) and nucleotide. The K68A form of the human HGPRTase was cocrystallized with 7-hydroxy [4,3-d] pyrazolo pyrimidine (HPP) and Mg PRPP, and the refined structure reported. The PRPP molecule built into the [(Fo − Fc)φcalc] electron density shows atomic interactions between the Mg PRPP and enzyme residues in the pyrophosphate binding domain as well as in a long flexible loop (residues Leu101 to Gly111) that closes over the active site. Loop closure reveals the functional roles for the conserved SY dipeptide of the loop as well as the molecular basis for one form of gouty arthritis (S103R). In addition, the closed loop conformation provides structural information relevant to the mechanism of catalysis in human HGPRTase.

At least four mRNAs for oat phytochrome A (phyA) are present in etiolated oat tissue. The complete amino acid sequences of two phyA isoforms (A3 and A4) and the N-terminal amino acid sequence of a third isoform (A5) were deduced from cDNA sequencing (Hershey et al., 1985). In the present study, heterogeneity of phyA on a protein level was studied by tryptic mapping using electrospray ionization mass-spectrometry (ESIMS). The total tryptic digest of iodoacetamide-modified phyA was fractionated by gel filtration chromatography followed by reversed-phase high-performance liquid chromatography. ESIMS was used to identify peptides. Amino acid sequences of the peptides were confirmed or determined by collision-induced dissociation mass spectrometry (CID MS), MS/MS, or by subdigestion of the tryptic peptides followed by ESIMS analysis. More than 97% of the phyA3 sequence (1,128 amino acid residues) was determined in the present study. Mass-spectrometric analysis of peptides unique to each form showed that phyA purified from etiolated oat seedling is represented by three isoforms A5, A3, and A4, with ratio 3.4:2.3:1.0. Possible light-induced changes in phytochrome in vivo phosphorylation site at Ser7 (Lapko VN et al., 1997, Biochemistry 36:10595–10599) as well at Ser17 and Ser598 (known as in vitro phosphorylation sites) were also analyzed. The extent of phosphorylation at Ser7 appears to be the same for phyA isolated from dark-grown and red-light illuminated seedlings. In addition to Ser7, Ser598 was identified as an in vivo phosphorylation site in oat phyA. Ser598 phosphorylation was found only in phyA from the red light-treated seedlings, suggesting that the protein phosphorylation plays a functional role in the phytochrome A-mediated light-signal transduction.

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

The “calorimetric criterion” is one of the important experimental approaches for determining whether protein folding is an “all-or-none” two-state transition (i.e., whether intermediates are present at equilibrium). The calorimetric criterion states that the equivalence of the “measured” calorimetric enthalpy change and the effective two-state van't Hoff enthalpy change demonstrates that there is a two-state transition. This paper addresses the essential question of whether the calorimetric criterion is a necessary and sufficient condition for a two-state process and shows that it is necessary but not sufficient by means of specific examples. Analysis of simple models indicates that the heat capacity curve, regardless of whether it originates from a two-state process or not, can always be decomposed in such a way that the calorimetric criterion is satisfied. Exact results for a three-state model and a homopolymer tetramer demonstrate that the deviation from the calorimetric criterion is not simply related to the population of intermediate states. Analysis of a three-helix bundle protein model, which has a two-state folding from a random coil to ordered (molten) globule, shows that the calorimetric criterion may not be satisfied if the standard linear interpolation of baselines (weighted or unweighted) is employed. A specific example also suggests that the more recently introduced deconvolution method is not necessarily better than the simple calorimetric criterion for distinguishing a two-state transition from a three-state transition. Although the calorimetric criterion is not a sufficient condition for a two-state process, it is likely to continue to be of practical utility, particularly when its results are shown to be consistent with those from other experimental 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.

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.

Recent advances in attempts to unravel the protein folding mechanism have indicated the need to identify the folding intermediates. Despite their transient nature, in a number of cases it has been possible to detect and characterize some of the equilibrium intermediates, for example, the molten globule (MG) state. The key features of the MG state are retention of substantial secondary structure of the native state, considerable loss of tertiary structure leading to increased hydrophobic exposure, and a compact structure. NMR, circular dichroism, and fluorescence spectroscopies have been most useful in characterizing such intermediates. We report here a new method for structural characterization of the MG state that involves probing the exposed hydrophobic sites with a hydrophobic photoactivable reagent—2[3H]diazofluorene. This carbene-based reagent binds to hydrophobic sites, and on photolysis covalently attaches itself to the neighboring amino acid side chains. The reagent photolabels α-lactalbumin as a function of pH (3–7.4), the labeling at neutral pH being negligible and maximal at pH 3. Chemical and proteolytic fragmentation of the photolabeled protein followed by peptide sequencing permitted identification of the labeled residues. The results obtained indicate that the sequence corresponding to B (23–34) and C (86–98) helix of the native structure are extensively labeled. The small β-domain (40–50) is poorly labeled, Val42 being the only residue that is significantly labeled. Our data, like NMR data, indicate that in the MG state of α-lactalbumin, the α-domain has a greater degree of persistent structure than the β-domain. However, unlike the NMR method, the photolabeling method is not limited by the size of the protein and can provide information on several new residues, for example, Leu115. The current method using DAF thus allows identification of stable and hydrophobic exposed regions in folding intermediates as the reagent binds and on photolysis covalently links to these regions.

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

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

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

A popular approach to the computational modeling of ligand/receptor interactions is to use an empirical free energy like model with adjustable parameters. Parameters are learned from one set of complexes, then used to predict another set. To improve these empirical methods requires an independent way to study their inherent errors. We introduce a toy model of ligand/receptor binding as a workbench for testing such errors. We study the errors incurred from the two state binding assumption—the assumption that a ligand is either bound in one orientation, or unbound. We find that the two state assumption can cause large errors in free energy predictions, but it does not affect rank order predictions significantly. We show that fitting parameters using data from high affinity ligands can reduce two state errors; so can using more physical models that do not use the two state assumption. We also find that when using two state models to predict free energies, errors are more severe on high affinity ligands than low affinity ligands. And we show that two state errors can be diagnosed by systematically adding new binding modes when predicting free energies: if predictions worsen as the modes are added, then the two state assumption in the fitting step may be at fault.