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

Recent advances in utilization of ferredoxins for biosynthesis of valuable compounds

Ferredoxin (Fd) is a small metalloprotein holding one or two Fe-S clusters in its inner shell. Like many other metalloproteins, Fd is redox active and involved in electron transfer during cellular metabolism. The electrons from reduced Fd are mostly used to regenerate NADPH under physiological conditions. Increasing number of attempts have been reported, however, where Fd delivers electrons to enable biosynthesis of valuable compounds. Various compounds ranging from H₂ to vitamin D₃ have been synthesized successfully using electrons mediated by Fd molecules. In this review, we provide an overview of the engineering studies utilizing Fd for biosynthesis of targeted molecules. The emphasis is on the role and activity of Fd as well as the methods used to improve the rate of electron transfer. Both microbial and electrochemical biosynthesis technologies are described and compared with respect to productivity and the compound being produced. In addition to the ferredoxins from the microbial organisms, artificially designed de novo types are described, highlighting the potential of the emerging computational methods used in metabolic and protein engineering. We believe that the recent advances in utilization of Fd for biosynthesis can result in breakthrough innovation across the biotechnology industry.

Crystal structure of a photosynthetic LH1-RC in complex with its electron donor HiPIP

Photosynthetic electron transfers occur through multiple components ranging from small soluble proteins to large integral membrane protein complexes. Co-crystallization of a bacterial photosynthetic electron transfer complex that employs weak hydrophobic interactions was achieved by using high-molar-ratio mixtures of a soluble donor protein (high-potential iron-sulfur protein, HiPIP) with a membrane-embedded acceptor protein (reaction center, RC) at acidic pH. The structure of the co-complex offers a snapshot of a transient bioenergetic event and revealed a molecular basis for thermodynamically unfavorable interprotein electron tunneling. HiPIP binds to the surface of the tetraheme cytochrome subunit in the light-harvesting (LH1) complex-associated RC in close proximity to the low-potential heme-1 group. The binding interface between the two proteins is primarily formed by uncharged residues and is characterized by hydrophobic features. This co-crystal structure provides a model for the detailed study of long-range trans-protein electron tunneling pathways in biological systems.

The high potential iron-sulfur (HiPIP) proteins are direct electron donors to the light-harvesting-reaction center complexes (LH1-RC) in photosynthetic β- and γ-Proteobacteria. Here, the authors present the 2.9 Å crystal structure of the HiPIP-bound LH1-RC complex from the thermophilic purple sulfur bacterium Thermochromatium tepidum and discuss mechanistic implications for the electron transfer pathway.

Osmotic pressure effects identify dehydration upon cytochrome c-cytochrome c oxidase complex formation contributing to a specific electron pathway formation

In the electron transfer (ET) reaction from cytochrome c (Cyt c) to cytochrome c oxidase (CcO), we determined the number and sites of the hydration water released from the protein surface upon the formation of the ET complex by evaluating the osmotic pressure dependence of kinetics for the ET from Cyt c to CcO. We identified that ∼20 water molecules were dehydrated in complex formation under turnover conditions, and systematic Cyt c mutations in the interaction site for CcO revealed that nearly half of the released hydration water during the complexation were located around Ile81, one of the hydrophobic amino acid residues near the exposed heme periphery of Cyt c. Such a dehydration dominantly compensates for the entropy decrease due to the association of Cyt c with CcO, resulting in the entropy-driven ET reaction. The energetic analysis of the interprotein interactions in the ET complex predicted by the docking simulation suggested the formation of hydrophobic interaction sites surrounding the exposed heme periphery of Cyt c in the Cyt c-CcO interface (a 'molecular breakwater'). Such sites would contribute to the formation of the hydrophobic ET pathway from Cyt c to CcO by blocking water access from the bulk water phase.

Binding Thermodynamics to Intrinsically Disordered Protein Domains

Many proteins are intrinsically disordered or contain one or more disordered domains. These domains can participate in binding interactions with other proteins or small ligands. Binding to intrinsically disordered protein domains requires the folding or structuring of those regions such that they can establish well-defined stoichiometric interactions. Since, in such a situation binding is coupled to folding, the energetics of those two events is reflected in the measured binding thermodynamics. In this protocol, we illustrate the thermodynamic differences between binding coupled to folding and binding independent of folding for the same protein. As an example, we use the HIV-1 envelope glycoprotein gp120 that contains structured as well as disordered domains. In the experiments presented, the binding of gp120 to molecules that bind to disordered regions and trigger structuring (CD4 or MAb 17b) and to molecules that bind to structured regions and do not induce conformational structuring (MAb b12) is discussed.

Thermodynamics of ferredoxin binding to cyanobacterial nitrate reductase

The role of the seven negatively charged amino acids of Synechocystis sp. PCC 6803 ferredoxin (Fd), i.e., Glu29, Glu30, Asp60, Asp65, Asp66, Glu92, and Glu93, predicted to form complex with nitrate reductase (NR), was investigated using site-directed mutagenesis and isothermal titration calorimetry (ITC). These experiments identified four Fd amino acids, i.e., Glu29, Asp60, Glu92, and Glu93, that are essential for the Fd binding and efficient electron transfer to the NR. ITC measurements showed that the most likely stoichiometry for the wild-type NR/wild-type Fd complex is 1:1, a K_d value 4.7 μM for the complex at low ionic strength residues and both the enthalpic and entropic components are associated with complex formation. ITC titrations of wild-type NR with four Fd variants, E29N, D60N, E92Q, and E93N demonstrated that the complex formation, although favorable, was less energetically favorable when compared to complex formation between the two wild-type proteins, suggesting that these negatively charged Fd residues at these positions are important for the effective and productive interaction with wild-type enzyme.

NADP(H) allosterically regulates the interaction between ferredoxin and ferredoxin-NADP+ reductase

Ferredoxin‐NADP⁺ reductase (FNR) in plants receives electrons from ferredoxin (Fd) at the end of the photosynthetic electron transfer chain and converts NADP⁺ to NADPH. The interaction between Fd and FNR in plants was previously shown to be attenuated by NADP(H). Here, we investigated the molecular mechanism of this phenomenon using maize FNR and Fd, as the three‐dimensional structure of this complex is available. NADPH, NADP⁺, and 2′5′‐ADP differentially affected the interaction, as revealed through kinetic and physical binding analyses. Site‐directed mutations of FNR which change the affinity for NADPH altered the affinity for Fd in the opposite direction to that for NADPH. We propose that the binding of NADP(H) causes a conformational change of FNR which is transferred to the Fd‐binding region through different domains of FNR, resulting in allosteric changes in the affinity for Fd.

The interaction between ferredoxin (Fd) and Fd‐NADP⁺ reductase (FNR) in plants is attenuated by NADP(H). Site‐directed mutations of FNR which change the affinity for NADPH altered the affinity for Fd in the opposite direction. We propose that the binding of NADP(H) leads to conformational changes of FNR, resulting in allosteric changes in the affinity for Fd.

Binding of ferredoxin NADP+ oxidoreductase (FNR) to plant photosystem I

