Interior and surface of monomeric proteins

Hydrophilic framework in proteins?

The spatial neighborhood composition of residues was determined in a 511-structure set by taking only side-chain atoms into account to generate a hydrophobicity scale. This scale is symmetrical and has been divided into seven functional groups. Hydrophobic (LIVFMCAWYG) and hydrophilic (PTHSQRNKED) residues obey an equipartition rule: not only are they found in equal proportions, but they play equivalent roles in many of their properties. The nearest neighbors of all residues are always hydrophilic. However, hydrophobic residues are mostly surrounded by other hydrophobic residues located at a peak at 3.9 A, while hydrophilic residues show three peaks at 5.0, 6.5, and 8.0 A, suggesting a hydrophilic structural framework. This leads us to question the importance of hydrophobic cores believed to be at the origin of protein folding.

Structure-based thermodynamic analysis of the dissociation of protein phosphatase-1 catalytic subunit and microcystin-LR docked complexes

The relationship between the structure of a free ligand in solution and the structure of its bound form in a complex is of great importance to the understanding of the energetics and mechanism of molecular recognition and complex formation. In this study, we use a structure-based thermodynamic approach to study the dissociation of the complex between the toxin microcystin-LR (MLR) and the catalytic domain of protein phosphatase-1 (PP-1c) for which the crystal structure of the complex is known. We have calculated the thermodynamic parameters (enthalpy, entropy, heat capacity, and free energy) for the dissociation of the complex from its X-ray structure and found the calculated dissociation constant (4.0 x 10(-11)) to be in excellent agreement with the reported inhibitory constant (3.9 x 10(-11)). We have also calculated the thermodynamic parameters for the dissociation of 47 PP-1c:MLR complexes generated by docking an ensemble of NMR solution structures of MLR onto the crystal structure of PP-1c. In general, we observe that the lower the root-mean-square deviation (RMSD) of the docked complex (compared to the X-ray complex) the closer its free energy of dissociation (deltaGd(o)) is to that calculated from the X-ray complex. On the other hand, we note a significant scatter between the deltaGd(o) and the RMSD of the docked complexes. We have identified a group of seven docked complexes with deltaGd(o) values very close to the one calculated from the X-ray complex but with significantly dissimilar structures. The analysis of the corresponding enthalpy and entropy of dissociation shows a compensation effect suggesting that MLR molecules with significant structural variability can bind PP-1c and that substantial conformational flexibility in the PP-1c:MLR complex may exist in solution.

Self-assembled complexes of oligopeptides and metalloporphyrins: measurements of the reorganization and electronic interaction energies for photoinduced electron-transfer reactions

