Electrostatic orientation of the electron-transfer complex between plastocyanin and cytochrome c

Universally Accessible Structural Data on Macromolecular Conformation, Assembly, and Dynamics by Small Angle X-Ray Scattering for DNA Repair Insights

Structures provide a critical breakthrough step for biological analyses, and small angle X-ray scattering (SAXS) is a powerful structural technique to study dynamic DNA repair proteins. As toxic and mutagenic repair intermediates need to be prevented from inadvertently harming the cell, DNA repair proteins often chaperone these intermediates through dynamic conformations, coordinated assemblies, and allosteric regulation. By measuring structural conformations in solution for both proteins, DNA, RNA, and their complexes, SAXS provides insight into initial DNA damage recognition, mechanisms for validation of their substrate, and pathway regulation. Here, we describe exemplary SAXS analyses of a DNA damage response protein spanning from what can be derived directly from the data to obtaining super resolution through the use of SAXS selection of atomic models. We outline strategies and tactics for practical SAXS data collection and analysis. Making these structural experiments in reach of any basic and clinical researchers who have protein, SAXS data can readily be collected at government-funded synchrotrons, typically at no cost for academic researchers. In addition to discussing how SAXS complements and enhances cryo-electron microscopy, X-ray crystallography, NMR, and computational modeling, we furthermore discuss taking advantage of recent advances in protein structure prediction in combination with SAXS analysis.

X-ray scattering reveals disordered linkers and dynamic interfaces in complexes and mechanisms for DNA double-strand break repair impacting cell and cancer biology

Evolutionary selection ensures specificity and efficiency in dynamic metastable macromolecular machines that repair DNA damage without releasing toxic and mutagenic intermediates. Here we examine non‐homologous end joining (NHEJ) as the primary conserved DNA double‐strand break (DSB) repair process in human cells. NHEJ has exemplary key roles in networks determining the development, outcome of cancer treatments by DSB‐inducing agents, generation of antibody and T‐cell receptor diversity, and innate immune response for RNA viruses. We determine mechanistic insights into NHEJ structural biochemistry focusing upon advanced small angle X‐ray scattering (SAXS) results combined with X‐ray crystallography (MX) and cryo‐electron microscopy (cryo‐EM). SAXS coupled to atomic structures enables integrated structural biology for objective quantitative assessment of conformational ensembles and assemblies in solution, intra‐molecular distances, structural similarity, functional disorder, conformational switching, and flexibility. Importantly, NHEJ complexes in solution undergo larger allosteric transitions than seen in their cryo‐EM or MX structures. In the long‐range synaptic complex, X‐ray repair cross‐complementing 4 (XRCC4) plus XRCC4‐like‐factor (XLF) form a flexible bridge and linchpin for DNA ends bound to KU heterodimer (Ku70/80) and DNA‐PKcs (DNA‐dependent protein kinase catalytic subunit). Upon binding two DNA ends, auto‐phosphorylation opens DNA‐PKcs dimer licensing NHEJ via concerted conformational transformations of XLF‐XRCC4, XLF–Ku80, and LigIV^BRCT–Ku70 interfaces. Integrated structures reveal multifunctional roles for disordered linkers and modular dynamic interfaces promoting DSB end processing and alignment into the short‐range complex for ligation by LigIV. Integrated findings define dynamic assemblies fundamental to designing separation‐of‐function mutants and allosteric inhibitors targeting conformational transitions in multifunctional complexes.

Emerging critical roles of Fe-S clusters in DNA replication and repair

Fe-S clusters are partners in the origin of life that predate cells, acetyl-CoA metabolism, DNA, and the RNA world. The double helix solved the mystery of DNA replication by base pairing for accurate copying. Yet, for genome stability necessary to life, the double helix has equally important implications for damage repair. Here we examine striking advances that uncover Fe-S cluster roles both in copying the genetic sequence by DNA polymerases and in crucial repair processes for genome maintenance, as mutational defects cause cancer and degenerative disease. Moreover, we examine an exciting, controversial role for Fe-S clusters in a third element required for life - the long-range coordination and regulation of replication and repair events. By their ability to delocalize electrons over both Fe and S centers, Fe-S clusters have unbeatable features for protein conformational control and charge transfer via double-stranded DNA that may fundamentally transform our understanding of life, replication, and repair. This article is part of a Special Issue entitled: Fe/S proteins: Analysis, structure, function, biogenesis and diseases.

Envisioning the dynamics and flexibility of Mre11-Rad50-Nbs1 complex to decipher its roles in DNA replication and repair

The Mre11-Rad50-Nbs1 (MRN) complex is a dynamic macromolecular machine that acts in the first steps of DNA double strand break repair, and each of its components has intrinsic dynamics and flexibility properties that are directly linked with their functions. As a result, deciphering the functional structural biology of the MRN complex is driving novel and integrated technologies to define the dynamic structural biology of protein machinery interacting with DNA. Rad50 promotes dramatic long-range allostery through its coiled-coil and zinc-hook domains. Its ATPase activity drives dynamic transitions between monomeric and dimeric forms that can be modulated with mutants modifying the ATPase rate to control end joining versus resection activities. The biological functions of Mre11's dual endo- and exonuclease activities in repair pathway choice were enigmatic until recently, when they were unveiled by the development of specific nuclease inhibitors. Mre11 dimer flexibility, which may be regulated in cells to control MRN function, suggests new inhibitor design strategies for cancer intervention. Nbs1 has FHA and BRCT domains to bind multiple interaction partners that further regulate MRN. One of them, CtIP, modulates the Mre11 excision activity for homologous recombination repair. Overall, these combined properties suggest novel therapeutic strategies. Furthermore, they collectively help to explain how MRN regulates DNA repair pathway choice with implications for improving the design and analysis of cancer clinical trials that employ DNA damaging agents or target the DNA damage response.

Structural and mechanistic insights into the electron flow through protein for cytochrome c-tethering copper nitrite reductase

Copper-containing nitrite reductases (CuNiRs), which catalyse the reversible one-electron reduction of nitrite to nitric oxide, are members of a large family of multi-copper enzymes that require an interprotein electron transfer (ET) reaction with redox partner proteins. Here, we show that the naturally fused type of CuNiR tethering a cytochrome c (Cyt c) at the C-terminus folds as a unique trimeric domain-swapped structure and has a self-sufficient electron flow system. The C-terminal Cyt c domain is located at the surface of the type 1 copper (T1Cu) site in the N-terminal CuNiR domain from the adjacent subunit, the heme-to-Cu distance (10.6 Å) of which is comparable to the transient ET complex of normal CuNiR with Cyt c. The structural aspects for the domain-domain interface and the ET kinetics indicate that the Cyt c-CuNiR domain interaction should be highly transient. The further electrochemical analysis of the interprotein ET reaction with a cognate redox partner protein suggested that an electron is directly transferred from the partner to the T1Cu. Structural and mechanistic comparisons of Cyt c-CuNiR with another cupredoxin-tethering CuNiR highlight the behaviours of extra domains on the fusion types of CuNiRs required for ET through proteins.

The disordered C-terminal domain of human DNA glycosylase NEIL1 contributes to its stability via intramolecular interactions

NEIL1 [Nei (endonuclease VIII)-like protein 1], one of the five mammalian DNA glycosylases that excise oxidized DNA base lesions in the human genome to initiate base excision repair, contains an intrinsically disordered C-terminal domain (CTD; ~100 residues), not conserved in its Escherichia coli prototype Nei. Although dispensable for NEIL1's lesion excision and AP lyase activities, this segment is required for efficient in vivo enzymatic activity and may provide an interaction interface for many of NEIL1's interactions with other base excision repair proteins. Here, we show that the CTD interacts with the folded domain in native NEIL1 containing 389 residues. The CTD is poised for local folding in an ordered structure that is induced in the purified fragment by osmolytes. Furthermore, deletion of the disordered tail lacking both Tyr and Trp residues causes a red shift in NEIL1's intrinsic Trp-specific fluorescence, indicating a more solvent-exposed environment for the Trp residues in the truncated protein, which also exhibits reduced stability compared to the native enzyme. These observations are consistent with stabilization of the native NEIL1 structure via intramolecular, mostly electrostatic, interactions that were disrupted by mutating a positively charged (Lys-rich) cluster of residues (amino acids 355-360) near the C-terminus. Small-angle X-ray scattering (SAXS) analysis confirms the flexibility and dynamic nature of NEIL1's CTD, a feature that may be critical to providing specificity for NEIL1's multiple, functional interactions.

