<<

1

2

3 Thecytochromef–complexasamodeltostudytransientinteractions

4 betweenproteins ⇑ 5 Q1 IsabelCruz-Gallardo,IreneDíaz-Moreno,AntonioDíaz-Quintana,MiguelA.DelaRosa

6 InstitutodeBioquímicaVegetalyFotosíntesis,CentrodeInvestigacionesCientíficas,IsladelaCartuja(cicCartuja),UniversidaddeSevilla-CSIC,Avda.AméricoVespucio49, 7 Sevilla41092,Spain

8

2910 11 Transientcomplexes,withalifetimerangingbetweenmicrosecondsandseconds,areessential30 for 12 biochemicalreactionsrequiringafastturnover.Thatisthecaseoftheinteractionsbetween31proteins 13 engagedinelectrontransferreactions,whichareinvolvedinrelevantphysiologicalprocesses32 such 14 asrespirationand.Inthelatter,theproteinplastocyaninactsas33asoluble 15 34 16 carriertransferringelectronsbetweenthetwomembrane-embeddedcomplexescytochromeb6f 17 andphotosystemI.Herewereviewthecombinationofexperimentaleffortsintheliterature35 to 18 unveilthefunctionalandstructuralfeaturesofthecomplexbetweencytochromefandplastocya-36 19 nin,whichhavewidelybeenusedasasuitablemodelforanalyzingtransientredoxQ2inter37actions. 20 Ó 2011PublishedbyElsevierB.V.onbehalfoftheFederationofEuropeanBiochemical3Societies.8 21

22 Keywords: 23 Cytochrome f 24 Electrontransfer 25 Plastocyanin 26 Transientcomplex 27 Photosynthesis 28 39 40 41 42 1.Introduction recognitionasanencounterthatcanbetransientlystabilized54to yieldaproductivecomplexasanoutcome. 55 43 Protein–proteininteractionsarekeyprocessesintheproperElectrontransfer(ET)reactionsareexcellentexamplesof56tran- 44 operationoflivingcells.Differentkindsofinteractionscansientbecomplexes.Infact,solubleproteinsmediateredoxexchange57 45 distinguishedonthebasisofproteinbindingaffinities.Complexesbetweenlargemembranecomplexesinthephotosynthetic58and 46 characterizedbylowbindingaffinitiesoccurintransientprotein–respiratoryelectrontransportchainsviashort-livedinteraction59s. 47 proteininteractions,whichdisplayhighdissociationequilibriumForinstance,inoxygenicphotosynthesisthereisasolublemetal-60

48 constants ðKDÞ,withintherangelM–mM,of orlifetimesintheloproteinthatshuttleselectronsfromcytochromef (Cf),whichis 61

49 orderofmilliseconds[1,2]. acomponentofthemembrane-embeddedcytochrome6f (Cb6bf) 62 50 Transientcomplexesaretypicalofphysiologicalprocessesthatcomplex,toP700,whichisthespecialchlorophyllpairof63Photo-

51 requireacompromisebetweenbindingspecificityandturnoversystemto I(PSI)[5–7].Eitherplastocyanin(Pc)orcytochrome6 64c 52 runatanappropriaterate[3].Onthebasisofatwo-stepmodel(Cc6)can[4],playsucharoleofelectroncarrierbetweenf and 65C 53 theformationofafinalcomplexfirstentailsaprimaryunspecificPSI.HigherplantsonlycontainPc,whereasmostcyanobacte66 ria

andgreensynthesizeeitherPcc6 ordependingC onthe 67 relativeavailabilityofcopperand,theirrespectivecofactor68 Abbreviations: BD,Browniandynamics;bfC,cytochromeb f;Cc ,cytochrome 6 6 6 metals[8,9]. 69 c6;Cf, cytochrome f;CSP, chemical-shiftperturbations;ET,electrontransfer; Cf isanchoredtothemembrane,withinCb6thef 70 EXAFS,extendedX-rayabsorptionfinestructurek2,bimolecular; rateconstant;KD, dissociationequilibriumconstant;MD,moleculardynamics;NIR,nitritereductascomplex,e; byaC-terminaltransmembranehelixleavinga7128- NMR,nuclearmagneticresonance;Pc,Plastocyanin;PCS,pseudocontactshifts;kDapI,N-terminalsolubleportionexposedtothelumenwith72 a isoelectricpoint;PRE,paramagneticrelaxationenhancement;PSI,PhotosystemclearI; two-domainstructure.Thelargedomainharborstheheme73 RMSD,rootmeansquaredeviation;WT,wildtype;XANE,X-rayabsorptionnear 74 edgestructure;XAS,X-rayabsorptionspectrometry group,andthesmalldomainpossessesapatchofchargedresi- ⇑ Correspondingauthor.Fax:+34954460065. dues.Cf isconsideredanunusualc-typecytochromebecause 75 E-mailaddress:[email protected](M.A.DelaRosa). ofitsb-sheet-basedstructure,elongatedformandparticular76

f 2 I.Cruz-Gallardoetal./FEBSLettersxxx(2011)xxx–xxx