The binding of FNR to PSI has been postulated long ago, however, a clear evidence is still missing. In this work, using isothermal titration calorimetry (ITC), we found that FNR binds to photosystem I with its light harvesting complex I (PSI-LHCI) from C. reinhardtii with a 1:1 stoichiometry, a Kd of ~0.8 μM and ∆H of -20.7 kcal/mol. Titrations at different temperatures were used to determine the heat capacity change, ∆CP, of the binding, through which the size of the interface area between the proteins was assessed as ~3000 Ų. In a different set of ITC experiments, introduction of various sucrose concentrations was used to estimate that ~95 water molecules are released to the solvent. These observations support the notion of a binding site shared by few of the photosystem I - light harvesting complex I (PSI-LHCI) subunits in addition to PsaE. Based on these results, a hypothetical model was built for the binding site of FNR at PSI, using known crystallographic structures of: cyanobacterial PSI in complex with ferredoxin (Fd), plant PSI-LHCI and Fd:FNR complex from cyanobacteria. FNR binding site location is proposed to be at the foot of the stromal ridge and above the inner LHCI belt. It is expected to form contacts with PsaE, PsaB, PsaF and at least one of the LHCI. In addition, a ~4.5-fold increased affinity between FNR and PSI-LHCI under crowded 1 M sucrose environment led us to conclude that in C. reinhardtii FNR also functions as a subunit of PSI-LHCI.

Polyethylene glycol promotes autoxidation of cytochrome c

Cytochrome c (Cyt c) was rapidly oxidized by molecular oxygen in the presence, but not absence of PEG. The redox potential of heme c was determined by the potentiometric titration to be +236 ± 3 mV in the absence of PEG, which was negatively shifted to +200 ± 4 mV in the presence of PEG. The underlying the rapid oxidation was explored by examining the structural changes in Cyt c in the presence of PEG using UV-visible absorption, circular dichroism, resonance Raman, and fluorescence spectroscopies. These spectroscopic analyses suggested that heme oxidation was induced by a modest tertiary structural change accompanied by a slight shift in the heme position (<1.0 Å) rather than by partial denaturation, as is observed in the presence of cardiolipin. The near-infrared spectra showed that PEG induced dehydration from Cyt c, which triggered heme displacement. The primary dehydration site was estimated to be around surface-exposed hydrophobic residues near the heme center: Ile81 and Val83. These findings and our previous studies, which showed that hydrated water molecules around Ile81 and Val83 are expelled when Cyt c forms a complex with CcO, proposed that dehydration of these residues is functionally significant to electron transfer from Cyt c to CcO.

Interaction and electron transfer between ferredoxin-NADP+ oxidoreductase and its partners: structural, functional, and physiological implications

Ferredoxin-NADP⁺ reductase (FNR) catalyzes the last step of linear electron transfer in photosynthetic light reactions. The FAD cofactor of FNR accepts two electrons from two independent reduced ferredoxin molecules (Fd) in two sequential steps, first producing neutral semiquinone and then the fully anionic reduced, or hydroquinone, form of the enzyme (FNR_hq). FNR_hq transfers then both electrons in a single hydride transfer step to NADP⁺. We are presenting the recent progress in studies focusing on Fd:FNR interaction and subsequent electron transfer processes as well as on interaction of FNR with NADP⁺/H followed by hydride transfer, both from the structural and functional point of views. We also present the current knowledge about the physiological role(s) of various FNR isoforms present in the chloroplasts of higher plants and the functional impact of subchloroplastic location of FNR. Moreover, open questions and current challenges about the structure, function, and physiology of FNR are discussed.

Aspects of Weak Interactions between Folate and Glycine Betaine

Folate, or vitamin B9, is an important compound in one-carbon metabolism. Previous studies have found weaker binding of dihydrofolate to dihydrofolate reductase in the presence of osmolytes. In other words, osmolytes are more difficult to remove from the dihydrofolate solvation shell than water; this shifts the equilibrium toward the free ligand and protein species. This study uses vapor-pressure osmometry to explore the interaction of folate with the model osmolyte, glycine betaine. This method yields a preferential interaction potential (μ23/RT value). This value is concentration-dependent as folate dimerizes. The μ23/RT value also tracks the deprotonation of folate's N3-O4 keto-enol group, yielding a pKa of 8.1. To determine which folate atoms interact most strongly with betaine, the interaction of heterocyclic aromatic compounds (as well as other small molecules) with betaine was monitored. Using an accessible surface area approach coupled with osmometry measurements, deconvolution of the μ23/RT values into α values for atom types was achieved. This allows prediction of μ23/RT values for larger molecules such as folate. Molecular dynamics simulations of folate show a variety of structures from extended to L-shaped. These conformers possess μ23/RT values from -0.18 to 0.09 m-1, where a negative value indicates a preference for solvation by betaine and a positive value indicates a preference for water. This range of values is consistent with values observed in osmometry and solubility experiments. As the average predicted folate μ23/RT value is near zero, this indicates folate interacts almost equally well with betaine and water. Specifically, the glutamate tail prefers to interact with water, while the aromatic rings prefer betaine. In general, the more protonated species in our small molecule survey interact better with betaine as they provide a source of hydrogens (betaine is not a hydrogen bond donor). Upon deprotonation of the small molecule, the preference swings toward water interaction because of its hydrogen bond donating capacities.

Thermodynamics of Nanobody Binding to Lactose Permease

Camelid nanobodies (Nbs) raised against the outward-facing conformer of a double-Trp mutant of the lactose permease of Escherichia coli (LacY) stabilize the permease in outward-facing conformations. Isothermal titration calorimetry is applied herein to dissect the binding thermodynamics of two Nbs, one that markedly improves access to the sugar-binding site and another that dramatically increases the affinity for galactoside. The findings presented here show that both enthalpy and entropy contribute favorably to binding of the Nbs to wild-type (WT) LacY and that binding of Nb to double-Trp mutant G46W/G262W is driven by a greater enthalpy at an entropic penalty. Thermodynamic analyses support the interpretation that WT LacY is stabilized in outward-facing conformations like the double-Trp mutant with closure of the water-filled cytoplasmic cavity through conformational selection. The LacY conformational transition required for ligand binding is reflected by a favorable entropy increase. Molecular dynamics simulations further suggest that the entropy increase likely stems from release of immobilized water molecules primarily from the cytoplasmic cavity upon closure.

Posttranslational modifications of FERREDOXIN-NADP+ OXIDOREDUCTASE in Arabidopsis chloroplasts

Rapid responses of chloroplast metabolism and adjustments to photosynthetic machinery are of utmost importance for plants' survival in a fluctuating environment. These changes may be achieved through posttranslational modifications of proteins, which are known to affect the activity, interactions, and localization of proteins. Recent studies have accumulated evidence about the crucial role of a multitude of modifications, including acetylation, methylation, and glycosylation, in the regulation of chloroplast proteins. Both of the Arabidopsis (Arabidopsis thaliana) leaf-type FERREDOXIN-NADP(+) OXIDOREDUCTASE (FNR) isoforms, the key enzymes linking the light reactions of photosynthesis to carbon assimilation, exist as two distinct forms with different isoelectric points. We show that both AtFNR isoforms contain multiple alternative amino termini and undergo light-responsive addition of an acetyl group to the α-amino group of the amino-terminal amino acid of proteins, which causes the change in isoelectric point. Both isoforms were also found to contain acetylation of a conserved lysine residue near the active site, while no evidence for in vivo phosphorylation or glycosylation was detected. The dynamic, multilayer regulation of AtFNR exemplifies the complex regulatory network systems controlling chloroplast proteins by a range of posttranslational modifications, which continues to emerge as a novel area within photosynthesis research.

Plant type ferredoxins and ferredoxin-dependent metabolism