Human flap endonuclease structures, DNA double-base flipping, and a unified understanding of the FEN1 superfamily

Flap endonuclease (FEN1), essential for DNA replication and repair, removes RNA and DNA 5' flaps. FEN1 5' nuclease superfamily members acting in nucleotide excision repair (XPG), mismatch repair (EXO1), and homologous recombination (GEN1) paradoxically incise structurally distinct bubbles, ends, or Holliday junctions, respectively. Here, structural and functional analyses of human FEN1:DNA complexes show structure-specific, sequence-independent recognition for nicked dsDNA bent 100° with unpaired 3' and 5' flaps. Above the active site, a helical cap over a gateway formed by two helices enforces ssDNA threading and specificity for free 5' ends. Crystallographic analyses of product and substrate complexes reveal that dsDNA binding and bending, the ssDNA gateway, and double-base unpairing flanking the scissile phosphate control precise flap incision by the two-metal-ion active site. Superfamily conserved motifs bind and open dsDNA; direct the target region into the helical gateway, permitting only nonbase-paired oligonucleotides active site access; and support a unified understanding of superfamily substrate specificity.

XPB and XPD helicases in TFIIH orchestrate DNA duplex opening and damage verification to coordinate repair with transcription and cell cycle via CAK kinase

Helicases must unwind DNA at the right place and time to maintain genomic integrity or gene expression. Biologically critical XPB and XPD helicases are key members of the human TFIIH complex; they anchor CAK kinase (cyclinH, MAT1, CDK7) to TFIIH and open DNA for transcription and for repair of duplex distorting damage by nucleotide excision repair (NER). NER is initiated by arrested RNA polymerase or damage recognition by XPC-RAD23B with or without DDB1/DDB2. XP helicases, named for their role in the extreme sun-mediated skin cancer predisposition xeroderma pigmentosum (XP), are then recruited to asymmetrically unwind dsDNA flanking the damage. XPB and XPD genetic defects can also cause premature aging with profound neurological defects without increased cancers: Cockayne syndrome (CS) and trichothiodystrophy (TTD). XP helicase patient phenotypes cannot be predicted from the mutation position along the linear gene sequence and adjacent mutations can cause different diseases. Here we consider the structural biology of DNA damage recognition by XPC-RAD23B, DDB1/DDB2, RNAPII, and ATL, and of helix unwinding by the XPB and XPD helicases plus the bacterial repair helicases UvrB and UvrD in complex with DNA. We then propose unified models for TFIIH assembly and roles in NER. Collective crystal structures with NMR and electron microscopy results reveal functional motifs, domains, and architectural elements that contribute to biological activities: damaged DNA binding, translocation, unwinding, and ATP driven changes plus TFIIH assembly and signaling. Coupled with mapping of patient mutations, these combined structural analyses provide a framework for integrating and unifying the rich biochemical and cellular information that has accumulated over forty years of study. This integration resolves puzzles regarding XP helicase functions and suggests that XP helicase positions and activities within TFIIH detect and verify damage, select the damaged strand for incision, and coordinate repair with transcription and cell cycle through CAK signaling.

Effect of the ordered interfacial water layer in protein complex formation: A nonlocal electrostatic approach

Using a nonlocal electrostatic approach that incorporates the short-range structure of the contacting media, we evaluated the electrostatic contribution to the energy of the complex formation of two model proteins. In this study, we have demonstrated that the existence of an ordered interfacial water layer at the protein-solvent interface reduces the charging energy of the proteins in the aqueous solvent, and consequently increases the electrostatic contribution to the protein binding (change in free energy upon the complex formation of two proteins). This is in contrast with the finding of the continuum electrostatic model, which suggests that electrostatic interactions are not strong enough to compensate for the unfavorable desolvation effects.

The structural biochemistry of the superoxide dismutases

The discovery of superoxide dismutases (SODs), which convert superoxide radicals to molecular oxygen and hydrogen peroxide, has been termed the most important discovery of modern biology never to win a Nobel Prize. Here, we review the reasons this discovery has been underappreciated, as well as discuss the robust results supporting its premier biological importance and utility for current research. We highlight our understanding of SOD function gained through structural biology analyses, which reveal important hydrogen-bonding schemes and metal-binding motifs. These structural features create remarkable enzymes that promote catalysis at faster than diffusion-limited rates by using electrostatic guidance. These architectures additionally alter the redox potential of the active site metal center to a range suitable for the superoxide disproportionation reaction and protect against inhibition of catalysis by molecules such as phosphate. SOD structures may also control their enzymatic activity through product inhibition; manipulation of these product inhibition levels has the potential to generate therapeutic forms of SOD. Markedly, structural destabilization of the SOD architecture can lead to disease, as mutations in Cu,ZnSOD may result in familial amyotrophic lateral sclerosis, a relatively common, rapidly progressing and fatal neurodegenerative disorder. We describe our current understanding of how these Cu,ZnSOD mutations may lead to aggregation/fibril formation, as a detailed understanding of these mechanisms provides new avenues for the development of therapeutics against this so far untreatable neurodegenerative pathology.

Investigating the mechanisms of photosynthetic proteins using continuum electrostatics

Computational methods based on continuum electrostatics are widely used in theoretical biochemistry to analyze the function of proteins. Continuum electrostatic methods in combination with quantum chemical and molecular mechanical methods can help to analyze even very complex biochemical systems. In this article, applications of these methods to proteins involved in photosynthesis are reviewed. After giving a short introduction to the basic concepts of the continuum electrostatic model based on the Poisson-Boltzmann equation, we describe the application of this approach to the docking of electron transfer proteins, to the comparison of isofunctional proteins, to the tuning of absorption spectra, to the analysis of the coupling of electron and proton transfer, to the analysis of the effect of membrane potentials on the energetics of membrane proteins, and to the kinetics of charge transfer reactions. Simulations as those reviewed in this article help to analyze molecular mechanisms on the basis of the structure of the protein, guide new experiments, and provide a better and deeper understanding of protein functions.

Evaluation of the influence of the internal aqueous solvent structure on electrostatic interactions at the protein-solvent interface by nonlocal continuum electrostatic approach

The dielectric properties of the polar solvent on the protein-solvent interface at small intercharge distances are still poorly explored. To deconvolute this problem and to evaluate the pair-wise electrostatic interaction (PEI) energies of the point charges located at the protein-solvent interface we used a nonlocal (NL) electrostatic approach along with a static NL dielectric response function of water. The influence of the aqueous solvent microstructure (determined by a strong nonelectrostatic correlation effect between water dipoles within the orientational Debye polarization mode) on electrostatic interactions at the interface was studied in our work. It was shown that the PEI energies can be significantly higher than the energies evaluated by the classical (local) consideration, treating water molecules as belonging to the bulk solvent with a high dielectric constant. Our analysis points to the existence of a rather extended, effective low-dielectric interfacial water shell on the protein surface. The main dielectric properties of this shell (effective thickness together with distance- and orientation-dependent dielectric permittivity function) were evaluated. The dramatic role of this shell was demonstrated when estimating the protein association rate constants.

Reduction of plastocyanin by tyrosine-containing oligopeptides

Oxidized plastocyanin (PC) was reduced with TyrTyrTyr and LysLysLysLysTyrTyrTyr (KKKKYYY) oligopeptides at neutral pH. The TyrTyrTyr site of the peptides provided an electron to the copper active site of PC, whereas the tetralysine site of KKKKYYY functioned as the recognition site for the negative patch of PC. The reciprocal initial rate constant (1/k(int)) increased linearly with the reciprocal TyrTyrTyr concentration and proton concentration, although the electron transfer rate decreased gradually with time. The results showed that PC was reduced by the deprotonated species of TyrTyrTyr. A linear increase of log k(int) with increase in the ionic strength was observed due to decrease in the electrostatic repulsion between negatively charged PC and deprotonated (TyrTyrTyr)(-). PC was reduced faster by an addition of KKKKYYY to the PC-TyrTyrTyr solution, although KKKKYYY could not reduce PC without TyrTyrTyr. The ESI-LCMS spectrum of the products from the reaction between PC and TyrTyrTyr showed molecular ion peaks at m/z 1015.7 and 1037.7, which suggested formation of a dimerized peptide that may be produced from the reaction of a tyrosyl radical. The results indicate that PC and the tyrosine-containing oligopeptides form an equilibrium, PC(ox)/(oligopeptide)(-)-->/<--PC(red)/(oligopeptide)(*). The equilibrium is usually shifted to the left, but could shift to the right when the produced oligopeptide radical reacts with unreacted peptides. For the reaction of PC with KKKKYYY in the absence of TyrTyrTyr, the produced KKKK(YYY)(*) radical peptide could not react with other KKKKYYY peptides, since they were positively charged. In the presence of both KKKKYYY and TyrTyrTyr, PC may interact effectively with KKKKYYY through its tetralysine site and receive an electron from its TyrTyrTyr site, where the produced KKKK(YYY)(*) may interact with TyrTyrTyr peptides.

