pH-dependent charge equilibria between tyrosine-D and the S states in photosystem II. Estimation of relative midpoint redox potentials

The effect of protonation events on the charge equilibrium between tyrosine-D and the water-oxidizing complex in photosystem II has been studied by time-resolved measurements of the EPR signal IIslow at room temperature. The flash-induced oxidation of YD by the water-oxidizing complex in the S2 state is a monophasic process above pH 6.5 and biphasic at lower pHs, showing a slow and a fast phase. The half-time of the slow phase increases from about 1 s at pH 8.0 to about 20 s at pH 5.0, whereas the half-time of the fast phase is pH independent (0.4-1 s). The dark reduction of YD+ was followed by measuring the decay of signal IIslow at room temperature. YD+ decays in a biphasic way on the tens of minutes to hours time scale. The minutes phase is due to the electron transfer to YD+ from the S0 state of the water-oxidizing complex. The half-time of this process increases from about 5 min at pH 8.0 to 40 min at pH 4.5. The hours phase of YD+ has a constant half-time of about 500 min between pH 4.7 and 7.2, which abruptly decreases above pH 7.2 and below pH 4.7. This phase reflects the reduction of YD+ either from the medium or by an unidentified redox component of PSII in those centers that are in the S1 state. The titration curve of the half-times for the oxidation of YD reveals a proton binding with a pK around 7.3-7.5 that retards the electron transfer from YD to the water-oxidizing complex. We propose that this monoprotic event reflects the protonation of an amino acid residue, probably histidine-190 on the D2 protein, to which YD is hydrogen bonded. The titration curves for the oxidation of YD and for the reduction of YD+ show a second proton binding with pK approximately 5.8-6.0 that accelerates the electron transfer from YD to the water-oxidizing complex and retards the process in the opposite direction. This protonation most probably affects the water-oxidizing complex. From the measured kinetic parameters, the lowest limits for the equilibrium constants between the S0YD+ and the S1YD as well as between the S1YD+ and S2YD states were estimated to be 5 and 750-1000, respectively.(ABSTRACT TRUNCATED AT 400 WORDS)     

pH-dependent redox potential: how to use it correctly in the activation energy analysis

The activation barrier (the activation free energy) for the reaction's elementary act proper does not depend on the presence of reactants outside the reaction complex. The barrier is determined directly by the concentration-independent configurational free energy. In the case of redox reactants with pH-dependent redox potential, only the pH-independent quantity, the configurational redox potential enters immediately into expression for activation energy. Some typical cases of such reactions have been discussed (e.g., simultaneous proton and electron detachment, acid dissociation followed by oxidation, dissociation after oxidation, and others). For these mechanisms, the algorithms for calculation of the configurational redox potential from the experimentally determined redox potentials have been described both for the data related to a dissolved reactant or to a prosthetic group of an enzyme. Some examples of pH-dependent enzymatic redox reactions, in particular for the Rieske iron-sulfur protein, have been discussed.     

pH-dependent redox potential: how to use it correctly in the activation energy analysis

The activation barrier (the activation free energy) for the reaction's elementary act proper does not depend on the presence of reactants outside the reaction complex. The barrier is determined directly by the concentration-independent configurational free energy. In the case of redox reactants with pH-dependent redox potential, only the pH-independent quantity, the configurational redox potential enters immediately into expression for activation energy. Some typical cases of such reactions have been discussed (e.g., simultaneous proton and electron detachment, acid dissociation followed by oxidation, dissociation after oxidation, and others). For these mechanisms, the algorithms for calculation of the configurational redox potential from the experimentally determined redox potentials have been described both for the data related to a dissolved reactant or to a prosthetic group of an enzyme. Some examples of pH-dependent enzymatic redox reactions, in particular for the Rieske iron-sulfur protein, have been discussed.     

pH-dependent redox activity and fluxionality of the copper site in amicyanin from Thiobacillus yersutus as studied by 300- and 600- MHz 1H NMR

The kinetics of the deuteronation of one of the copper ligand histidines of the reduced Type I blue-copper protein amicyanin from Thiobacillus versutus was studied as a function of temperature by 300- and 600- MHz 1H NMR. The NMR data were analyzed with the help of a three site exchange model. Deuteron exchange between the histidine ligand and the solution appears to be catalyzed by phosphate. After deuteronation the histidine can occur in two conformations. The electron self-exchange rate of amicyanin was determined as a function of temperature and ionic strength. At 298 K, pD = 8.6, I = 0.05 M, the ese rate amounts to 1.3 x 10(5) M-1 S-1. The activation parameters amount to delta H not equal to = (52 +/- 3) kJ/mol and delta S not equal to = (26 +/- 9) J/mol.K. The dependence of the ese rate on ionic strength is small. The deuteronated amicyanin appears to be redox-inactive. The experimental findings clearly distinguish amicyanin from other classes of blue-copper proteins like the azurins and the pseudo-azurins.     