Ferredoxin (Fd) is a small [2Fe-2S] cluster-containing protein found in all organisms performing oxygenic photosynthesis. Fd is the first soluble acceptor of electrons on the stromal side of the chloroplast electron transport chain, and as such is pivotal to determining the distribution of these electrons to different metabolic reactions. In chloroplasts, the principle sink for electrons is in the production of NADPH, which is mostly consumed during the assimilation of CO2 . In addition to this primary function in photosynthesis, Fds are also involved in a number of other essential metabolic reactions, including biosynthesis of chlorophyll, phytochrome and fatty acids, several steps in the assimilation of sulphur and nitrogen, as well as redox signalling and maintenance of redox balance via the thioredoxin system and Halliwell-Asada cycle. This makes Fds crucial determinants of the electron transfer between the thylakoid membrane and a variety of soluble enzymes dependent on these electrons. In this article, we will first describe the current knowledge on the structure and function of the various Fd isoforms present in chloroplasts of higher plants and then discuss the processes involved in oxidation of Fd, introducing the corresponding enzymes and discussing what is known about their relative interaction with Fd.

Biochemical characterization of ferredoxin-NADP(+) reductase interaction with flavodoxin in Pseudomonas putida

Flavodoxin (Fld) has been demonstrated to bind to ferredoxin- NADP(+) reductase A (FprA) in Pseudomonas putida. Two residues (Phe(256), Lys(259)) of FprA are likely to be important for interacting with Fld based on homology modeling. Sitedirected mutagenesis and pH-dependent enzyme kinetics were performed to further examine the role of these residues. The catalytic efficiencies of FprA-Ala(259) and FprA-Asp(259) proteins were two-fold lower than those of the wild-type FprA. Homology modeling also strongly suggested that these two residues are important for electron transfer. Thermodynamic properties such as entropy, enthalpy, and heat capacity changes of FprA-Ala(259) and FprA-Asp(259) were examined by isothermal titration calorimetry. We demonstrated, for the first time, that Phe(256) and Lys(259) are critical residues for the interaction between FprA and Fld. Van der Waals interactions and hydrogen bonding were also more important than ionic interactions for forming the FprA-Fld complex.

A new structural insight into differential interaction of cyanobacterial and plant ferredoxins with nitrite reductase as revealed by NMR and X-ray crystallographic studies

Ferredoxin (Fd), which plays a pivotal role in photosynthesis as an electron carrier, forms a transient complex with various Fd-dependent enzymes, such as nitrite reductase (NiR), to achieve efficient intermolecular electron transfer. We studied the protein-protein interaction of Fd and NiR by NMR spectroscopy and determined three acidic regions of Fd to be major sites for the interaction with NiR, indicating that the complex is stabilized through electrostatic interaction. During this study, we found Fds from higher plant and cyanobacterium, in spite of their high structural similarities including the above acidic regions, differ remarkably in the interaction with cyanobacterial NiR. In activity assay of NiR, K(m) value for maize Fd (74.6 µM) was 9.6 times larger than that for Leptolyngbya boryana Fd (7.8 µM). The two Fds also showed a similar difference in binding assay to NiR-immobilized resin. Comparative site-specific mutagenesis of two Fds revealed that their discriminative ability for the interaction with NiR is attributed mainly to non-charged residues in the peripheral region of [2Fe-2S] cluster. These non-charged residues are conserved separately between Fds of plant and cyanobacterial origins. Our data highlight that intermolecular force(s) other than electrostatic attraction is(are) also crucial for the molecular interaction between Fd and partner enzyme.

Evidence from thermodynamics that DNA photolyase recognizes a solvent-exposed CPD lesion

Binding of a cis,syn-cyclobutane pyrimidine dimer (CPD) to Escherichia coli DNA photolyase was examined as a function of temperature, enzyme oxidation state, salt, and substrate conformation using isothermal titration calorimetry. While the overall ΔG° of binding was relatively insensitive to most of the conditions examined, the enthalpic and entropic terms that make up the free energy of binding are sensitive to the conditions of the experiment. Substrate binding to DNA photolyase is generally driven by a negative change in enthalpy. Electrostatic interactions and protonation are affected by the oxidation state of the required FAD cofactor and substrate conformation. The fully reduced enzyme appears to bind approximately two additional water molecules as part of substrate binding. More significantly, the experimental change in heat capacity strongly suggests that the CPD lesion must be flipped out of the intrahelical base stacking prior to binding to the protein; the DNA repair enzyme appears to recognize a solvent-exposed CPD as part of its damage recognition mechanism.

Analysis of the binding forces driving the tight interactions between beta-lactamase inhibitory protein-II (BLIP-II) and class A beta-lactamases

β-Lactamases hydrolyze β-lactam antibiotics to provide drug resistance to bacteria. β-Lactamase inhibitory protein-II (BLIP-II) is a potent proteinaceous inhibitor that exhibits low picomolar affinity for class A β-lactamases. This study examines the driving forces for binding between BLIP-II and β-lactamases using a combination of presteady state kinetics, isothermal titration calorimetry, and x-ray crystallography. The measured dissociation rate constants for BLIP-II and various β-lactamases ranged from 10(-4) to 10(-7) s(-1) and are comparable with those found in some of the tightest known protein-protein interactions. The crystal structures of BLIP-II alone and in complex with Bacillus anthracis Bla1 β-lactamase revealed no significant side-chain movement in BLIP-II in the complex versus the monomer. The structural rigidity of BLIP-II minimizes the loss of the entropy upon complex formation and, as indicated by thermodynamics experiments, may be a key determinant of the observed potent inhibition of β-lactamases.

Isothermal titration calorimetry for measuring macromolecule-ligand affinity

Isothermal titration calorimetry (ITC) is a useful tool for understanding the complete thermodynamic picture of a binding reaction. In biological sciences, macromolecular interactions are essential in understanding the machinery of the cell. Experimental conditions, such as buffer and temperature, can be tailored to the particular binding system being studied. However, careful planning is needed since certain ligand and macromolecule concentration ranges are necessary to obtain useful data. Concentrations of the macromolecule and ligand need to be accurately determined for reliable results. Care also needs to be taken when preparing the samples as impurities can significantly affect the experiment. When ITC experiments, along with controls, are performed properly, useful binding information, such as the stoichiometry, affinity and enthalpy, are obtained. By running additional experiments under different buffer or temperature conditions, more detailed information can be obtained about the system. A protocol for the basic setup of an ITC experiment is given.

Thermodynamics and solvent effects on substrate and cofactor binding in Escherichia coli chromosomal dihydrofolate reductase

Chromosomal dihydrofolate reductase from Escherichia coli catalyzes the reduction of dihydrofolate to tetrahydrofolate using NADPH as a cofactor. The thermodynamics of ligand binding were examined using an isothermal titration calorimetry approach. Using buffers with different heats of ionization, zero to a small, fractional proton release was observed for dihydrofolate binding, while a proton was released upon NADP(+) binding. The role of water in binding was additionally monitored using a number of different osmolytes. Binding of NADP(+) is accompanied by the net release of ∼5-24 water molecules, with a dependence on the identity of the osmolyte. In contrast, binding of dihydrofolate is weakened in the presence of osmolytes, consistent with "water uptake". Different effects are observed depending on the identity of the osmolyte. The net uptake of water upon dihydrofolate binding was previously observed in the nonhomologous R67-encoded dihydrofolate reductase (dfrB or type II enzyme) [Chopra, S., et al. (2008) J. Biol. Chem. 283, 4690-4698]. As R67 dihydrofolate reductase possesses a nonhomologous sequence and forms a tetrameric structure with a single active site pore, the observation of weaker DHF binding in the presence of osmolytes in both enzymes implicates cosolvent effects on free dihydrofolate. Consistent with this analysis, stopped flow experiments find betaine mostly affects DHF binding via changes in k(on), while betaine mostly affects NADPH binding via changes in k(off). Finally, nonadditive enthalpy terms when binary and ternary cofactor binding events are compared suggest the presence of long-lived conformational transitions that are not included in a simple thermodynamic cycle.

