Models for the complexes formed between cytochrome b5 and the subunits of methemoglobin

Computer graphics-generated models for the electron transfer complexes formed between cytochrome b5 and the subunits of methemoglobin are proposed. For both complexes, the orientation allowing optimal hydrogen bonding involves interaction between negatively charged residues on cytochrome b5 and positively charged residues on methemoglobin. In each complex, the heme groups of the interacting species are coplanar with the edges of the heme groups separated by 7-8 A and with the iron atoms 16 A apart. For the alpha-chain X cytochrome b5 complex, alpha-chain residues 56 (Lys), 60 (Lys), and 90 (Lys) interact with cytochrome b5 residues 44 (Glu), 43 (Glu), and 60 (Asp) respectively. A fourth hydrogen bond involves alpha-61 (Lys) bridging between a heme propionate from cytochrome b5 and a heme propionate from the alpha-chain. The contacts present in the beta-chain X cytochrome b5 complex involve hydrogen-bonding between beta-chain lysyl residues 59, 61, 65, and 95, and cytochrome b5 residues 48 (Glu), 44 (Glu), 43 (Glu), and 60 (Asp) respectively. An additional hydrogen bond can be formed by bridging of the epsilon-amino group of beta-66 (Lys) between a heme propionate from cytochrome b5 and a beta-chain heme propionate. In each complex, two nonionic interactions, one on each side of the heme groups, are also suggested. These interactions appear to effectively exclude external water molecules from the center of the protein-protein interaction domain. Comparison of the proposed binding loci for cytochrome b5 on the methemoglobin subunits with those proposed on cytochrome c reveals considerable structural homology between the cytochrome b5 binding sites.     

Characterization of the EGF-like module pair 3-4 from vitamin K-dependent protein S using NMR spectroscopy reveals dynamics on three separate time scales and extensive effects from calcium binding

Protein S, a cofactor of anticoagulant activated protein C, exhibits three high-affinity Ca(2+)-binding sites in a region comprising four EGF modules. The EGF 3-4 module pair constitutes the smallest fragment that retains one high-affinity Ca(2+)-binding site and is therefore useful for investigation of the structural basis of the unusually high-affinity Ca(2+) binding compared to other EGF-containing proteins characterized so far. Extensive chemical shift effects caused by Ca(2+) binding to the EGF 3-4 module pair are observed, particularly from Ca(2+) binding to the high-affinity site in EGF 4. Ca(2+) binding to the high-affinity site in EGF 4 and the low-affinity site in EGF 3 is associated with slow and fast exchange on the NMR time-scale, respectively. We show the presence of two isoforms, characterized by a cis or trans Lys 167-Pro 168 peptide bond, that do not convert on time scales that were accessible to the experiments (k(ex) < 0.2 s(-1)). Both conformers have similar Ca(2+) affinities and backbone dynamics. Further, broadening of (1)H resonances involving residues in the major beta-sheet of EGF 3 and (15)N exchange terms, primarily in the N-terminal part of the protein, indicate the presence of slow exchange on a microsecond to millisecond time scale. (15)N spin relaxation data suggest that the module pair has a well-defined relative orientation between EGF modules 3 and 4 and has a significantly anisotropic rotational diffusion tensor in solution.     

Mutating protein kinase cAMP-binding sites into cGMP-binding sites. Mechanism of cGMP selectivity.

The cAMP-dependent protein kinase contains two different cAMP-binding sites referred to as the slow and fast sites. Mutation of Ala-334 to a threonine in the slow site of the bovine type I regulatory subunit created a site with marked increase in cGMP affinity without changing cAMP affinity (Shabb, J. B., Ng. L., Corbin, J. D. (1990) J. Biol. Chem. 265, 16031-16034). The corresponding fast site residue (Ala-210) was changed to a threonine by oligonucleotide-directed mutagenesis, and a double mutant containing a threonine in each site was also made. Holoenzymes were formed from native catalytic subunit and each recombinant regulatory subunit. The fast site mutant holoenzyme exhibited an improved cGMP activation constant and an impaired cAMP activation constant. The double mutant cGMP/cAMP selectivity was 200-fold greater than that of wild-type holoenzyme, making it as responsive to cGMP as native cGMP-dependent protein kinase. The increased intrinsic binding energies of mutated sites for cGMP were 2.7-3.0 kcal mol-1, consistent with the presence of an extra hydrogen bond. Cyclic nucleotide analog studies implied that this hydrogen bond was between the threonine hydroxyl and the 2-amino of cGMP. Comparisons of amino acid sequences and cyclic nucleotide specificities suggested that the Ala/Thr difference may also impart cAMP/cGMP binding selectivity to related proteins such as cyclic nucleotide-gated ion channels.     