Cationic porphyrins form ground state electrostatically associated complexes with anionic oligo-electrolytes such as those formed by a series of glutamic acid (E) residues. Temperature dependencies were measured of the rate constants for intra-complex electron transfer to the triplet state of Pd(II)TMPyP4+ from a tyrosine (tyr, Y) or tryptophan (trp, W) moiety connected to a glutamic acid tetramer. In complexes such as YE4, E2YE2, YE4G10E (G, glycine), and WE4 these data were used to estimate the reorganization energy (lambda) and electronic interaction energy (HDA) relevant to the process. For all tyr-peptide complexes, lambda values were found to be large (lambda approximately 1.60 +/- 0.06 eV), reflecting a relatively high medium polarity in the vicinity of tyr residues. It further indicates that the tyr residues in all oligo-peptides are exposed to the aqueous medium in a similar way irrespective of the position of the aromatic moiety in the peptide chain. A significantly lower lambda value (lambda = 1.08 eV) was derived for the tryptophan-containing peptide complex, indicating a relatively higher hydrophobic character of trp compared to tyr. The electronic coupling matrix elements (HDA) derived for tyr-peptide complexes (5.1 meV for YE4, 5.4 meV for YE4G10E and 7.5 meV for E2YE2) were larger than that found for WE4 (1.1 meV). Molecular dynamics calculations were employed to obtain structural features of the porphyrin-peptide complexes. These showed average distances between the center of mass (COM) of the porphyrin ring and the center of mass of the amino acid aromatic ring of 816 +/- 140 pm (YE4), 800 +/- 80 pm (E2YE2), 900 +/- 130 pm (YE4G10E) and 970 +/- 160 pm (WE4). The molecular dynamics calculations were shown to be in good agreement with the experimentally determined electronic interaction energies, strongly suggesting that HDA is primarily responsible for the dependence of the electron-transfer rate constant (KET) on the donor-acceptor separation distance and relative orientation. The higher HDA (7.55 meV) derived for tyr incorporated into the middle of the peptide backbone (E2YE2) was presumed to be associated with a higher degree of orbital overlap due to a more favorable ring-ring orientation. Overlap parameters (beta derived for all peptide-porphyrin complexes were similar (approximately 0.95 +/- 0.06 A-1), being in good agreement with most literature values for similar systems. Finally, the intra-complex electron-transfer ratio (ktrp/ktyr) derived from flash photolysis experiments and the corresponding ratio derived from Marcus' theory combined with experimental data from the temperature-dependence investigations and electrochemical measurements were found to be in excellent agreement. This same consistency was found for the couple E4Y and E2YE2. The empirical expression (Moser and Dutton) governing the intraprotein electron-transfer rate constant in native systems combined with our experimental data (kET, lambda, delta G0) yielded tunneling pathway distances in excellent agreement with those arising from the molecular modeling studies. The exception was for the long peptide YE4G10E, for which the Quenched Molecular Dynamic (QMD) sampling technique was complicated and is probably inadequate.

kPROT: a knowledge-based scale for the propensity of residue orientation in transmembrane segments. Application to membrane protein structure prediction

Modeling of integral membrane proteins and the prediction of their functional sites requires the identification of transmembrane (TM) segments and the determination of their angular orientations. Hydrophobicity scales predict accurately the location of TM helices, but are less accurate in computing angular disposition. Estimating lipid-exposure propensities of the residues from statistics of solved membrane protein structures has the disadvantage of relying on relatively few proteins. As an alternative, we propose here a scale of knowledge-based Propensities for Residue Orientation in Transmembrane segments (kPROT), derived from the analysis of more than 5000 non-redundant protein sequences. We assume that residues that tend to be exposed to the membrane are more frequent in TM segments of single-span proteins, while residues that prefer to be buried in the transmembrane bundle interior are present mainly in multi-span TMs. The kPROT value for each residue is thus defined as the logarithm of the ratio of its proportions in single and multiple TM spans. The scale is refined further by defining it for three discrete sections of the TM segment; namely, extracellular, central, and intracellular. The capacity of the kPROT scale to predict angular helical orientation was compared to that of alternative methods in a benchmark test, using a diversity of multi-span alpha-helical transmembrane proteins with a solved 3D structure. kPROT yielded an average angular error of 41 degrees, significantly lower than that of alternative scales (62 degrees -68 degrees ). The new scale thus provides a useful general tool for modeling and prediction of functional residues in membrane proteins. A WWW server (http://bioinfo.weizmann.ac.il/kPROT) is available for automatic helix orientation prediction with kPROT.

De novo protein design. II. Plasticity in sequence space

It is generally accepted that many different protein sequences have similar folded structures, and that there is a relatively high probability that a new sequence possesses a previously observed fold. An indirect consequence of this is that protein design should define the sequence space accessible to a given structure, rather than providing a single optimized sequence. We have recently developed a new approach for protein sequence design, which optimizes the complete sequence of a protein based on the knowledge of its backbone structure, its amino acid composition and a physical energy function including van der Waals interactions, electrostatics, and environment free energy. The specificity of the designed sequence for its template backbone is imposed by keeping the amino acid composition fixed. Here, we show that our procedure converges in sequence space, albeit not to the native sequence of the protein. We observe that while polar residues are well conserved in our designed sequences, non-polar amino acids at the surface of a protein are often replaced by polar residues. The designed sequences provide a multiple alignment of sequences that all adopt the same three-dimensional fold. This alignment is used to derive a profile matrix for chicken triose phosphate isomerase, TIM. The matrix is found to recognize significantly the native sequence for TIM, as well as closely related sequences. Possible application of this approach to protein fold recognition is discussed.

Synthesis of a stable form of tertiapin: a high-affinity inhibitor for inward-rectifier K+ channels

Tertiapin (TPN), a small protein derived from honey bee venom, inhibits the GIRK1/4 and ROMK1 channels with nanomolar affinities. Methionine residue 13 in TPN interacts with residue F148 in the channel, located just outside of the narrow region of the ROMK1 pore. The methionine residue in TPN can be oxidized by air, which significantly hinders TPN binding to the channels. To overcome the reduction in TPN affinity due to oxidation of M13, we replaced M13 in TPN with fourteen different residues. Out of the fourteen derivatives, only the one in which M13 was replaced by glutamine, TPNQ, binds to the channel with a Ki value very similar to that of native TPN. Since TPNQ is stable and functionally resembles native TPN, it will be a very useful molecular probe for studying the inward-rectifier K+ channels.

Evolution of the genetic code

Comparative path lengths in amino acid biosynthesis and other molecular indicators of the timing of codon assignment were examined to reconstruct the main stages of code evolution. The codon tree obtained was rooted in the 4 N-fixing amino acids (Asp, Glu, Asn, Gln) and 16 triplets of the NAN set. This small, locally phased (commaless) code evidently arose from ambiguous translation on a poly(A) collector strand, in a surface reaction network. Copolymerisation of these amino acids yields polyanionic peptide chains, which could anchor uncharged amide residues to a positively charged mineral surface. From RNA virus structure and replication in vitro, the first genes seemed to be RNA segments spliced into tRNA. Expansion of the code reduced the risk of mutation to an unreadable codon. This step was conditional on initiation at the 5'-codon of a translated sequence. Incorporation of increasingly hydrophobic amino acids accompanied expansion. As codons of the NUN set were assigned most slowly, they received the most nonpolar amino acids. The origin of ferredoxin and Gln synthetase was traced to mid-expansion phase. Surface metabolism ceased by the end of code expansion, as cells bounded by a proteo-phospholipid membrane, with a protoATPase, had emerged. Incorporation of positively charged and aromatic amino acids followed. They entered the post-expansion code by codon capture. Synthesis of efficient enzymes with acid-base catalysis was then possible. Both types of aminoacyl-tRNA synthetases were attributed to this stage. tRNA sequence diversity and error rates in RNA replication indicate the code evolved within 20 million yr in the preIsuan era. These findings on the genetic code provide empirical evidence, from a contemporaneous source, that a surface reaction network, centred on C-fixing autocatalytic cycles, rapidly led to cellular life on Earth.

Design of highly stable functional GroEL minichaperones

GroEL minichaperones have potential in the biotechnology industry for the refolding of recombinant proteins. With the aim of enhancing and widening their use, we have created two highly stable functional variants of minichaperone GroEL(193-345). A sequence alignment of 130 members of the chaperonin 60 (Cpn60) family was used to design 37 single mutations. Two small-to-large mutations, A223T, A223V and one similar-size mutation, M233L, all located in the hydrophobic core were found to stabilize the protein by more than 1 kcal mol(-1) each. Six stabilizing mutations were combined, yielding two multiple mutants that were 6.99 and 6.15 kcal mol(-1) more stable than wild-type protein. Even though some of the substituted residue pairs are close to each other in the protein structure, the energetic effects of mutation are approximately additive. In particular, the stabilizing substitution A223T is unexpected and would have been missed by purely structural analysis. In the light of previously reported successes employing similar methods with several other proteins, our results show that a homology based approach is a simple and efficient method of increasing the stability of a protein.

An immunoglobulin-like fold in a major plant allergen: the solution structure of Phl p 2 from timothy grass pollen

BACKGROUND

Grass pollen allergens are the most important and widespread elicitors of pollen allergy. One of the major plant allergens which millions of people worldwide are sensitized to is Phl p 2, a small protein from timothy grass pollen. Phl p 2 is representative of the large family of cross-reacting plant allergens classified as group 2/3. Recombinant Phl p 2 has been demonstrated by immunological cross-reactivity studies to be immunologically equivalent to the natural protein.

RESULTS

We have solved the solution structure of recombinant Phl p 2 by means of nuclear magnetic resonance techniques. The three-dimensional structure of Phl p 2 consists of an all-beta fold with nine antiparallel beta strands that form a beta sandwich. The topology is that of an immunoglobulin-like fold with the addition of a C-terminal strand, as found in the C2 domain superfamily. Lack of functional and sequence similarity with these two families, however, suggests an independent evolution of Phl p 2 and other homologous plant allergens.

CONCLUSIONS

Because of the high homology with other plant allergens of groups 1 and 2/3, the structure of Phl p 2 can be used to rationalize some of the immunological properties of the whole family. On the basis of the structure, we suggest possible sites of interaction with IgE antibodies. Knowledge of the Phl p 2 structure may assist the rational structure-based design of synthetic vaccines against grass pollen allergy.

X-ray structure of the FimC-FimH chaperone-adhesin complex from uropathogenic Escherichia coli

Type 1 pili-adhesive fibers expressed in most members of the Enterobacteriaceae family-mediate binding to mannose receptors on host cells through the FimH adhesin. Pilus biogenesis proceeds by way of the chaperone/usher pathway. The x-ray structure of the FimC-FimH chaperone-adhesin complex from uropathogenic Escherichia coli at 2.5 angstrom resolution reveals the basis for carbohydrate recognition and for pilus assembly. The carboxyl-terminal pilin domain of FimH has an immunoglobulin-like fold, except that the seventh strand is missing, leaving part of the hydrophobic core exposed. A donor strand complementation mechanism in which the chaperone donates a strand to complete the pilin domain explains the basis for both chaperone function and pilus biogenesis.

A new set of peptide-based group heat capacities for use in protein stability calculations

The partial molar heat capacities of the tripeptides of the sequence glycyl-X-glycine, where X is one of the amino acids leucine, threonine, glutamine, phenylalanine, histidine, cysteine, proline, glutamic acid or arginine, and of the two tetrapeptides tetraglycine and glycyltryptophanylglycylglycine in aqueous solution over the temperature range 10-100 degrees C have been determined using high sensitivity scanning microcalorimetry. These results were used to derive the partial molar heat capacities of the various amino acid side-chains. This report completes our programme to derive reliable side-chain heat capacities for all 20 amino acids of proteins over a wide temperature range using the tripeptides Gly-X-Gly as realistic model compounds. Included in the study is a summary of the partial molar heat capacities of all 20 amino acid side-chains. These results, along with the heat capacity of the peptide backbone group, were used to calculate the partial molar heat capacities of some oligopeptides and of the random coil form of some unfolded proteins in water. The calculated heat capacities of the proteins obtained using this new set of heat capacities for the constituent groups are consistent with the heat capacities of the denatured state determined experimentally.

Distinguishing between sequential and nonsequentially folded proteins: implications for folding and misfolding

We describe here an algorithm for distinguishing sequential from nonsequentially folding proteins. Several experiments have recently suggested that most of the proteins that are synthesized in the eukaryotic cell may fold sequentially. This proposed folding mechanism in vivo is particularly advantageous to the organism. In the absence of chaperones, the probability that a sequentially folding protein will misfold is reduced significantly. The problem we address here is devising a procedure that would differentiate between the two types of folding patterns. Footprints of sequential folding may be found in structures where consecutive fragments of the chain interact with each other. In such cases, the folding complexity may be viewed as being lower. On the other hand, higher folding complexity suggests that at least a portion of the polypeptide backbone folds back upon itself to form three-dimensional (3D) interactions with noncontiguous portion(s) of the chain. Hence, we look at the mechanism of folding of the molecule via analysis of its complexity, that is, through the 3D interactions formed by contiguous segments on the polypeptide chain. To computationally splice the structure into consecutively interacting fragments, we either cut it into compact hydrophobic folding units or into a set of hypothetical, transient, highly populated, contiguous fragments ("building blocks" of the structure). In sequential folding, successive building blocks interact with each other from the amino to the carboxy terminus of the polypeptide chain. Consequently, the results of the parsing differentiate between sequentially vs. nonsequentially folded chains. The automated assessment of the folding complexity provides insight into both the likelihood of misfolding and the kinetic folding rate of the given protein. In terms of the funnel free energy landscape theory, a protein that truly follows the mechanism of sequential folding, in principle, encounters smoother free energy barriers. A simple sequentially folded protein should, therefore, be less error prone and fold faster than a protein with a complex folding pattern.

Modularity and homology: modelling of the type II module family from titin

We report the homology modelling of the structures of the 162 type II modules from the giant multi-domain protein titin (also known as connectin). The package MODELLER was used and implemented in an automated fashion using four experimentally determined structures as templates. Validation of the models was assessed in terms of divergence from the templates and consensus of the alignments. The homology within the whole family of type II modules as well as with the templates is relatively high (20-35% identity and ca 50% similarity). Comparison between the models of domains for which an NMR structure has been solved and the experimental solution gives an estimate of the quality of the modelling. Our results allow us to distinguish between a set of structurally relevant residues, which are conserved throughout the whole family and buried in the hydrophobic core, from the residues that are conserved and exposed. These latter residues are potentially functionally important. Comparison of exposed conserved patches for modules in different regions of the titin molecule suggests potential interaction surfaces. Our results may be tested directly for those modules whose binding partner is known.

Crystal structure of the vesicular transport protein Sec17: implications for SNAP function in SNARE complex disassembly

SNAP proteins play an essential role in membrane trafficking in eukaryotic cells. They activate and recycle SNARE proteins by serving as adaptors between SNAREs and the cytosolic chaperone NSF. We have determined the crystal structure of Sec17, the yeast homolog of alpha-SNAP, to 2.9 A resolution. Sec17 is composed of an N-terminal twisted sheet of alpha-helical hairpins and a C-terminal alpha-helical bundle. The N-terminal sheet has local similarity to the tetratricopeptide repeats from protein phosphatase 5 but has a different overall twist. Sec17 also shares structural features with HEAT and clathrin heavy chain repeats. Possible models of SNAP:SNARE binding suggest that SNAPs may function as lever arms, transmitting forces generated by conformational changes in NSF/Sec18 to drive disassembly of SNARE complexes.

Ultraviolet resonance Raman examination of horse apomyoglobin acid unfolding intermediates

We have used UV resonance Raman spectroscopy to study the acid-induced denaturation of horse apomyoglobin (apoMb) between pH 7. 0 and 1.8. The 206.5 nm excited Raman spectra are dominated by amide vibrations, which are used to quantitatively determine the apoMb secondary structure. The 229 nm excited Raman spectra are dominated by the Tyr and Trp Raman bands, which are analyzed to examine changes of Tyr and Trp environments and solvent exposures. We observe two partially unfolded apoMb intermediates at pH 4 and pH 2, while we observe only one partially unfolded holoMb intermediate at 2, in which the G and H helices are mainly intact, while the rest of protein is unfolded. This partially unfolded holoMb intermediate at pH 2 is essentially identical to the pH 2 apoMb intermediate. The partially unfolded pH 4 apoMb intermediate is composed of the three folded A, G, and H helices and contains 38% helical structure. The changes in the Trp Raman cross sections during the acid-induced denaturation indicates that Trp 7 is likely to be fully exposed in the apoMb pH 4 intermediate and that the A helix melts with a pKa approximately 3.5.

Position 127 amino acid substitutions affect the formation of CRP:cAMP:lacP complexes but not CRP:cAMP:RNA polymerase complexes at lacP

The lacP DNA binding and activation characteristics of CRP having amino acid substitutions at position 127 were investigated. Wild-type (WT) and T127C CRP footprinted lacP DNA in the presence of DNase I in a cAMP-dependent manner. The T127G, T127I, and T127S forms of CRP failed to footprint lacP both in the absence and in the presence of cAMP. Consistent with these data, WT and T127C CRP:cAMP complexes exhibited high affinity for the lacP CRP site whereas T127G, T127I, or T127S CRP:cAMP complexes exhibited low affinity for the lacP CRP site. CRP:cAMP:RNA polymerase (RNAP) complexes formed at lacP in reactions that contained WT, T127C, T127G, T127I, or T127S CRP. These results demonstrate that allosteric changes important for cAMP-mediated CRP activation are differentially affected by amino acid substitution at position 127. Proper cAMP-mediated reorientation of the DNA binding helices required either threonine or cysteine at position 127. However, cAMP-dependent interaction of CRP with RNAP was accomplished regardless of the amino acid at position 127. RNAP:lacP complexes that supported high-level lac RNA synthesis formed rapidly in reactions that contained WT or T127C CRP whereas RNAP:lacP complexes that supported only low-level lac RNA synthesis formed at slower rates in reactions that contained T127I or T127S CRP. The T127G CRP:cAMP:RNAP:lacP complex failed to activate lacP. The results of this study lead us to conclude that threonine 127 plays an important role in transduction of the signal from the CRP cyclic nucleotide binding pocket that promotes proper orientation of the DNA binding helices and only a minor, if any, role in the functional exposure of the CRP RNAP interaction domain.

Hydration-coupled dynamics in proteins studied by neutron scattering and NMR: the case of the typical EF-hand calcium-binding parvalbumin

The influence of hydration on the internal dynamics of a typical EF-hand calciprotein, parvalbumin, was investigated by incoherent quasi-elastic neutron scattering (IQNS) and solid-state 13C-NMR spectroscopy using the powdered protein at different hydration levels. Both approaches establish an increase in protein dynamics upon progressive hydration above a threshold that only corresponds to partial coverage of the protein surface by the water molecules. Selective motions are apparent by NMR in the 10-ns time scale at the level of the polar lysyl side chains (externally located), as well as of more internally located side chains (from Ala and Ile), whereas IQNS monitors diffusive motions of hydrogen atoms in the protein at time scales up to 20 ps. Hydration-induced dynamics at the level of the abundant lysyl residues mainly involve the ammonium extremity of the side chain, as shown by NMR. The combined results suggest that peripheral water-protein interactions influence the protein dynamics in a global manner. There is a progressive induction of mobility at increasing hydration from the periphery toward the protein interior. This study gives a microscopic view of the structural and dynamic events following the hydration of a globular protein.

Interaction of DnaK with native proteins and membrane proteins correlates with their accessible hydrophobicity

Molecular chaperones are involved in protein folding, protein targeting to membranes, and protein renaturation after stress. They interact specifically with hydrophobic sequences that are exposed in unfolded proteins, and buried in native proteins. We have studied the interaction of DnaK with native water-soluble proteins and membrane proteins. DnaK-native protein interactions are characterized by dissociation constants between 1 and 50 microM (compared with 0.01-1 microM for unfolded proteins). This affinity is within the range of most intracellular protein concentrations, suggesting that DnaK interacts with a greater number of native proteins than previously suspected. We found a correlation between the affinity of native proteins for DnaK and their affinity for hydrophobic-interaction chromatography adsorbents, suggesting that DnaK interacts with exposed hydrophobic groups in native proteins. The interaction between DnaK and membrane proteins is characterized by DnaK's high affinity for detergent-solubilized membrane proteins, and its lower affinity for membrane proteins inserted in lipid bilayers, suggesting that the chaperone can interact with the hydrophobic sequences of the former, while it cannot penetrate the hydrophobic core of lipid bilayers. Thus, the specificity of DnaK for hydrophobic sequences is involved in its interaction with not only unfolded proteins, but also native water-soluble proteins and membrane proteins. All proteins interact with DnaK according to their exposed hydrophobicity.

Residue 2 of TIMP-1 is a major determinant of affinity and specificity for matrix metalloproteinases but effects of substitutions do not correlate with those of the corresponding P1' residue of substrate

The unregulated activities of matrix metalloproteinases (MMPs) are implicated in disease processes including arthritis and tumor cell invasion and metastasis. MMP activities are controlled by four homologous endogenous protein inhibitors, tissue inhibitors of metalloproteinases (TIMPs), yet different TIMPs show little specificity for individual MMPs. The large interaction interface in the TIMP-1.MMP-3 complex includes a contiguous region of TIMP-1 around the disulfide bond between Cys1 and Cys70 that inserts into the active site of MMP-3. The effects of fifteen different substitutions for threonine 2 of this region reveal that this residue makes a large contribution to the stability of complexes with MMPs and has a dominant influence on the specificity for different MMPs. The size, charge, and hydrophobicity of residue 2 are key factors in the specificity of TIMP. Threonine 2 of TIMP-1 interacts with the S1' specificity pocket of MMP-3, which is a key to substrate specificity, but the structural requirements in TIMP-1 residue 2 for MMP binding differ greatly from those for the corresponding residue of a peptide substrate. These results demonstrate that TIMP variants with substitutions for Thr2 represent suitable starting points for generating more targeted TIMPs for investigation and for intervention in MMP-related diseases.

Exploring protein interiors: the role of a buried histidine in the KH module fold

The K-homology (KH) module is a novel RNA-binding motif. The structures of a representative KH motif from vigilin (vig-KH6) and of the first KH domain of fmr1 have been recently solved by nuclear magnetic resonance (NMR) and automated assignment-refinement techniques (ARIA). While a hydrophobic residue is found at position 21 in most of the KH modules, a buried His is conserved in all the 15 KH repeats of vigilin. This position must therefore have a key structural role in stabilizing the hydrophobic core. In the present work, we have addressed the following questions in order to obtain a detailed description of the role of His 21: i) what is the exact role of the histidine in the hydrophobic core of vig-KH6? ii) can we define the interactions that allow a conserved buried position to be occupied by a histidine both in vig-KH6 and in the whole vigilin KH sub-family? iii) how is the structure and stability of vig-KH6 influenced by the state of protonation of this histidine? To answer these questions, we have carried out an extensive refinement of the vig-KH6 structure using both an improved ARIA protocol starting from different initial structures and successively running restrained and unrestrained trajectories in water. An analysis of the stability of secondary structural elements, solvent accessibility, and hydrogen bonding patterns allows hypothesis on the structural role of residue His 21 and on the interactions that this residue forms with the environment. The importance of the protonation state of His 21 on the stability of the KH fold was addressed and validated by experimental results.

