Identification of the barstar binding site of barnase by NMR spectroscopy and hydrogen-deuterium exchange

Design for Solubility May Reveal Induction of Amide Hydrogen/Deuterium Exchange by Protein Self-Association

Structural heterogeneity often constrains the characterization of aggregating proteins to indirect or low-resolution methods, obscuring mechanistic details of association. Here, we report progress in understanding the aggregation of Adnectins, engineered binding proteins with an immunoglobulin-like fold. We rationally design Adnectin solubility and measure amide hydrogen/deuterium exchange (HDX) under conditions that permit transient protein self-association. Protein-protein binding commonly slows rates of HDX; in contrast, we find that Adnectin association may induce faster HDX for certain amides, particularly in the C-terminal β-strand. In aggregation-prone proteins, we identify a pattern of very different rates of amide HDX for residues linked by reciprocal hydrogen bonds in the native structure. These results may be explained by local loss of native structure and formation of an inter-protein interface. Amide HDX induced by self-association, detected here by deliberate modulation of propensity for such interactions, may be a general phenomenon with the potential to expose mechanisms of aggregation by diverse proteins.

Barnase-Barstar Pair: Contemporary Application in Cancer Research and Nanotechnology

Barnase is an extracellular ribonuclease secreted by Bacillus amyloliquefaciens that was originally studied as a small stable enzyme with robust folding. The identification of barnase intracellular inhibitor barstar led to the discovery of an incredibly strong protein-protein interaction. Together, barnase and barstar provide a fully genetically encoded toxin-antitoxin pair having an extremely low dissociation constant. Moreover, compared to other dimerization systems, the barnase-barstar module provides the exact one-to-one ratio of the complex components and possesses high stability of each component in a complex and high solubility in aqueous solutions without self-aggregation. The unique properties of barnase and barstar allow the application of this pair for the engineering of different variants of targeted anticancer compounds and cytotoxic supramolecular complexes. Using barnase in suicide gene therapy has also found its niche in anticancer therapy. The application of barnase and barstar in contemporary experimental cancer therapy is reflected in the review.

CID fragmentation, H/D exchange and supermetallization of Barnase-Barstar complex

The barnase-barstar complex is one of the most stable protein-protein complexes and has a very wide range of possible applications. Here we report the use of top-down mass spectrometry for the investigation of the structure of this complex, its ionization via ESI, isolation and fragmentation. It was found that the asymmetry of the resulting charge state distributions of the protein monomer product ions increased as the charge state of the precursor ions increased. For the investigation of the 3D structure of the complex, the gas phase H/D exchange reaction was used. In addition, supermetallized ions of the complex with Zn were produced and investigated. It was observed that an increase in the number of metals bound to the complex results in a change in complex stability and the charge distribution between protein fragment. Analysis of the fragmentation pattern of the supermetallized complex [bn-b* + 5Zn]¹⁰⁺ indicated that this ion is present in different conformations with different charges and Zn distributions. Since Zn cannot migrate, such structures must be formed during ionization.

Development of a Model Protein Interaction Pair as a Benchmarking Tool for the Quantitative Analysis of 2-Site Protein-Protein Interactions

A significant challenge in the molecular interaction field is to accurately determine the stoichiometry and stepwise binding affinity constants for macromolecules having >1 binding site. The mission of the Molecular Interactions Research Group (MIRG) of the Association of Biomolecular Resource Facilities (ABRF) is to show how biophysical technologies are used to quantitatively characterize molecular interactions, and to educate the ABRF members and scientific community on the utility and limitations of core technologies [such as biosensor, microcalorimetry, or analytic ultracentrifugation (AUC)]. In the present work, the MIRG has developed a robust model protein interaction pair consisting of a bivalent variant of the Bacillus amyloliquefaciens extracellular RNase barnase and a variant of its natural monovalent intracellular inhibitor protein barstar. It is demonstrated that this system can serve as a benchmarking tool for the quantitative analysis of 2-site protein-protein interactions. The protein interaction pair enables determination of precise binding constants for the barstar protein binding to 2 distinct sites on the bivalent barnase binding partner (termed binase), where the 2 binding sites were engineered to possess affinities that differed by 2 orders of magnitude. Multiple MIRG laboratories characterized the interaction using isothermal titration calorimetry (ITC), AUC, and surface plasmon resonance (SPR) methods to evaluate the feasibility of the system as a benchmarking model. Although general agreement was seen for the binding constants measured using solution-based ITC and AUC approaches, weaker affinity was seen for surface-based method SPR, with protein immobilization likely affecting affinity. An analysis of the results from multiple MIRG laboratories suggests that the bivalent barnase-barstar system is a suitable model for benchmarking new approaches for the quantitative characterization of complex biomolecular interactions.