Crystal structures of aspartate carbamoyltransferase ligated with phosphonoacetamide, malonate, and CTP or ATP at 2.8-A resolution and neutral pH

The R-state structures of the ATP and CTP complexes of aspartate carbamoyltransferase ligated with phosphonoacetamide and malonate have been determined at 2.8-A resolution and neutral pH. These structures were solved by the method of molecular replacement and were refined to crystallographic residuals between 0.167 and 0.182. The triphosphate, the ribose, and the purine and pyrimidine moieties of ATP and CTP interact with similar regions of the allosteric domain of the regulatory dimer. ATP and CTP relatively increase and decrease the size of the allosteric site in the vicinity of the base, respectively. For both CTP and ATP at pH 7, the gamma-phosphates are bound to His20 and are also near Lys94, while the alpha-phosphates interact exclusively with Lys94. The 2'-hydroxyls of both CTP and ATP are near the amino group of Lys60. The pyrimidine ring of CTP makes specific hydrogen bonds at the allosteric site: the NH2 group donates hydrogen bonds to the main-chain carbonyls of Ile12 and Tyr89 and the pyrimidine ring carbonyl oxygen accepts a hydrogen bond from the amino group of Lys60; the nitrogen at position 3 in the pyrimidine ring is hydrogen bonded to a main-chain NH group of Ile12. The purine ring of ATP also makes numerous interactions with residues at the allosteric site: the purine NH2 (analogous to the amino group of CTP) donates a hydrogen bond to the main-chain carbonyl oxygen of Ile12, the N3 nitrogen interacts with the amino group of Lys60, and the N1 nitrogen hydrogen bonds to the NH group of Ile12. The binding of CTP and ATP to the allosteric site in the presence of phosphonoacetamide and malonate does not dramatically alter the structure of the allosteric binding site or of the allosteric domain. Nonetheless, in the CTP-ligated structure, the average separation between the catalytic trimers decreases by approximately 0.5 A, indicating a small shift of the quaternary structure toward the T state. In the CTP- and ATP-ligated R-state structures, the binding and occupancy of phosphonoacetamide and malonate are similar and the structures of the active sites are similar at the current resolution of 2.8 A.     

Pathways of cytochrome c oxidation by soluble and membrane-bound cytochrome aa3