Photo Processes on Self-Associated Cationic Porphyrins and Plastocyanin Complexes 1. Ligation of Plastocyanin Tyrosine 83 onto Metalloporphyrins and Electron-Transfer Fluorescence Quenching

The spectroscopic properties of the self-associated complexes formed between the anionic surface docking site of spinach plastocyanin and the cationic metalloporphyrins, in which the tyrosine 83 (Y83) moiety is placed just below the docking site, tetrakis(N-methyl-4-pyridyl)porphyrin (Pd(II)TMPyP(4+) and Zn(II)TMPyP(4+)), have been studied and reported herein. The fluorescence quenching phenomenon of the self-assembled complex of Zn(II)TMPyP(4+)/plastocyanin has also been discovered. The observed red-shifting of the Soret and Q-bands of the UV-visible spectra, ca. 9 nm for Pd(II)TMPyP(4+)/plastocyanin and ca. 6 nm for the Zn(II)TMPyP(4+)/plastocyanin complexes, was explained in terms of exciton theory coupled with the Gouterman model. Thus, the hydroxyphenyl terminus of the Y83 residue of the self-associated plastocyanin/cationic porphyrin complexes was implicated in the charge-transfer ligation with the central metal atoms of these metalloporphyrins. Moreover, ground-state spectrometric-binding studies between Pd(II)TMPyP(4+) and the Y83 mutant plastocyanin (Y83F-PC) system proved that Y83 moiety of plastocyanin played a critical role in the formation of such ion-pair complexes. Difference absorption spectra and the Job's plots showed that the electrostatic attractions between the cationic porphyrins and the anionic patch of plastocyanin, bearing the nearby Y83 residue, led to the predominant formation of a self-associated 1:1 complex in the ground-state with significantly high binding constants (K = (8.0 +/- 1.1) x 10(5) M(-1) and (2.7 +/- 0.8) x 10(6) M(-1) for Pd(II)TMPyP(4+) and zinc variant, respectively) in low ionic strength buffer, 1 mM KCl and 1 mM phosphate buffer (pH 7.4). Molecular modeling calculations supported the formation of a 1:1 self-associated complex between the porphyrin and plastocyanin with an average distance of ca. 9 A between the centers of mass of the porphyrin and Y83 positioned just behind the anionic surface docking site on the protein surface. The photoexcited singlet state of Zn(II)TMPyP(4+) was quenched by the Y83 residue of the self-associated plastocyanin in a static mechanism as evidenced by steady-state and time-resolved fluorescence experiments. Even when all the porphyrin was complexed (more than 97%), significant residual fluorescence from the complex was observed such that the amplitude of quenching of the singlet state of uncomplexed species was enormously obscured.

Application of new multi-resolution methods for the comparison of biomolecular electrostatic properties in the absence of global structural similarity

In this paper we present a method for the multi-resolution comparison of biomolecular electrostatic potentials without the need for global structural alignment of the biomolecules. The underlying computational geometry algorithm uses multi-resolution attributed contour trees (MACTs) to compare the topological features of volumetric scalar fields. We apply the MACTs to compute electrostatic similarity metrics for a large set of protein chains with varying degrees of sequence, structure, and function similarity. For calibration, we also compute similarity metrics for these chains by a more traditional approach based upon 3D structural alignment and analysis of Carbo similarity indices. Moreover, because the MACT approach does not rely upon pairwise structural alignment, its accuracy and efficiency promises to perform well on future large-scale classification efforts across groups of structurally-diverse proteins. The MACT method discriminates between protein chains at a level comparable to the Carbo similarity index method; i.e., it is able to accurately cluster proteins into functionally-relevant groups which demonstrate strong dependence on ligand binding sites. The results of the analyses are available from the linked web databases http://ccvweb.cres.utexas.edu/MolSignature/ and http://agave.wustl.edu/similarity/. The MACT analysis tools are available as part of the public domain library of the Topological Analysis and Quantitative Tools (TAQT) from the Center of Computational Visualization, at the University of Texas at Austin (http://ccvweb.csres.utexas.edu/software). The Carbo software is available for download with the open-source APBS software package at http://apbs.sf.net/.

Role of the Electrostatic Interactions in Pre-orientation of Subunits in the Formation of Protein-protein Complexes

Theoretical estimation of contribution of the electrostatic interactions to pre-orientation of ribonuclease subunits in process of complex formation was carried out. The subunit was considered as a multipole consisting of partial charges of all atoms of the molecule. The object of investigation was a system of two subunits with their centers of gravity fixed at some distance in vacuum. It was proposed that each subunit independently could rotate freely around its fixed center of gravity. The relative orientation states of the subunits in such system were searched at which the system has electrostatic energy minima (equilibrium states). In first approximation the equilibrium states were found using especially designed approximate method for electrostatic interaction energy calculation, which permitted to calculate and compare the energies of the system in 24(5) (approximately 8 10(6)) states with different mutual orientation of subunits. The angular coordinates of the found equilibrium states were further specified by calculation with gradient sliding method. Angular coordinates of the equilibrium states and the shapes of energy surface cuts along each coordinate angle were calculated also for the intersubunits distances diminished down to 50 angstroms. The dispersions of the angular coordinates of equilibrium states caused by heat movement (at T=300 degrees) and their changes with shortening the distance between centers of gravity of subunits were estimated. Mutual orientation of subunits in the equilibrium states of the system under consideration was found to be similar to their mutual orientations in complex. Also it was found that relaxation time of the system, caused by electrostatic interaction of subunits, after removing the system from an equilibrium state, is much less in vacuum than the mean time between their Brownian collisions at room temperature. It follows from these results that in the case of ribonuclease in vacuum the electrostatic interactions of its subunits must be strong enough to realize the effective pre-orientation of subunits during their Brownian approach from distances of the order 100 angstroms. Preliminary consideration taking into account the effect of surrounding water molecules on the electrostatic interactions of ribonuclease subunits showed that weakening of the interaction must be much less than in the case when one uses in its calculation the macroscopic dielectric permeability value equal to 80. So the results obtained for vacuum seem to be true for water solution also. More strict theoretical analysis of this problem will be carried out in the following publication.

Interaction of plastocyanin with oligopeptides: effect of lysine distribution within the peptide

We synthesized and purified four oligopeptides containing four lysines (KKKK, GKKGGKK, KKGGGKK, and KGKGKGK) as models for the plastocyanin (PC) interacting site of cytochrome f. These peptides competitively inhibited electron transfer between cytochrome c and PC. The inhibitory effect increased as the peptide concentrations were increased. The association constants between PC and the peptides did not differ significantly (3500-5100 M(-1)), although the association constant of PC-KGKGKGK was a little larger than the constants between PC and other peptides. Changes in the absorption spectrum of PC were observed when the peptides were added to the PC solution: peaks and troughs were detected at about 460 and 630 nm and at about 560 and 700 nm, respectively, in the difference absorption spectra between the spectra with and without peptides. These changes were attributed to the structural change at the copper site of PC by interaction with the peptides. The structural change was most significant when tetralysine was used. These results show that binding of the oligopeptide to PC is slightly more efficient when lysines are distributed uniformly within the peptide, whereas the structural change of PC becomes larger when the lysines are close to each other within the peptide.

Electrostatic analysis and Brownian dynamics simulation of the association of plastocyanin and cytochrome f

The oxidation of cytochrome f by the soluble cupredoxin plastocyanin is a central reaction in the photosynthetic electron transfer chain of all oxygenic organisms. Here, two different computational approaches are used to gain new insights into the role of molecular recognition and protein-protein association processes in this redox reaction. First, a comparative analysis of the computed molecular electrostatic potentials of seven single and multiple point mutants of spinach plastocyanin (D42N, E43K, E43N, E43Q/D44N, E59K/E60Q, E59K/E60Q/E43N, Q88E) and the wt protein was carried out. The experimentally determined relative rates (k(2)) for the set of plastocyanin mutants are found to correlate well (r(2) = 0.90 - 0.97) with the computed measure of the similarity of the plastocyanin electrostatic potentials. Second, the effects on the plastocyanin/cytochrome f association rate of these mutations in the plastocyanin "eastern site" were evaluated by simulating the association of the wild type and mutant plastocyanins with cytochrome f by Brownian dynamics. Good agreement between the computed and experimental relative rates (k(2)) (r(2) = 0.89 - 0.92) was achieved for the plastocyanin mutants. The results obtained by applying both computational techniques provide support for the fundamental role of the acidic residues at the plastocyanin eastern site in the association with cytochrome f and in the overall electron-transfer process.