77 hemeaxialcoordinationwiththeN-terminusTyr1[10,11].Pc,withinspinachandpeaPc[21].Thosechargedresiduesatthe122basic 78 itsturn,isan11-kDacupredoxinwithaP-barrelstructureCf patcharecrucialfortheinvitrointeractionwithPc.123Lysines 79 formedbyeightP-strandsandasmalla-helix,alongwith58,65a and187fofdirectlyC interactwithacidicresiduesatsite124 80 coppercentrecoordinatedbytwo,onemethionine2ofPc.Theelectrostaticattractionbetweenthetwopartners125 81 andone[12]. resultsinabell-shapedreactionratedependenceonionicstrength126 82 Theoverallstructuresofbothf andCPcarehighlyconserved[16](Fig.1,upperpanel).Thisbehaviorhasbeenexplained127by 83 fromtohigherplants[13–15],butstrikingdiffer-assumingthatthecomplexgets‘locked’inanon-produ128ctive 84 encesoccurinsomeoftheirphysicalproperties.Thesurfaceelectrosnear taticorientationatlowionicstrength.Astheionicstrength129 85 thehememoietyinf isC mainlyhydrophobicforallorganisms.increases,therearrangementofbothproteinsintothecomplex130 86 However,aremarkablebasicridgeisfoundinthesmalldomaintakesplaceinordertoachieveawell-orientedandproductive131 87 forhigherplantsandgreenalgaeturningtoacidicincyanobacomplex,cte- whichcorrespondstothemaximumofthebell.Further132 88 ria.Ontheotherhand,Pcreliesontwofunctionalsites:aincreasehydro-inionicstrengthleadstocomplexdissociation[6].133Ithas 89 phobicpatchsurroundingthesolvent-accessiblecopperrecentlybeenproposedthattheexpression‘lockedcomplex’134 90 ligand,ortheso-called‘‘site1’’,andanelectrostaticallychargedshouldbeavoidedassuchanon-productivestateatlow135ionic 91 surfacearea,theso-called‘‘site2’’,whosenaturevariesfromstrengthone couldrathercorrespondtoasetofnon-productive136orien- 92 organismtoanother.Actually,site2ismainlyacidicintationsplantsslightlydifferentinenergy[4]. 137 93 andgreenalgae,whereasitrangesfromacidictobasicincyano-Invitrokineticanalysesoftheinteractionbetweenf andPcCin 138 94 bacteria.Thisfeatureisakeyforthedynamicsofthecomplexcyanobacteriahavemainlybeenaddressedinthemesophi139lic 95 andisthushereinreviewed. Nostoc [22,24]andthermophiliPhormidiumc [23,25,26]species. 140 96 TheCf–PccomplexhasextensivelybeenstudiedintheThelastelectrontransferrateconstantbetweentheNostoctwoproteins 141 97 yearsasamodeltounderstandthenatureofprotein–proteinmonotointer- nicallydecreaseswithincreasingionicstrength(Fig.1421, 98 actionsinETchains.Themaingoalofthisarticleistomiddlebrieflypanel).Thechargemutationsatsite2reducethe143ability 99 reviewnotonlyourcurrentunderstandingofthemechanismofPcoftooxidizef [22C],whereasthemutationsatthehydrophobic144 100 ETintransientcomplexes[16],butalsothedifferenttechniquespatchdecreaseitsreactivitytowardsthehemeprotein.Itis145worth 101 usedtoanalyzethestructuralfeaturesofsuchcomplexes[17]tomentionin thatthechargemutationofArg93–aresidue146ofPc 102 anamplesetofprokaryoticandeukaryoticorganisms. locatedattheinterfacebetweensites1and2–drastically147dimin- ishesthereactionrate,sorevealingthecrucialroleofsuch148amino 103 2.Kineticswithinthe–PcCf complex acidinETwithinthecomplex(Table1).Incontrast,the149charge replacementsatthesmalldomainf hardlyofC affecttheinterac- 150 104 KineticanalysesprovideinformationabouttheETreactiontionwithPc[24](Table1),therebysuggestingthatthespecificity151 105 mechanism,namelythelimitingstepandthenatureofthepartnerinthef–PcC interactionismostlydeterminedbythecupredoxin.152 106 interactions.Someyearsago,ourgroupproposedthreedifferentThishasbeencorroboratedbyNMRstudiesofthef–Pcmixed153C 107 kineticmechanismstoanalyzethewell-relatedinteractionscomplexesof betweenNostoc and Phormidium cyanobacterialpro- 154 155 108 PSIwithPcandc6 fromC awiderangeofphotosyntheticorgan-teins[27]. 109 isms[18–20].ThesekineticmodelscanalsobeappliedtotheredoxTheETrateconstantbetweenthewild-type(WT)forms156ofthe 157 110 interactionsbetweenf CandthetwosolublecarriersPcc6andtwoC Phormidium proteinsslightlydependsonionicstrength 111 [21–23].Someofthecurrentlyavailableexperimentaldata(Fig.for1,lowerpanel),butsite-directedmutagenesisofcertain158 112 theinteractionbetweenf andC PcaresummarizedinTablecharged1, residuesofbothproteinsrevealedaclearinfluence159over 113 wherethevaluesforthefollowingconstantsarepresented:bimo-thereactionrate(Table1)[23,25].Actually,adeeperanalysis160 of 161 114 lecularrateconstantforcomplexformationðk2Þ,equilibriumcon- theCf–Pccomplexbycontinuumelectrostaticsshowsthatnet 162 115 stantforassociationbetweenpartnersðKAÞ,andeffectiveelectron coulombicforcesbetweenproteinchargesareslightlyrepulsive 0 [28].Assumingadiffusioncontrolledreaction,theelectrostatic163 116 transferrateconstantðketÞ.Theinteractionbetweenf andC Cc6 117 cannotbeanalyzedbecauseofspectraloverlappingoftheforcestwodoinfluenceencountercomplexformation.Theprocess164also 118 hemeproteins. involvesspecifichydrophobicinteractionsofaromaticresidues165in 119 Electrostaticsdeterminesthekineticsoffthe–PcCreaction.In theN-terminalpeptideoff [26C]. 166 120 plants,sucheffectontheETratewasmainlyestablishedbyanalyz-Altogether,themajorcontributionk2toinhigherplantsand 167 121 ingtheinteractionsoflysinemutantsatthebasicridgef ofturnipNostocCcyanobacterium is electrostatics, which leads to the 168

Table1 KineticdatafortheredoxreactionbetweenWTandmutantf andformsPc.ofC