Probability analysis of variational crystallization and its application to gp120, the exterior envelope glycoprotein of type 1 human immunodeficiency virus (HIV-1)

The extensive glycosylation and conformational mobility of gp120, the envelope glycoprotein of type 1 human immunodeficiency virus (HIV-1), pose formidable barriers for crystallization. To surmount these difficulties, we used probability analysis to determine the most effective crystallization approach and derive equations which show that a strategy, which we term variational crystallization, substantially enhances the overall probability of crystallization for gp120. Variational crystallization focuses on protein modification as opposed to crystallization screening. Multiple variants of gp120 were analyzed with an iterative cycle involving a limited set of crystallization conditions and biochemical feedback on protease sensitivity, glycosylation status, and monoclonal antibody binding. Sources of likely conformational heterogeneity such as N-linked carbohydrates, flexible or mobile N and C termini, and variable internal loops were reduced or eliminated, and ligands such as CD4 and antigen-binding fragments (Fabs) of monoclonal antibodies were used to restrict conformational mobility as well as to alter the crystallization surface. Through successive cycles of manipulation involving 18 different variants, we succeeded in growing six different types of gp120 crystals. One of these, a ternary complex composed of gp120, its receptor CD4, and the Fab of the human neutralizing monoclonal antibody 17b, diffracts to a minimum Bragg spacing of at least 2.2 A and is suitable for structural analysis.