VASP-E: specificity annotation with a volumetric analysis of electrostatic isopotentials

Algorithms for comparing protein structure are frequently used for function annotation. By searching for subtle similarities among very different proteins, these algorithms can identify remote homologs with similar biological functions. In contrast, few comparison algorithms focus on specificity annotation, where the identification of subtle differences among very similar proteins can assist in finding small structural variations that create differences in binding specificity. Few specificity annotation methods consider electrostatic fields, which play a critical role in molecular recognition. To fill this gap, this paper describes VASP-E (Volumetric Analysis of Surface Properties with Electrostatics), a novel volumetric comparison tool based on the electrostatic comparison of protein-ligand and protein-protein binding sites. VASP-E exploits the central observation that three dimensional solids can be used to fully represent and compare both electrostatic isopotentials and molecular surfaces. With this integrated representation, VASP-E is able to dissect the electrostatic environments of protein-ligand and protein-protein binding interfaces, identifying individual amino acids that have an electrostatic influence on binding specificity. VASP-E was used to examine a nonredundant subset of the serine and cysteine proteases as well as the barnase-barstar and Rap1a-raf complexes. Based on amino acids established by various experimental studies to have an electrostatic influence on binding specificity, VASP-E identified electrostatically influential amino acids with 100% precision and 83.3% recall. We also show that VASP-E can accurately classify closely related ligand binding cavities into groups with different binding preferences. These results suggest that VASP-E should prove a useful tool for the characterization of specific binding and the engineering of binding preferences in proteins.

Proteins, the ubiquitous worker molecules of the cell, are a diverse class of molecules that perform very specific tasks. Understanding how proteins achieve specificity is a critical step towards understanding biological systems and a key prerequisite for rationally engineering new proteins. To examine electrostatic influences on specificity in proteins, this paper presents VASP-E, a software tool that generates solid representations of the electrostatic potential fields that surround proteins. VASP-E compares solids with constructive solid geometry, a class of techniques developed first for modeling complex machine parts. We observed that solid representations could quantify the degree of charge complementarity in protein-protein interactions and identify key residues that strengthen or weaken them. VASP-E correctly identified amino acids with established experimental influences on protein-protein binding specificity. We also observed that solid representations of electrostatic fields could identify electrostatic conservations and variations that relate to similarities and differences in binding specificity between proteins and small molecules.

Molecular simulations of mutually exclusive folding in a two-domain protein switch

A major challenge with testing designs of protein conformational switches is the need for experimental probes that can independently monitor their individual protein domains. One way to circumvent this issue is to use a molecular simulation approach in which each domain can be directly observed. Here we report what we believe to be the first molecular simulations of mutually exclusive folding in an engineered two-domain protein switch, providing a direct view of how folding of one protein drives unfolding of the other in a barnase-ubiquitin fusion protein. These simulations successfully capture the experimental effects of interdomain linker length and ligand binding on the extent of unfolding in the less stable domain. In addition, the effect of linker length on the potential for oligomerization, which eliminates switch activity, is in qualitative agreement with analytical ultracentrifugation experiments. We also perform what we believe to be the first study of protein unfolding via progressive localized compression. Finally, we are able to explore the kinetics of mutually exclusive folding by determining the effect of linker length on rates of unfolding and refolding of each protein domain. Our results demonstrate that molecular simulations can provide seemingly novel biological insights on the behavior of individual protein domains, thereby aiding in the rational design of bifunctional switches.

The Kinetic Pathway of Folding of Barnase

To search for folding intermediates, we have examined the folding and unfolding kinetics of wild-type barnase and four representative mutants under a wide range of conditions that span two-state and multi-state kinetics. The choice of mutants and conditions provided in-built controls for artifacts that might distort the interpretation of kinetics, such as the non-linearity of kinetic and equilibrium data with concentration of denaturant. We measured unfolding rate constants over a complete range of denaturant concentration by using by 1H/2H-exchange kinetics under conditions that favour folding, conventional stopped-flow methods at higher denaturant concentrations and continuous flow. Under conditions that favour multi-state kinetics, plots of the rate constants for unfolding against denaturant concentration fitted quantitatively to the equation for three-state kinetics, with a sigmoid component for a change of rate determining step, as did the refolding kinetics. The position of the transition state on the reaction pathway, as measured by solvent exposure (the Tanford beta value) also moved with denaturant concentration, fitting quantitatively to the same equations with a change of rate determining step. The sigmoid behaviour disappeared under conditions that favoured two-state kinetics. Those data combined with direct structural observations and simulation support a minimal reaction pathway for the folding of barnase that involves two detectable folding intermediates. The first intermediate, I(1), is the denatured state under physiological conditions, D(Phys), which has native-like topology, is lower in energy than the random-flight denatured state U and is suggested by molecular dynamics simulation of unfolding to be on-pathway. The second intermediate, I(2), is high energy, and is proven by the change in rate determining step in the unfolding kinetics to be on-pathway. The change in rate determining step in unfolding with structure or environment reflects the change in partitioning of this intermediate to products or starting materials.