In media of low ionic strength, membraneous cytochrome c oxidase, isolated cytochrome c oxidase, and proteoliposomal cytochrome c oxidase each bind cytochrome c at two sites, one of low affinity (1 microM greater than Kd' greater than 0.2 microM) and readily reversible and the other of high affinity (0.01 microM greater than Kd) and weakly reversible. When cytochrome c occupies both sites, including the low affinity site, the maximal turnover measured polarographically with ascorbate and N,N,N',N'-tetramethyl-p-phenylenediamine (TMPD) is independent of TMPD concentration, and lies between 250 and 400 s-1 (30 degrees C, pH 7.4) for fully activated systems. The apparent affinity of the enzyme for cytochrome c is, however, TMPD dependent. When cytochrome c occupies only the high-affinity site, the maximal turnover is closely dependent upon the concentration of TMPD, which, unlike ascorbate, can reduce bound cytochrome c. As TMPD concentration is increased, the maximal turnover approaches that seen when both sites as occupied. The lower activity of isolated cytochrome aa3 is due to the presence of inactive or inaccessible enzyme molecules. Incorporation of isolated enzyme into phospholipid vesicles restores full activity to all the subsequently accessible cytochrome aa3 molecules. Negatively charged (asolectin) vesicles show a higher cytochrome c affinity at the low-affinity sites than do the other enzyme preparations. A model for the cytochrome c-cytochrome aa3 complexes is put forward in which both sites, when occupied, are fully catalytically competent, but in which occupation of the "tight" site by a catalytically functional cytochrome c molecule is required for overall oxidation of cytochrome c via the "loose" site.     

Nucleotide binding of an ADP analog to cooperating sites of chloroplast F1-ATPase (CF1).

Pre-steady state nucleotide binding to the chloroplast F1-ATPase (CF1) was measured in a stopped-flow apparatus by monitoring the change of fluorescence intensity of TNP-ADP upon binding. The analysis of the time courses led to the proposal of a mechanism of nucleotide binding with the following characteristics. (a) It involves three sites binding nucleotides in a concerted manner. (b) Each binding site is able to undergo a conformational change from a loose binding state into a state refraining from any direct release of the bound nucleotide into the medium. Only the reverse reaction via the loose binding state enables release out of the tight binding state. (c) Due to very strong negative cooperativity, a maximum of two of the three sites can be found in the state of tight binding. (d) Cooperativity between the three sites leads to a slower nucleotide binding of the second nucleotide compared to the first nucleotide. Furthermore, the conformational change from the loose binding state to the tight binding state is slowed down if one of the other sites already is in the tight binding state. Three-sites mechanisms in which rotation leads to an exchange of the properties of the binding sites failed to simulate the observed time courses of nucleotide binding. However, as the experimental set up was designed to prevent catalysis taking place, our results entirely agree with the current finding that rotation requires catalytic turnover of the enzyme.     

Intermolecular aggregations are responsible for the slow kinetics observed in the folding of cytochrome c at neutral pH.

Folding of equine cytochrome c at a low protein concentration (26 microM) eliminated a slow kinetic phase (time constant three seconds) that was observed in the previous hydrogen exchange pulse-labeling experiments at pH 6.2 and 10 degrees C. It was demonstrated that this slow folding phase was caused by intermolecular aggregations. Because heterogeneous kinetics is a very general feature in the folding of proteins characterized by pulsed hydrogen exchange coupled with two-dimensional NMR, our experimental results suggest aggregations might also be responsible for the complex folding kinetics of other proteins. This is possible since these experiments were performed at relatively high protein concentrations.     

Desulfovibrio desulfuricans G20 tetraheme cytochrome structure at 1.5 Angstrom and cytochrome interaction with metal complexes

The structure of the type I tetraheme cytochrome c(3) from Desulfovibrio desulfuricans G20 was determined to 1.5 Angstrom by X-ray crystallography. In addition to the oxidized form, the structure of the molybdate-bound form of the protein was determined from oxidized crystals soaked in sodium molybdate. Only small structural shifts were obtained with metal binding, consistent with the remarkable structural stability of this protein. In vitro experiments with pure cytochrome showed that molybdate could oxidize the reduced cytochrome, although not as rapidly as U(VI) present as uranyl acetate. Alterations in the overall conformation and thermostability of the metal-oxidized protein were investigated by circular dichroism studies. Again, only small changes in protein structure were documented. The location of the molybdate ion near heme IV in the crystal structure suggested heme IV as the site of electron exit from the reduced cytochrome and implicated Lys14 and Lys56 in binding. Analysis of structurally conserved water molecules in type I cytochrome c(3) crystal structures identified interactions predicted to be important for protein stability and possibly for intramolecular electron transfer among heme molecules.     

Refinement of the structure of human basic fibroblast growth factor at 1.6 A resolution and analysis of presumed heparin binding sites by selenate substitution.

The three-dimensional structure of human basic fibroblast growth factor has been refined to a crystallographic residual of 16.1% at 1.6 A resolution. The structure has a Kunitz-type fold and is composed of 12 antiparallel beta-strands, 6 of which form a beta-barrel. One bound sulfate ion has been identified in the model, hydrogen bonded to the side chains of Asn 27, Arg 120, and Lys 125. The side chain of Arg 120 has two conformations, both of which permit hydrogen bonds to the sulfate. This sulfate binding site has been suggested as the binding site for heparin (Eriksson, A.E., Cousens, L.S., Weaver, L.H., & Matthews, B.W., 1991, Proc. Natl. Acad. Sci. USA 88, 3441-3445). Two beta-mercaptoethanol (BME) molecules are also included in the model, each forming a disulfide bond to the S gamma atoms of Cys 69 and Cys 92, respectively. The side chain of Cys 92 has two conformations of which only one can bind BME. Therefore the BME molecule is half occupied at this site. The locations of possible sulfate binding sites on the protein were examined by replacing the ammonium sulfate in the crystallization medium with ammonium selenate. Diffraction data were measured to 2.2 A resolution and the structure refined to an R-factor of 13.8%. The binding of the more electron-dense selenate ion was identified at two positions. One position was identical to the sulfate binding site identified previously. The second selenate binding site, which is of lower occupancy, is situated 5.6 A from the first. This ion is hydrogen bonded by the side chain of Lys 135 and Arg 120. Thus the side chain of Arg 120 binds two selenate ions simultaneously. It is suggested that the observed second selenate binding site should also be considered as a possible binding site for heparin, or that both selenate binding sites might simultaneously contribute to the binding of heparin.     

Identification of a Cd2+- and Zn2+-binding site in cytochrome c using FTIR coupled to an ATR microdialysis setup and NMR spectroscopy

Fourier transform infrared (FTIR) difference spectroscopy allows the study of molecular changes occurring at active sites in proteins with high sensitivity. Reactions are triggered by light, potential, or temperature steps and more recently by the diffusion of buffers containing effectors above membrane proteins deposited as films on ATR crystals. We have adapted a microdialysis system to an ATR, to study metal sites in soluble proteins. In this study, we identified a Cd(2+)- or Zn(2+)-binding site in cytochrome c with dissociation constants of 17 and 42 microM, respectively, which affects the oxidation rate of ferrocytochrome c by hydrogen peroxide. Using the microdialysis ATR-FTIR setup, we determined that a histidine and the carboxylate group of a glutamate are involved in Zn(2+) binding. The implication of His 33 and Glu 104 in the binding site was deduced from the comparison of FTIR data recorded with horse heart and the variant tuna cytochrome c lacking these two amino acids. A two-dimensional NMR analysis of the Zn(2+)-binding site in horse heart cytochrome c confirmed that His 33 and residues close to the C terminus are sensitive to Zn(2+) binding. This study demonstrates that the microdialysis ATR-FTIR setup is promising for the analysis of metal sites in proteins. From H(2)O/(2)H(2)O exchange experiments, we concluded that the impact of Zn(2+) and Cd(2+) binding on the oxidation kinetics of ferrocytochrome c by H(2)O(2) is associated to the perturbation of a hydrogen-bonding network involving His 33 that is sensitive to the redox state of cytochrome c.     

Catalytic effect of ferricyanide on the rate of electron transfer between myoglobin and cytochrome c]