The atomic structure of protein-protein recognition sites

The non-covalent assembly of proteins that fold separately is central to many biological processes, and differs from the permanent macromolecular assembly of protein subunits in oligomeric proteins. We performed an analysis of the atomic structure of the recognition sites seen in 75 protein-protein complexes of known three-dimensional structure: 24 protease-inhibitor, 19 antibody-antigen and 32 other complexes, including nine enzyme-inhibitor and 11 that are involved in signal transduction.The size of the recognition site is related to the conformational changes that occur upon association. Of the 75 complexes, 52 have "standard-size" interfaces in which the total area buried by the components in the recognition site is 1600 (+/-400) A2. In these complexes, association involves only small changes of conformation. Twenty complexes have "large" interfaces burying 2000 to 4660 A2, and large conformational changes are seen to occur in those cases where we can compare the structure of complexed and free components. The average interface has approximately the same non-polar character as the protein surface as a whole, and carries somewhat fewer charged groups. However, some interfaces are significantly more polar and others more non-polar than the average. Of the atoms that lose accessibility upon association, half make contacts across the interface and one-third become fully inaccessible to the solvent. In the latter case, the Voronoi volume was calculated and compared with that of atoms buried inside proteins. The ratio of the two volumes was 1.01 (+/-0.03) in all but 11 complexes, which shows that atoms buried at protein-protein interfaces are close-packed like the protein interior. This conclusion could be extended to the majority of interface atoms by including solvent positions determined in high-resolution X-ray structures in the calculation of Voronoi volumes. Thus, water molecules contribute to the close-packing of atoms that insure complementarity between the two protein surfaces, as well as providing polar interactions between the two proteins.