Binding thermodynamics of ferredoxin:NADP+ reductase: two different protein substrates and one energetics

The thermodynamics of the formation of binary and ternary complexes between Anabaena PCC 7119 FNR and its substrates, NADP+ and Fd, or Fld, has been studied by ITC. Despite structural dissimilarities, the main difference between Fd and Fld binding to FNR relates to hydrophobicity, reflected in different binding heat capacity and number of water molecules released from the interface. At pH 8, the formation of the binary complexes is both enthalpically and entropically driven, accompanied by the protonation of at least one ionizable group. His299 FNR has been identified as the main responsible for the proton exchange observed. However, at pH 10, where no protonation occurs and intrinsic binding parameters can be obtained, the formation of the binary complexes is entropically driven, with negligible enthalpic contribution. Absence of the FMN cofactor in Fld does not alter significantly the strength of the interaction, but considerably modifies the enthalpic and entropic contributions, suggesting a different binding mode. Ternary complexes show negative cooperativity (6-fold and 11-fold reduction in binding affinity, respectively), and an increase in the enthalpic contribution (more favorable) and a decrease in the entropic contribution (less favorable), with regard to the binary complexes energetics.

Thermodynamics of Protein–Ligand Interactions: History, Presence, and Future Aspects

The understanding of molecular recognition processes of small ligands and biological macromolecules requires a complete characterization of the binding energetics and correlation of thermodynamic data with interacting structures involved. A quantitative description of the forces that govern molecular associations requires determination of changes of all thermodynamic parameters, including free energy of binding (deltaG), enthalpy (deltaH), and entropy (deltaS) of binding and the heat capacity change (deltaCp). A close insight into the binding process is of significant and practical interest, since it provides the fundamental know-how for development of structure-based molecular design-strategies. The only direct method to measure the heat change during complex formation at constant temperature is provided by isothermal titration calorimetry (ITC). With this method one binding partner is titrated into a solution containing the interaction partner, thereby generating or absorbing heat. This heat is the direct observable that can be quantified by the calorimeter. The use of ITC has been limited due to the lack of sensitivity, but recent developments in instrument design permit to measure heat effects generated by nanomol (typically 10-100) amounts of reactants. ITC has emerged as the primary tool for characterizing interactions in terms of thermodynamic parameters. Because heat changes occur in almost all chemical and biochemical processes, ITC can be used for numerous applications, e.g., binding studies of antibody-antigen, protein-peptide, protein-protein, enzyme-inhibitor or enzyme-substrate, carbohydrate-protein, DNA-protein (and many more) interactions as well as enzyme kinetics. Under appropriate conditions data analysis from a single experiment yields deltaH, K(B), the stoichiometry (n), deltaG and deltaS of binding. Moreover, ITC experiments performed at different temperatures yield the heat capacity change (deltaCp). The informational content of thermodynamic data is large, and it has been shown that it plays an important role in the elucidation of binding mechanisms and, through the link to structural data, also in rational drug design. In this review we will present a comprehensive overview to ITC by giving some historical background to calorimetry, outline some critical experimental and data analysis aspects, discuss the latest developments, and give three recent examples of studies published with respect to macromolecule-ligand interactions that have utilized ITC technology.

Identification, characterization and activation mechanism of a tyrosine kinase of Bacillus anthracis

Bacillus subtilis has three active tyrosine kinases, PtkA, PtkB and McsB, which play an important role in the physiology of the bacterium. Genome sequence analysis and biochemical experiments indicated that the ortholog of McsB, BAS0080, is the only active tyrosine kinase present in Bacillus anthracis. The autophosphorylation of McsB of B. anthracis was enhanced in the presence of an activator protein McsA (BAS0079), a property similar to that reported for B. subtilis. However, the process of enhanced phosphorylation of McsB in the presence of McsA remains elusive. To understand the activation mechanism of McsB, we carried out spectroscopic and calorimetric experiments with McsB and McsA. The spectroscopic results suggest that the binding affinity of Mg-ATP for McsB increased by one order from 10(3) to 10(4) in the presence of McsA. The calorimetric experiments revealed that the interaction between McsB and McsA is endothermic in nature, with unfavourable positive enthalpy (DeltaH) and favourable entropy (DeltaS) changes leading to an overall favourable free energy change (DeltaG). Kinetics of binding of both ATP and McsA with McsB showed low association rates (k(a)) and fast dissociation rates (k(d)). These results suggest that enhanced phosphorylation of McsB in the presence of McsA is due to increased affinity of ATP for McsB.

Kinetic analysis of interaction of BRCA1 tandem breast cancer c-terminal domains with phosphorylated peptides reveals two binding conformations

Tandem breast cancer C-terminal (BRCT) domains, present in many DNA repair and cell cycle checkpoint signaling proteins, are phosphoprotein binding modules. The best-characterized tandem BRCT domains to date are from the protein BRCA1 (BRCA1-BRCT), an E3 ubiquitin ligase that has been linked to breast and ovarian cancer. While X-ray crystallography and NMR spectroscopy studies have uncovered the structural determinants of specificity of BRCA1-BRCT for phosphorylated peptides, a detailed kinetic and thermodynamic characterization of the interaction is also required to understand how structure and dynamics are connected and therefore better probe the mechanism of phosphopeptide recognition by BRCT domains. Through a global analysis of binding kinetics data obtained from surface plasmon resonance (SPR) and stopped-flow fluorescence spectroscopy, we show that the recognition mechanism is complex and best modeled by two equilibrium conformations of BRCA1-BRCT in the free state that both interact with a phosphopeptide, with dissociation constants ( K d) in the micromolar range. We show that the apparent global dissociation constant derived from this kinetic analysis is similar to the K d values measured using steady-state SPR, isothermal titration calorimetry, and fluorescence anisotropy. The dynamic nature of BRCA1-BRCT may facilitate the binding of BRCA1 to different phosphorylated protein targets.

Structural and functional diversity of ferredoxin-NADP(+) reductases

Although all ferredoxin-NADP(+) reductases (FNRs) catalyze the same reaction, i.e. the transfer of reducing equivalents between NADP(H) and ferredoxin, they belong to two unrelated families of proteins: the plant-type and the glutathione reductase-type of FNRs. Aim of this review is to provide a general classification scheme for these enzymes, to be used as a framework for the comparison of their properties. Furthermore, we report on some recent findings, which significantly increased the understanding of the structure-function relationships of FNRs, i.e. the ability of adrenodoxin reductase and its homologs to catalyze the oxidation of NADP(+) to its 4-oxo derivative, and the properties of plant-type FNRs from non-photosynthetic organisms. Plant-type FNRs from bacteria and Apicomplexan parasites provide examples of novel ways of FAD- and NADP(H)-binding. The recent characterization of an FNR from Plasmodium falciparum brings these enzymes into the field of drug design.

A balancing act between net uptake of water during dihydrofolate binding and net release of water upon NADPH binding in R67 dihydrofolate reductase