The effect of redox state on the folding free energy of azurin

 The unfolding of oxidized and reduced azurin by guanidine hydrochloride has been monitored by circular dichroism. Dilution experiments showed the unfolding to be reversible, and the equilibrium data have been interpreted in terms of a two-state model. The protein is stabilized by the strong metal binding in the native state, so that the folding free energy is as high as –52.2 kJ mol^–1 for the oxidized protein. The reduced protein is less stable, with a folding free energy of –40.0 kJ mol^–1. A thermodynamic cycle shows, as a consequence, that unfolded azurin has a reduction potential 0.13 V above that of the folded protein. This is explained by the bipyramidal site in the native fold stabilizing Cu(II) by a rack mechanism, with the same geometry being maintained in the Cu(I) form. In the unfolded protein, on the other hand, the coordination geometries are expected to differ for the two oxidation states, Cu(I) being stabilized by the cysteine thiol group in a linear or trigonal symmetry, whereas Cu(II) prefers oxygen ligands in a tetragonal geometry. Received: 15 January 1997 / Accepted: 3 April 1997     

Conformational equilibria and free energy profiles for the allosteric transition of the ribose-binding protein

The ribose-binding protein (RBP) is a sugar-binding bacterial periplasmic protein whose function is associated with a large allosteric conformational change from an open to a closed conformation upon binding to ribose. The crystal structures of RBP in open and closed conformations have been solved. It has been hypothesized that the open and closed conformations exist in a dynamic equilibrium in solution, and that sugar binding shifts the population from open conformations to closed conformations. Here, we study by computer simulations the thermodynamic changes that accompany this conformational change, and model the structural changes that accompany the allosteric transition, using umbrella sampling molecular dynamics and the weighted histogram analysis method. The open state is comprised of a diverse ensemble of conformations; the open ribose-free X-ray crystal conformations being representative of this ensemble. The unligated open form of RBP is stabilized by conformational entropy. The simulations predict detectable populations of closed ribose-free conformations in solution. Additional interdomain hydrogen bonds stabilize this state. The predicted shift in equilibrium from the open to the closed state on binding to ribose is in agreement with experiments. This is driven by the energetic stabilization of the closed conformation due to ribose-protein interactions. We also observe a significant population of a hitherto unobserved ribose-bound partially open state. We believe that this state is the one that has been suggested to play a role in the transfer of ribose to the membrane-bound permease complex.     

Shape of the free energy barriers for protein folding probed by multiple perturbation analysis

The characterization of the free energy barriers has been a major goal in studies on the mechanism of protein folding. Testing the effect of mutations or denaturants on protein folding reactions revealed that transition state movement is rare, suggesting that folding barriers are robust and narrow maxima on the free energy landscape. Here we demonstrate that the application of multiple perturbations allows the observation of small transition state movements that escape detection in single perturbation experiments. We used tendamistat as a model protein to test the broadness of the free energy barriers. Tendamistat folds over two consecutive transition states and through a high-energy intermediate. Measuring the combined effect of temperature and denaturant on the position of the transition state in the wild-type protein and in several mutants revealed that the early transition state shows significant transition state movement. Its accessible surface area state becomes more native-like with destabilization of the native state by temperature. To the same extent, the entropy of the early transition state becomes more native-like with increasing denaturant concentration, in accordance with Hammond behavior. The position of the late transition state, in contrast, is much less sensitive to the applied perturbations. These results suggest that the barriers in protein folding become increasingly narrow as the folding polypeptide chain approaches the native state.     

Low free energy cost of very long loop insertions in proteins