Tolerance of a protein to multiple polar-to-hydrophobic surface substitutions

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

Protein structure prediction and design

Proteins have a unique native conformation, which can be proven in many instances to be determined by the amino acid sequence alone. The folding problem, that is the understanding of how the amino acid sequence directs folding, is still unsolved, despite more than 30 years of effort. However, many new methods have appeared in the past few years. This chapter describes the different principles underlying them and tries to give an overview of their successes and pitfalls.

Knowledge-based structure prediction of MHC class I bound peptides: a study of 23 complexes

BACKGROUND

The binding of T-cell antigenic peptides to MHC molecules is a prerequisite for their immunogenicity. The ability to identify binding peptides based on the protein sequence is of great importance to the rational design of peptide vaccines. As the requirements for peptide binding cannot be fully explained by the peptide sequence per se, structural considerations should be taken into account and are expected to improve predictive algorithms. The first step in such an algorithm requires accurate and fast modeling of the peptide structure in the MHC-binding groove.

RESULTS

We have used 23 solved peptide-MHC class I complexes as a source of structural information in the development of a modeling algorithm. The peptide backbones and MHC structures were used as the templates for prediction. Sidechain conformations were built based on a rotamer library, using the 'dead end elimination' approach. A simple energy function selects the favorable combination of rotamers for a given sequence. It further selects the correct backbone structure from a limited library. The influence of different parameters on the prediction quality was assessed. With a specific rotamer library that incorporates information from the peptide sidechains in the solved complexes, the algorithm correctly identifies 85% (92%) of all (buried) sidechains and selects the correct backbones. Under cross-validation, 70% (78%) of all (buried) residues are correctly predicted and most of all backbones. The interaction between peptide sidechains has a negligible effect on the prediction quality.

CONCLUSIONS

The structure of the peptide sidechains follows from the interactions with the MHC and the peptide backbone, as the prediction is hardly influenced by sidechain interactions. The proposed methodology was able to select the correct backbone from a limited set. The impairment in performance under cross-validation suggests that, currently, the specific rotamer library is not satisfactorily representative. The predictions might improve with an increase in the data.

Self-consistent estimation of inter-residue protein contact energies based on an equilibrium mixture approximation of residues

Pairwise contact energies for 20 types of residues are estimated self-consistently from the actual observed frequencies of contacts with regression coefficients that are obtained by comparing "input" and predicted values with the Bethe approximation for the equilibrium mixtures of residues interacting. This is premised on the fact that correlations between the "input" and the predicted values are sufficiently high although the regression coefficients themselves can depend to some extent on protein structures as well as interaction strengths. Residue coordination numbers are optimized to obtain the best correlation between "input" and predicted values for the partition energies. The contact energies self-consistently estimated this way indicate that the partition energies predicted with the Bethe approximation should be reduced by a factor of about 0.3 and the intrinsic pairwise energies by a factor of about 0.6. The observed distribution of contacts can be approximated with a small relative error of only about 0.08 as an equilibrium mixture of residues, if many proteins were employed to collect more than 20,000 contacts. Including repulsive packing interactions and secondary structure interactions further reduces the relative errors. These new contact energies are demonstrated by threading to have improved their ability to discriminate native structures from other non-native folds.

Identification of regions in which positive selection may operate in S-RNase of Rosaceae: implication for S-allele-specific recognition sites in S-RNase