Blue copper proteins: a comparative analysis of their molecular interaction properties

Blue copper proteins are type-I copper-containing redox proteins whose role is to shuttle electrons from an electron donor to an electron acceptor in bacteria and plants. A large amount of experimental data is available on blue copper proteins; however, their functional characterization is hindered by the complexity of redox processes in biological systems. We describe here the application of a semiquantitative method based on a comparative analysis of molecular interaction fields to gain insights into the recognition properties of blue copper proteins. Molecular electrostatic and hydrophobic potentials were computed and compared for a set of 33 experimentally-determined structures of proteins from seven blue copper subfamilies, and the results were quantified by means of similarity indices. The analysis provides a classification of the blue copper proteins and shows that (I) comparison of the molecular electrostatic potentials provides useful information complementary to that highlighted by sequence analysis; (2) similarities in recognition properties can be detected for proteins belonging to different subfamilies, such as amicyanins and pseudoazurins, that may be isofunctional proteins; (3) dissimilarities in interaction properties, consistent with experimentally different binding specificities, may be observed between proteins belonging to the same subfamily, such as cyanobacterial and eukaryotic plastocyanins; (4) proteins with low sequence identity, such as azurins and pseudoazurins, can have sufficient similarity to bind to similar electron donors and acceptors while having different binding specificity profiles.

Definition of the interaction domain for cytochrome c on cytochrome c oxidase. III. Prediction of the docked complex by a complete, systematic search

The electron transfer complex between bovine cytochrome c oxidase and horse cytochrome c has been predicted with the docking program DOT, which performs a complete, systematic search over all six rotational and translational degrees of freedom. Energies for over 36 billion configurations were calculated, providing a free-energy landscape showing guidance of positively charged cytochrome c to the negative region on the cytochrome c oxidase surface formed by subunit II. In a representative configuration, the solvent-exposed cytochrome c heme edge is within 4 A of the indole ring of subunit II residue Trp(104), indicating a likely electron transfer path. These two groups are surrounded by a small, hydrophobic contact region, which is surrounded by electrostatically complementary hydrophilic interactions. Cytochrome c/cytochrome c oxidase interactions of Lys(13) with Asp(119) and Lys(72) with Gln(103) and Asp(158) are the most critical polar interactions due to their proximity to the hydrophobic region and exclusion from bulk solvent. The predicted complex matches previous mutagenesis, binding, and time-resolved kinetics studies that implicate Trp(104) in electron transfer and show the importance of specific charged residues to protein affinity. Electrostatic forces not only enhance long range protein/protein association; they also predominate in short range alignment, creating the transient interaction needed for rapid turnover.

Structural and kinetic studies of the pyruvate-ferredoxin oxidoreductase/ferredoxin complex from Desulfovibrio africanus

The pyruvate-ferredoxin oxidoreductase (PFOR)/ferredoxin (Fd) system of Desulfovibrio africanus has been investigated with the aim of understanding more fully protein-protein interaction and the kinetic characteristics of electron transfer between the two redox partners. D. africanus contains three Fds (Fd I, Fd II and Fd III) able to function as electron acceptors for PFOR. The complete amino acid sequence of Fd II was determined by automatic Edman degradation. It revealed a striking similarity to that of Fd I. The protein consists of 64 residues and its amino acid sequence is in agreement with a molecular mass of 6822.5 Da as measured by electrospray MS. Fd II contains five cysteine residues of which the first four (Cys11, Cys14, Cys17 and Cys54) are likely ligands for the single [4Fe-4S] cluster. A covalently cross-linked complex between PFOR and Fd I or Fd II was obtained by using a water soluble carbodiimide. This complex exhibited a stoichiometry of one ferredoxin for one PFOR subunit and is dependent on the ionic strength. The second-order rate constants for electron transfer between PFOR and Fds determined electrochemically using cyclic voltammetry are 7 x 107 M-1.s-1 for Fd I and 2 x 107 M-1.s-1 for Fd II and Fd III. The Km values of PFOR for Fd I and Fd II measured both by the electrochemical and the spectrophotometric method have been found to be 3 microM and 5 microM, respectively. The three-dimensional modelling of Fd II and surface analysis of Fd I, Fd II and PFOR suggest that a protein-protein complex is likely to be formed between aspartic acid/glutamic acid invariant residues of Fds and lysine residues surrounding the distal [4Fe-4S] cluster of PFOR. All of these studies are indicative of the involvement of electrostatic interactions between the two redox partners.

Comparing the rates and the activation parameters for the forward reaction between the triplet state of zinc cytochrome c and cupriplastocyanin and the back reaction between the zinc cytochrome c cation radical and cuproplastocyanin

This is a comparative study of the photoinduced (so-called forward) electron-transfer reaction 3Zncyt/pc(II) --> Zncyt+/pc(I), between the triplet state of zinc cytochrome c (3Zncyt) and cupriplastocyanin [pc(II)], and the thermal (so-called back) electron-transfer reaction Zncyt+/pc(I) --> Zncyt/pc(II), between the cation (radical) of zinc cytochrome c (Zncyt+) and cuproplastocyanin [pc(I)], which follows it. Both reactions occur between associated (docked) reactants, and the respective unimolecular rate constants are kF and kB. Our previous studies showed that the forward reaction is gated by a rearrangement of the diprotein complex. Now we examine the back reaction and complare the two. We study the effects of temperature (in the range 273.3-302.9 K) and viscosity (in the range 1.00-17.4 cP) on the rate constants and determine enthalpies (DeltaH), entropies (DeltaS), and free energies (DeltaG) of activation. We compare wild-type spinach plastocyanin, the single mutants Tyr83Leu and Glu59Lys, and the double mutant Glu59Lys/Glu60Gln. The rate constant kB for wild-type spinach plastocyanin and its mutants markedly depends on viscosity, an indication that the back reaction is also gated. The activation parameters DeltaH and DeltaS show that the forward and back reactions have similar mechanisms, involving a rearrangement of the diprotein complex from the initial binding configuration to the reactive configuration. The rearrangements of the complexes 3Zncyt/pc(II) and Zncyt+/pc(I) that gate their respective reactions are similar but not identical. Since the back reaction of all plastocyanin variants is faster than the forward reaction, the difference in free energy between the docking and the reactive configuration is smaller for the back reaction than for the forward reaction. This difference is explained by the change in the electrostatic potential on the plastocyanin surface as Cu(II) is reduced to Cu(I). It is the smaller DeltaH that makes DeltaG smaller for the back reaction than for the forward reaction.

Effects of single and double mutations in plastocyanin on the rate constant and activation parameters for the rearrangement gating the electron-transfer reaction between the triplet state of zinc cytochrome c and cupriplastocyanin

The unimolecular rate constant for the photoinduced electron-transfer reaction 3Zncyt/pc(II) --> Zncyt+/pc(I) within the electrostatic complex of zinc cytochrome c and spinach cupriplastocyanin is kF. We report the effects on kF of the following factors, all at pH 7.0: 12 single mutations on the plastocyanin surface (Leu12Asn, Leu12Glu, Leu12Lys, Asp42Asn, Asp42Lys, Glu43Asn, Glu59Gln, Glu59Lys, Glu60Gln, Glu60Lys, Gln88Glu, and Gln88Lys), the double mutation Glu59Lys/Glu60Gln, temperature (in the range 273.3-302.9 K), and solution viscosity (in the range 1. 00-116.0 cP) at 283.2 and 293.2 K. We also report the effects of the plastocyanin mutations on the association constant (Ka) and the corresponding free energy of association (DeltaGa) with zinc cytochrome c at 298.2 K. Dependence of kF on temperature yielded the activation parameters DeltaH, DeltaS, and DeltaG. Dependence of kF on solution viscosity yielded the protein friction and confirmed the DeltaG values determined from the temperature dependence. The aforementioned intracomplex reaction is not a simple electron-transfer reaction because donor-acceptor electronic coupling (HAB) and reorganizational energy (lambda), obtained by fitting of the temperature dependence of kF to the Marcus equation, deviate from the expectations based on precedents and because kF greatly depends on viscosity. This last dependence and the fact that certain mutations affect Ka but not kF are two lines of evidence against the mechanism in which the electron-transfer step is coupled with the faster, but thermodynamically unfavorable, rearrangement step. The electron-transfer reaction is gated by the slower, and thus rate determining, structural rearrangement of the diprotein complex; the rate constant kF corresponds to this rearrangement. Isokinetic correlation of DeltaH and DeltaS parameters and Coulombic energies of the various configurations of the Zncyt/pc(II) complex consistently show that the rearrangement is a facile configurational fluctuation of the associated proteins, qualitatively the same process regardless of the mutations in plastocyanin. Correlation of kF with the orientation of the cupriplastocyanin dipole moment indicates that the reactive configuration of the diprotein complex involves the area near the residue 59, between the upper acidic cluster and the hydrophobic patch. Kinetic effects and noneffects of plastocyanin mutations show that the rearrangement from the initial (docking) configuration, which involves both acidic clusters, to the reactive configuration does not involve the lower acidic cluster and the hydrophobic patch but involves the upper acidic cluster and the area near the residue 88.