Long insertions into a loop of a folded host protein are expected to have destabilizing effects because of the entropic cost associated with loop closure unless the inserted sequence adopts a folded structure with amino- and carboxy-termini in close proximity. A loop entropy reduction screen based on this concept was used in an attempt to retrieve folded sequences from random sequence libraries. A library of long random sequences was inserted into a loop of the SH2 domain, displayed on the surface of M13 phage, and the inserted sequences that did not disrupt SH2 function were retrieved by panning using beads coated with a phosphotyrosine containing SH2 peptide ligand. Two sequences of a library of 2 x 10(8) sequences were isolated after multiple rounds of panning, and were found to have recovery levels similar to the wild-type SH2 domain and to be relatively intolerant to further mutation in PCR mutagenesis experiments. Surprisingly, although these inserted sequences exhibited little nonrandom structure, they do not significantly destabilize the host SH2 domain. Additional insertion variants recovered at lower levels in the panning experiments were also found to have a minimal effect on the stability and peptide-binding function of the SH2 domain. The additional level of selection present in the panning experiments is likely to involve in vivo folding and assembly, as there was a rough correlation between recovery levels in the phage-panning experiments and protein solubility. The finding that loop insertions of 60-80 amino acids have minimal effects on SH2 domain stability suggests that the free energy cost of inserting long loops may be considerably less than polymer theory estimates based on the entropic cost of loop closure, and, hence, that loop insertion may have provided an evolutionary route to multidomain protein structures.     

Standardizing the free energy change of transmembrane helix-helix interactions

Side-to-side associations of transmembrane alpha-helices are integral components of the structure and function of helical membrane proteins. A fundamental unknown in the understanding of the chemical principles driving the lateral interactions between transmembrane alpha-helices is the balance of forces arising from the polypeptide sequence versus the hydrophobic solvent. To begin to address this question, a consideration of basic thermodynamic principles has been applied to assess the experimental free energy change associated with transmembrane helix dimerization in micelles. This analysis demonstrates the ability to partition the apparent free energy of transmembrane helix-helix association into two components. The first component is a statistical energy term, which arises from the fact that there are an unequal number of reactants and products. The second component is a standard state free energy change, which informs on the molecular details of the transmembrane helix self-association reaction. The advantage of separating these two energy terms arises from the fact that extrapolation to the standard state free energy change normalizes the statistical energy term so that it applies equivalently in all experimental systems. Accompanying experimental results for the glycophorin A transmembrane alpha-helix dimer measured in micelles are well described by these theoretical components assuming an ideal-dilute solution.     

NQO1-dependent redox cycling of idebenone: effects on cellular redox potential and energy levels

Short-chain quinones are described as potent antioxidants and in the case of idebenone have already been under clinical investigation for the treatment of neuromuscular disorders. Due to their analogy to coenzyme Q10 (CoQ10), a long-chain quinone, they are widely regarded as a substitute for CoQ10. However, apart from their antioxidant function, this provides no clear rationale for their use in disorders with normal CoQ10 levels. Using recombinant NAD(P)H:quinone oxidoreductase (NQO) enzymes, we observed that contrary to CoQ10 short-chain quinones such as idebenone are good substrates for both NQO1 and NQO2. Furthermore, the reduction of short-chain quinones by NQOs enabled an antimycin A-sensitive transfer of electrons from cytosolic NAD(P)H to the mitochondrial respiratory chain in both human hepatoma cells (HepG2) and freshly isolated mouse hepatocytes. Consistent with the substrate selectivity of NQOs, both idebenone and CoQ1, but not CoQ10, partially restored cellular ATP levels under conditions of impaired complex I function. The observed cytosolic-mitochondrial shuttling of idebenone and CoQ1 was also associated with reduced lactate production by cybrid cells from mitochondrial encephalomyopathy, lactic acidosis and stroke-like episodes (MELAS) patients. Thus, the observed activities separate the effectiveness of short-chain quinones from the related long-chain CoQ10 and provide the rationale for the use of short-chain quinones such as idebenone for the treatment of mitochondrial disorders.     

Nitroxides as redox probes of melanins: dark-induced and photoinduced changes in redox equilibria

The interaction of nitroxide free radicals and their reduced products (hydroxylamines) with synthetic and natural melanins has been studied. Electron spin resonance spectroscopy was used to measure changes in radical concentration in the dark and during irradiation with visible or uv light. Some reduction of nitroxide occurs in the dark, and is reversible: the nitroxide can be completely regenerated by the one-electron oxidant ferricyanide. The kinetics of the process depend strongly on radical charge and pH. For positively charged nitroxides the rate is much faster than for either neutral or anionic radicals. At pH 10 the rate is about 20 times faster than at pH 5. Oxidation of hydroxylamine also can occur so that a redox equilibrium is established. The equilibrium constant has been estimated for the reaction between a nitroxide and melanin from autoxidation of 3,4-dihydroxyphenylalanine. Results are also dependent upon the type of melanin used and chemical modification (oxidation or reduction) of the melanin. Redox equilibria are altered during irradiation with either visible or uv light. Rapid oxidation of hydroxylamine to nitroxide is apparent, together with a slower reduction of nitroxide. Action spectra for these processes are related to those for melanin radical production and oxygen consumption in nitroxide-free melanin systems. Reduction of nitroxide is inhibited by oxygen, suggesting a competition between nitroxide and oxygen for photoinduced reducing equivalents.     