R67 dihydrofolate reductase (DHFR) catalyzes the reduction of dihydrofolate (DHF) to tetrahydrofolate using NADPH as a cofactor. This enzyme is a homotetramer possessing 222 symmetry, and a single active site pore traverses the length of the protein. A promiscuous binding surface can accommodate either DHF or NADPH, thus two nonproductive complexes can form (2NADPH or 2DHF) as well as a productive complex (NADPH.DHF). The role of water in binding was monitored using a number of different osmolytes. From isothermal titration calorimetry (ITC) studies, binding of NADPH is accompanied by the net release of 38 water molecules. In contrast, from both steady state kinetics and ITC studies, binding of DHF is accompanied by the net uptake of water. Although different osmolytes have similar effects on NADPH binding, variable results are observed when DHF binding is probed. Sensitivity to water activity can also be probed by an in vivo selection using the antibacterial drug, trimethoprim, where the water content of the media is decreased by increasing concentrations of sorbitol. The ability of wild type and mutant clones of R67 DHFR to allow host Escherichia coli to grow in the presence of trimethoprim plus added sorbitol parallels the catalytic efficiency of the DHFR clones, indicating water content strongly correlates with the in vivo function of R67 DHFR.

Atomic resolution modeling of the ferredoxin:[FeFe] hydrogenase complex from Chlamydomonas reinhardtii

The [FeFe] hydrogenases HydA1 and HydA2 in the green alga Chlamydomonas reinhardtii catalyze the final reaction in a remarkable metabolic pathway allowing this photosynthetic organism to produce H(2) from water in the chloroplast. A [2Fe-2S] ferredoxin is a critical branch point in electron flow from Photosystem I toward a variety of metabolic fates, including proton reduction by hydrogenases. To better understand the binding determinants involved in ferredoxin:hydrogenase interactions, we have modeled Chlamydomonas PetF1 and HydA2 based on amino-acid sequence homology, and produced two promising electron-transfer model complexes by computational docking. To characterize these models, quantitative free energy calculations at atomic resolution were carried out, and detailed analysis of the interprotein interactions undertaken. The protein complex model we propose for ferredoxin:HydA2 interaction is energetically favored over the alternative candidate by 20 kcal/mol. This proposed model of the electron-transfer complex between PetF1 and HydA2 permits a more detailed view of the molecular events leading up to H(2) evolution, and suggests potential mutagenic strategies to modulate electron flow to HydA2.

Calorimetric analysis of binding of two consecutive DNA strands to RecA protein illuminates mechanism for recognition of homology

RecA protein recognises two complementary DNA strands for homologous recombination. To gain insight into the molecular mechanism, the thermodynamic parameters of the DNA binding have been characterised by isothermal calorimetry. Specifically, conformational changes of protein and DNA were searched for by measuring variations in enthalpy change (DeltaH) with temperature (heat capacity change, DeltaC(p)). In the presence of the ATP analogue ATPgammaS, the DeltaH for the binding of the first DNA strand depends upon temperature (large DeltaC(p)) and the type of buffer, in a way that is consistent with the organisation of disordered parts and the protonation of RecA upon complex formation. In contrast, the binding of the second DNA strand occurs without any pronounced DeltaC(p), indicating the absence of further reorganisation of the RecA-DNA filament. In agreement with these findings, a significant change in the CD spectrum of RecA was observed only upon the binding of the first DNA strand. In the absence of nucleotide cofactor, the DeltaH of DNA binding is almost independent of temperature, indicating a requirement for ATP in the reorganisation of RecA. When the second DNA strand is complementary to the first, the DeltaH is larger than that for non-complementary DNA strand, but less than the DeltaH of the annealing of the complementary DNA without RecA. This small DeltaH could reflect a weak binding that may facilitate the dissociation of only partly complementary DNA and thus speed the search for complementary DNA. The DeltaH of binding DNA sequences displaying strong base-base stacking is small for both the first and second binding DNA strand, suggesting that the second is also stretched upon interaction with RecA. These results support the proposal that the RecA protein restructures DNA, preparing it for the recognition of a complementary second DNA strand, and that the recognition is due mainly to direct base-base contacts between DNA strands.

Thermodynamics of tryptophan-mediated activation of the trp RNA-binding attenuation protein

The trp RNA-binding attenuation protein (TRAP) functions in many bacilli to control the expression of the tryptophan biosynthesis genes. Transcription of the trp operon is controlled by TRAP through an attenuation mechanism, in which competition between two alternative secondary-structural elements in the 5' leader sequence of the nascent mRNA is influenced by tryptophan-dependent binding of TRAP to the RNA. Previously, NMR studies of the undecamer (11-mer) suggested that tryptophan-dependent control of RNA binding by TRAP is accomplished through ligand-induced changes in protein dynamics. We now present further insights into this ligand-coupled event from hydrogen/deuterium (H/D) exchange analysis, differential scanning calorimetry (DSC), and isothermal titration calorimetry (ITC). Scanning calorimetry showed tryptophan dissociation to be independent of global protein unfolding, while analysis of the temperature dependence of the binding enthalpy by ITC revealed a negative heat capacity change larger than expected from surface burial, a hallmark of binding-coupled processes. Analysis of this excess heat capacity change using parameters derived from protein folding studies corresponds to the ordering of 17-24 residues per monomer of TRAP upon tryptophan binding. This result is in agreement with qualitative analysis of residue-specific broadening observed in TROSY NMR spectra of the 91 kDa oligomer. Implications for the mechanism of ligand-mediated TRAP activation through a shift in a preexisting conformational equilibrium and an induced-fit conformational change are discussed.

Exact analysis of heterotropic interactions in proteins: Characterization of cooperative ligand binding by isothermal titration calorimetry

Intramolecular interaction networks in proteins are responsible for heterotropic ligand binding cooperativity, a biologically important, widespread phenomenon in nature (e.g., signaling transduction cascades, enzymatic cofactors, enzymatic allosteric activators or inhibitors, gene transcription, or repression). The cooperative binding of two (or more) different ligands to a macromolecule is the underlying principle. To date, heterotropic effects have been studied mainly kinetically in enzymatic systems. Until now, approximate approaches have been employed for studying equilibrium heterotropic ligand binding effects, except in two special cases in which an exact analysis was developed: independent binding (no cooperativity) and competitive binding (maximal negative cooperativity). The exact analysis and methodology for characterizing ligand binding cooperativity interactions in the general case (any degree of cooperativity) using isothermal titration calorimetry are presented in this work. Intramolecular interaction pathways within the allosteric macromolecule can be identified and characterized using this methodology. As an example, the thermodynamic characterization of the binding interaction between ferredoxin-NADP+ reductase and its three substrates, NADP+, ferredoxin, and flavodoxin, as well as the characterization of their binding cooperativity interaction, is presented.

Structural analysis of interactions for complex formation between Ferredoxin-NADP+ reductase and its protein partners

The three-dimensional structures of K72E, K75R, K75S, K75Q, and K75E Anabaena Ferredoxin-NADP+ reductase (FNR) mutants have been solved, and particular structural details of these mutants have been used to assess the role played by residues 72 and 75 in optimal complex formation and electron transfer (ET) between FNR and its protein redox partners Ferredoxin (Fd) and Flavodoxin (Fld). Additionally, because there is no structural information available on the interaction between FNR and Fld, a model for the FNR:Fld complex has also been produced based on the previously reported crystal structures and on that of the rat Cytochrome P450 reductase (CPR), onto which FNR and Fld have been structurally aligned, and those reported for the Anabaena and maize FNR:Fd complexes. The model suggests putative electrostatic and hydrophobic interactions between residues on the FNR and Fld surfaces at the complex interface and provides an adequate orientation and distance between the FAD and FMN redox centers for efficient ET without the presence of any other molecule as electron carrier. Thus, the models now available for the FNR:Fd and FNR:Fld interactions and the structures presented here for the mutants at K72 and K75 in Anabaena FNR have been evaluated in light of previous biochemical data. These structures confirm the key participation of residue K75 and K72 in complex formation with both Fd and Fld. The drastic effect in FNR activity produced by replacement of K75 by Glu in the K75E FNR variant is explained not only by the observed changes in the charge distribution on the surface of the K75E FNR mutant, but also by the formation of a salt bridge interaction between E75 and K72 that simultaneously "neutralizes" two essential positive charged side chains for Fld/Fd recognition.