7 1 1 0 3 1 3 1 Cf Pc k2 (10 M s ) ket (10 s ) KA (10 M ) References TurnipWT SpinachWT 17.6±2.2 3a – [21,22] TurnipK187E SpinachWT 2.5±0.3 – – [21] TurnipWT PeaWT 17.5±0.3 – 6.9±0.18 [21] TurnipK65Q PeaWT 3.5±0.1 – 0.55±0.38 [21] Phormidium WT Phormidium WT 4.7b – 0.3 [23,39] Phormidium WT Phormidium D44A 6.0b ––[23] Phormidium WT Phormidium R93E 1.0b ––[23] Phormidium D63A Phormidium WT 3.1b ––[25] Nostoc WT Nostoc WT – 13.4 26±1 [22,27] Nostoc WT Nostoc D54K – 25.5 – [22] Nostoc WT Nostoc R93E – 1.8 – [22] Nostoc D64A Nostoc WT – 13.3 – [24] Phormidium WT Nostoc WT – – 12±1 [27] Prochlorothrix WT Prochlorothrix WT 20 – 25±2 [40]

0 k2,bimolecularrateconstantforcomplexformation;KA,equilibriumconstantforassociationbetweenpartners;ket,effectiveelectrontransferrateconstant. a FrenchbeanPcwasemployed. b Theoverallerrorkof2 withallPhormidium variantsisestimatedto6be5%.

Pleasecitethisarticleinpressas:Cruz-Gallardo,I.,etal.Thef–plastocytochrcyaninome complexasamodeltostudytransientinteractionsbetweenredox proteins.FEBSLett.(2011),doi:10.1016/j.febslet.2011.08.035 FEBS 34903 No. of Pages 8, Model 5G 30 August 2011

I. Cruz-Gallardo et al. / FEBS Letters xxx (2011) xxx–xxx 3

shifts (PCS) experienced by the Pc signals as restraints for 192 rigid-body docking calculations. Diamagnetic, chemical-shift per- 193 turbations (CSP) are caused by changes in the chemical environ- 194 ment suffered by the atoms located at the binding interface. 195 Structural restrains derived from CSP do not provide precise geo- 196 metrical information about angles and distances between the 197 two partners. However, measuring PCS renders geometrical 198 restraints that may be sufficient to define the relative orientation 199 of both proteins. PCS restraints are provided by the intrinsic para- 200 magnetic probe of Cf, which supplies the oxidized iron atom (Fe3+) 201 from the heme phorphyrin ring. PCS depend on the axial and rhom- 202 bic components of the magnetic susceptibility tensor and the an- 203 gles that the vector joining the heme Fe atom of Cf to the target 204 amide proton in Pc forms with them. Moreover, PCS are inversely 205 proportional to the cubed distance between any particular nucleus 206 of Pc and the paramagnetic center. Structure elucidation of the 207 Cf–Pc complex for different organisms shown in this section fits 208 well with the orientation of Cf within the Cbf complex as deter- 209 mined by X-ray crystallography [11,33]. 210 Some similar structural features, like interface areas of ca. 211 600–850 Å2 per protein, are shared by the Cf–Pc complexes from 212 prokaryotic and eukaryotic organisms. All the cases studied show 213 that Cu-ligand His87 in Pc is close to iron-coordinating Tyr1 in Cf 214 thereby providing an efficient pathway [34]. 215 However, remarkable differences are found not only when compar- 216 ing the structures of plant and cyanobacterial complexes, but also 217 when the comparisons are made between complexes from cyano- 218 bacteria (Fig. 1). 219 Plant heterologous complexes between Cf from turnip and Pc 220 from several other sources like spinach, parsley or poplar have 221 been solved [35–37]. Their overall structures show a binding mode 222 named side-on (Fig. 1, upper panel) in which the hydrophobic 223

Fig. 1. Comparison of Cf–Pc complexes from different organisms. Right, Superpo- patches of the two proteins keep in close contact with each other 224 sition of the average structure, with the lowest energy value, calculated from the 10 and the complementary charged regions are juxtaposed. In other 225 best NMR structures of each complex. Left, profile of the ionic strength dependence words, both Pc sites 1 and 2 are involved in the binding interface. 226 of the observed rate constant in each case. Upper, poplar [36]; middle, Nostoc The Fe–Cu distance varies from 10.9 Å in spinach to 13.9 Å in 227 [22,38]; and lower, Phormidium [25,39].Cf and Pc ribbons are in red and blue, poplar and 13.0 Å in parsley. Actually, poplar and parsley Pc are 228 respectively. The heme porphyrin ring is depicted in green sticks, whereas the copper and iron atoms are grey spheres. slightly tilted with regard to the position of spinach Pc [36]. Within 229 the poplar Cf–Pc complex, the approach of the side-chain of His87 230 from Pc to the heme ring at Cf is restricted by the loop containing 231 this copper ligand. In contrast, the parsley Cf–Pc model shows a 232 169 encounter between proteins and to the formation of a well-defined rotation of Pc, with respect to the spinach, in a direction opposite 233 170 functional complex. The opposite is valid for the thermophilic to that in poplar complex [37]. 234 171 Phormidium cyanobacterium, in which non-polar interactions play The cyanobacterial complexes have been analyzed using both Pc 235 172 the main role because of the small net contribution of electrostat- and Cf from the same source. The structure of the Cf–Pc complex 236 173 ics [23,26,28]. Nevertheless, the hydrophobic interactions are from Nostoc [38] shows a relative orientation that resembles the 237 174 essential for ET in all complexes, since the reorganization of water plant side-on fashion, but with the charges reversed. Actually, the 238 175 molecules at the interface should impair the charge transfer pro- positive charged residues at site 2 of Nostoc Pc – which indeed are 239 176 cess. In fact, the exclusion of water molecules yields a productive determinant in the isoelectric point (pI) of 8.4 for the whole protein 240 177 complex within which the ET takes place. Intriguingly, the reaction – are keeping salt bridges with negative residues at the hinge and 241 178 rates reach similar values under physiological salt concentration small domain of Cf. However, the pI of plant Pc is ca. 4.0 and its neg- 242 8 1 1 179 and temperature ðk2 10 M s Þ despite the apparent differ- atively charged site 2 interacts with the positive charged residues of 243 180 ences among organisms [23]. Cf. The hydrophobic interactions involving the two metal centers of 244 181 In vivo approaches are quite scarce. In the alga Chlamydomonas both proteins are also present (Fig. 1, middle panel). In contrast, Pc 245 182 reinhardtii [29,30], neutralization of Cf basic ridge has a negligible in the Phormidium complex is oriented in a head-on manner relative 246 183 effect on the rate of protein re-oxidation, suggesting that electro- to Cf, with only site 1 (hydrophobic patch) making contact to Cf 247 184 static interactions within the Cf–Pc complex are not so relevant (Fig. 1, lower panel). In Phormidium,Cf lacks the typical basic ridge, 248 185 in vivo in this organism. Experimental conditions mimicking vis- whereas Pc has fewer charged residues at site 2 [39] than Pc from 249 186 cosity, macromolecular crowding and temperature have also been plants and shows a dipole moment aligned with the main molecule 250 187 explored, but no relevant conclusions can be reached [31]. axis. This results in the attraction of the copper site towards Cf and 251 repulsion of site 2. The resulting net coulombic term between the 252 188 3. Defining the relative orientation between Cf and Pc two proteins is repulsive, but is mostly compensated by solvent 253 polarization phenomena to yield a small ionic strength effect [28]. 254 189 The structure of the Cf–Pc complex, as well as that of the Cf–Cc6 The Phormidium complex shows a highly dynamic nature that could 255 190 one, has been analyzed by nuclear magnetic resonance (NMR) explain the lower precision of the solved structure – with an RMSD 256 191 spectroscopy [17,32] by using the chemical and pseudo-contact of 3.7 Å – compared with other Cf–Pc models. It could also be 257