Lack of definable nucleation sites in the rate-limiting transition state of barnase under native conditions

It has been shown that the burst-phase (submillisecond) intermediate of barnase, if it exists, can be only marginally more stable than the fully unfolded state at pH 6.3 and 25 degrees C. In the study reported here, no stable burst-phase intermediate could be detected, even in the presence of stabilizing salt (0.4 M Na(2)SO(4)). These results suggest that a burst-phase intermediate with even marginal stability does not exist. The absence of such an intermediate in turn suggests the need for re-examination of the rate-limiting transition state (RLTS) under native conditions, which was previously characterized by using a three-state model with a stable intermediate and protein engineering. Surprisingly, mutations throughout the structure of barnase do not significantly affect the folding rate, suggesting a lack of specific favorable interactions among the side-chains in the RLTS. This RLTS is clearly different from that previously characterized under denaturing conditions, indicating that changes take place in the RLTS under native and denaturing conditions. The occurrence of such changes is further supported by the observation that the unfolding rate constants of barnase and its mutants were divergent or convergent as a function of denaturant concentrations. Consistent with changes in the RLTS, a re-analysis of data from native-state hydrogen exchange studies has shown that the logarithm of the unfolding rate constant inflects down under low concentrations of denaturant. Here, we discuss in detail the question of whether changes in the RLTS involve a kinetically silent intermediate that occurs after the initial RLTS.

A novel approach for assessing macromolecular complexes combining soft-docking calculations with NMR data

We present a novel and efficient approach for assessing protein-protein complex formation, which combines ab initio docking calculations performed with the protein docking algorithm BiGGER and chemical shift perturbation data collected with heteronuclear single quantum coherence (HSQC) or TROSY nuclear magnetic resonance (NMR) spectroscopy. This method, termed "restrained soft-docking," is validated for several known protein complexes. These data demonstrate that restrained soft-docking extends the size limitations of NMR spectroscopy and provides an alternative method for investigating macromolecular protein complexes that requires less experimental time, effort, and resources. The potential utility of this novel NMR and simulated docking approach in current structural genomic initiatives is discussed.

Can non-mechanical proteins withstand force? Stretching barnase by atomic force microscopy and molecular dynamics simulation

Atomic force microscopy (AFM) experiments have provided intriguing insights into the mechanical unfolding of proteins such as titin I27 from muscle, but will the same be possible for proteins that are not physiologically required to resist force? We report the results of AFM experiments on the forced unfolding of barnase in a chimeric construct with I27. Both modules are independently folded and stable in this construct and have the same thermodynamic and kinetic properties as the isolated proteins. I27 can be identified in the AFM traces based on its previous characterization, and distinct, irregular low-force peaks are observed for barnase. Molecular dynamics simulations of barnase unfolding also show that it unfolds at lower forces than proteins with mechanical function. The unfolding pathway involves the unraveling of the protein from the termini, with much more native-like secondary and tertiary structure being retained in the transition state than is observed in simulations of thermal unfolding or experimentally, using chemical denaturant. Our results suggest that proteins that are not selected for tensile strength may not resist force in the same way as those that are, and that proteins with similar unfolding rates in solution need not have comparable unfolding properties under force.

Calculation of z-coordinates and orientational restraints using a metal binding tag