Herbicide binding and thermal stability of photosystem II isolated from Thermosynechococcus elongatus

Binding of herbicides to photosystem II inhibits the electron transfer from Q(A) to Q(B) due to competition of herbicides with plastoquinone bound at the Q(B) site. We investigated herbicide binding to monomeric and dimeric photosystem II core complexes (PSIIcc) isolated from Thermosynechococcus elongatus by a combination of different methods (isothermal titration and differential scanning calorimetry, CD spectroscopy and measurements of the oxygen evolution) yielding binding constants, enthalpies and stoichiometries for various herbicides as well as information regarding stabilization/destabilization of the complex. Herbicide binding to detergent-solubilized PSIIcc can be described by a model of single independent binding sites present on this important membrane protein. Interestingly, binding stoichiometries herbicide:PSIIcc are lower than 1:1 and vary depending on the herbicide under study. Strong binding herbicides such as terbutryn stabilize PSIIcc in thermal unfolding experiments and endothermically binding herbicides like ioxynil probably cause large structural changes accompanied with the binding process as shown by differential scanning calorimetry experiments of the unfolding reaction of PSIIcc monomer in the presence of ioxynil. In addition we studied the occupancy of the Q(B) sites with plastoquinone (PQ9) by measuring flash induced fluorescence relaxation yielding a possible explanation for the deviations of herbicide binding from a 1:1 herbicide/binding site model.

Glutamate synthase: structural, mechanistic and regulatory properties, and role in the amino acid metabolism

Ammonium ion assimilation constitutes a central metabolic pathway in many organisms, and glutamate synthase, in concert with glutamine synthetase (GS, EC 6.3.1.2), plays the primary role of ammonium ion incorporation into glutamine and glutamate. Glutamate synthase occurs in three forms that can be distinguished based on whether they use NADPH (NADPH-GOGAT, EC 1.4.1.13), NADH (NADH-GOGAT, EC 1.4.1.14) or reduced ferredoxin (Fd-GOGAT, EC 1.4.7.1) as the electron donor for the (two-electron) conversion of L-glutamine plus 2-oxoglutarate to L-glutamate. The distribution of these three forms of glutamate synthase in different tissues is quite specific to the organism in question. Gene structures have been determined for Fd-, NADH- and NADPH-dependent glutamate synthases from different organisms, as shown by searches in nucleic acid sequence data banks. Fd-glutamate synthase contains two electron-carrying prosthetic groups, the redox properties of which are discussed. A description of the ferredoxin binding by Fd-glutamate synthase is also presented. In plants, including nitrogen-fixing legumes, Fd-glutamate synthase and NADH-glutamate synthase supply glutamate during the nitrogen assimilation and translocation. The biological functions of Fd-glutamate synthase and NADH-glutamate synthase, which show a highly tissue-specific distribution pattern, are tightly related to the regulation by the light and metabolite sensing systems. Analysis of mutants and transgenic studies have provided insights into the primary individual functions of Fd-glutamate synthase and NADH-glutamate synthase. These studies also provided evidence that glutamate dehydrogenase (NADH-GDH, EC 1.4.1.2) does not represent a significant alternate route for glutamate formation in plants. Taken together, biochemical analysis and genetic and molecular data imply that Fd-glutamate synthase incorporates photorespiratory and non-photorespiratory ammonium and provides nitrogen for transport to maintain nitrogen status in plants. Fd-glutamate synthase also plays a role that is redundant, in several important aspects, to that played by NADH-glutamate synthase in ammonium assimilation and nitrogen transport.

Positive contribution of hydration on DNA binding by E2c protein from papillomavirus

Protein-nucleic acid interactions are responsible for the regulation of key biological events such as genomic transcription and recombination and viral replication. However, the recognition mechanisms involved in these processes are not completely understood. Here, we investigate the dominant forces involved in protein-protein and protein-DNA interactions for the 80-amino-acid C-terminal domain of the E2 protein (E2c) from human papillomavirus (HPV-16). The E2c protein is a homodimer that specifically binds to double-stranded DNA containing the consensus sequence ACCG-N(4)-CGGT, where N is any nucleotide. DNA binding affinity is reduced by lowering water chemical potential, accompanied by an increase in cooperativity. Wyman linkage relations between affinity and water chemical potential indicate that 11 additional water molecules are bound in the formation of the complex between E2c and DNA. Salt dissociation isotherms showed that 10 counterions are released upon association, even at low water activity, indicating that this latter variable does not change the electrostatic component of the interaction. Further analysis demonstrates a strong dependence of cooperativity of binding on the protein concentration. Altogether, these results reveal a novel binding pathway in which the consolidated complex may achieve its final form via a monomer-DNA intermediate, which favors the binding of a second monomer. This molecular mechanism reveals the contributions of multiple conformers in a tight virus genome modulation that seems to be important in the cell infection scenario.

Thermodynamic characterization of the interaction of human cyclophilin 18 with cyclosporin A

Isothermal titration calorimetry (ITC) was used to investigate thermodynamic parameters of the cyclosporin A (CsA)-cyclophilin 18 (hCyp18) association reaction. We have calculated the thermodynamic parameters (enthalpy, entropy, heat capacity, and free energy of binding) of the CsA/hCyp18 complexation. All but two methods described in the literature underestimate the affinity to hCyp18 of CsA. We found that the association constant (1.1.10(8) M(-1) at 10 degrees C) of CsA to hCyp18 is in close agreement with the reciprocal of the reported inhibitory constant of the peptidylprolyl cis/trans isomerase activity of hCyp18. Interpretation of the thermodynamic parameters in buffered solution of water, 30% glycerol and D(2)O leads to the conclusion that the highly specific binding of CsA to hCyp18 is mainly mediated through hydrogen bonding and to a lesser degree through hydrophobic interaction. Furthermore, the pH dependence of the association constant was determined and analyzed according to a single proton linkage model, resulting in a pK(a) value of 5.7 in free hCyp18 and below 4.5 in the CsA complexed form. Titration experiments using different single component buffers possessing different heats of ionization allowed us to estimate that statistically half a proton is transferred upon CsA binding from the binding interface of hCyp18 to the buffer at pH 5.5. No proton transfer was detected at pH 7.5. The thermodynamic results are discussed in relation to the published X-ray and NMR structure of the free and CsA complexed hCyp18.

Structure-function relationships in Anabaena ferredoxin/ferredoxin:NADP(+) reductase electron transfer: insights from site-directed mutagenesis, transient absorption spectroscopy and X-ray crystallography

The interaction between reduced Anabaena ferredoxin and oxidized ferredoxin:NADP(+) reductase (FNR), which occurs during photosynthetic electron transfer (ET), has been investigated extensively in the authors' laboratories using transient and steady-state kinetic measurements and X-ray crystallography. The effect of a large number of site-specific mutations in both proteins has been assessed. Many of the mutations had little or no effect on ET kinetics. However, non-conservative mutations at three highly conserved surface sites in ferredoxin (F65, E94 and S47) caused ET rate constants to decrease by four orders of magnitude, and non-conservative mutations at three highly conserved surface sites in FNR (L76, K75 and E301) caused ET rate constants to decrease by factors of 25-150. These residues were deemed to be critical for ET. Similar mutations at several other conserved sites in the two proteins (D67 in Fd; E139, L78, K72, and R16 in FNR) caused smaller but still appreciable effects on ET rate constants. A strong correlation exists between these results and the X-ray crystal structure of an Anabaena ferredoxin/FNR complex. Thus, mutations at sites that are within the protein-protein interface or are directly involved in interprotein contacts generally show the largest kinetic effects. The implications of these results for the ET mechanism are discussed.