Please cite this article in press as: Cruz-Gallardo, I., et al. The cytochrome f–plastocyanin complex as a model to study transient interactions between redox proteins. FEBS Lett. (2011), doi:10.1016/j.febslet.2011.08.035 FEBS 34903 No. of Pages 8, Model 5G 30 August 2011

4 I. Cruz-Gallardo et al. / FEBS Letters xxx (2011) xxx–xxx

258 explained by differences in the dynamic behavior of the complexes, 259 a fact that is especially significant in those organisms that inhabit 260 niches at high temperatures. 261 Crossed complexes between proteins from Nostoc and Phormidi- 262 um revealed that it must be a low net electrostatic contribution in 263 the latter complex [27], in agreement with results on ET kinetics. 264 Although early NMR data [39] suggested the absence of electrostatic 265 interactions, several transient conformations involving electrostat- 266 ics could remain invisible by NMR, at least in the absence of para- 267 magnetic relaxation enhancement measurements (PRE). Therefore, 268 only the long-life conformation into the complex is detectable. In 269 addition, studies of such mixed complexes show that the differences 270 in interactions are mainly attributable to the surface properties of 271 Pc [27]. The NMR analyses (PCS measurements) show that the 272 Cf–Pc complex from Prochlorothrix hollandica, unlike the Nostoc 273 Cf–Pc complex, exhibits a side-on orientation, in which Pc site 2 is 274 not interacting with the small domain of Cf, thereby suggesting a 275 dynamic nature that resembles the Phormidium complex [40]. 276 An overall view of structural features for the well-defined ori- 277 ented Cf–Pc complexes postulates a predominant side-on confor- 278 mation highly dependent on electrostatics. As the electrostatic 279 contribution becomes lower, the orientation of Pc relative to Cf 280 shifts from side-on to head-on. All these observations are also in 281 agreement with the ET proposed kinetic mechanisms.

282 4. Metal cofactors within the Cf–Pc complex

283 Although there are a lot of experimental data on the interaction 284 between the two redox partners, the information on the effect of 285 binding on the metal cofactors and their properties is rather scarce. 286 In fact, the redox potential of Pc decreases in 30 mV upon binding 287 to Cf [41]. Noteworthy, the way on how metal sites can adapt to 288 changes in the protein matrix and modulate the ET reaction highly 289 contributes to understand the transient complexes involved in this 290 kind of processes. 291 The cyanobacterial Nostoc Cf–Pc complex has been analyzed by 292 X-ray absorption spectroscopy (XAS). Both Fe and Cu K-edge XAS 293 measurements of free and bound redox proteins have been studied 294 in solution, using either the oxidized or reduced species [42,43]. 295 In Cf the Fe atom of the heme group is axially coordinated by

296 two nitrogen atoms belonging to Tyr1 (N) and His26 (N2). Such 297 atypical coordination geometry involves the N-terminal , 298 which also takes part in hydrophobic interactions with Pc, accord- 299 ing to the structural elucidations of this transient complex [38]. 300 XAS measurements reveals that the Fe coordination geometry 301 remains unaltered upon binding to Pc. Getting a deeper insight into 302 the data, a slight distortion of the metal (Fe2þ) center geometry 303 seems to happen when reduced Cf binds to reduced Pc. The result- 304 ing geometry is closer to that of oxidized Cf, either free or bound to 305 oxidized Pc. Because of the smaller size of Fe3þ compared to Fe2þ, 306 the first fits better than the second into the heme group as Fe2þ is Fig. 2. Distortion of the copper geometry of Pc upon Cf binding in Phormidium. 307 slightly out of the plane [44]. However, Fe2þ could be driven back Upper, XAFS spectra at the Cu K-edge of oxidized Pc upon binding to oxidized Cf (black) and fitting in the k space either considering that the contribution of 308 to the ring upon complex formation. Iron-N-terminus nitrogen 2 2 Sd(Met97) is negligible (red, v : 16.97) or taking it into account (blue, v : 11.3). 309 bond in Cf is strong enough to prevent distortions over the coordi- Middle, Fourier Transform modules of free Pc (blue) and Pc-bound to Cf (red), both 310 nation geometry of the metal center, so facilitating a stable binding in the reduced state. The best fits for the data are depicted in dotted lines following 311 site to Pc and enhancing ET. the same color-code. The contribution of the Sd(Met97) atom is considered for the fit in R space. Lower, XANES region of the Cu K-edge XAS spectra for free Pc (blue) 312 The copper atom in Pc is coordinated by two nitrogen atoms and Pc-bound to Cf, both in the oxidized state (red). 313 from His39 and His92, as well as by two sulphur atoms from 314 Cys89 and Met97. All of them are well-conserved residues placed 315 at a protein loop that is at the interface with Cf [35,38,39]. Cu–K 316 edge XAS shows a remarkable distortion in the trigonal pyramidal (EXAFS) wave, indicating that the mobility of this side-chain is sub- 321 317 geometry of the copper coordination sphere upon binding to Cf, stantially restrained upon binding of Pc to Cf. Actually, the result- 322 318 regardless of its redox state (Fig. 2, upper and middle panel). The ing tetrahedral structure of the copper center within the Cf–Pc 323