The influence of small amounts of low-molecular electron acceptor, potassium ferricyanide, 1 to 20% relative to the cytohrome c concentration, on the rate of electron transfer in the sperm whale oxymyoglobin--horse heart cytochrome c and deoxymyoglobin--cytochrome c systems (under aerobic and anaerobic conditions, respectively) was studied. At low ionic strength, the redox reaction rate was found to increase proportionally to the concentration of ferricyanide in both redox systems. The effect depends on pH in the pH range 5-8, increasing sharply at pH < 6. It was shown that the enhancing of electron transfer is caused by the complexing of [Fe(CN)6]3- with cytohrome c in the Lys72 region, where one of the two strong binding sites for this anion is determined by NMR. Both the high ionic strength and the chemical modification of Lys72 residue inhibit this effect at low ionic strength, markedly decreasing the rate of reaction with myoglobin. Under the same conditions, the effect of ferricyanide in the reaction of oxy-Mb with yeast cytohrome c, which is isopotential to animal cytochromes c but possesses trimethylated Lys72, was several times smaller. In turn, the chemical modification of His residues in myoglobin and the complexing of zinc ion to His119(GH1) almost completely inhibit electron transfer in the systems. Thus, electron transfer between the proteins must proceed through the formation of the Mb.[Fe(CN)6]3-.Cyt c ternary complex, the contacting sites being localized in the His119(GH1) region of myoglobin and near Lys72 of cytohrome c. The increased electron transfer rate in the presence of [Fe(CN)6]3- can be explained by that its binding near Lys72, firstly, provides better electrostatic interactions in the electron transfer complex and, besides, decreases significantly (about 2-fold) the tunneling distance between the two hemes (two lengths of 1.7 and 1.2 nm instead of one of 2.9 nm).     

Structures of HO(2)-Co(III)bleomycin A(2) bound to d(GAGCTC)(2) and d(GGAAGCTTCC)(2): structure-reactivity relationships of Co and Fe bleomycins