We introduce a new simple methodology allowing the measurement of (1)H-(15)N residual dipolar couplings, dipolar shifts, and unpaired electron-amide proton distances. This method utilizes a zinc finger tag fused at either the N- or the C-terminus of a protein. We have demonstrated this fusion strategy by incorporating the zinc finger of the retroviral gag protein onto the C-terminus of barnase, a ribonuclease produced by Bacillus amiloliquifaciance. We show that this tag can be substituted with cobalt and manganese. Binding of cobalt to the gag zinc finger-barnase fusion protein introduced sufficient anisotropic paramagnetic susceptibility for orientation of the molecule in the magnetic field. Partial alignment permitted measurement of (1)J(HN) scalar couplings along with dipolar couplings. Replacement of bound cobalt with diamagnetic zinc removes the paramagnetic-induced orientation of barnase, permitting the measurement of only (1)J(HN) scalar couplings. Dipolar couplings, ranging from -0.9 to 0.6 Hz, were easily measured from the difference in splitting frequencies in the presence of cobalt and zinc. The observed paramagnetic anisotropy induced by cobalt binding to the metal binding tag also permitted measurement of dipolar shifts. Substitution of manganese into the metal binding tag permitted the measurement of unpaired electron-amide proton distances using paramagnetic relaxation enhancement methodology. The availability of both amide proton dipolar shifts and unpaired electron to amide proton distances permitted the direct calculation of z-coordinates for individual amide protons. This approach is robust and will prove powerful for global fold determination of proteins identified in genome initiatives.

Backbone dynamics of barstar: a (15)N NMR relaxation study

Backbone dynamics of uniformly (15)N-labeled barstar have been studied at 32 degrees C, pH 6.7, by using (15)N relaxation data obtained from proton-detected 2D (1)H-(15)N NMR spectroscopy. (15)N spin-lattice relaxation rate constants (R(1)), spin-spin relaxation rate constants (R(2)), and steady-state heteronuclear (1)H-(15)N NOEs have been determined for 69 of the 86 (excluding two prolines and the N-terminal residue) backbone amide (15)N at a magnetic field strength of 14.1 Tesla. The primary relaxation data have been analyzed by using the model-free formalism of molecular dynamics, using both isotropic and axially symmetric diffusion of the molecule, to determine the overall rotational correlation time (tau(m)), the generalized order parameter (S(2)), the effective correlation time for internal motions (tau(e)), and NH exchange broadening contributions (R(ex)) for each residue. As per the axially symmetric diffusion, the ratio of diffusion rates about the unique and perpendicular axes (D( parallel)/D( perpendicular)) is 0.82 +/- 0.03. The two results have only marginal differences. The relaxation data have also been used to map reduced spectral densities for the NH vectors of these residues at three frequencies: 0, omega(H), and omega(N), where omega(H),(N) are proton and nitrogen Larmor frequencies. The value of tau(m) obtained from model-free analysis of the relaxation data is 5.2 ns. The reduced spectral density analysis, however, yields a value of 5.7 ns. The tau(m) determined here is different from that calculated previously from time-resolved fluorescence data (4.1 ns). The order parameter ranges from 0.68 to 0.98, with an average value of 0.85 +/- 0.02. A comparison of the order parameters with the X-ray B-factors for the backbone nitrogens of wild-type barstar does not show any considerable correlation. Model-free analysis of the relaxation data for seven residues required the inclusion of an exchange broadening term, the magnitude of which ranges from 2 to 9.1 s(-1), indicating the presence of conformational averaging motions only for a small subset of residues.

Absence of stable intermediates on the folding pathway of barnase

Barnase is one of the few protein models that has been studied extensively for protein folding. Previous studies led to the conclusion that barnase folds through a very stable submillisecond intermediate ( approximately 3 kcal/mol). The structure of this intermediate was characterized intensively by using a protein engineering approach. This intermediate has now been reexamined with three direct and independent methods. (i) Hydrogen exchange experiments show very small protection factors ( approximately 2) for the putative intermediate, indicating a stability of approximately 0.0 kcal/mol. (ii) Denaturant-dependent unfolding of the putative intermediate is noncooperative and indicates a stability less than 0.0 kcal/mol. (iii) The logarithm of the unfolding rate constant of native barnase vs. denaturant concentrations is not linear. Together with the measured rate ("I" to N), this nonlinear behavior accounts for almost all of the protein stability, leaving only about 0.3 kcal/mol that could be attributed to the rapidly formed intermediate. Other observations previously interpreted to support the presence of an intermediate are now known to have alternative explanations. These results cast doubts on the previous conclusions on the nature of the early folding state in barnase and therefore should have important implications in understanding the early folding events of barnase and other proteins in general.

Protein global fold determination using site-directed spin and isotope labeling