The effects of osmotic and hydrostatic pressures on macromolecular systems

Osmotic pressure and hydrostatic pressure can be used effectively to probe the behavior of biologically important macromolecules and their complexes. Using the two techniques requires a theoretical framework as well as knowledge of the more common pitfalls. Both are discussed in this review in the context of several examples.

Rapid and selective binding to the synaptic SNARE complex suggests a modulatory role of complexins in neuroexocytosis

The Ca(2+)-triggered release of neurotransmitters is mediated by fusion of synaptic vesicles with the plasma membrane. The molecular machinery that translates the Ca(2+) signal into exocytosis is only beginning to emerge. The soluble N-ethylmaleimide-sensitive factor attachment protein receptor (SNARE) proteins syntaxin, SNAP-25, and synaptobrevin are central components of the fusion apparatus. Assembly of a membrane-bridging ternary SNARE complex is thought to initiate membrane merger, but the roles of other factors are less understood. Complexins are two highly conserved proteins that modulate the Ca(2+) responsiveness of neurotransmitter release. In vitro, they bind in a 1:1 stoichiometry to the assembled synaptic SNARE complex, making complexins attractive candidates for controlling the exocytotic fusion apparatus. We have now performed a detailed structural, kinetic, and thermodynamic analysis of complexin binding to the SNARE complex. We found that no major conformational changes occur upon binding and that the complexin helix is aligned antiparallel to the four-helix bundle of the SNARE complex. Complexins bound rapidly (approximately 5 x 10(7) m(-1) s(-1)) and with high affinity (approximately 10 nm), making it one of the fastest protein-protein interactions characterized so far in membrane trafficking. Interestingly, neither affinity nor binding kinetics was substantially altered by Ca(2+) ions. No interaction of complexins was detectable either with individual SNARE proteins or with the binary syntaxin x SNAP-25 complex. Furthermore, complexin did not promote the formation of SNARE complex oligomers. Together, our data suggest that complexins modulate neuroexocytosis after assembly of membrane-bridging SNARE complexes.

Binding specificity and the ligand dissociation process in the E. coli biotin holoenzyme synthetase

The binding of the Escherichia coli biotin holoenzyme synthetase to the two ligands, biotin and bio-5'-AMP, is coupled to disorder-to-order transitions in the protein. In the structure of the biotin complex, a "glycine-rich" loop that is disordered in the apo-enzyme is folded over the ligand. Mutations in three residues in this loop result in significant changes in the affinity of the enzyme for both biotin and bio-5'-AMP. The kinetic basis of these losses in the affinity resides primarily in changes in the unimolecular rates of dissociation of the complexes. In this work, isothermal titration calorimetry has been employed to examine the detailed thermodynamics of binding of three loop mutants to biotin and bio-5'-AMP. The energetic features of dissociation of the protein*ligand complexes also have been probed by measuring the temperature dependencies of the unimolecular dissociation rates. Analysis of the data using the Eyring formalism yielded entropic and enthalpic contributions to the energetic barrier to dissociation. The thermodynamic results coupled with the known structures of the apo-enzyme and biotin complex have been used to formulate a model for progression from the ground-state complex to the transition state in biotin dissociation. In this model, the transition-state is characterized by both partial disruption of noncovalent bonds and acquisition of some of the disorder that characterizes the glycine-rich loop in the absence of ligand.

Thermodynamic analysis of the binding of glutathione to glutathione S-transferase over a range of temperatures

The binding properties of a glutathione S-transferase (EC 2.5.1.18) from Schistosoma japonicum to substrate glutathione (GSH) has been investigated by intrinsic fluorescence and isothermal titration calorimetry (ITC) at pH 6.5 over a temperature range of 15-30 degrees C. Calorimetric measurements in various buffer systems with different ionization heats suggest that protons are released during the binding of GSH at pH 6.5. We have also studied the effect of pH on the thermodynamics of GSH-GST interaction. The behaviour shown at different pHs indicates that at least three groups must participate in the exchange of protons. Fluorimetric and calorimetric measurements indicate that GSH binds to two sites in the dimer of 26-kDa glutathione S-transferase from Schistosoma japonicum (SjGST). On the other hand, noncooperativity for substrate binding to SjGST was detected over a temperature range of 15-30 degrees C. Among thermodynamic parameters, whereas DeltaG degrees remains practically invariant as a function of temperature, DeltaH and DeltaS degrees both decrease with an increase in temperature. While the binding is enthalpically favorable at all temperatures studied, at temperatures below 25 degrees C, DeltaG degrees is also favoured by entropic contributions. As the temperature increases, the entropic contributions progressively decrease, attaining a value of zero at 24.3 degrees C, and then becoming unfavorable. During this transition, the enthalpic contributions become progressively favorable, resulting in an enthalpy-entropy compensation. The temperature dependence of the enthalpy change yields the heat capacity change (DeltaCp degrees ) of -0.238 +/- 0.04 kcal per K per mol of GSH bound.

Role of a cluster of hydrophobic residues near the FAD cofactor in Anabaena PCC 7119 ferredoxin-NADP+ reductase for optimal complex formation and electron transfer to ferredoxin

In the ferredoxin-NADP(+) reductase (FNR)/ferredoxin (Fd) system, an aromatic amino acid residue on the surface of Anabaena Fd, Phe-65, has been shown to be essential for the electron transfer (ET) reaction. We have investigated further the role of hydrophobic interactions in complex stabilization and ET between these proteins by replacing three hydrophobic residues, Leu-76, Leu-78, and Val-136, situated on the FNR surface in the vicinity of its FAD cofactor. Whereas neither the ability of FNR to accept electrons from NADPH nor its structure appears to be affected by the introduced mutations, different behaviors with Fd are observed. Thus, the ET interaction with Fd is almost completely lost upon introduction of negatively charged side chains. In contrast, only subtle changes are observed upon conservative replacement. Introduction of Ser residues produces relatively sizable alterations of the FAD redox potential, which can explain the modified behavior of these mutants. The introduction of bulky aromatic side chains appears to produce rearrangements of the side chains at the FNR/Fd interaction surface. Thus, subtle changes in the hydrophobic patch influence the rates of ET to and from Fd by altering the binding constants and the FAD redox potentials, indicating that these residues are especially important in the binding and orientation of Fd for efficient ET. These results are consistent with the structure reported for the Anabaena FNR.Fd complex.

The role of water molecules in the association of cytochrome P450cam with putidaredoxin. An osmotic pressure study

We have investigated the osmotic pressure dependence of the association between ferric cytochrome P450cam and putidaredoxin (Pdx) to gain an insight into the role of water molecules in the P450cam-reduced Pdx complexation amenable to physiological electron transfer. The association constant was evaluated from the electron transfer rates from reduced Pdx to P450cam. The natural logarithm of the association constant K(a) was linearly reduced by the osmotic pressure, and osmotic stress yields uptake of 25 waters upon association. In contrast, uptake of only 13 waters is observed from the osmotic pressure dependence of the association in the nonphysiological redox partners P450cam and oxidized Pdx. Although general protein-protein associations proceed through dehydration around the complex interface, the interfacial waters could mediate hydrogen-bonding interactions. Therefore, about 10 more interfacial waters imply an additional water-mediated hydrogen-bonding network in the P450cam.reduced Pdx complex, which does not exist in the complex with oxidized Pdx. It is also possible that the water-mediated hydrogen-bonding interactions support a high P450cam affinity for reduced (K(a) = 0.83 microm(-1)) relative to oxidized (K(a) = 0.058 microm(-1)) Pdx. This study points to a novel role of solvents in assisting redox state-dependent interaction between P450cam and Pdx.