HO(2)-Co(III)bleomycin is a model for HO(2)-Fe(III)bleomycin, which initiates single and double strand cleavage of DNA. In order to enlarge the understanding of its structure and reactivity, three-dimensional structures of HO(2)-Co(III)bleomycin bound to two DNA oligomers, d(GAGCTC)(2) (I) and d(GGAAGCTTCC)(2) (II), that have 5'-GC-3' binding sites, have been determined by nuclear magnetic resonance (NMR) methods. Besides previously recognized determinants of binding selectivity, a probable hydrogen bond was detected between the pyrimidinyl acetamido NH(2) and the carbonyl of cytosine base paired to G at the recognition site. Another hydrogen bond between the NH of the dimethylsulfonium R group and N7 of guanine opposite cytosine at the GC site may contribute to specification of the pyrimidine. Substitution of G with inosine shifted HO(2)-Co(III)Blm A(2)[bond]I and Fe(III)Blm[bond]I into fast exchange on the NMR time scale, supporting the role of the 2-amino group in site specification for each molecule. The conformationally stable metal-domain linker established a close-packed adduct with the minor groove in which the hydroperoxide ligand occupies a sterically constrained pocket that is isolated from the solvent. The hydroperoxide group is directed toward one of the two cytosine H4' hydrogens but is sterically blocked from access to the other by the drug. These findings enlarge the structural understanding of selective binding of Co(III)/Fe(III)Blm species at G-pyrimidine sites. They also rationalize the instability of a number of ligands bound to Co(III)/Fe(III)Blm at specific binding sequences and the relative unreactivity of Fe(III)Blm[bond]I with ascorbate as well as its lack of interaction with spin labels.     

Experimental and computational investigations of Ser10 and Lys13 in the binding and cleavage of DNA substrates by Escherichia coli DNA topoisomerase I

Ser10 and Lys13 found near the active site tyrosine of Escherichia coli DNA topoisomerase I are conserved among the type IA topoisomerases. Site-directed mutagenesis of these two residues to Ala reduced the relaxation and DNA cleavage activity, with a more severe effect from the Lys13 mutation. Changing Ser10 to Thr or Lys13 to Arg also resulted in loss of DNA cleavage and relaxation activity of the enzyme. In simulations of the open form of the topoisomerase–DNA complex, Lys13 interacts directly with Glu9 (proposed to be important in the catalytic mechanism). This interaction is removed in the K13A mutant, suggesting the importance of lysine as either a proton donor or a stabilizing cation during strand cleavage, while the Lys to Arg mutation significantly distorts catalytic residues. Ser10 forms a direct hydrogen bond with a phosphate group near the active site and is involved in direct binding of the DNA substrate; this interaction is disturbed in the S10A and S10T mutants. This combination of a lysine and a serine residue conserved in the active site of type IA topoisomerases may be required for correct positioning of the scissile phosphate and coordination of catalytic residues relative to each other so that DNA cleavage and subsequent strand passage can take place.     

Electron transfer complexes of cytochrome c peroxidase from Paracoccus denitrificans containing more than one cytochrome

According to the model proposed in previous papers [Pettigrew, G. W., Prazeres, S., Costa, C., Palma, N., Krippahl, L., and Moura, J. J. (1999) The structure of an electron-transfer complex containing a cytochrome c and a peroxidase, J. Biol. Chem. 274, 11383-11389; Pettigrew, G. W., Goodhew, C. F., Cooper, A., Nutley, M., Jumel, K., and Harding, S. E. (2003) Electron transfer complexes of cytochrome c peroxidase from Paracoccus denitrificans, Biochemistry 42, 2046-2055], cytochrome c peroxidase of Paracoccus denitrificans can accommodate horse cytochrome c and Paracoccus cytochrome c(550) at different sites on its molecular surface. Here we use (1)H NMR spectroscopy, analytical ultracentrifugation, molecular docking simulation, and microcalorimetry to investigate whether these small cytochromes can be accommodated simultaneously in the formation of a ternary complex. The pattern of perturbation of heme methyl and methionine methyl resonances in binary and ternary solutions shows that a ternary complex can be formed, and this is confirmed by the increase in the sedimentation coefficient upon addition of horse cytochrome c to a solution in which cytochrome c(550) fully occupies its binding site on cytochrome c peroxidase. Docking experiments in which favored binary solutions of cytochrome c(550) bound to cytochrome c peroxidase act as targets for horse cytochrome c and the reciprocal experiments in which favored binary solutions of horse cytochrome c bound to cytochrome c peroxidase act as targets for cytochrome c(550) show that the enzyme can accommodate both cytochromes at the same time on adjacent sites. Microcalorimetric titrations are difficult to interpret but are consistent with a weakened binding of horse cytochrome c to a binary complex of cytochrome c peroxidase and cytochrome c(550) and binding of cytochrome c(550) to the cytochrome c peroxidase that is affected little by the presence of horse cytochrome c in the other site. The presence of a substantial capture surface for small cytochromes on the cytochrome c peroxidase has implications for rate enhancement mechanisms which ensure that the two electrons required for re-reduction of the enzyme after reaction with hydrogen peroxide are delivered efficiently.     