We describe a simple experimental approach for the rapid determination of protein global folds. This strategy utilizes site-directed spin labeling (SDSL) in combination with isotope enrichment to determine long-range distance restraints between amide protons and the unpaired electron of a nitroxide spin label using the paramagnetic effect on relaxation rates. The precision and accuracy of calculating a protein global fold from only paramagnetic effects have been demonstrated on barnase, a well-characterized protein. Two monocysteine derivatives of barnase, (H102C) and (H102A/Q15C), were 15N enriched, and the paramagnetic nitroxide spin label, MTSSL, attached to the single Cys residue of each. Measurement of amide 1H longitudinal relaxation times, in both the oxidized and reduced states, allowed the determination of the paramagnetic contribution to the relaxation processes. Correlation times were obtained from the frequency dependence of these relaxation processes at 800, 600, and 500 MHz. Distances in the range of 8 to 35 A were calculated from the magnitude of the paramagnetic contribution to the relaxation processes and individual amide 1H correlation times. Distance restraints from the nitroxide spin to amide protons were used as restraints in structure calculations. Using nitroxide to amide 1H distances as long-range restraints and known secondary structure restraints, barnase global folds were calculated having backbone RMSDs <3 A from the crystal structure. This approach makes it possible to rapidly obtain the overall topology of a protein using a limited number of paramagnetic distance restraints.

Amide proton hydrogen exchange rates for sperm whale myoglobin obtained from 15N-1H NMR spectra

The hydrogen exchange behavior of exchangeable protons in proteins can provide important information for understanding the principles of protein structure and function. The positions and exchange rates of the slowly-exchanging amide protons in sperm whale myoglobin have been mapped using 15N-1H NMR spectroscopy. The slowest-exchanging amide protons are those that are hydrogen bonded in the longest helices, including members of the B, E, and H helices. Significant protection factors were observed also in the A, C, and G helices, and for a few residues in the D and F helices. Knowledge of the identity of slowly-exchanging amide protons forms the basis for the extensive quench-flow kinetic folding experiments that have been performed for myoglobin, and gives insights into the tertiary interactions and dynamics in the protein.

Relationship between the native-state hydrogen exchange and the folding pathways of barnase

The previous native-state hydrogen exchange experiment with barnase failed to detect any partially unfolded intermediate state which was contrary to the experimental results from kinetic deuterium hydrogen exchange pulse labeling and protein engineering studies. This has been taken to suggest that the native-state hydrogen exchange method cannot be used alone as an analytical tool to study the folding pathways of proteins. Here, we revisited the pulse labeling experiment with barnase and detected no stable folding intermediate. This finding allows a reconciliation of the native-state HX data and the folding pathway of barnase. Along with alternative theoretical interpretations for a curved chevron plot of protein folding, these data suggest that further investigation of the nature of the intermediate of barnase is needed.

An NMR analysis of ubiquitin recognition by yeast ubiquitin hydrolase: evidence for novel substrate recognition by a cysteine protease

Yeast ubiquitin hydrolase 1 (YUH1), a cysteine protease that catalyzes the removal of ubiquitin C-terminal adducts, is important for the generation of monomeric ubiquitin. Heteronuclear NMR spectroscopy has been utilized to map the YUH1 binding surface on ubiquitin. When YUH1 was titrated into a sample of ubiquitin, approximately 50% of the (1)H-(15)N correlation peaks of ubiquitin were affected to some degree, as a result of binding to YUH1. It is noteworthy that the amide resonances of the basic residues (Arg42, Lys48, Arg72, and Lys74) were highly perturbed. These positively charged basic residues may be involved in direct interactions with the negatively charged acidic residues on YUH1. In addition to the electrostatic surface, the hydrophobic surfaces on ubiquitin (Leu8, Ile44, Phe45, Val70, Leu71, and Leu73) and YUH1 are also likely to contribute to the binding interaction. Furthermore, the amide resonances of Ile13, Leu43, Leu50, and Leu69, the side chains of which are not on the surface, were also highly perturbed, indicating substrate-induced changes in the environments of these residues as well. These large changes, observed from residues located throughout the five-stranded beta-sheet surface and the C-terminus, suggest that substrate recognition by YUH1 involves a wider area on ubiquitin.

A mini-protein designed by removing a module from barnase: molecular modeling and NMR measurements of the conformation