A stylar S-RNase is associated with gametophytic self-incompatibility in the Rosaceae, Solanaceae, and Scrophulariaceae. This S-RNase is responsible for S-allele-specific recognition in the self-incompatible reaction, but how it functions in specific discrimination is not clear. Window analysis of the numbers of synonymous (dS) and non-synonymous (dN) substitutions in rosaceous S-RNases detected four regions with an excess of dN over dS in which positive selection may operate (PS regions). The topology of the secondary structure of the S-RNases predicted by the PHD method is very similar to that of fungal RNase Rh whose tertiary structure is known. When the sequences of S-RNases are aligned with the sequence of RNase Rh based on the predicted secondary structures, the four PS regions correspond to two surface sites on the tertiary structure of RNase Rh. These findings suggest that in S-RNases the PS regions also form two sites and are candidates for the recognition sites for S-allele-specific discrimination.

Surface properties of adipocyte lipid-binding protein: Response to lipid binding, and comparison with homologous proteins

Adipocyte lipid-binding protein (ALBP) is one of a family of intracellular lipid-binding proteins (iLBPs) that bind fatty acids, retinoids, and other hydrophobic ligands. The different members of this family exhibit a highly conserved three-dimensional structure; and where structures have been determined both with (holo) and without (apo) bound lipid, observed conformational changes are extremely small (Banaszak, et al., 1994, Adv. Prot. Chem. 45, 89; Bernlohr, et al., 1997, Annu. Rev. Nutr. 17, 277). We have examined the electrostatic, hydrophobic, and water accessible surfaces of ALBP in the apo form and of holo forms with a variety of bound ligands. These calculations reveal a number of previously unrecognized changes between apo and holo ALBP, including: 1) an increase in the overall protein surface area when ligand binds, 2) expansion of the binding cavity when ligand is bound, 3) clustering of individual residue exposure increases in the area surrounding the proposed ligand entry portal, and 4) ligand-binding dependent variation in the topology of the electrostatic potential in the area surrounding the ligand entry portal. These focused analyses of the crystallographic structures thus reveal a number of subtle but consistent conformational and surface changes that might serve as markers for differential targeting of protein-lipid complexes within the cell. Most changes are consistent from ligand to ligand, however there are some ligand-specific changes. Comparable calculations with intestinal fatty-acid-binding protein and other vertebrate iLBPs show differences in the electrostatic topology, hydrophobic topology, and in localized changes in solvent exposure near the ligand entry portal. These results provide a basis toward understanding the functional and mechanistic differences among these highly structurally homologous proteins. Further, they suggest that iLBPs from different tissues exhibit one of two predominant end-state structural distributions of the ligand entry portal.

Protein engineering the surface of enzymes

The protein surface is the interface through which a protein senses the external world. Its composition of charged, polar and hydrophobic residues is crucial for the stability and activity of the protein. The charge state of seven of the twenty naturally occurring amino acids is pH dependent. A total of 95% of all titratable residues are located on the surface of soluble proteins. In evolutionary related families of proteins such residues are particularly prone to substitutions, insertions and deletions. We present here an analysis of the residue composition of 4038 proteins, selected from 125 protein families with < 25% identity between core members of each family. Whereas only 16.8% of the residues were truly buried, 40.7% were > 30% exposed on the surface and the remainder were < 30% exposed. The individual residue types show distinct differences. The data presented provides an important new approach to protein engineering of protein surfaces. Guidelines for the optimization of solvent exposure for a given residue are given. The cutinase family of enzymes has been investigated. The stability of native cutinase has been studied as a function of pH, and has been compared with the cutinase activity towards tributyrin. Whereas the onset of enzymatic activity is linked with the deprotonation of the active site HIS188, destabilization of the 3D structure as determined by differential scanning calorimetry is coupled with the loss of activity at very basic pH values. A modeling investigation of the pH dependence of the electrostatic potentials reveals that the activity range is accompanied by the development of a highly significant negative potential in the active site cleft. The 3D structures of three mutants of the Fusarium solani pisi cutinase have been solved to high resolution using X-ray diffraction analysis. Preliminary X-ray data are presented.

Protein sidechain conformer prediction: a test of the energy function

BACKGROUND

Homology modeling is an important technique for making use of the rapidly increasing number of protein sequences in the absence of structural information. The major problems in such modeling, once the alignment has been made, concern the positions of loops and the orientations of sidechains. Although progress has been made in recent years for sidechain prediction, current methods appear to have a limit on the order of 70% in their accuracy. It is important to have an understanding of this limitation, which for energy-based methods could arise from inaccuracies of the potential function.

RESULTS

A test of the CHARMM function for sidechain prediction was performed. To eliminate the multiple-residue search problem, the minimum energy positions of individual sidechains in ten proteins were calculated in the presence of all other sidechains in their crystal orientations. This test provides a necessary condition that any energy function useful for sidechain placement must satisfy. For chi1 x chi2 rotations, the accuracies were 77.4% and 89.5%, respectively, and in the presence of crystal waters were 86.5% and 94.9%, respectively. If there was an error, the crystal structure usually corresponded to an alternative local minimum on the calculated energy map. Prediction accuracy correlated with the size of the energy gap between primary and secondary minima.

CONCLUSIONS

The results indicate that the errors in current sidechain prediction schemes cannot be attributed to the potential energy function per se. The test used here establishes a necessary condition that any proposed energy-based sidechain prediction method, as well as many statistically based methods, must satisfy.

The novel substrate recognition mechanism utilized by aspartate aminotransferase of the extreme thermophile Thermus thermophilus HB8

Aspartate aminotransferase (AspAT) is a unique enzyme that can react with two types of substrate with quite different properties, acidic substrates, such as aspartate and glutamate, and neutral substrates, although the catalytic group Lys-258 acts on both types of substrate. The dynamic properties of the substrate-binding site are indispensable to the interaction with hydrophobic substrates (Kawaguchi, S., Nobe, Y., Yasuoka, J., Wakamiya, T., Kusumoto, S., and Kuramitsu, S. (1997) J. Biochem. (Tokyo) 122, 55-63). AspATs from various organisms are classified into two subgroups, Ia and Ib. The former includes AspATs from Escherichia coli and higher eukaryotes, whereas the latter includes those from Thermus thermophilus and many prokaryotes. The AspATs belonging to subgroup Ia each have an Arg-292 residue, which interacts with the distal carboxyl groups of dicarboxylic (acidic) substrates, but the functionally similar residue of subgroup Ib AspATs has not been identified. In view of the x-ray crystallographic structure of T. thermophilus AspAT, we expected Lys-109 to be this residue in the subgroup Ib AspATs and constructed K109V and K109S mutants. Replacing Lys-109 with Val or Ser resulted in loss of activity toward acidic substrates but increased that toward the neutral substrate, alanine, considerably. These results indicate that Lys-109 is a major determinant of the acidic substrate specificity of subgroup Ib AspATs. Kinetic analysis of the interactions with neutral substrates indicated that T. thermophilus AspAT is subject to less steric hindrance and its substrate-binding pocket has a more flexible conformation than E. coli AspAT. A flexible active site in the rigid T. thermophilus AspAT molecule may explain its high activity even at room temperature.