A redox-dependent interaction between two electron-transfer partners involved in photosynthesis

Ferredoxin:NADP+:reductase (FNR) catalyzes one terminal step of the conversion of light energy into chemical energy during photosynthesis. FNR uses two high energy electrons photoproduced by photosystem I (PSI) and conveyed, one by one, by a ferredoxin (Fd), to reduce NADP+ to NADPH. The reducing power of NADPH is finally involved in carbon assimilation. The interaction between oxidized FNR and Fd was studied by crystallography at 2.4 A resolution leading to a three-dimensional picture of an Fd-FNR biologically relevant complex. This complex suggests that FNR and Fd specifically interact prior to each electron transfer and disassemble upon a redox-linked conformational change of the Fd.

Electrostatic forces involved in orienting Anabaena ferredoxin during binding to Anabaena ferredoxin:NADP+ reductase: site-specific mutagenesis, transient kinetic measurements, and electrostatic surface potentials

Transient absorbance measurements following laser flash photolysis have been used to measure the rate constants for electron transfer (et) from reduced Anabaena ferredoxin (Fd) to wild-type and seven site-specific charge-reversal mutants of Anabaena ferredoxin:NADP+ reductase (FNR). These mutations have been designed to probe the importance of specific positively charged amino acid residues on the surface of the FNR molecule near the exposed edge of the FAD cofactor in the protein-protein interaction during et with Fd. The mutant proteins fall into two groups: overall, the K75E, R16E, and K72E mutants are most severely impaired in et, and the K138E, R264E, K290E, and K294E mutants are impaired to a lesser extent, although the degree of impairment varies with ionic strength. Binding constants for complex formation between the oxidized proteins and for the transient et complexes show that the severity of the alterations in et kinetics for the mutants correlate with decreased stabilities of the protein-protein complexes. Those mutated residues, which show the largest effects, are located in a region of the protein in which positive charge predominates, and charge reversals have large effects on the calculated local surface electrostatic potential. In contrast, K138, R264, K290, and K294 are located within or close to regions of intense negative potential, and therefore the introduction of additional negative charges have considerably smaller effects on the calculated surface potential. We attribute the relative changes in et kinetics and complex binding constants for these mutants to these characteristics of the surface charge distribution in FNR and conclude that the positively charged region of the FNR surface located in the vicinity of K75, R16, and K72 is especially important in the binding and orientation of Fd during electron transfer.

Mapping EBNA-1 domains involved in binding to metaphase chromosomes

The Epstein-Barr virus (EBV) genome can persist in dividing human B cells as multicopy circular episomes. Viral episomes replicate in synchrony with host cell DNA and are maintained at a relatively constant copy number for a long time. Only two viral elements, the replication origin OriP and the EBNA-1 protein, are required for the persistence of viral genomes during latency. EBNA-1 activates OriP during the S phase and may also contribute to the partition and/or retention of viral genomes during mitosis. Indeed, EBNA-1 has been shown to interact with mitotic chromatin. Moreover, viral genomes are noncovalently associated with metaphase chromosomes. This suggests that EBNA-1 may facilitate the anchorage of viral genomes on cellular chromosomes, thus ensuring proper partition and retention. In the present paper, we have investigated the chromosome-binding activity of EBV EBNA-1, herpesvirus papio (HVP) EBNA-1, and various derivatives of EBV EBNA-1, fused to a variant of the green fluorescent protein. The results show that binding to metaphase chromosomes is a common property of EBV and HVP EBNA-1. Further studies indicated that at least three independent domains (CBS-1, -2, and -3) mediate EBNA-1 binding to metaphase chromosomes. In agreement with the anchorage model, two of these domains mapped to a region that has been previously demonstrated to be required for the long-term persistence of OriP-containing plasmids.

Isothermal titration calorimetry and differential scanning calorimetry as complementary tools to investigate the energetics of biomolecular recognition

The principles of isothermal titration calorimetry (ITC) and differential scanning calorimetry (DSC) are reviewed together with the basic thermodynamic formalism on which the two techniques are based. Although ITC is particularly suitable to follow the energetics of an association reaction between biomolecules, the combination of ITC and DSC provides a more comprehensive description of the thermodynamics of an associating system. The reason is that the parameters DeltaG, DeltaH, DeltaS, and DeltaCp obtained from ITC are global properties of the system under study. They may be composed to varying degrees of contributions from the binding reaction proper, from conformational changes of the component molecules during association, and from changes in molecule/solvent interactions and in the state of protonation.

Protein-protein interaction in electron transfer reactions: the ferredoxin/flavodoxin/ferredoxin:NADP+ reductase system from Anabaena

Electron transfer reactions involving protein-protein interactions require the formation of a transient complex which brings together the two redox centres exchanging electrons. This is the case for the flavoprotein ferredoxin:NADP+ reductase (FNR) from the cyanobacterium Anabaena, an enzyme which interacts with ferredoxin in the photosynthetic pathway to receive the electrons required for NADP+ reduction. The reductase shows a concave cavity in its structure into which small proteins such as ferredoxin can fit. Flavodoxin, an FMN-containing protein that is synthesised in cyanobacteria under iron-deficient conditions, plays the same role as ferredoxin in its interaction with FNR in spite of its different structure, size and redox cofactor. There are a number of negatively charged amino acid residues on the surface of ferredoxin and flavodoxin that play a role in the electron transfer reaction with the reductase. Thus far, in only one case has charge replacement of one of the acidic residues produced an increase in the rate of electron transfer, whereas in several other cases a decrease in the rate is observed. In the most dramatic example, replacement of Glu at position 94 of Anabaena ferredoxin results in virtually the complete loss of ability to transfer electrons. Charge-reversal of positively charged amino acid residues in the reductase also produces strong effects on the rate of electron transfer. Several degrees of impairment have been observed, the most significant involving a positively charged Lys at position 75 which appears to be essential for the stability of the complex between the reductase and ferredoxin. The results presented in this paper provide a clear demonstration of the importance of electrostatic interactions on the stability of the transient complex formed during electron transfer by the proteins presently under study.

Roles of negatively charged surface residues of putidaredoxin in interactions with redox partners in p450cam monooxygenase system

To investigate the interaction of putidaredoxin (Pdx) with its redox partners in the cytochrome P450cam system, we focused on the role of negatively charged surface amino acid residues. The amino acid residues we examined in this mutational study are Asp-58, Glu-65, Glu-72, and Glu-77, which are located on the alpha-helical segment to form a negatively charged region on the surface of Pdx and have been supposed to play key roles in the association with the redox partners, NADH-putidaredoxin reductase (PdR) and P450cam. The neutralization of the single negative charge on these amino acid residues did not significantly inhibit the electron-transfer reaction with the redox partners, except for the mutation at Glu-72. Together with the previous results, we can conclude that the negatively charged cluster on the alpha-helical segment is not so crucial for the electron transfer of the Pdx/PdR complex, and, instead of the negative charges, the steric hindrance is essential for the binding of Pdx with PdR. In the electron transfer from Pdx to P450cam, the alpha-helical region would not be included in the binding site with P450cam and some specific hydrogen bonds on the surface loop near the Fe-S center contribute to the electron transfer to P450cam. Such different binding sites and interactions for Pdx will shed light on the electron-transfer mechanism mediated by Pdx, the shuttle mechanism.