A globular domain can be decomposed into compact modules consisting of contiguous 10-30 amino acid residues. The correlation between modules and exons observed in different proteins suggests that each module was encoded by an ancestral exon and that modules were combined into globular domains by exon fusion. Barnase is a single domain RNase consisting of 110 amino acid residues and was decomposed into six modules. We designed a mini-protein by removing the second module, M2, from barnase in order to gain an insight into the structural and functional roles of the module. In the molecular modeling of the mini-protein, we evaluated thermodynamic stability and aqueous solubility together with mechanical stability of the model. We chemically synthesized a mini-barnase with (15)N-labeling at 10 residues, whose corresponding residues in barnase are all found in the region around the hydrophobic core. Circular dichroism and NMR measurements revealed that mini-barnase takes a non-random specific conformation that has a similar hydrophobic core structure to that of barnase. This result, that a module could be deleted without altering the structure of core region of barnase, supports the view that modules act as the building blocks of protein design.

Superslow backbone protein dynamics as studied by 1D solid-state MAS exchange NMR spectroscopy

Superslow backbone dynamics of the protein barstar and the polypeptide polyglycine was studied by means of a solid-state MAS 1D exchange NMR method (time-reverse ODESSA) that can detect reorientation of nuclei carrying anisotropic chemical shift tensors. Experiments were performed on carbonyl 13C in polyglycine (natural abundance) and backbone 15N nuclei in uniformly 15N-enriched barstar within a wide range of temperatures in dry and wet powders for both samples. Two exchange processes were observed in the experiments: molecular reorientation and spin diffusion. Experimental conditions that are necessary to separate these two processes are discussed on a quantitative level. It was revealed that the wet protein undergoes molecular motion in the millisecond range of correlation times, whereas in dry protein and polyglycine molecular reorientations could not be detected. The correlation time of the motion in the wet barstar at room temperature is 50-100 ms; the activation energy is about 80 kJ/mol. Previously, protein motions with such a long correlation time could be observed only by methods detecting chemical exchange in solution (e.g., hydrogen exchange). The application of solid-state MAS exchange spectroscopy provides new opportunities in studying slow biomolecular dynamics that is important for the biological function of proteins.

Three-dimensional structure of binase in solution

We present the spatial structure of binase, a small extracellular ribonuclease, derived from 1H-NMR* data in aqueous solution. The total of 20 structures were obtained via torsion angle dynamics using DYANA program with experimental NOE and hydrogen bond distance constraints and phi and chi1 dihedral angle constraints. The final structures were energy minimised with ECEPP/2 potential in FANTOM program. Binase consists of three alpha-helices in N-terminal part (residues 6-16, 26-32 and 41-44), five-stranded antiparallel beta-sheet in C-terminal part (residues 51-55, 70-75, 86-90, 94-99 and 104-108) and two-stranded parallel beta-sheet (residues 22-24 and 49-51). Three loops (residues 36-39, 56-67 and 76-83), which play significant role in biological functioning of binase, are flexible in solution. The differences between binase and barnase spatial structures in solution explain the differences in thermostability of binase, barnase and their hybrids.

N epsilon,N epsilon-dimethyl-lysine cytochrome c as an NMR probe for lysine involvement in protein-protein complex formation

The reductively dimethylated derivatives of horse and yeast iso-1-ferricytochromes c have been prepared and characterized for use as NMR probes of the complexes formed by cytochrome c with bovine liver cytochrome b5 and yeast cytochrome c peroxidase. The electrostatic properties and structures of the derivatized cytochromes are not significantly perturbed by the modifications; neither are the electrostatics of protein-protein complex formation or rates of interprotein electron transfer. Two-dimensional 1H-13C NMR spectroscopy of the complexes formed by the derivatized cytochromes with cytochrome b5 and cytochrome c peroxidase has been used to investigate the number and identity of lysine residues of cytochrome c that are involved in interprotein interactions of cytochrome c. The NMR data are incompatible with simple static models proposed previously for the complexes formed by these proteins, but are consistent with the presence of multiple, interconverting complexes of comparable stability, consistent with studies employing Brownian dynamics to model the complexes. The NMR characteristics of the Nepsilon,Nepsilon-dimethyl-lysine groups, their chemical shift dispersion, oxidation state and temperature dependences and the possibility of chemical exchange phenomena are discussed with relevance to the utility of Nepsilon, Nepsilon-dimethyl-lysine's being a generally useful derivative for characterizing protein-protein complexes.

Discrepancies between the NMR and X-ray structures of uncomplexed barstar: analysis suggests that packing densities of protein structures determined by NMR are unreliable

