Coupling of Protonation, Reduction, and Conformational Change in azurin from Pseudomonas aeruginosa Investigated with Free Energy Measures of Cooperativity

The Influence of pH on Long-Range Electron Transfer and Proton-Coupled Electron Transfer in Ruthenium-Modified Azurin

Long-range electron transfer (ET) is an essential component of all biological systems. Reactions of metalloproteins are important in this context. Recent work on protein "charge ladders" has revealed how the redox state of embedded metal ions can influence the ionization of amino acid residues at protein surface sites. Inspired by these observations, we carried out a variable pH investigation of intramolecular ET reactions in a ruthenium-modified protein system built on azurin from Pseudomonas aeruginosa. We also generate a Pourbaix diagram that describes the variable pH redox behavior of a Ru model complex, Ru(2,2'-bipyridyl)2(imidazole)2(PF6)2. The intramolecular ET rate constants for the oxidation of azurin-Cu+ by flash-quench-generated Ru3+ oxidants do not follow the predictions of the semi-classical ET rate expression with fixed values of reorganization energy (λ) and electronic coupling (HDA). Based on the pH dependence of the Ru3+/2+ redox couple, we propose a model where pure ET is operative at acidic pH values (≤ 7) and the mechanism changes to proton-coupled electron transfer at pH ≥ 7.5. The implications of this mechanistic proposal are discussed in the context of biological redox reactions and with respect to other examples of intramolecular ET reactions in the literature.

Shift from Entropic Cu2+ Binding to Enthalpic Cu+ Binding Determines the Reduction Thermodynamics of Blue Copper Proteins

The enthalpic and entropic components of Cu²⁺ and Cu⁺ binding to the blue copper protein azurin have been quantified with isothermal titration calorimetry (ITC) measurements and analysis, providing the first such experimental values for Cu⁺ binding to a protein. The high affinity of azurin for Cu²⁺ is entirely due to a very favorable binding entropy, while its even higher affinity for Cu⁺ is due to a favorable binding enthalpy and entropy. The binding thermodynamics provide insight into bond enthalpies at the blue copper site and entropic contributions from desolvation and proton displacement. These values were used in thermodynamic cycles to determine the enthalpic and entropic contributions to the free energy of reduction and thus the reduction potential. The reduction thermodynamics obtained with this method are in good agreement with previous results from temperature-dependent electrochemical measurements. The calorimetry method, however, provides new insight into contributions from the initial (oxidized) and final (reduced) states of the reduction. Since ITC measurements quantify the protons that are displaced upon metal binding, the proton transfer that is coupled with electron transfer is also determined with this method. Preliminary results for Cu²⁺ and Cu⁺ binding to the Phe114Pro variant of azurin demonstrate the insight about protein tuning of the reduction potential that is provided by the binding thermodynamics of each metal oxidation state.

DEPC modification of the CuA protein from Thermus thermophilus

The Cu_A center is the initial electron acceptor in cytochrome c oxidase, and it consists of two copper ions bridged by two cysteines and ligated by two histidines, a methionine, and a carbonyl in the peptide backbone of a nearby glutamine. The two ligating histidines are of particular interest as they may influence the electronic and redox properties of the metal center. To test for the presence of reactive ligating histidines, a portion of cytochrome c oxidase from the bacteria Thermus thermophilus that contains the Cu_A site (the TtCu_A protein) was treated with the chemical modifier diethyl pyrocarbonate (DEPC) and the reaction followed through UV-visible, circular dichroism, and electron paramagnetic resonance spectroscopies at pH 5.0-9.0. A mutant protein (H40A/H117A) with the non-ligating histidines removed was similarly tested. Introduction of an electron-withdrawing DEPC-modification onto the ligating histidine 157 of TtCu_A increased the reduction potential by over 70 mV, as assessed by cyclic voltammetry. Results from both proteins indicate that DEPC reacts with one of the two ligating histidines, modification of a ligating histidine raises the reduction potential of the Cu_A site, and formation of the DEPC adduct is reversible at room temperature. The existence of the reactive ligating histidine suggests that this residue may play a role in modulating the electronic and redox properties of TtCu_A through kinetically-controlled proton exchange with the solvent. Lack of reactivity by the metalloproteins Sco and azurin, both of which contain a mononuclear copper center, indicate that reactivity toward DEPC is not a characteristic of all ligating histidines.

Charge-Neutral Constant pH Molecular Dynamics Simulations Using a Parsimonious Proton Buffer

In constant pH molecular dynamics simulations, the protonation states of titratable sites can respond to changes of the pH and of their electrostatic environment. Consequently, the number of protons bound to the biomolecule, and therefore the overall charge of the system, fluctuates during the simulation. To avoid artifacts associated with a non-neutral simulation system, we introduce an approach to maintain neutrality of the simulation box in constant pH molecular dynamics simulations, while maintaining an accurate description of all protonation fluctuations. Specifically, we introduce a proton buffer that, like a buffer in experiment, can exchange protons with the biomolecule enabling its charge to fluctuate. To keep the total charge of the system constant, the uptake and release of protons by the buffer are coupled to the titration of the biomolecule with a constraint. We find that, because the fluctuation of the total charge (number of protons) of a typical biomolecule is much smaller than the number of titratable sites of the biomolecule, the number of buffer sites required to maintain overall charge neutrality without compromising the charge fluctuations of the biomolecule, is typically much smaller than the number of titratable sites, implying markedly enhanced simulation and sampling efficiency.

Cooperative binding: a multiple personality

Cooperative binding has been described in many publications and has been related to or defined by several different properties of the binding behavior of the ligand to the target molecule. In addition to the commonly used Hill coefficient, other characteristics such as a sigmoidal shape of the overall titration curve in a linear plot, a change of ligand affinity of the other binding sites when a site of the target molecule becomes occupied, or complex roots of the binding polynomial have been used to define or to quantify cooperative binding. In this work, we analyze how the different properties are related in the most general model for binding curves based on the grand canonical partition function and present several examples which highlight differences between the cooperativity characterizing properties which are discussed. Our results mainly show that among the presented definitions there are not two which fully coincide. Moreover, this work poses the question whether it can make sense to distinguish between positive and negative cooperativity based on the macroscopic binding isotherm only. This article shall emphasize that scientists who investigate cooperative effects in biological systems could help avoiding misunderstandings by stating clearly which kind of cooperativity they discuss.

Converting the yeast arginine can1 permease to a lysine permease

Amino acid uptake in yeast cells is mediated by about 16 plasma membrane permeases, most of which belong to the amino acid-polyamine-organocation (APC) transporter family. These proteins display various substrate specificity ranges. For instance, the general amino acid permease Gap1 transports all amino acids, whereas Can1 and Lyp1 catalyze specific uptake of arginine and lysine, respectively. Although Can1 and Lyp1 have different narrow substrate specificities, they are close homologs. Here we investigated the molecular rules determining the substrate specificity of the H(+)-driven arginine-specific permease Can1. Using a Can1-Lyp1 sequence alignment as a guideline and a three-dimensional Can1 structural model based on the crystal structure of the bacterial APC family arginine/agmatine antiporter, we introduced amino acid substitutions liable to alter Can1 substrate specificity. We show that the single substitution T456S results in a Can1 variant transporting lysine in addition to arginine and that the combined substitutions T456S and S176N convert Can1 to a Lyp1-like permease. Replacement of a highly conserved glutamate in the Can1 binding site leads to variants (E184Q and E184A) incapable of any amino acid transport, pointing to a potential role for this glutamate in H(+) coupling. Measurements of the kinetic parameters of arginine and lysine uptake by the wild-type and mutant Can1 permeases, together with docking calculations for each amino acid in their binding site, suggest a model in which residues at positions 176 and 456 confer substrate selectivity at the ligand-binding stage and/or in the course of conformational changes required for transport.

Thermodynamics of transport through the ammonium transporter Amt-1 investigated with free energy calculations

Amt-1 from Archaeoglobus fulgidus (AfAmt-1) belongs to the Amt/Rh family of ammonium/ammonia transporting membrane proteins. The transport mode and the precise microscopic permeation mechanism utilized by these proteins are intensely debated. Open questions concern the identity of the transported substrate (ammonia and/or ammonium) and whether the transport is passive or active. To address these questions, we studied the overall thermodynamics of the different transport modes as a function of the environmental conditions. Then, we investigated the thermodynamics of the underlying microscopic transport mechanisms with free energy calculations within a continuum electrostatics model. The formalism developed for this purpose is of general utility in the calculation of binding free energies for ligands with multiple protonation forms or other binding forms. The results of our calculations are compared to the available experimental and theoretical data on Amt/Rh proteins and discussed in light of the current knowledge on the physiological conditions experienced by microorganisms and plants. We found that microscopic models of electroneutral and electrogenic transport modes are in principle thermodynamically viable. However, only the electrogenic variants have a net thermodynamic driving force under the physiological conditions experienced by microorganisms and plants. Thus, the transport mechanism of AfAmt-1 is most likely electrogenic.

A mathematical view on the decoupled sites representation

The decoupled sites representation (DSR) is a theoretical instrument which allows to regard complex pH titration curves of biomolecules with several interacting proton binding sites as composition of isolated, non-interacting sites, each with a standard Henderson-Hasselbalch titration curve. In this work, we present the mathematical framework in which the DSR is embedded and give mathematical proofs for several statements in the periphery of the DSR. These proofs also identify exceptions. To apply the DSR to any molecule, it is necessary to extend the set of binding energies from R to a stripe within C. An important observation in this context is that even positive interaction energies (repulsion) between the binding sites will not guarantee real binding energies in the decoupled system, at least if the molecule has more than four proton binding sites. Moreover, we show that for a given overall titration curve it is not only possible to find a corresponding system with an interaction energy of zero but with any arbitrary fix interaction energy. This result also effects practical work as it shows that for any given titration curve, there is an infinite number of corresponding hypothetical molecules. Furthermore, this implies that--using a common definition of cooperative binding on the level of interaction energies--a meaningful measure of cooperativity between the binding sites cannot be defined solely on the basis of the overall titration. Consequently, all measures of cooperativity based on the overall binding curve do not measure the type of cooperativity commonly defined on the basis of interaction energies. Understanding the DSR mathematically provides the basis of transferring the DSR to biomolecules with different types of interacting ligands, such as protons and electrons, which play an important role within electron transport chains like in photosynthesis.

GMCT : a Monte Carlo simulation package for macromolecular receptors

Generalized Monte Carlo titration (GMCT) is a versatile suite of computer programs for the efficient simulation of complex macromolecular receptor systems as for example proteins. The computational model of the system is based on a microstate description of the receptor and an average description of its surroundings in terms of chemical potentials. The receptor can be modeled in great detail including conformational flexibility and many binding sites with multiple different forms that can bind different ligand types. Membrane embedded systems can be modeled including electrochemical potential gradients. Overall properties of the receptor as well as properties of individual sites can be studied with a variety of different Monte Carlo (MC) simulation methods. Metropolis MC, Wang-Landau MC and efficient free energy calculation methods are included. GMCT is distributed as free open source software at www.bisb.uni-bayreuth.de under the terms of the GNU Affero General Public License.