Cytochrome c at charged interfaces. 1. Conformational and redox equilibria at the electrode/electrolyte interface probed by surface-enhanced resonance Raman spectroscopy

The structure and the electron-transfer properties of cytochrome c (cyt c) absorbed on a silver electrode were analyzed by surface-enhanced resonance Raman spectroscopy. It was found that the absorbed cyt c exists in various conformational states depending on the electrode potential. In state I the native structure of the heme protein is fully preserved and the redox potential (+0.02 V vs saturated calomel electrode) is close to the value for cyt c in solution. In state II the heme iron exists in a mixture of five-coordinated high-spin and six-coordinated low-spin configurations. It had been shown that these configurations form a thermal equilibrium [Hildebrandt, P., & Stockburger, M. (1986) J. Phys. Chem. 90,6017]. It is demonstrated that these equilibria strongly depend on the charge distribution within the electrical double layer of the silver electrode/electrolyte interface, indicating that the changes in the coordination shell are induced by electrostatic interactions. The structural alterations in state II are apparently restricted to the heme crevice, which assumes an open conformation compared to the close structure in state I. This leads to a strong decrease of the redox potentials, which were determined to be -0.31 and -0.41 V for the five-coordinated high-spin and six-coordinated low-spin configurations, respectively. On the other hand, gross distortions of the protein structure can be excluded since the reversible proton-induced conformational change of cyt c as found in solution at low pH also takes place in state II of the absorbed cyt c. The linkage of cyt c molecules to the surface is mediated by charged amino acid groups, and it depends on the potential which groups are thermodynamically favored to form such a molecular binding site. The conformational states I and II, which are in potential-dependent equilibrium, therefore refer to two different molecular binding sites. At potentials below zero charge (less than approximately -0.6 V) a rapid denaturation of the absorbed cyt c is noted, which is reflected by drastic and irreversible changes in the surface-enhanced resonance Raman spectrum. Our results are discussed on the background of previous electrochemical studies of cyt c at electrodes.     

Ligand escape pathways and (un)binding free energy calculations for the hexameric insulin-phenol complex

Cooperative binding of phenolic species to insulin hexamers is known to stabilize pharmaceutical preparations of the hormone. Phenol dissociation is rapid on hexamer dissolution timescales, and phenol unbinding upon dilution is likely the first step in the conversion of (pharmaceutical) hexameric insulin to the active monomeric form upon injection. However, a clear understanding of the determinants of the rates of phenol unbinding remains obscure, chiefly because residues implicated in phenol exchange as determined by NMR are not all associated with likely unbinding routes suggested by the best-resolved hexamer structures. We apply random acceleration molecular dynamics simulation to identify potential escape routes of phenol from hydrophobic cavities in the hexameric insulin-phenol complex. We find three major pathways, which provide new insights into (un)binding mechanisms for phenol. We identify several residues directly participating in escape events that serve to resolve ambiguities from recent NMR experiments. Reaction coordinates for dissociation of phenol are developed based on these exit pathways. Potentials of mean force along the reaction coordinate for each pathway are resolved using multiple independent steered molecular dynamics simulations with second-order cumulant expansion of Jarzynski's equality. Our results for ΔF agree reasonably well within the range of known experimental and previous simulation magnitudes of this quantity. Based on structural analysis and energetic barriers for each pathway, we suggest a plausible preferred mechanism of phenolic exchange that differs from previous mechanisms. Several weakly-bound metastable states are also observed for the first time in the phenol dissociation reaction.     

Conservation of the free energy change of the alkaline isomerization in mitochondrial and bacterial cytochromes c