319 main evidence for such a distortion is the contribution of the Sd complex exhibits a shorter Cu–Sd(Met97) distance with respect 324 320 atom from Met97 to the extended X-ray absorption fine structure to the crystallographic structure of free Pc [45]. Noteworthy, the 325

Please cite this article in press as: Cruz-Gallardo, I., et al. The cytochrome f–plastocyanin complex as a model to study transient interactions between redox proteins. FEBS Lett. (2011), doi:10.1016/j.febslet.2011.08.035 FEBS 34903 No. of Pages 8, Model 5G 30 August 2011

I. Cruz-Gallardo et al. / FEBS Letters xxx (2011) xxx–xxx 5

326 Cu–Sd(Met97) distance of Cf–bound Pc by EXAFS resembles that of The Cf–Pc complexes of a wide range of organisms, from cyano- 343 327 the crystal structure of oxidized nitrite reductase (NIR) [46], which bacteria to plants, have been studied by Brownian Dynamics (BD) 344 328 belongs to the so-called ‘‘perturbed’’ copper centers. In addition, simulations using crystal or NMR structures. The role of electro- 345 329 the electronic density around the copper atom increases when Pc static interactions in higher plant and algal complexes has been 346 330 binds to Cf in their oxidized states, according to the data from confirmed by BD simulations that are in good agreement with the 347

331 the X-ray absorption near edge structure (XANES) region (Fig. 2, k2 values for ET reactions previously reported [21,23,25]. The com- 348 332 lower panel). In fact, the redox potential of Pc becomes more neg- plex between turnip Cf and spinach Pc [47] and that between Cf and 349 333 ative upon binding to Cf [41]. As a result, the driving force for ET Pc from the green alga Chlamydomonas reinhardtii [48] have both 350 334 between both is significantly decreased by the been studied by performing systematic mutations of charged resi- 351 335 Cf–Pc interaction. In addition, the observed geometrical changes dues at the protein interfaces. These data support the model of an 352 336 in the first coordination sphere of copper within the Cf–Pc complex electrostatically driven encounter complex as a previous step of a 353 337 modulate the electron coupling along the different ET pathways final productive complex wherein efficient ET occurs. BD analyses 354 338 and hence the redox reaction. with Phormidium Cf and Pc from different sources [49] show consid- 355 erable variations in electrostatic forces, a finding that highlights the 356 339 5. Theoretical approaches key role of a few charged residues close to His92 (His87 in plants 357 and algae) for the complex interface. However, the computational 358 340 The resulting data from analyzing the weak Cf–Pc complex by studies show that Pc is highly mobile in all these complexes, oppo- 359 341 means of theoretical methods add useful information to our site to its behavior in plant complexes. Still, hydrophobic interac- 360 342 current knowledge of transient interactions. tions remain highly significant in all organisms. 361

Fig. 3. of the Phormidium Cf–Pc complex. (A) Drift of Pc around Cf along the MD trajectory in Ref. [28]. The dihedral first angle projections of the trajectories of iron from Cf (brown dots), copper (blue dots) and center of mass of sPc (grey dots) are represented with respect to the main axes and mass center of the large domain of Cf. (B) Interaction of the copper site and Arg93 in Pc with negative residues of Cf upon approaching of the two partners. Pc and copper center are in purple, and Cf in gold. The side chains of charged residues and the heme porphyrin ring are depicted in sticks following the same color-pattern. (C) Coordinate covariance matrix of a trajectory corresponding to the Cf–Pc complex obtained after aligning main chain atoms of the large domain of Cf. Residues 1–105 correspond to Pc, and 106–354 to Cf Positive values are shown in blue, whereas negative are in red. (D) Analysis of continuum electrostatics averaged along different time intervals of a MD computation of the Cf–Pc complex.

The insert summarizes the thermodynamic cycle used in the computation of the different terms in which DGele-bind stands for the sum of free energies that account for

distorting the free monomers to get the conformation they adopt in the complex, DGele-strain, and their subsequent rigid-body association, DGele-rigid. DGele-bind is represented by

white bars; DGele-strain by red-filled bars; DGele-rigid bars are colored in grey.

Please cite this article in press as: Cruz-Gallardo, I., et al. The cytochrome f–plastocyanin complex as a model to study transient interactions between redox proteins. FEBS Lett. (2011), doi:10.1016/j.febslet.2011.08.035 FEBS 34903 No. of Pages 8, Model 5G 30 August 2011

6 I. Cruz-Gallardo et al. / FEBS Letters xxx (2011) xxx–xxx