Protein thermostability in extremophiles

Thermostability of a protein is a property which cannot be attributed to the presence of a particular amino acid or to a post synthetic modification. Thermostability seems to be a property acquired by a protein through many small structural modifications obtained with the exchange of some amino acids and the modulation of the canonical forces found in all proteins such as electrostatic (hydrogen bonds and ion-pairs) and hydrophobic interactions. Proteins produced by thermo and hyperthermophilic microorganisms, growing between 45 and 110 degrees C are in general more resistant to thermal and chemical denaturation than their mesophilic counterparts. The observed structural resistance may reflect a restriction on the flexibility of these proteins, which, while allowing them to be functionally competent at elevated temperatures, renders them unusually rigid at mesophilic temperatures (10-45 degrees C). The increased rigidity at mesophilic temperatures may find a structural determinant in increased compactness. In thermophilic proteins a number of amino acids are often exchanged. These exchanges with some strategic placement of proline in beta-turns give rise to a stabilization of the protein. Mutagenesis experiments have confirmed this statement. From the comparative analysis of the X-ray structures available for several families of proteins, including at least one thermophilic structure in each case, it appears that thermal stabilization is accompanied by an increase in hydrogen bonds and salt bridges. Thermostability appears also related to a better packing within buried regions. Despite these generalisations, no universal rules can be found in these proteins to achieve thermostability.

Structural determinants of cytochrome P450 substrate specificity, binding affinity and catalytic rate

The structural characteristics of cytochrome P450 substrates are summarised, showing that molecular descriptors can discriminate between chemicals of differing P450 isozyme specificity. Procedures for the estimation of P450 substrate binding interaction energies and rates of metabolism are described, providing specific examples in both individual compounds binding to P450s, including those of known crystal structure, and within series of structurally related chemicals. It is demonstrated that binding energy components are primarily hydrophobic/desolvation and electrostatic/hydrogen-bonded in nature, whereas electronic factors are of importance in determining variations in reaction rates. It is thus shown that the prediction of P450 substrate binding affinities and catalytic rates may be feasible, provided that sufficient structural information is available for the relevant enzyme-substrate complex.

Thermodynamic analysis of biomolecules: a volumetric approach

Fundamental thermodynamic relationships reveal that volumetric studies on molecules of interest can yield useful new information. In particular, appropriately designed volumetric studies can characterize the properties of molecules as a function of solution conditions, including the role of solvation. Until recently, such studies on biologically interesting molecules have been limited because of the lack of readily available instrumentation with the requisite sensitivity; however, during the past decade, advances in the development of highly sensitive, small-volume densimetric, acoustic and high-pressure spectroscopic instrumentation have enabled biological molecules to be subjected to a wide range of volumetric studies. In fact, the volumetric methods used in these studies have already provided unique insights into the molecular origins of the intramolecular and intermolecular recognition events that modulate biomolecular processes. Of particular note are recent volumetric studies on globular proteins and nucleic acid duplexes. These studies have provided unique insights into the role of hydration in modulating the stabilities of these biopolymers, as well as their conformational transitions and ligand-binding properties.

Proteins from hyperthermophiles: stability and enzymatic catalysis close to the boiling point of water

It has become clear since about a decade ago, that the biosphere contains a variety of microorganisms that can live and grow in extreme environments. Hyperthermophilic microorganisms, present among Archaea and Bacteria, proliferate at temperatures of around 80-100 degrees C. The majority of the genera known to date are of marine origin, however, some of them have been found in continental hot springs and solfataric fields. Metabolic processes and specific biological functions of these organisms are mediated by enzymes and proteins that function optimally under these extreme conditions. We are now only starting to understand the structural, thermodynamic and kinetic basis for function and stability under conditions of high temperature, salt and extremes of pH. Insights gained from the study of such macromolecules help to extend our understanding of protein biochemistry and -biophysics and are becoming increasingly important for the investigation of fundamental problems in structure biology such as protein stability and protein folding. Extreme conditions in the biosphere require either the adaptation of the amino acid sequence of a protein by mutations, the optimization of weak interactions within the protein and at the protein-solvent boundary, the influence of extrinsic factors such as metabolites, cofactors, compatible solutes. Furthermore folding catalysts, known as chaperones, that assist the folding of proteins may be involved or increased protein protein synthesis in order to compensate for destruction by extreme conditions. The comparison of structure and stability of homologous proteins from mesophiles and hyperthermophiles has revealed important determinants of thermal stability of proteins. Rather than being the consequence of one dominant type of interactions or of a general stabilization strategy, it appears that the adaptation to high temperatures reflects a number of subtle interactions, often characteristic for each protein species, that minimize the surface energy and the hydration of apolar surface groups while burying hydrophobic residues and maximizing packing of the core as well as the energy due to charge-charge interactions and hydrogen bonds. In this article, mechanisms of intrinsic stabilization of proteins are reviewed. These mechanisms are found on different levels of structural organization. Among the extrinsic stabilization factors, emphasis is put on archaea chaperonins and their still strongly debated function. It will be shown, that optimization of weak protein-protein and protein-solvent interactions plays a key role in gaining thermostability. The difficulties in correlating suitable optimization criteria with real thermodynamic stability measures are due to experimental difficulties in measuring stabilization energies in large proteins or protein oligomers and will be discussed. Thus small single domain proteins or isolated domains of larger proteins may serve as model systems for large or multidomain proteins which due to the complexity of their thermal unfolding transitions cannot be analyzed by equilibrium thermodynamics. The analysis of the energetics of the thermal unfolding of a small, hyperthermostable DNA binding protein from Sulfolobus has revealed that a high melting temperature is not synonymous with a larger maximum thermodynamic stability. Finally, it is now well documented, that many thermophilic and hyperthermophilic proteins show a statistically increased number of salt bridges and salt bridge networks. However their contribution to thermodynamic and functional stability is still obscure.

Hydrophobic core substitutions in calbindin D9k: effects on stability and structure