Plasticity of the TSG-6 HA-binding loop and mobility in the TSG-6-HA complex revealed by NMR and X-ray crystallography

Tumour necrosis factor-stimulated gene-6 (TSG-6) is a glycosaminoglycan-binding protein expressed during inflammatory and inflammation-like processes. Previously NMR structures were calculated for the Link module of TSG-6 (Link_TSG6) in its free state and when bound to an octasaccharide of hyaluronan (HA(8)). Heparin was found to compete for HA binding even though it interacts at a site that is distinct from the HA-binding surface. Here we present crystallography data on the free protein, and (15)N NMR relaxation data for the uncomplexed and HA(8)-bound forms of Link_TSG6. Although the Link module is comparatively rigid overall, the free protein shows a high degree of mobility in the beta4/beta5 loop and at the Cys47-Cys68 disulfide bond, both of which are regions involved in HA binding. When bound to HA(8), this dynamic behaviour is dampened, but not eliminated, suggesting a degree of dynamic matching between the protein and sugar that may decrease the entropic penalty of complex formation. A further highly dynamic residue is Lys54, which is distant from the HA-binding site, but was previously shown to be involved in heparin binding. When HA is bound, Lys54 becomes less mobile, providing evidence for an allosteric effect linking the HA and heparin-binding sites. A mechanism is suggested involving the beta2-strand and alpha2-helix. The crystal structure of free Link_TSG6 contains five molecules in the asymmetric unit that are highly similar to the NMR structure and support the dynamic behaviour seen near the HA-binding site: they show little or no electron density for the beta4/beta5 loop and display multiple conformations for the Cys47-Cys68 disulfide bond. The crystal structures were used in docking calculations with heparin. An extended interface between a Link_TSG6 dimer and heparin 11-mer was identified that is in excellent agreement with previous mutagenesis and calorimetric data, providing the basis for further investigation of this interaction.     

Molecular mechanics PBSA ligand binding energy and interaction of Efavirenz derivatives with HIV-1 reverse transcriptase

In order to evaluate the properties of several HIV-1 reverse transcripase(RT) inhibitors, Efavirenz (SUSTIVA) and a set of its derivatives (benzoxazinones) have been placed into the nonnucleoside analogue binding site of the enzyme by molecular docking. The resulting geometries were used for a molecular dynamics simulation and binding energy calculations. The enzyme-inhibitor binding energies were estimated from experimental inhibitory activities (IC90). The correlation of the predicted and experimental binding energies were satisfactory acceptable as indicated by r2 = 0.865. Based on MD simulations, the obtained results indicate that the tight association of the ligand to the HIV-1 RT binding pocket was based on hydrogen bonding between Efavirenz's N1 and the oxygen of the backbone of Lys 101, with an estimated average distance of 1.88 A. Moreover, electrostatic interaction was mainly contributed by two amino acid residues in the binding site; Lys 101 and His 235. MD simulations open the possibility to study the reaction of the flexible enzyme to those substances as well as the overall affinity.     

Auto-FACE: An NMR Based Binding Site Mapping Program for Fast Chemical Exchange Protein-Ligand Systems

Background

Nuclear Magnetic Resonance (NMR) spectroscopy offers a variety of experiments to study protein-ligand interactions at atomic resolution. Among these experiments, N Heteronuclear Single Quantum Correlation (HSQC) experiment is simple, less time consuming and highly informative in mapping the binding site of the ligand. The interpretation of N HSQC becomes ambiguous when the chemical shift perturbations are caused by non-specific interactions like allosteric changes and local structural rearrangement. Under such cases, detailed chemical exchange analysis based on chemical shift perturbation will assist in locating the binding site accurately.

Methodology/Principal Findings