The crystal structure of the C82A mutant of barstar, the intracellular inhibitor of the Bacillus amyloliquefaciens ribonuclease barnase, has been solved to a resolution of 2.8 A. The molecule crystallizes in the space group I41 with a dimer in the asymmetric unit. An identical barstar dimer is also found in the crystal structure of the barnase-barstar complex. This structure of uncomplexed barstar is compared to the structure of barstar bound to barnase and also to the structure of barstar solved using NMR. The free structure is similar to the bound state, and there are no significant main-chain differences in the 27-44 region involved in barstar binding to barnase. The C82A structure shows significant differences from the average NMR structure, both overall and in the binding region. In contrast to the crystal structure, the NMR structure shows an unusually high packing value based on the occluded surface algorithm, indicating errors in the packing of the structure. We show that the NMR structures of homologous proteins generally show large differences in packing value, while the crystal structures of such proteins have very similar packing values, suggesting that protein packing density is not well determined by NMR.

A new method to characterize hydrophobic organization of proteins: application to rational protein engineering of barnase

We present a new algorithm for characterization of protein spatial structure basing on the molecular hydrophobicity potential approach. The method is illustrated by the analysis of three-dimensional structure of barnase and barnase-barstar complex. Current approach enables identification of amino acid residues situated in unfavorable environment (these residues may be "active" for binding), and to map quantitatively hydrophobic, hydrophilic and unfavorable hydrophobic-hydrophilic intra- and inter-molecular contacts involving backbone and side-chain segments of amino acid residues. Calculation of individual contributions of amino acid residues to such contacts permits identification of structurally-important residues. The contact plots obtained with molecular hydrophobicity potential calculations, provide easy rules to choose sites for mutations, which can increase a strength of intra- or inter-molecular hydrophobic interactions. The unfavorable hydrophobic-hydrophilic contact can be mutated to favorable hydrophobic, and already existing weak hydrophobic contact can be strengthen by increasing hydrophobicity of residues in contact. Basing on the analysis of the contact plots, we suggest several mutations of barnase which are supposed to increase intramolecular hydrophobic interactions, and thus might lead to increased stability of the protein. Part of these mutations was studied previously experimentally, and indeed stabilized barnase. The other of predicted mutations were not studied experimentally yet. Several new mutations of barnase and barstar are also proposed to enhance the hydrophobic interactions on their binding interface.

Hydrogen exchange and protein folding

Amide hydrogen-deuterium exchange is a sensitive probe of the structure, stability and dynamics of proteins. The significant increase in the number of small, model proteins that have been studied has allowed a better understanding of the structural fluctuations that lead to hydrogen exchange. Recent technical advances enable the methodology to be applied to the study of protein-protein interactions in much larger, more complex systems.

A phage display technique for a fast, sensitive, and systematic investigation of protein-protein interactions

Phage display is a technique in which a foreign protein or peptide is presented at the surface of a (filamentous) bacteriophage. This system, developed by Smith [(1985), Science 228, 1315-1317], was originally used to create large libraries of antibodies for the purpose of selecting those that strongly bound a particular antigen. More recently it was also employed to present peptides, domains of proteins, or intact proteins at the surface of phages, again to identify high-affinity interactions with ligands. Here we want to illustrate the use of phage display, in combination with PCR saturation mutagenesis, for the study of protein-protein interactions. Rather than selecting for mutants having high affinity, we systematically investigate the binding of every variant with its natural ligand. Via a modified ELISA we can calculate a relative affinity. As a model system we chose to display thymosin beta 4 on the phage surface in order to study its interaction with actin.

Identification of critical residues in the colicin E9 DNase binding region of the Im9 protein

1H-15N NMR studies, in conjunction with mutagenesis experiments, have been used to delineate the DNase-binding surface of the colicin E9 inhibitor protein Im9 (where Im stands for immunity protein). Complexes were formed between the 15 kDa unlabelled E9 DNase domain and the 9.5 kDa Im9 protein uniformly labelled with 15N. Approx. 90% of the amide resonances of the bound Im9 were assigned and spectral parameters obtained from 1H-15N heteronuclear single quantum coherence (HSQC) spectra were compared with those for the free Im9 assigned previously. Many of the amide resonances were shifted on complex formation, some by more than 2 p.p.m. in the 15N dimension and more than 0.5 p.p.m. in the 1H dimension. Most of the strongly shifted amides are located on the surfaces of two of the four helices, helix II and helix III. Whereas helix II had already been identified through genetic and biochemical investigations as an important determinant of biological specificity, helix III had not previously been implicated in binding to the DNase. To test the robustness of the NMR-delineated DNase-binding site, a selection of Im9 alanine mutants were constructed and their dissociation rate constants from E9 DNase-immunity protein complexes quantified by radioactive subunit exchange kinetics. Their off-rates correlated well with the NMR perturbation analysis; for example, residues that were highly perturbed in HSQC experiments, such as residues 34 (helix II) and 54 (helix III), had a marked effect on the DNase-immunity protein dissociation rate when replaced by alanine. The NMR and mutagenesis data are consistent with a DNase-binding region on Im9 composed of invariant residues in helix III and variable residues in helix II. The relationship of this binding site model to the wide range of affinities (Kd values in the range 10(-4) to 10(-16)M) that have been measured for cognate and non-cognate colicin DNase-immunity protein interactions is discussed.