The effects of hydrophobic core mutations on the stability and structure of the four-helix calcium-binding protein, calbindin D9k, have been investigated. Eleven mutations involving eight residues distributed within the hydrophobic core of calbindin D9k were examined. Stabilities were measured by denaturant and thermal induced unfolding monitored by circular dichroism spectroscopy. The mutations were found to exert large effects on the stability with midpoints in the urea induced unfolding varying from 1.8 M for Leu23 --> Gly up to 6.6 M for Val70 --> Leu and free energies of unfolding in the absence of denaturant ranging from 6.6 to 27.4 kJ/mol for the Phe66 --> Ala mutant and the wild-type, respectively. A significant correlation was found between the difference in free energy of unfolding (Delta Delta GNU) and the change in the surface area of the side chain caused by the mutation, in agreement with other studies. Notably, both increases and decreases in side-chain surface area caused quantitatively equivalent effects on the stability. In other words, a correlation between the absolute value of the change in the surface of the side chain and Delta DeltaGNU was observed with a value of approximately 0.14 kJ M-1 A-2. The generality of this observation is discussed. Significant effects on the cooperativity of the unfolding reaction were also observed. However, a correlation between the cooperativity and Delta Delta GNU, which has been reported in other systems as an indication of effects of mutations on the unfolded state, was not observed for calbindin D9k. Despite the large effects on Delta Delta GNU and cooperativity, the structures of the mutants in the native form remained intact as indicated by circular dichroism, NMR, and fluorescence measurements. The structural response to calcium-binding was also conserved. The following paper in this issue [Kragelund, B. B., et al. (1998) Biochemistry 37, 8926-8937] examines the effects of these mutations on the calcium binding properties of calbindin D9k.

CD2 and the nature of protein interactions mediating cell-cell recognition

Rapid progress has recently been made in characterising the structures of leukocyte cell-surface molecules. Detailed analyses of the structure and interactions of CD2 were the first involving a molecule that has not been directly linked to antigen recognition in the manner of antigen receptors or co-receptors. It seems highly likely that the properties of ligand binding by CD2 are relevant to the general mechanisms of cell-cell recognition. As an example of biological recognition, the defining characteristic of cell-cell contact is that it involves the simultaneous interaction of hundreds, if not thousands, of molecules. Affinity and kinetic analyses of ligand binding by CD2 indicated that the protein interactions mediating cell-cell contact, whilst highly specific, are much weaker than initially anticipated, probably due to the requirement that such contacts be easily reversible. Simultaneously, in addressing the mechanism of this mode of recognition, structural and mutational studies focussed on the role of charged residues clustered in the ligand-binding face of CD2, yielding the concept that electrostatic complementarity, rather than surface-shape complementarity, is the dominant feature of specific, low-affinity protein recognition at the cell surface by CD2. The crystallographic analysis of the CD2-binding domain of CD58 strongly supports this concept.

A search for single substitutions that eliminate enzymatic function in a bacterial ribonuclease

Exhaustive-substitution studies, where many amino acid replacements are individually tested at all positions in a natural protein, have proven to be very valuable in probing the relationship between sequence and function. The broad picture that has emerged from studies of this sort is one of functional tolerance of substitution. We have applied this approach to barnase, a 110-residue bacterial ribonuclease. Because the selection system used to score barnase mutants as active or inactive detects activity down to a level that can be approached by nonenzyme catalysts, mutants that test inactive are essentially devoid of enzymatic function. Of the 109 barnase positions subjected to substitution, only 15 (14%) are vulnerable to this extreme level of inactivation, and only 2 could not be substituted without such inactivation. A total of 33 substitutions (amounting to 5% of the explored substitutions) were found to render barnase wholly inactive. The profoundly disruptive effects of all of these inactivating substitutions appear to result from either (1) replacement of a side chain that is directly involved in substrate binding or catalysis, (2) replacement of a substantially buried side chain, (3) introduction of a proline residue, or (4) replacement of a glycine residue. Although substitutions of these types are functionally tolerated more often than not, the system used here indicates that only these sorts of substitution are capable of single-handedly reducing catalytic function to, or nearly to, levels that can be achieved by nonenzyme catalysts.

The role of charged residues mediating low affinity protein-protein recognition at the cell surface by CD2

Insights into the structural basis of protein-protein recognition have come principally from the analysis of proteins such as antibodies, hormone receptors, and proteases that bind their ligands with relatively high affinity (Ka approximately 10(9) M-1). In contrast, few studies have been done on the very low affinity interactions mediating cell adhesion and cell-cell recognition. As a site of protein-protein recognition, the ligand binding face of the T lymphocyte cell-cell recognition molecule, CD2, which binds its ligands 10(4)- to 10(5)-fold more weakly than do antibodies and proteases, is unusual in being both very flat and highly charged. An analysis of the effect of mutations and ionic strength on CD2 binding to its ligand, CD48, indicates that these charged residues contribute little, if any, binding energy to this interaction. However, the loss of these charged residues is shown to markedly reduce ligand-binding specificity. Thus, the charged residues increase the specificity of CD2 binding without increasing the affinity. This phenomenon is likely to result from a requirement for electrostatic complementarity between charged binding surfaces to compensate for the removal, upon binding, of water interacting with the charged residues. It is proposed that this mode of recognition is highly suited to biological interactions requiring a low affinity because it uncouples increases in specificity from increases in affinity.

Protein hydration and unfolding--insights from experimental partial specific volumes and unfolded protein models

BACKGROUND

The partial specific volume of a protein is an experimental quantity containing information about solute-solvent interactions and protein hydration. We use a hydration-shell model to partition the partial specific volume into an intrinsic volume occupied by the protein and a change in the volume occupied by the solvent resulting from the solvent interactions with the protein. We seek to extract microscopic information about protein hydration and unfolding from experimental volume measurements without using computer simulations. We employ the idea that the protein-solvent interaction will be proportional to the surface area of the protein.

RESULTS

A linear relationship is obtained when the difference between the experimental protein partial specific volume and its intrinsic volume is plotted as a function of the protein solvent-accessible surface area. The effect of using different protein volume definitions on the analysis of protein volumetric properties is discussed. Volumetric data are used to test a model for the unfolded state of proteins and to make predictions about the denatured state.

CONCLUSIONS

The linear relationship between hydration-shell volume change and accessible surface area reflects the similar surface properties (fractional composition of nonpolar, polar and charged surface) among a diverse set of proteins. This linear relationship is found to be independent of how the solution is partitioned into solute and solvent components. The interpretation of hydration shell versus bulk water properties is found to be very model dependent, however. The maximally exposed unfolded protein model is found to be inconsistent with experimental volume changes of unfolding.