The thermodynamic parameters of the alkaline transition for oxidized native yeast iso-1 cytochrome c and Rhodopseudomonas palustris cytochrome c(2) (cytc(2)) have been determined through direct electrochemistry experiments carried out at variable pH and temperature and compared to those for horse and beef heart cytochromes c. We have found that both transition enthalpy and entropy are remarkably species dependent, following the order R. palustris cytc(2) > beef (horse) heart cytc>yeast iso-1 cytc. Considering the high homology at the heme-protein interface in the native species, this variability is likely to be mainly determined by differences in the structural and solvation properties and the relative abundance of the various alkaline conformers. Notably, changes in transition enthalpy and entropy among these cytochromes c are compensative and result in small variations in the free energy change of the process (which amounts approximately to +50 kJ mol(-1)) and consequently in the apparent pK(a) value. This compensation indicates that solvent reorganization effects play an important role in the thermodynamics of the transition. This mechanism is functional to ensure a relatively high pK(a) value for the alkaline transition, which is needed to preserve His,Met ligation to the heme iron in cytochrome c at physiological pH and temperature, hence the E(o) value required for the biological function.     

pH-dependent effects of chlorpromazine on liposomes and erythrocyte membranes

Chlorpromazine (CPZ) is an amphipathic antipsychotic drug that binds to erythrocytes reaching in this way the central nervous system. CPZ is a basic molecule with pK=8.6. This paper reports on CPZ-induced lysis of red blood cells and liposomes. Haemolysis was tested under hypotonic conditions, in the pH range 5.0-10.0. Cell sensitivity towards CPZ increased with increasing pH. Increasing pH caused also a decrease in the critical micellar concentrations of CPZ. These results are interpreted in terms of a competition between repulsive electrostatic forces and attractive hydrophobic forces, that would act both in pure CPZ and in mixed CPZ-phospholipid micelles. In order to eliminate possible pH effects mediated by red blood cell proteins, experiments were carried out in which CPZ induced release of a fluorescent dye from liposomes (large unilamellar vesicles). The latter observations confirmed that membrane sensitivity towards CPZ was increased at higher pH.     

Detection of non-coding RNAs on the basis of predicted secondary structure formation free energy change

Background  

Non-coding RNAs (ncRNAs) have a multitude of roles in the cell, many of which remain to be discovered. However, it is difficult to detect novel ncRNAs in biochemical screens. To advance biological knowledge, computational methods that can accurately detect ncRNAs in sequenced genomes are therefore desirable. The increasing number of genomic sequences provides a rich dataset for computational comparative sequence analysis and detection of novel ncRNAs.     

A bulk water-dependent desolvation energy model for analyzing the effects of secondary solutes on biological equilibria

A new phenomenological model for interpreting the effects of solutes on biological equilibria is presented. The model attributes changes in equilibria to differences in the desolvation energy of the reacting species that, in turn, reflect changes in the free energy of the bulk water upon addition of secondary solutes. The desolvation approach differs notably from that of other solute models by treating the free energy of bulk water as a variable and by not ascribing the observed shifts in reaction equilibria to accumulation or depletion of solutes next to the surfaces of the reacting species. On the contrary, the partitioning of solutes is viewed as a manifestation of the different subpopulations of water that arise in response to the surface boundary conditions. A thermodynamic framework consistent with the proposed model is used to derive a relationship for a specific reaction, an aqueous solubility equilibrium, in two or more solutions. The resulting equation reconciles some potential issues with the transfer free energy model of Tanford. Application of the desolvation energy model to the analysis of a two-state protein folding equilibrium is discussed and contrasted to the application of two other solute models developed by Timasheff and by Parsegian. Future tabulation of solvation energies and bulk water energies may allow biophysical chemists to confirm the mechanism by which secondary solutes influence binding and conformational equilibria and may provide a common ground on which experimentalists and theoreticians can compare and evaluate their results.     

Prediction of the mutation-induced change in thermodynamic stabilities of membrane proteins from free energy simulations

Comparative protein structure modeling and free energy perturbation simulation have been applied in a consecutive manner to investigate the mutation-induced stabilization of membrane proteins (MPs) in aqueous solution without knowledge of their three-dimensional structures. The calculated difference in protein solvation free energy between the wild type and a mutant compares well with their relative thermodynamic stabilities in solution. For monomeric MPs, a mutant reveals a higher stability than the wild type if the calculated solvation free energy indicates a favorable change. On the contrary, for oligomeric MPs the stability of a mutant increases as the solvation free energy of a mutated monomer becomes less favorable, indicating that the oligomeric MP mutant would be stabilized in solution due to the reduced desolvation cost for oligomerization. The present computational strategy is expected to find its way as a useful tool for assessing the relative stability of a mutant MP with respect to its wild type in solution.