Photoinduced electron transfer from carboxymethylated cytochrome c to plastocyanin

Photoinduced electron transfer from cytochrome c to plastocyanin was investigated using a novel method. Reduced carboxymethylated cytochrome c (CmCyt c), with carbon monoxide bound to the heme iron, and oxidized plastocyanin were mixed. At 1 mM CO the reduced state of CmCyt c is stabilized by about 350 meV. After flash photolysis of CO the apparent redox potential of CmCyt c drops resulting in electron transfer to plastocyanin. The electron transfer characteristics were investigated at approximately 30 different wavelengths in the range 390-460 nm. A global fit of the data showed that the electron transfer rate is 960+/-30 s-1 at pH 7.

Effects of viscosity and temperature on the kinetics of the electron-transfer reaction between the triplet state of zinc cytochrome c and cupriplastocyanin

This is a study of the effects of viscosity (in the range of 0.8-790 cP), of temperature (in the range of 260.7-307.7 K), and of ionic strength (in the range of 2.5-20.0 mM) on the kinetics of photoinduced electron-transfer reaction 3Zncyt/pc(II) --> Zncyt+/pc(I) within the electrostatic complex of zinc cytochrome c and cupriplastocyanin at pH 7.0. The unimolecular rate constant is kF. The apparent activation parameters DeltaH*, DeltaS*, and DeltaG* for this reaction were obtained in experiments with aqueous glycerol solutions having a constant composition. The interpolation of kF values obtained at the constant composition into the dependence of kF on temperature at constant viscosity gave the proper activation parameters, which agree with those obtained in experiments with solutions having a constant viscosity. This agreement validates the latter method, which is more efficient than the former, for determining activation parameters of processes that are modulated by viscosity. The smooth change in kF is governed by the change in viscosity, not in other properties of the solvent, and it does not depend on the choice of the viscosigen. Donor/acceptor electronic coupling (HAB) and reorganizational energy (lambda), obtained by fitting of the temperature dependence of kF to the Marcus equation, are consistent with true electron transfer and with electron transfer that is coupled to, or gated by, a preceding structural rearrangement of the diprotein complex 3Zncyt/pc(II). The fact that at very high viscosity kF approaches zero shows that the reaction is probably gated throughout the investigated range of viscosity. Kinetic effects and noneffects of ionic strength, viscosity, and thermodynamic driving force indicate, but do not prove, that the reaction under consideration is gated. The kinetic effect of viscosity is analyzed in terms of two models. Because ln kF is a nonlinear function of ln eta, protein friction has to be considered in the analysis of viscosity effects on kinetics.

Electrostatic complementarity at protein/protein interfaces

Calculation of the electrostatic potential of protein-protein complexes has led to the general assertion that protein-protein interfaces display "charge complementarity" and "electrostatic complementarity". In this study, quantitative measures for these two terms are developed and used to investigate protein-protein interfaces in a rigorous manner. Charge complementarity (CC) was defined using the correlation of charges on nearest neighbour atoms at the interface. All 12 protein-protein interfaces studied had insignificantly small CC values. Therefore, the term charge complementarity is not appropriate for the description of protein-protein interfaces when used in the sense measured by CC. Electrostatic complementarity (EC) was defined using the correlation of surface electrostatic potential at protein-protein interfaces. All twelve protein-protein interfaces studied had significant EC values, and thus the assertion that protein-protein association involves surfaces with complementary electrostatic potential was substantially confirmed. The term electrostatic complementarity can therefore be used to describe protein-protein interfaces when used in the sense measured by EC. Taken together, the results for CC and EC demonstrate the relevance of the long-range effects of charges, as described by the electrostatic potential at the binding interface. The EC value did not partition the complexes by type such as antigen-antibody and proteinase-inhibitor, as measures of the geometrical complementarity at protein-protein interfaces have done. The EC value was also not directly related to the number of salt bridges in the interface, and neutralisation of these salt bridges showed that other charges also contributed significantly to electrostatic complementarity and electrostatic interactions between the proteins. Electrostatic complementarity as defined by EC was extended to investigate the electrostatic similarity at the surface of influenza virus neuraminidase where the epitopes of two monoclonal antibodies, NC10 and NC41, overlap. Although NC10 and NC41 both have quite high values of EC for their interaction with neuraminidase, the similarity in electrostatic potential generated by the two on the overlapping region of the epitopes is insignificant. Thus, it is possible for two antibodies to recognise the electrostatic surface of a protein in dissimilar ways.

Temperature and viscosity dependence of the electron-transfer reaction between plastocyanin and cytochrome c labeled with a ruthenium(II) bipyridine complex