362 BD analyses are quite sensitive to small changes in conforma- Acknowledgements 425 363 tion of the proteins that are being modeled. For instance, the Cf– 364 Pc complex with spinach Pc adopted two different conformations The authors wish to thank the Spanish Ministry of Science and 426 365 when using distinct structural data sources [50]. On the basis of Innovation (BFU2009-07190) and the Andalusian Government 427 366 such finding, it was suggested that Pc may assume distinct confor- (BIO198 and P08-CVI-3876) for financial support. 428 367 mations in solution so as to yield a productive ET complex. 368 BD and other docking approaches treat proteins as rigid bodies References 429 369 [40], whereas Molecular Dynamics (MD) calculations provide addi- 370 tional information about internal protein motions. Despite all the [1] Perkins, J.R., Diboun, I., Dessailly, B.H., Lees, J.C. and Orengo, C. (2010) Transient 430 protein–protein interactions: structural, functional, and network properties. 431 371 kinetic, structural and theoretical studies performed on the Pho- Structure 18, 1233–1243. 432 372 rmidium Cf–Pc complex, there is certain controversy regarding [2] Janin, J. (2000) Kinetics and thermodynamics of protein–protein interactions 433 373 NMR and functional data. Hence, MD simulations combined to con- in: Protein–protein recognition (Kleanthous, C., Ed.), pp. 1–32, Oxford 434 University Press, Oxford. 435 374 tinuum electrostatic calculations on this particular interaction [28] [3] Bashir, Q., Scanu, S. and Ubbink, M. (2011) Dynamics in electron transfer 436 375 try to harmonize the two sets of data. The MD trajectories reveal protein complexes. FEBS J. 278, 1391–1400. 437 376 that Pc tilts towards the small domain of Cf approaching site 2 to [4] Ubbink, M. (2009) The courtship of proteins: Understanding the encounter 438 439 377 the heme protein (Fig. 3A). This involves the interaction of the complex. FEBS Lett. 583, 1060–1066. [5] Allen, J.F. (2004) Cytochrome b6: structure for signalling and vectorial 440 378 positive charges of the copper site and Arg93 with the negatively metabolism. Trends Plant Sci. 9, 130–137. 441 379 charged residues at the loop of Cf surrounding the heme group [6] Hervás, M., Navarro, J.A. and De la Rosa, M.A. (2003) Electron transfer between 442 443 380 (Fig. 3B). Thus, the relative orientation of Pc respect to Cf in the membrane complexes and soluble proteins in photosynthesis. Acc. Chem. Res. 36, 798–805. 444 381 Phormidium complex can be redefined as a tilted head-on confor- [7] Sun, J., Xu, W., Navarro, J.A., Hervás, M., De la Rosa, M.A. and Chitnis, P.R. 445 382 mation (see Fig. 1, lower panel). However, the repulsive forces (1999) Oxidizing side of the cyanobacterial : evidence for 446 383 make both proteins swing in an opposite but concerted manner interaction between the electron donor proteins and a lumenal surface helix of 447 the PsaB subunit. J. Biol. Chem. 274, 19048–19054. 448 384 within the transient Cf–Pc complex, as inferred from the negative [8] De la Rosa, M.A., Navarro, J.A., Hervás, M., De la Cerda, B., Molina-Heredia, F.P., 449 385 covariance between motions of Pc and those of the small domain Balme, A., Murdoch, P.S., Díaz-Moreno, I., Durán, R.V. and Hervás, M. (2002) An 450 386 of Cf (Fig. 3C). Interestingly, the conformation of Cf is strained upon evolutionary analysis of the reaction mechanisms of photosystem I reduction 451 by cytochrome c and plastocyanin. Bioelectrochemistry 55, 41–45. 452 387 6 binding of its partner and relaxes upon release. Although there are [9] Durán, R., Hervás, M., De la Rosa, M.A. and Navarro, J.A. (2004) The efficient 453 388 no direct contacts between Pc and the small domain of Cf, contin- functioning of photosynthesis and respiration in Synechocystis sp. PCC 6803 454 389 uum electrostatic calculations indicate that the long-range electro- strictly requires the presence of either cytochrome c6 or plastocyanin. J. Biol. 455 Chem. 279, 7229–7233. 456 390 static interactions between them are responsible for straining the [10] Martinez, S.E., Huang, D., Szczepaniak, A., Cramer, W.A. and Smith, J.L. (1994) 457 391 conformation of the partners (Fig. 3D). The thermodynamic cycle Crystal structure of cytochrome f reveals a novel cytochrome fold 458 459 392 shown in Fig. 3D reveals that DGele-bind can be split in two terms and unexpected haem ligation. Structure 2, 95–105. 460 393 if any one of the partners changes within the complex. The first [11] Stroebel, D., Choquet, Y., Popot, J.L. and Picot, D. (2003) An atypical haem in the . Nature 426, 413–418. 461 394 term is the ‘‘strain energy’’ ðDGele-strainÞ, which accounts for the cost [12] Coleman, P.M., Freeman, H.C., Guss, J.M., Murata, M., Norris, V.A., Ramshaw, 462 395 of structural changes needed for docking; the second one is the J.A.M. and Venkatappa, M.P. (1978) X-ray crystal structure analysis of 463 plastocyanin at 2.7 Aresolution. Nature 272, 319–324. 464 396 ‘‘rigid binding’’ ðDGele-rigidÞ component, which represents the [13] Carrel, C.J., Schlarb, B.G., Bendall, D.S., Howe, C.J., Cramer, W.A. and Smith, J.L. 465 397 binding energy of the strained partners. Fig. 3D also shows that (1999) Structure of the soluble domain of cytochrome f from cyanobacterium 466 398 the time intervals corresponding to the two major ensembles of Phormidium laminosum. Biochemistry 38, 9590–9599. 467 399 conformations (from 0 to 3.7 ns, and from 7 to 13 ns) are character- [14] Bond, C.S., Bendall, D.S., Freeman, H.C., Guss, J.M., Howe, C.J., Wagner, M.J. and 468 Wilce, M.C.J. (1999) The structure of plastocyanin from the cyanobacterium 469 400 ized by a positive DGele-strain term and by large strain energy values Phormidium laminosum. Acta Cryst. D55, 414–421. 470 401 (DGele-bind). However, DGele-bind becomes negative, with a concomi- [15] Chi, Y.I., Huang, L.S., Zhang, Z., Fernández-Velasco, J.G. and Berry, E.A. (2000) X- 471 ray structure of a truncated form of cytochrome f from Chlamydomonas 472 402 tant drop in DGele-strain, in the time interval from 4 to 6 ns corre- reindhardtii. Biochemistry 39, 7689–7701. 473 403 sponding to the transition between these two major populations. [16] Hope, A.B. (2000) Electron transfers amongst cytochrome f, plastocyanin and 474 404 An explanation for this binding energy decrease comes from the photosystem I: kinetics and mechanisms. Biochim. Biophys. Acta 1456, 5–26. 475 405 reaction field term that accounts for solvent polarization phenom- [17] Ubbink, M. (2004) Complexes of photosynthetic redox proteins studied by 476 NMR. Photosynth. Res. 81, 277–287. 477 406 ena. In addition, the electrostatic strain may play a key role in [18] Hervás, M., Navarro, J.A., Diaz, A., Bottin, H. and De la Rosa, M.A. (1995) Laser- 478 407 breaking the complex off after the attraction between the copper flash kinetic analysis of the fast electron transfer from plastocyanin and 479 408 surroundings and Cf is weakened upon charge transfer. cytochrome c6 to photosystem I. Experimental evidence on the evolution of 480 the reaction mechanism. Biochemistry 34, 11321–11326. 481 [19] Hervas, M., Navarro, J.A., Diaz, A. and De la Rosa, M.A. (1996) A comparative 482 409 6. Conclusion and outlook thermodynamic analysis by laser-flash absorption spectroscopy of 483 photosystem I reduction by plastocyanin and cytochrome c6 in Anabaena 484 PCC 7119, Synechocystis PCC 6803, and spinach. Biochemistry 35, 485 410 A multidisciplinary effort has been made in the last years to 2693–2698. 486 411 understand the behavior of transient ET complexes and, in partic- [20] De la Cerda, B., Navarro, J.A., Hervas, M. and De la Rosa, M.A. (1997) Changes in 487 412 ular, of the short-lived Cf–Pc complex from different organisms. In the reaction mechanism of electron transfer from plastocyanin to 488 photosystem I in the cyanobacterium Synechocystis sp. PCC 6803 as induced 489 413 most cases, the functional and structural approaches coincide in by site-directed mutagenesis of the . Biochemistry 36, 490 414 postulating the relevance of long-range electrostatic interactions 10125–10130. 491 415 in vitro. The existence of an encounter complex may indeed play [21] Gong, X.S., Wen, J.Q., Fisher, N.E., Young, S., Howe, C.J., Bendall, D.S. and Gray, 492 J.C. (2000) The role of individual lysine residues in the basic patch on turnip 493 416 an important role in the Cf–Pc complexes, a finding that should cytochrome f for electrostatic interactions with plastocyanin in vitro. Eur. J. 494 417 be further explored in detail. Some controversy still exists over Biochem. 267, 3461–3468. 495 418 experimental data interpretation in the cyanobacterial model sys- [22] Albarran, C., Navarro, J.A., Molina-Heredia, F.P., Murdoch, P.S., De la Rosa, M.A. 496 497 419 tem from Phormidium, but theoretical simulations are providing and Hervas, M. (2005) Laser flash-induced kinetic analysis of cytochrome f oxidation by wild-type and mutant plastocyanin from the cyanobacterium 498 420 valuable information to conciliate them. Actually, recent modeling Nostoc sp. PCC 7119. Biochemistry 44, 11601–11607. 499 421 developments could help to widen the scope of this subject. The [23] Schlarb-Ridley, B.G., Bendall, D.S. and Howe, C.J. (2002) Role of electrostatics in 500 501 422 bottleneck is to unveil the relevance of these data in vivo. Specifi- the interaction between cytochrome f and plastocyanin of the cyanobacterium Phormidium laminosum. Biochemistry 41, 3279–3285. 502 423 cally, the role played by electrostatics on the ET process inside [24] Albarran, C., Navarro, J.A., De la Rosa, M.A. and Hervas, M. (2007) The 503 424 crowded cells needs further research. specificity in the interaction between cytochrome f and plastocyanin from the 504