We have automated the mapping of binding sites for fast chemical exchange systems using information obtained from N HSQC spectra of protein serially titrated with ligand of increasing concentrations. The automated program Auto-FACE (Auto-FAst Chemical Exchange analyzer) determines the parameters, e.g. rate of change of perturbation, binding equilibrium constant and magnitude of chemical shift perturbation to map the binding site residues. Interestingly, the rate of change of perturbation at lower ligand concentration is highly sensitive in differentiating the binding site residues from the non-binding site residues. To validate this program, the interaction between the protein and the ligand BH3I-1 was studied. Residues in the hydrophobic BH3 binding groove of were easily identified to be crucial for interaction with BH3I-1 from other residues that also exhibited perturbation. The geometrically averaged equilibrium constant () calculated for the residues present at the identified binding site is consistent with the values obtained by other techniques like isothermal calorimetry and fluorescence polarization assays (). Adjacent to the primary site, an additional binding site was identified which had an affinity of 3.8 times weaker than the former one. Further NMR based model fitting for individual residues suggest single site model for residues present at these binding sites and two site model for residues present between these sites. This implies that chemical shift perturbation can represent the local binding event much more accurately than the global binding event.

Conclusion/Significance

Detail NMR chemical shift perturbation analysis enabled binding site residues to be distinguished from non-binding site residues for accurate mapping of interaction site in complex fast exchange system between small molecule and protein. The methodology is automated and implemented in a program called “Auto-FACE”, which also allowed quantitative information of each interaction site and elucidation of binding mechanism.     

Preferential binding of equine ferricytochrome c to the bacterial photosynthetic reaction center from Rhodobacter sphaeroides

Redox titration of horse heart cytochrome c (cyt c), in the presence of varying concentrations of detergent-solubilized photosynthetic reaction center (RC) from Rhodobacter sphaeroides, revealed an RC concentration-dependent decrease in the measured cyt c midpoint potential that is indicative of a 3.6 +/- 0.2-fold stronger binding affinity of oxidized cytochrome to a single binding site. This effect was correlated with preferential binding in the functional complex by redox titration of the fraction of RCs exhibiting microsecond, first-order, special pair reduction by cytochrome. A binding affinity ratio of 3.1 +/- 0.4 was determined by this second technique, confirming the result. Redox titration of flash-induced intracomplex electron transfer also showed the association in the electron transfer-active complex to be strong, with a dissociation constant of 0.17 +/- 0.03 microM. The tight binding is associated with a slow off-rate which, in the case of the oxidized form, can influence the kinetics of P(+) reduction. The pitfalls of the common use of xenon flashlamps to photoexcite fast electron-transfer reactions are discussed with relation to the first electron transfer from primary to secondary RC quinone acceptors. The results shed some light on the diversity of kinetic behavior reported for the cytochrome to RC electron-transfer reaction.     

A two-process model describes the hydrogen exchange behavior of cytochrome c in the molten globule state with various extents of acetylation

Acetylation of Lys residues of horse cytochrome c steadily stabilizes the molten globule state in 18 mM HCl as more Lys residues are acetylated [Goto and Nishikiori (1991) J. Mol. Biol. 222, 679-686]. The dynamic features of the molten globule state were characterized by hydrogen/deuterium exchange of amide protons, monitored by mass spectrometry as each deuteration increased the protein mass by 1 Da. Electrospray mass spectrometry enabled us to monitor simultaneously the exchange kinetics of more than seven species with a different number of acetyl groups. One to four Lys residue-acetylated cytochrome c showed almost no protection of the amide protons from rapid exchange. The transition from the unprotected to the protected state occurred between five and eight Lys residue-acetylated species. For species with more than nine acetylated Lys residues, the exchange kinetics were independent of the extent of acetylation, and 26 amide protons were protected at 60 min of exchange, indicating the formation of a rigid hydrophobic core with hydrogen-bonded secondary structures. The apparent transition to the protected state required a higher degree of acetylation than the conformational transition measured by circular dichroism, which had a midpoint at about four acetylated residues. This difference in the transitions suggested a two-process model in which the exchange occurs either from the protected folded state or from the unprotected unfolded state through global unfolding. On the basis of a two-process model and with the reported values of the exchange and stability parameters, we simulated the exchange kinetics of a series of acetylated cytochrome c species. The simulated kinetics reproduced the observed kinetics well, indicating validity of this model for hydrogen exchange of the molten globule state.