The temperature and viscosity dependence of the photo-induced electron-transfer reaction between plastocyanin and cytochrome c labeled at Lys13 with Ru(4,4'-dicarboxybipyridine)(bipyridine)(2+)2 have been investigated. In these studies, a short pulse of 450 nm light was used to excite the ruthenium complex which was oxidatively quenched by the iron center of cytochrome c. The resulting Fe(II) cytochrome c was then rapidly reoxidized by plastocyanin. The reactions were investigated over a temperature range of 3.5 to 37 degrees C under low ionic strength conditions such that protein/protein complex formation was favored. The enthalpy of activation was 7 kcal mol-1 and the entropy of activation was -20 cal mol-1 K-1. Increasing the viscosity by the addition of sucrose up to 70% resulted in a 4-fold decrease in the rate constant for electron transfer. The overall results suggest a rate-limiting step that involves either dissociation of the dominant protein/protein complex or surface diffusion of the associated proteins.

Effects of mutations in plastocyanin on the kinetics of the protein rearrangement gating the electron-transfer reaction with zinc cytochrome c. Analysis of the rearrangement pathway

We study, by flash kinetic spectrophotometry on the microsecond time scale, the effects of ionic strength and viscosity on the kinetics of oxidative quenching of the triplet state of zinc cytochrome c (3Zncyt) by the wild-type form and the following nine mutants of cupriplastocyanin: Leu12Glu, Leu12Asn, Phe35Tyr, Gln88Glu, Tyr83Phe, Tyr83His, Asp42Asn, Glu43Asn, and the double mutant Glu59Lys/Glu60Gln. The unimolecular rate constants for the quenching reactions within the persistent diprotein complex, which predominates at low ionic strengths, and within the transient diprotein complex, which is involved at higher ionic strengths, are equal irrespective of the mutation. Evidently, the two complexes are the same. In both reactions, the rate-limiting step is rearrangement of the diprotein complex from a configuration optimal for docking to the one optimal for the subsequent electron-transfer step, which is fast. We investigate the effects of plastocyanin mutations on this rearrangement, which gates the overall electron-transfer reaction. Conversion of the carboxylate anions into amide groups in the lower acidic cluster (residues 42 and 43), replacement of Tyr83 with other aromatic residues, and mutations in the hydrophobic patch in plastocyanin do not significantly affect the rearrangement. Conversion of a pair of carboxylate anions into a cationic and a neutral residue in the upper acidic cluster (residues 59 and 60) impedes the rearrangement. Creation of an anion at position 88, between the upper acidic cluster and the hydrophobic patch, facilitates the rearrangement. The rate constant for the rearrangement smoothly decreases as the solution viscosity increases, irrespective of the mutation. Fittings of this dependence to the modified Kramers's equation and to an empirical equation show that zinc cytochrome c follows the same trajectory on the surfaces of all the plastocyanin mutants but that the obstacles along the way vary as mutations alter the electrostatic potential. Mutations that affect protein association (i.e., change the binding constant) do not necessarily affect the reaction between the associated proteins (i.e., the rate constant) and vice versa. All of the kinetic and thermodynamic effects and noneffects of mutations consistently indicate that in the protein rearrangement the basic patch of zinc cytochrome c moves from a position between the two acidic clusters to a position at or near the upper acidic cluster.

Effects of temperature on the kinetics of the gated electron-transfer reaction between zinc cytochrome c and plastocyanin. Analysis of configurational fluctuation of the diprotein complex

This is a study of the effects of temperature (in the range 273.3-307.7 K) and of ionic strength (in the range 2.5-100 mM) on the kinetics of photoinduced electron-transfer reaction 3Zncyt/pc(II)--> Zncyt+/pc(I) within the electrostatic complex of zinc cytochrome c and cupriplastocyanin at pH 7.0. In order to separate direct and indirect effects of temperature on the rate constants, viscosity of the solutions was fixed, at different values, by additions of sucrose. The activation parameters for the reaction within the preformed complex, at the low ionic strength, are delta H++ = 13 +/- 2 kJ/mol and delta S++ = -97 +/- 4 J/K mol. The activation parameters for the reaction within the encounter complex, at the higher ionic strength, are delta H++ = 13 +/- 1 kJ/mol and delta S++ = -96 +/- 3 J/K mol. Evidently, the two complexes are the same. The proteins associate similarly in the persistent and the transient complex, i.e., at different ionic strengths. In both complexes, however, electron transfer is gated by a rearrangement, as previous studies from this laboratory showed. Changes in the solution viscosity modulate this rearrangement by affecting delta H++, not delta S++. The activation parameters are analyzed by empirical methods. The thermodynamic parameters delta H and delta S for the formation of the complex Zncyt/pc(II) are determined and related to changes in hydrophilic and hydrophobic surfaces upon protein association in three configurations. A difference between the values of delta H for the configuration providing optimal electronic coupling between the redox sites and the configuration providing optimal docking equals the experimental value delta H++ = 13 kJ/mol for the rearrangement of the latter configuration into the former. Enthalpy of activation may reflect a change in the character of the exposed surface as the diprotein complex rearranges. Entropy of activation may reflect tightening of the contact between the associated proteins.

Voltammetry of Plastocyanin at a Graphite Electrode: Effects of Structure, Charge, and Electrolyte

Comparative voltammetric studies on Anabaena variabilis plastocyanin (positively charged at neutral pH) and spinach and poplar plastocyanins (negatively charged at neutral pH) have been undertaken at an edge-plane graphite electrode as a function of ionic strength, pH, and Mg(2+) concentration at 3 degrees C. The aim was to provide a more detailed understanding of the influence of the electrode-protein (solution) interfacial characteristics, as well as the variation of the formal potential with both the nature of the plastocyanin species and the pH. As might be expected, some of the interfacial properties associated with the positive charge on A. variabilis plastocyanin are the opposite of those observed with the negatively charged plastocyanins. For example, the linear diffusion component of the mass transport process for A. variabilis plastocyanin under the conditions of cyclic voltammetry is decreased and the radial diffusion component is increased by the addition of Mg(2+), whereas the reverse occurs with poplar and spinach plastocyanins. The voltammetrically determined reversible potentials for A. variabilis plastocyanin are considerably less positive than those for spinach and poplar plastocyanins, in agreement with values calculated from chemically based redox studies. Ionic strength effects, as determined by addition of NaClO(4) over the concentration range 0.005-0.20 M, are negligible for all three proteins. The addition of Mg(2+) causes a significant shift in the reversible potential toward more positive values for spinach and poplar plastocyanin but only a small positive shift for A. variabilis plastocyanin. The difference is attributed to a specific binding effect. The addition of Mg(2+) also dramatically alters the pH dependence of the reversible potential, indicating that the equilibrium between the protonated and unprotonated forms of reduced plastocyanin is modified by binding of Mg(2+) to the protein. It is concluded that the effects of biologically relevant redox-inactive cations such as Mg(2+) or Ca(2+) have to be considered carefully in studies of the redox chemistry of metalloproteins.

Analysis of the 1H-NMR chemical shifts of Cu(I)-, Cu(II)- and Cd-substituted pea plastocyanin. Metal-dependent differences in the hydrogen-bond network around the copper site

To compare cadmium-substituted plastocyanin with copper plastocyanin, the 1H-NMR spectra of CuI-, CuII- and Cd-plastocyanin from pea have been analyzed. Full assignments of the spectra of CuI- and Cd-plastocyanin indicate chemical shift differences up to 1 ppm. The affected protons are located in the four loops that surround the Cu site. The largest differences were found for protons in the hydrogen bond network which stabilizes this part of the protein. This suggests that the chemical shift differences are caused by very small but extensive structural changes in the network upon replacement of CuI by Cd. For CuII-plastocyanin the resonances of 72% of the protons observed in the CuI form have been identified. Protons within approximately 0.9 nm of the CuII were not observed due to fast paramagnetic relaxation. The protons between 0.9-1.7 nm from the CuII showed chemical shift differences up to 0.4 ppm compared to both CuI- and Cd-plastocyanin. These differences can be predicted assuming that they represent pseudocontact shifts. When corrected for the pseudocontact shift contribution, the CuII-plastocyanin chemical shifts were nearly all identical within error to those of the Cd form, but not of the CuI-plastocyanin, indicating that the CuII-plastocyanin structure, in as far as it can be observed, resembles Cd-rather than CuI-plastocyanin. In a single stretch of residues (64-69) chemical shift differences remained between all three forms after correction. The fact that pseudocontact shifts were observed for protons which were not broadened may be attributable to the weaker distance dependence of the pseudocontact shift effect compared to paramagnetic relaxation. This results in two shells around the Cu atom, an inner paramagnetic shell (0-0.9 nm), in which protons are not observed due to broadening, and an outer paramagnetic shell (0.9-1.7 nm), in which protons can be observed and show pseudocontact shifts. It is concluded that Cd-plastocyanin is a suitable redox-inactive substitute for Cu-plastocyanin.

SOME NEW STRUCTURAL ASPECTS AND OLD CONTROVERSIES CONCERNING THE CYTOCHROME b6f COMPLEX OF OXYGENIC PHOTOSYNTHESIS

The cytochrome b6f complex functions in oxygenic photosynthetic membranes as the redox link between the photosynthetic reaction center complexes II and I and also functions in proton translocation. It is an ideal integral membrane protein complex in which to study structure and function because of the existence of a large amount of primary sequence data, purified complex, the emergence of structures, and the ability of flash kinetic spectroscopy to assay function in a readily accessible ms-100 mus time domain. The redox active polypeptides are cytochromes f and b6 (organelle encoded) and the Rieske iron-sulfur protein (nuclear encoded) in a mol wt = 210,000 dimeric complex that is believed to contain 22-24 transmembrane helices. The high resolution structure of the lumen-side domain of cytochrome f shows it to be an elongate (75 A long) mostly beta-strand, two-domain protein, with the N-terminal alpha-amino group as orthogonal heme ligand and an internal linear 11-A bound water chain. An unusual electron transfer event, the oxidant-induced reduction of a significant fraction of the p (lumen)-side cytochrome b heme by plastosemiquinone indicates that the electron transfer pathway in the b6f complex can be described by a version of the Q-cycle mechanism, originally proposed to describe similar processes in the mitochondrial and bacterial bc1 complexes.

Enzymatic and electron transfer activities in crystalline protein complexes

Enzymatic and electron transfer activities have been studied by polarized absorption spectroscopy in single crystals of both binary and ternary complexes of methylamine dehydrogenase (MADH) with its redox partners. Within the crystals, MADH oxidizes methylamine, and the electrons are passed from the reduced tryptophan tryptophylquinone (TTQ) cofactor to the copper of amicyanin and to the heme of cytochrome c551i via amicyanin. The equilibrium distribution of electrons among the cofactors, and the rate of heme reduction after reaction with substrate, are both dependent on pH. The presence of copper in the ternary complex is not absolutely required for electron transfer from TTQ to heme, but its presence greatly enhances the rate of electron flow to the heme.

Enforced interaction of one molecule of plastocyanin with two molecules of cytochrome c and an electron-transfer reaction involving the hydrophobic patch on the plastocyanin surface

Laser flash photolysis is used to study the photoinduced electron-transfer reaction cyt(III)//pc(II) + 3Zncyt --> cyt(III)//pc(I) + Zincyt+ at pH 7.0 and 25 degrees. In the covalent (symbol//) complex cyt(III)//pc(II) the acidic patch in cupriplastocyanin is directly cross-linked to the basic patch in ferricytochrome c. The triplet state of zinc cytochrome c reduces the pc(II) moiety, not the cyt(III) moiety, of the covalent complex. The reaction is strictly bimolecular in the entire range of ionic strength studied, from 1.25 mM to 1.00 M. The two reactants interact only transiently, in a collisional complex, and do not form a persistent complex cyt(III)//pc(II)/Zncyt. Because noncovalent (symbol/) association of three separate protein molecules is far less probable than association of the covalent complex and another protein molecule, we conclude that, without the aid of covalent cross-links, one molecule of plastocyanin will not form a ternary complex with two molecules of cytochrome c, cyt/pc/cyt. Dependence of the rate constant on ionic strength is analyzed in terms of van Leeuwen theory of electrostatic interactions, which recognizes the importance of dipole moments of the proteins. This analysis shows that 3Zncyt reacts with the hydrophobic patch in the pc(II) moiety of the covalent complex cyt(III)//pc(II). At high ionic strength, at which electrostatic interactions are practically abolished, the blue copper site is reduced with approximately equal rates via the hydrophobic patch in the pc(II) moiety of the complex and via the acidic patch in free pc(II). This is evidence that the two distinct patches on the plastocyanin surface are comparable in their intrinsic "conductivity" for electrons coming to the copper site. Positively charged and electroneutral redox partners tend to react at the acidic patch (although not necessarily at the initial docking site in this broad patch) for electrostatic, not electronic, reasons. Earlier theorectical studies disagreed about the relative electronic conductivities of the two patches. This experimental study corroborates very recent theoretical studies that found the two patches to be comparable in the efficiency of electron transfer.

Electrostatic potential distribution of the gene V protein from Ff phage facilitates cooperative DNA binding: a model of the GVP-ssDNA complex

The crystal structure of the gene V protein (GVP) from the Ff filamentous phages (M13, fl, fd) has been solved for the wild-type and two mutant (Y41F and Y41H) proteins at high resolution. The Y41H mutant crystal structure revealed crystal packing interactions, which suggested a plausible scheme for constructing the polymeric protein shell of the GVP-single-stranded DNA (ssDNA) complex (Guan Y, et al., 1994, Biochemistry 33:7768-7778). The electrostatic potentials of the isolated and the cooperatively formed protein shell have been calculated using the program GRASP and they revealed a highly asymmetric pattern of the electrostatic charge distribution. The inner surface of the putative DNA-binding channel is positively charged, whereas the opposite outer surface is nearly neutral. The electrostatic calculation further demonstrated that the formation of the helical protein shell enhanced the asymmetry of the electrostatic distribution. A model of the GVP-ssDNA complex with the n = 4 DNA-binding mode could be built with only minor conformational perturbation to the GVP protein shell. The model is consistent with existing biochemical and biophysical data and provides clues to the properties of GVP, including the high cooperatively of the protein binding to ssDNA. The two antiparallel ssDNA strands form a helical ribbon with the sugar-phosphate backbones at the middle and the bases pointing away from each other. The bases are stacked and the Phe 73 residue is intercalated between two bases. The optimum binding to a tetranucleotide unit requires the participation of four GVP dimers, which may explain the cooperativity of the GVP binding to DNA.

Structure-function relationships in the ferredoxin/ferredoxin: NADP+ reductase system from Anabaena

We have used a combination of laser flash photolysis time-resolved spectrophotometry and site-specific mutagenesis of surface amino acid residues to investigate the structural factors which influence electron transfer from Anabaena ferredoxin to its physiological partner ferredoxin-NADP+ reductase. Two ferredoxin residues (E94 and F65) are found to be highly critical interaction sites, whereas other nearby residues are found to be either inconsequential or to have only moderate effects. Basic residues near the N-terminus of the reductase are also found to exert a significant influence on interprotein electron transfer. The mechanistic implications of these results are discussed.

Thermodynamics of ferredoxin binding to ferredoxin:NADP+ reductase and the role of water at the complex interface

The association of ferredoxin with ferredoxin:NADP+ reductase (both proteins from spinach chloroplasts) was characterized by isothermal titration calorimetry and fluorescence quenching titration. The formation of the complex is mainly driven by a positive entropy change (delta S = 125 +/- 8 J mol-1 K-1). The calorimetric enthalpy of binding is small between 10 and 37 degrees C and either negative or positive, with an inversion temperature near 25 degrees C. The pH dependence of the association constant [Batie, C. J., & Kamin, H. (1981) J. Biol. Chem. 256, 7756-7763] was shown to correlate with the uptake of a single proton by a group exhibiting a heat of protonation of -26 kJ mol-1. This value agrees with the protonation of an imidazole group. Possible residues to become protonated in the complex are His-19 or His-90 of ferredoxin:NADP+ reductase. The temperature dependence of the free energy of binding, delta G, is weak because of the enthalpy-entropy compensation caused by a heat capacity change, delta Cp, of -680 +/- 44 J mol-1 K-1. The favorable binding entropy and the negative delta Cp indicate a large contribution to binding from hydrophobic effects, which seem to originate from dehydration of the protein-protein interface. Dehydration was demonstrated by osmotic stress experiments in which the association constant was found to increase by 2-4-fold in the presence of 52% (w/w) glycerol. The increase in the association constant with osmotic pressure points to the release of several water molecules from the complex interface.

A "parallel plate" electrostatic model for bimolecular rate constants applied to electron transfer proteins

A "parallel plate" model describing the electrostatic potential energy of protein-protein interactions is presented that provides an analytical representation of the effect of ionic strength on a biomolecular rate constant. The model takes into account the asymmetric distribution of charge on the surface of the protein and localized charges at the site of electron transfer that are modeled as elements of a parallel plate condenser. Both monopolar and dipolar interactions are included. Examples of simple (monophasic) and complex (biphasic) ionic strength dependencies obtained from experiments with several electron transfer protein systems are presented, all of which can be accommodated by the model. The simple cases do not require the use of both monopolar and dipolar terms (i.e., they can be fit well by either alone). The biphasic dependencies can be fit only by using dipolar and monopolar terms of opposite sign, which is physically unreasonable for the molecules considered. Alternatively, the high ionic strength portion of the complex dependencies can be fit using either the monopolar term alone or the complete equation; this assumes a model in which such behavior is a consequence of electron transfer mechanisms involving changes in orientation or site of reaction as the ionic strength is varied. Based on these analyses, we conclude that the principal applications of the model presented here are to provide information about the structural properties of intermediate electron transfer complexes and to quantify comparisons between related proteins or site-specific mutants. We also conclude that the relative contributions of monopolar and dipolar effects to protein electron transfer kinetics cannot be evaluated from experimental data by present approximations.

Electron-transfer from cytochrome c to ascorbate oxidase and its type 2 copper-depleted derivatives

Rate constants have been determined for the electron-transfer reactions between reduced horse heart cytochrome c and resting cucumber ascorbate oxidase as functions of pH, ionic strength, and temperature. The second-order rate constant for the oxidation of reduced cytochrome c was determined to be k = 820 M-1 s-1 in 0.2 M phosphate buffer at pH 6.0 and 25 degrees C. The activation parameters were estimated to be delta H++ = 5 kJ mol-1 and delta S++ = -188 Jmol-1 K-1. The rate constants increased with decreasing buffer concentration, indicating that the electron-transfer from cytochrome c to ascorbate oxidase is realized by the local electrostatic interaction between them in spite of the reaction between positively charged proteins. Reactions of type 2 copper-depleted ascorbate oxidase whose type 3 coppers were in the reduced or oxidized form indicated that the type 1 copper site accepts an electron from cytochrome c. The reaction rate was remarkably increased with decreasing pH for both the native enzyme and derivatives. Further, on addition of hexametaphosphate anion the rate of the electron-transfer decreased because the association of both proteins to realize the electron-transfer was inhibited due to a change in distribution of the local charge on the protein surface(s).

Receptor binding properties of four-helix-bundle growth factors deduced from electrostatic analysis

Hormones of the hematopoietin class mediate signal transduction by binding to specific transmembrane receptors. Structural data show that the human growth hormone (hGH) forms a complex with a homodimeric receptor and that hGH is a member of a class of hematopoietins possessing an antiparallel 4-alpha-helix bundle fold. Mutagenesis experiments suggest that electrostatic interactions may have an important influence on hormone-receptor recognition. In order to examine the specificity of hormone-receptor complexation, an analysis was made of the electrostatic potentials of hGH, interleukin-2 (IL-2), interleukin-4 (IL-4), granulocyte colony-stimulating factor (G-CSF), granulocyte-macrophage colony-stimulating factor (GM-CSF), and the hGH and IL-4 receptors. The binding surfaces of hGH and its receptor, and of IL-4 and its receptor, show complementary electrostatic potentials. The potentials of the hGH and its receptor display approximately 2-fold rotational symmetry because the receptor subunits are identical. In contrast, the potentials of GM-CSF and IL-2 lack such symmetry, consistent with their known high affinity for hetero-oligomeric receptors. Analysis of the electrostatic potentials supports a recently proposed hetero-oligomeric model for a high-affinity IL-4 receptor and suggests a possible new receptor binding mode for G-CSF; it also provides valuable information for guiding structural and mutagenesis studies of signal-transducing proteins and their receptors.

Kinetics of the reduction of wild-type and mutant cytochrome c-550 by methylamine dehydrogenase and amicyanin from Thiobacillus versutus

To elucidate the kinetic properties of the methylamine dehydrogenase (MADH) redox chain of Thiobacillus versutus the reduction of cytochrome c-550 by MADH and amicyanin has been studied. Under steady state conditions, the rate constants of the reactions have been determined as a function of the ionic strength, both for wild type cytochrome c-550 and for mutants in which the conserved residue Lys14 has been replaced as follows: Lys14-->Gln (mutant [K14Q]cytochrome c-550) and Lys14-->Glu (mutant [K14E]cytochrome c-550). The second-order rate constant of the reduction of cytochrome c-550 by MADH shows a biphasic ionic-strength dependence. At low ionic strength the rate constant remains unchanged (wild type) or increases ([K14Q]cytochrome c-550) with increasing ionic strength, while at high salt concentrations the rate constant decreases monotonically as the ionic strength increases. It is suggested that conformational freedom exists in the association complex and that this is favourable for electron transfer. [K14Q]cytochrome c-550 and [K14E]cytochrome c-550 are reduced at rates 20-fold and 500-fold slower than wild-type cytochrome c-550 by MADH, due to a lower association constant. It is concluded that MADH possesses a negative patch with which cytochrome c-550 associates. Lys14 plays an important role in the formation of the reaction complex. The midpoint potentials of wild-type and mutant cytochrome c-550 have been determined by using cyclic voltammetry. [K14Q]cytochrome c-550 and [K14E]cytochrome c-550 show an increase in E0 of only 2 mV and 8 mV, respectively, compared to wild-type cytochrome c-550 (241 mV at pH 8.1). [K14Q]cytochrome c-550 and [K14E]cytochrome c-550 cytochrome c-550 are reduced by amicyanin at rates that are only slightly faster than for wild-type cytochrome c-550. The difference is partly attributable to the change in E0. High ionic strength results in a threefold increase in the rate in all three cases. These results indicate that charge interactions do not play a major role in the formation of the amicyanin/cytochrome c-550 reaction complex, suggesting an interaction at the hydrophobic patch of amicyanin. The reduction of cytochrome c-550 by MADH can be inhibited by Zn(2+)-substituted amicyanin. Ag(+)-amicyanin, however, has little effect on the reduction rate. These results suggest that MADH has a much higher affinity for Cu(2+)-amicyanin (substrate) than for Cu(+)-amicyanin (product). On the basis of these findings the roles of the components of the MADH redox chain are discussed.

Plastocyanin: structural and functional analysis

Plastocyanin is one of the best characterized of the photosynthetic electron transfer proteins. Since the determination of the structure of poplar plastocyanin in 1978, the structure of algal (Scenedesmus, Enteromorpha, Chlamydomonas) and plant (French bean) plastocyanins has been determined either by crystallographic or NMR methods, and the poplar structure has been refined to 1.33 A resolution. Despite the sequence divergence among plastocyanins of algae and vascular plants (e.g., 62% sequence identity between the Chlamydomonas and poplar proteins), the three-dimensional structures are remarkably conserved (e.g., 0.76 A rms deviation in the C alpha positions between the Chlamydomonas and poplar proteins). Structural features include a distorted tetrahedral copper binding site at one end of an eight-stranded antiparallel beta-barrel, a pronounced negative patch, and a flat hydrophobic surface. The copper site is optimized for its electron transfer function, and the negative and hydrophobic patches are proposed to be involved in recognition of physiological reaction partners. Chemical modification, cross-linking, and site-directed mutagenesis experiments have confirmed the importance of the negative and hydrophobic patches in binding interactions with cytochrome f and Photosystem I, and validated the model of two functionally significant electron transfer paths in plastocyanin. One putative electron transfer path is relatively short (approximately 4 A) and involves the solvent-exposed copper ligand His-87 in the hydrophobic patch, while the other is more lengthy (approximately 12-15 A) and involves the nearly conserved residue Tyr-83 in the negative patch.

Protein interaction sites obtained via sequence homology. The site of complexation of electron transfer partners of cytochrome c revealed by mapping amino acid substitutions onto three-dimensional protein surfaces

Amino acid substitutions in all but the most divergent of cytochromes c have been categorized as being conservative or radical and mapped onto the three-dimensional structure of yeast cytochrome c. Color-coded, space-filling representations reveal a large 24 A diameter surface area which is invariant or conservatively substituted on the front left face of the cytochrome c molecule. Chemical modifications and mutations which inhibit complex formation and electron transfer with reaction partners also map to this surface. In sharp contrast, the back side of the protein is randomly substituted with both conservative and radical replacements. The invariant/conservatively substituted surface on the front of cytochrome c thus defines the site of interaction with redox partners and provides a measure of its dimensions. Further, this analysis strongly suggests that there is only a single site of oxidation and reduction on cytochrome c for all of its physiological reactions. The same analysis applied to bacterial cytochrome c2 shows that its conserved surface is similar in size and location to that of cytochrome c. Analyses of native and model reaction partners of cytochromes c and c2, such as cytochrome b5, plastocyanin, and bacterial photosynthetic reaction centers, also reveal probable active site surfaces for complexation and electron transfer, which are complementary in size to that of the c-type cytochromes. The availability of a three-dimensional structure and of several closely related amino acid sequences for a given functional class of protein is the only limitation on this type of analysis, which can then serve as a basis for designing site-directed mutagenesis experiments.

Rigid-body docking with mutant constraints of influenza hemagglutinin with antibody HC19

An automatic docking algorithm has been applied to the modeling of the complex between hemagglutinin from influenza virus and the Fab fragment of a monoclonal antibody raised against this antigen. We have introduced here the use of biochemical information provided by mutants of hemagglutinin. The docking procedure finds a small number of candidate solutions where three sites of escape mutations are buried and form hydrogen bonds in the interface. The localization of the epitope is improved by additional biochemical data about mutants that do not affect antibody binding. Five candidate solutions with low energy, reasonably well-packed interfaces, and six to ten hydrogen bonds are compatible with mutant information. One of the five stands out as generally better than the others from these points of views.

Mutagenesis of the conserved lysine 14 of cytochrome c-550 from Thiobacillus versutus affects the protein structure and the electron self-exchange rate

The lysine residue K14 of cytochrome c-550 of Thiobacillus versutus has been mutated to a glutamine (Q) and a glutamate (E) residue. These mutations have a minimal effect on the pKa for replacement of the methionine ligand (the "alkaline transition"), indicating that a presumptive salt bridge between K14 and E11 does not help stabilize the native form. This is in contrast with mitochondrial cytochrome c, where the homologous K13 forms a structurally important salt bridge with glutamate 90. The NMR signals of protons close to the heme iron in wild-type and mutant ferricytochrome c-550 shift considerably with increasing ionic strength. These effects resemble those seen in mitochondrial cytochrome c upon addition of salt and upon complex formation with redox partners. It is likely that electrostatic screening of positive charges near the heme crevice leads to a slight redistribution of the electron density in the heme. At low ionic strength the NMR spectrum of wild-type cytochrome c-550 shows broad peaks. Line widths decrease upon addition of salt up to 200 mM. In K14Q and K14E cytochrome c-550 the line widths are much smaller at low ionic strength. Wild-type cytochrome c-550 may exist in two exchanging conformations, one of which may represent a more open (non-native) form, in analogy with cytochrome c. However, in the case of cytochrome c-550 this non-native form does not show ligand replacement. The electron self-exchange rates of wild type and mutants have been determined as a function of the ionic strength.(ABSTRACT TRUNCATED AT 250 WORDS)