Please cite this article in press as: Cruz-Gallardo, I., et al. The cytochrome f–plastocyanin complex as a model to study transient interactions between redox proteins. FEBS Lett. (2011), doi:10.1016/j.febslet.2011.08.035 FEBS 34903 No. of Pages 8, Model 5G 30 August 2011

I. Cruz-Gallardo et al. / FEBS Letters xxx (2011) xxx–xxx 7

505 cyanobacterium Nostoc sp. PCC 7119 is mainly determined by the copper [37] Crowley, P.B., Hunter, D.M., Sato, K., McFarlane, W. and Dennison, C. (2004) 551 506 protein. Biochemistry 46, 997–1003. The parsley plastocyanin-turnip cytochrome f complex: a structurally 552 507 [25] Hart, S.E., Schlarb-Ridley, B.G., Delon, C., Bendall, D.S. and Howe, C.J. (2003) distorted but kinetically functional acidic patch. Biochem. J. 378, 45–51. 553 508 Role of charges on cytochrome f from the cyanobacterium Phormidium [38] Diaz-Moreno, I., Diaz-Quintana, A., De la Rosa, M.A. and Ubbink, M. (2005) 554 509 laminosum in its interaction with plastocyanin. Biochemistry 42, 4829–4836. Structure of the complex between plastocyanin and cytochrome f from the 555 510 [26] Schlarb-Ridley, B.G., Bendall, D.S. and Howe, C.J. (2003) Relation between cyanobacterium Nostoc sp. PCC 7119 as determined by paramagnetic NMR. J. 556 511 interface properties and kinetics of electron transfer in the interaction of Biol. Chem. 280, 18908–18915. 557 512 cytochrome f and plastocyanin from plants and the cyanobacterium [39] Crowley, P.B., Otting, G., Schlarb-Ridley, B.G., Canters, G.W. and Ubbink, M. 558 513 Phormidium laminosum. Biochemistry 42, 4057–4063. (2001) Hydrophobic interactions in a cyanobacterial plastocyanin-cytochrome 559 514 [27] Diaz-Moreno, I., Diaz-Quintana, A., De la Rosa, M.A., Crowley, P.B. and Ubbink, f complex. J. Am. Chem. Soc. 123, 10444–10453. 560 515 M. (2005) Different modes of interaction in cyanobacterial complexes of [40] Hulsker, R., Baranova, M.V., Bullerjahn, G.S. and Ubbink, M. (2008) Dynamics 561 516 plastocyanin and cytochrome. Biochemistry 44, 3176–3183. in the transient complex of plastocyanin and cytochrome f from 562 517 [28] Diaz-Moreno, I., Munoz- Lopez, F.J., Frutos-Beltran, E., De la Rosa, M.A. and Prochlorothrix hollandica. J. Am. Chem. Soc. 130, 1985–1991. 563 518 Diaz-Quintana, A. (2009) Electrostatic strain and concerted motions in the [41] Malkin, R., Knaff, D.B. and Bearden, A.J. (1973) The oxidation–reduction 564 519 transient complex between plastocyanin and cytochrome f from the potential of membrane-bound chloroplast plastocyanin and cytochrome f. 565 520 cyanobacterium Phormidium laminosum. Bioelectrochemistry 77, 43–52. Biochim. Biophys. Acta 305, 675–678. 566 521 [29] Soriano, G.M., Ponamarev, M.V., Piskorowski, R.A. and Cramer, W.A. (1998) [42] Diaz-Moreno, I., Diaz-Moreno, S., Subias, G., De la Rosa, M.A. and Diaz- 567 522 Identification of the basic residues of cytochrome f responsible for Quintana, A. (2006) The atypical iron-coordination geometry of cytochrome f 568 523 electrostatic docking interactions with plastocyanin in vivo: relevance to the remains unchanged upon binding to plastocyanin, as inferred by XAS. 569 524 electron transfer in vivo. Biochemistry 37, 15120–15128. Photosynth. Res. 90, 23–28. 570 525 [30] Soriano, G.M., Ponamarev, M.V., Tae, G.S. and Cramer, W.A. (1996) Effect of [43] Diaz-Moreno, I., Diaz-Quintana, A., Diaz-Moreno, S., Subias, G. and De la Rosa, 571 526 interdomain basic region of cytochrome f on its redox reactions in vivo. M.A. (2006) Transient binding of plastocyanin to its physiological redox 572 527 Biochemistry 35, 14590–14598. partners modifies the copper site geometry. FEBS Lett. 580, 6187–6194. 573 528 [31] Schlarb-Ridley, B.G., Mi, H., Teale, W.D., Meyer, V.S., Howe, C.J. and Bendall, [44] Schnackenberg, J., Than, M.E., Mann, K., Wiegand, G., Huber, R. and Reuter, W. 574 529 D.S. (2005) Implications of the effects of viscosity, macromolecular crowding, (1999) Aminoacid sequence, crystallization and structure determination of 575 530 and temperature for the transient interactions between cytochrome f and reduced and oxidized cytochrome c6 from the green alga Scenedesmus 576 531 plastocyanin from the cyanobacterium Phormidium laminosum. Biochemistry obliquus. J. Mol. Biol. 290, 1019–1030. 577 532 44, 6232–6238. [45] Schmidt, L., Christensen, H.E.M. and Harris, P. (2006) Structure of plastocyanin 578 533 [32] Crowley, P.B., Díaz-Quintana, A., Molina-Heredia, F.P., Nieto, P., Sutter, M., from the cyanobacterium Anabaena variabilis. Acta Cryst. D 62, 1022–1029. 579 534 Haehnel, W., De la Rosa, M.A. and Ubbink, M. (2002) The interactions of [46] Adman, E.T., Godden, J.W. and Turley, S. (1995) The structure of copper-nitrite 580 535 cyanobacterial cytochrome c6 and cytochrome f, characterized by NMR. J. Biol. reductase from Achromobacter cycloclastes at five pH values, with NO2 bound 581 536 Chem. 277, 48685–48689. and with type II copper depleted. J. Biol. Chem. 270, 27458–27474. 582 537 [33] Kurisu, G., Zhang, H., Smith, J.L. and Cramer, W.A. (2003) Structure of the [47] De Rienzo, F., Gabdoulline, R.R., Menziani, M.C., De Benedetti, P.G. and Wade, 583 538 cytochrome b6f complex of oxygenic photosynthesis: tuning the cavity. R.C. (2001) Electrostatic analysis and Brownian dynamics simulation of the 584 539 Science 302, 1009–1014. association of plastocyanin and cytochrome f. Biophys. J. 81, 3090–3104. 585 540 [34] Diaz-Quintana, A., Hervas, M., Navarro, J.A. and De la Rosa, M.A. (2008) [48] Haddadian, E.J. and Gross, E.L. (2005) Brownian dynamics study of cytochrome 586 541 ;;Plastocyanin and cytochrome c6: the soluble electron carriers between the f interactions with cytochrome c6 and plastocyanin in Chlamydomonas 587 542 cytochrome b6f complex and photosystem I in: Structure of photosynthetic reinhardtii: plastocyanin and cytochrome c6 mutants. Biophys. J. 88, 2823– 588 543 proteins (Fromme, P., Ed.), pp. 181–200, Wiley-VCH, Weinheim, Germany. 2839. 589 544 [35] Ubbink, M., Ejdeback, M., Karlsson, B.G. and Bendall, D.S. (1998) The structure [49] Gross, E.L. and Rosenberg, A. (2006) A Brownian dynamics study of the 590 545 of the complex of plastocyanin and cytochrome f, determined by interaction of Phormidium cytochrome f with various cyanobacterial 591 546 paramagnetic NMR and restrained rigid-body molecular dynamics. Structure . Biophys. J. 90, 366–380. 592 547 6, 323–335. [50] Gross, E.L. (2007) A Brownian dynamics computational study of the 593 548 [36] Lange, C., Cornvik, T., Diaz-Moreno, I. and Ubbink, M. (2005) The transient interaction of spinach plastocyanin with turnip cytochrome f: the 594 549 complex of poplar plastocyanin with cytochrome f: effects of ionic strength importance of plastocyanin conformational changes. Photosynth. Res. 94, 595 550 and pH. Biochim. Biophys. Acta 1707, 179–188. 411–422. 596 597

Please cite this article in press as: Cruz-Gallardo, I., et al. The cytochrome f–plastocyanin complex as a model to study transient interactions between redox proteins. FEBS Lett. (2011), doi:10.1016/j.febslet.2011.08.035