Hydrogen exchange: the modern legacy of Linderstrøm-Lang

This discussion, prepared for the Protein Society's symposium honoring the 100th anniversary of Kaj Linderstrøm-Lang, shows how hydrogen exchange approaches initially conceived and implemented by Lang and his colleagues some 50 years ago are contributing to current progress in structural biology. Examples are chosen from the active protein folding field. Hydrogen exchange methods now make it possible to define the structure of protein folding intermediates in various contexts: as tenuous molten globule forms at equilibrium under destabilizing conditions, in kinetic intermediates that exist for less than one second, and as infinitesimally populated excited state forms under native conditions. More generally, similar methods now find broad application in studies of protein structure, energetics, and interactions. This article considers the rise of these capabilities from their inception at the Carlsberg Labs to their contemporary role as a significant tool of modern structural biology.

Protein complexes studied by NMR spectroscopy

Recent advances in NMR methods now allow protein complexes to be studied in great detail in a wide range of solution conditions. Isotope-enrichment strategies, resonance-assignment approaches and structural-determination methods have evolved to the point where almost any type of complex involving proteins of reasonable size may be studied in a straightforward way. A variety of isotope editing and filtering strategies underlie these powerful methodologies. Approaches to the characterization of the dynamics of protein complexes have also matured to the point where detailed studies of the effects of complexation on dynamics can be studied over a wide range of timescales.

Mechanisms and uses of hydrogen exchange

Recent work has largely completed our understanding of the hydrogen-exchange chemistry of unstructured proteins and nucleic acids. Some of the high-energy structural fluctuations that determine the hydrogen-exchange behavior of native macromolecules have been explained; others remain elusive. A growing number of applications are exploiting hydrogen-exchange behavior to study difficult molecular systems and elicit otherwise inaccessible information on protein structure, dynamics and energetics.

A calorimetric study of the thermal stability of barstar and its interaction with barnase

The temperature-induced unfolding of single, double, and triple mutants of barstar, the specific intracellular protein inhibitor of barnase from Bacillus amyloliquefaciens, has been studied by high-sensitivity differential scanning calorimetry. The thermal unfolding of barstar mutants, where at least one of the two cysteine residues in the molecule had been replaced by alanine, follows a two-state mechanism at neutral and alkaline pH. The unfolding enthalpy and heat capacity changes are slightly lower than those accepted for highly compact, small, globular proteins. We have found that at pH 2.5, where barstar seems to be in a molten globule state, the protein has a heat capacity between that of the native and the unfolded states and shows some tendency for association. Scanning calorimetry experiments were also extended to the barstar--barnase complex in the neutral and alkaline pH range. The binding constants obtained from DSC studies are similar to those already obtained from other (kinetic) studies. The interaction of barstar and barnase was also investigated by isothermal calorimetry in various buffers within the pH range 6.0-10.0 and a temperature range of 15-35 degrees C. The favorable enthalpy contribution to the binding is about 4 times higher than the entropic one at 25 degrees C. The overall data analysis of the combined calorimetric results has led to the thermodynamic characterization of barstar unfolding and the interaction of barstar and barnase over a wide range of temperatures.

Toward solving the folding pathway of barnase: the complete backbone 13C, 15N, and 1H NMR assignments of its pH-denatured state

The structures of the major folding intermediate, the transition state for folding, and the folded state of barnase have been previously characterized. We now add a further step toward a complete picture of the folding of barnase by reporting the backbone 15N, 13C, and 1H NMR assignments for barnase unfolded at pH 1.8 and 30 degrees C. These assignments, which were obtained from a combination of heteronuclear magnetization transfer and backbone triple-resonance NMR experiments, constitute the first stage in the structural characterization of this denatured state by NMR. Interresidue nuclear Overhauser effect contacts and deviations from 1H random-coil chemical shifts provide evidence for stable residual structure. The structured regions span residues in the native protein that contain its major alpha-helix and central strands of the beta-sheet. Earlier experiments have shown that these regions are predominantly intact in the major folding intermediate and that their docking is partly rate determining in folding.