Electrostatic environment of hemes in proteins: pK(a)s of hydroxyl ligands

Azole-Based Indoleamine 2,3-Dioxygenase 1 (IDO1) Inhibitors

The heme enzyme indoleamine 2,3-dioxygenase 1 (IDO1) plays an essential role in immunity, neuronal function, and aging through catalysis of the rate-limiting step in the kynurenine pathway of tryptophan metabolism. Many IDO1 inhibitors with different chemotypes have been developed, mainly targeted for use in anti-cancer immunotherapy. Lead optimization of direct heme iron-binding inhibitors has proven difficult due to the remarkable selectivity and sensitivity of the heme-ligand interactions. Here, we present experimental data for a set of closely related small azole compounds with more than 4 orders of magnitude differences in their inhibitory activities, ranging from millimolar to nanomolar levels. We investigate and rationalize their activities based on structural data, molecular dynamics simulations, and density functional theory calculations. Our results not only expand the presently known four confirmed chemotypes of sub-micromolar heme binding IDO1 inhibitors by two additional scaffolds but also provide a model to predict the activities of novel scaffolds.

Identifying the proton loading site cluster in the ba3 cytochrome c oxidase that loads and traps protons

Cytochrome c Oxidase (CcO) is the terminal electron acceptor in aerobic respiratory chain, reducing O₂ to water. The released free energy is stored by pumping protons through the protein, maintaining the transmembrane electrochemical gradient. Protons are held transiently in a proton loading site (PLS) that binds and releases protons driven by the electron transfer reaction cycle. Multi-Conformation Continuum Electrostatics (MCCE) was applied to crystal structures and Molecular Dynamics snapshots of the B-type Thermus thermophilus CcO. Six residues are identified as the PLS, binding and releasing protons as the charges on heme b and the binuclear center are changed: the heme a₃ propionic acids, Asp287, Asp372, His376 and Glu126B. The unloaded state has one proton and the loaded state two protons on these six residues. Different input structures, modifying the PLS conformation, show different proton distributions and result in different proton pumping behaviors. One loaded and one unloaded protonation states have the loaded/unloaded states close in energy so the PLS binds and releases a proton through the reaction cycle. The alternative proton distributions have state energies too far apart to be shifted by the electron transfers so are locked in loaded or unloaded states. Here the protein can use active states to load and unload protons, but has nearby trapped states, which stabilize PLS protonation state, providing new ideas about the CcO proton pumping mechanism. The distance between the PLS residues Asp287 and His376 correlates with the energy difference between loaded and unloaded states.

Ruffling drives coproheme decarboxylation by facilitating PCET: a theoretical investigation of ChdC

Coproheme decarboxylase (ChdC) is an essential enzyme in the coproporphyrin-dependent heme synthesis pathway, which catalyzes oxidative decarboxylation of coproheme at the positions p2 and p4 to generate heme b under the action of hydrogen peroxide.

Tracing the Pathways of Waters and Protons in Photosystem II and Cytochrome c Oxidase

Photosystem II (PSII) uses water as the terminal electron donor, producing oxygen in the Mn4CaO5 oxygen evolving complex (OEC), while cytochrome c oxidase (CcO) reduces O2 to water in its heme–Cu binuclear center (BNC). Each protein is oriented in the membrane to add to the proton gradient. The OEC, which releases protons, is located near the P-side (positive, at low-pH) of the membrane. In contrast, the BNC is in the middle of CcO, so the protons needed for O2 reduction must be transferred from the N-side (negative, at high pH). In addition, CcO pumps protons from N- to P-side, coupled to the O2 reduction chemistry, to store additional energy. Thus, proton transfers are directly coupled to the OEC and BNC redox chemistry, as well as needed for CcO proton pumping. The simulations that study the changes in proton affinity of the redox active sites and the surrounding protein at different states of the reaction cycle, as well as the changes in hydration that modulate proton transfer paths, are described.

AMOEBA Polarizable Force Field Parameters of the Heme Cofactor in Its Ferrous and Ferric Forms

We report the first parameters of the heme redox cofactors for the polarizable AMOEBA force field in both the ferric and ferrous forms. We consider two types of complexes, one with two histidine side chains as axial ligands and one with a histidine and a methionine side chain as ligands. We have derived permanent multipoles from second-order Møller-Plesset perturbation theory (MP2). The sets of parameters have been validated in a first step by comparison of AMOEBA interaction energies of heme and a collection of biologically relevant molecules with MP2 and Density Functional Theory (DFT) calculations. In a second validation step, we consider interaction energies with large aggregates comprising around 80 H₂O molecules. These calculations are repeated for 30 structures extracted from semiempirical PM7 DM simulations. Very encouraging agreement is found between DFT and the AMOEBA force field, which results from an accurate treatment of electrostatic interactions. We finally report long (10 ns) MD simulations of cytochromes in two redox states with AMOEBA testing both the 2003 and 2014 AMOEBA water models. These simulations have been carried out with the TINKER-HP (High Performance) program. In conclusion, owing to their ubiquity in biology, we think the present work opens a wide array of applications of the polarizable AMOEBA force field on hemeproteins.

Analyzing the electrogenicity of cytochrome c oxidase

Measurements of voltage changes in response to charge separation within membrane proteins can offer fundamental information on spectroscopically "invisible" steps. For example, results from studies of voltage changes associated with electron and proton transfer in cytochrome c oxidase could, in principle, be used to discriminate between different theoretical models describing the molecular mechanism of proton pumping. Earlier analyses of data from these measurements have been based on macroscopic considerations that may not allow for exploring the actual molecular mechanisms. Here, we have used a coarse-grained model describing the relation between observed voltage changes and specific charge-transfer reactions, which includes an explicit description of the membrane, the electrolytes, and the electrodes. The results from these calculations offer mechanistic insights at the molecular level. Our main conclusion is that previously assumed mechanistic evidence that was based on electrogenic measurements is not unique. However, the ability of our calculations to obtain reliable voltage changes means that we have a tool that can be used to describe a wide range of electrogenic charge transfers in channels and transporters, by combining voltage measurements with other experiments and simulations to analyze new mechanistic proposals.

Designing Light-Activated Charge-Separating Proteins with a Naphthoquinone Amino Acid

The first principles design of manmade redox-protein maquettes is used to clarify the physical/chemical engineering supporting the mechanisms of natural enzymes with a view to recapitulate and surpass natural performance. Herein, we use intein-based protein semisynthesis to pair a synthetic naphthoquinone amino acid (Naq) with histidine-ligated photoactive metal-tetrapyrrole cofactors, creating a 100 μs photochemical charge separation unit akin to photosynthetic reaction centers. By using propargyl groups to protect the redox-active para-quinone during synthesis and assembly while permitting selective activation, we gain the ability to employ the quinone amino acid redox cofactor with the full set of natural amino acids in protein design. Direct anchoring of quinone to the protein backbone permits secure and adaptable control of intraprotein electron-tunneling distances and rates.

Ferryl protonation in oxoiron(IV) porphyrins and its role in oxygen transfer

Ferryl porphyrins, P-Fe(IV)═O, are central reactive intermediates in the catalytic cycles of numerous heme proteins and a variety of model systems. There has been considerable interest in elucidating factors, such as terminal oxo basicity, that may control ferryl reactivity. Here, the sulfonated, water-soluble ferryl porphyrin complexes tetramesitylporphyrin, oxoFe(IV)TMPS (FeTMPS-II), its 2,6-dichlorophenyl analogue, oxoFe(IV)TDClPS (FeTDClPS-II), and two other analogues are shown to be protonated under turnover conditions to produce the corresponding bis-aqua-iron(III) porphyrin cation radicals. The results reveal a novel internal electromeric equilibrium, P-Fe(IV)═O ⇆ P(+)-Fe(III)(OH2)2. Reversible pKa values in the range of 4-6.3 have been measured for this process by pH-jump, UV-vis spectroscopy. Ferryl protonation has important ramifications for C-H bond cleavage reactions mediated by oxoiron(IV) porphyrin cation radicals in protic media. Both solvent O-H and substrate C-H deuterium kinetic isotope effects are observed for these reactions, indicating that hydrocarbon oxidation by these oxoiron(IV) porphyrin cation radicals occurs via a solvent proton-coupled hydrogen atom transfer from the substrate that has not been previously described. The effective FeO-H bond dissociation energies for FeTMPS-II and FeTDClPS-II were estimated from similar kinetic reactivities of the corresponding oxoFe(IV)TMPS(+) and oxoFe(IV)TDClPS(+) species to be ∼92-94 kcal/mol. Similar values were calculated from the two-proton P(+)-Fe(III)(OH2)2 pKa(obs) and the porphyrin oxidation potentials, despite a 230 mV range for the iron porphyrins examined. Thus, the iron porphyrin with the lower ring oxidation potential has a compensating higher basicity of the ferryl oxygen. The solvent-derived proton adds significantly to the driving force for C-H bond scission.

Molecular mechanisms for generating transmembrane proton gradients

Membrane proteins use the energy of light or high energy substrates to build a transmembrane proton gradient through a series of reactions leading to proton release into the lower pH compartment (P-side) and proton uptake from the higher pH compartment (N-side). This review considers how the proton affinity of the substrates, cofactors and amino acids are modified in four proteins to drive proton transfers. Bacterial reaction centers (RCs) and photosystem II (PSII) carry out redox chemistry with the species to be oxidized on the P-side while reduction occurs on the N-side of the membrane. Terminal redox cofactors are used which have pKas that are strongly dependent on their redox state, so that protons are lost on oxidation and gained on reduction. Bacteriorhodopsin is a true proton pump. Light activation triggers trans to cis isomerization of a bound retinal. Strong electrostatic interactions within clusters of amino acids are modified by the conformational changes initiated by retinal motion leading to changes in proton affinity, driving transmembrane proton transfer. Cytochrome c oxidase (CcO) catalyzes the reduction of O2 to water. The protons needed for chemistry are bound from the N-side. The reduction chemistry also drives proton pumping from N- to P-side. Overall, in CcO the uptake of 4 electrons to reduce O2 transports 8 charges across the membrane, with each reduction fully coupled to removal of two protons from the N-side, the delivery of one for chemistry and transport of the other to the P-side.

Coupled electron and proton transfer reactions during the O→E transition in bovine cytochrome c oxidase

A combined DFT/electrostatic approach is employed to study the coupling of proton and electron transfer reactions in cytochrome c oxidase (CcO) and its proton pumping mechanism. The coupling of the chemical proton to the internal electron transfer within the binuclear center is examined for the O→E transition. The novel features of the His291 pumping model are proposed, which involve timely well-synchronized sequence of the proton-coupled electron transfer reactions. The obtained pK(a)s and E(m)s of the key ionizable and redox-active groups at the different stages of the O→E transition are consistent with available experimental data. The PT step from E242 to H291 is examined in detail for various redox states of the hemes and various conformations of E242 side-chain. Redox potential calculations of the successive steps in the reaction cycle during the O→E transition are able to explain a cascade of equilibria between the different intermediate states and electron redistribution between the metal centers during the course of the catalytic activity. All four electrometric phases are discussed in the light of the obtained results, providing a robust support for the His291 model of proton pumping in CcO.

Slow dissociation of a charged ligand: analysis of the primary quinone Q(A) site of photosynthetic bacterial reaction centers

Reaction centers (RCs) are integral membrane proteins that undergo a series of electron transfer reactions during the process of photosynthesis. In the Q(A) site of RCs from Rhodobacter sphaeroides, ubiquinone-10 is reduced, by a single electron transfer, to its semiquinone. The neutral quinone and anionic semiquinone have similar affinities, which is required for correct in situ reaction thermodynamics. A previous study showed that despite similar affinities, anionic quinones associate and dissociate from the Q(A) site at rates ≈10(4) times slower than neutral quinones indicating that anionic quinones encounter larger binding barriers (Madeo, J.; Gunner, M. R. Modeling binding kinetics at the Q(A) site in bacterial reaction centers. Biochemistry 2005, 44, 10994-11004). The present study investigates these barriers computationally, using steered molecular dynamics (SMD) to model the unbinding of neutral ground state ubiquinone (UQ) and its reduced anionic semiquinone (SQ(-)) from the Q(A) site. In agreement with experiment, the SMD unbinding barrier for SQ(-) is larger than for UQ. Multi Conformational Continuum Electrostatics (MCCE), used here to calculate the binding energy, shows that SQ(-) and UQ have comparable affinities. In the Q(A) site, there are stronger binding interactions for SQ(-) compared to UQ, especially electrostatic attraction to a bound non-heme Fe(2+). These interactions compensate for the higher SQ(-) desolvation penalty, allowing both redox states to have similar affinities. These additional interactions also increase the dissociation barrier for SQ(-) relative to UQ. Thus, the slower SQ(-) dissociation rate is a direct physical consequence of the additional binding interactions required to achieve a Q(A) site affinity similar to that of UQ. By a similar mechanism, the slower association rate is caused by stronger interactions between SQ(-) and the polar solvent. Thus, stronger interactions for both the unbound and bound states of charged and highly polar ligands can slow their binding kinetics without a conformational gate. Implications of this for other systems are discussed.

Exploration of the cytochrome c oxidase pathway puzzle and examination of the origin of elusive mutational effects

Gaining detailed understanding of the energetics of the proton-pumping process in cytochrome c oxidase (CcO) is a problem of great current interest. Despite promising mechanistic proposals, so far, a physically consistent model that would reproduce all the relevant barriers needed to create a working pump has not been presented. In addition, there are major problems in elucidating the origin of key mutational effects and in understanding the nature of the apparent pK(a) values associated with the pH dependencies of specific proton transfer (PT) reactions in CcO. This work takes a key step in resolving the above problems, by considering mutations, such as the Asn139Asp replacement, that blocks proton pumping without affecting PT to the catalytic site. We first introduce a formulation that makes it possible to relate the apparent pK(a) of Glu286 to different conformational states of this residue. We then use the new formulation along with the calculated pK(a) values of Glu286 at these different conformations to reproduce the experimentally observed apparent pK(a) of the residue. Next, we take the X-ray structures of the native and Asn139Asp mutant of the Paracoccus denitrificans CcO (N131D in this system) and reproduce for the first time the change in the primary PT pathways (and other key features) based on simulations that start with the observed structural changes. We also consider the competition between proton transport to the catalytic site and the pump site, as a function of the bulk pH, as well as the H/D isotope effect, and use this information to explore the relative height of the two barriers. The paper emphasizes the crucial role of energy-based considerations that include the PT process, and the delicate control of PT in CcO.

Multiconformation continuum electrostatics analysis of the effects of a buried Asp introduced near heme a in Rhodobacter sphaeroides cytochrome c oxidase

Cytochrome c oxidase (CcO) reduces O(2) to water via a series of proton-coupled electron transfers, generating a transmembrane electrochemical gradient. Coupling electron and proton transfer requires changing the pK(a) values of buried residues at each stage in the reaction cycle. Heme a is a key cofactor in the CcO electron transfer chain. Mutation of Ser44 to Asp has been reported [Mills, D. A., et al. (2008) Biochemistry 47, 11499-11509], changing the hydrogen bond acceptor from His102, the heme a axial ligand in Rhodobactor sphaeroides CcO. This adds an acidic residue to the CcO interior. The electrochemical behavior of heme a in wild-type and S44D CcO is compared using the continuum electrostatics program MCCE. The introduced, deeply buried Asp remains ionized at physiological pH only when the nearby heme is oxidized. Heme a reduction is now calculated to be strongly coupled to Asp proton binding, while with Ser44, it is weakly coupled to small protonation shifts at multiple sites, increasing the pH dependence in the mutant. At pH 7, the partially ionized Asp 44 is calculated to lower the heme redox potential by 50 mV as expected given the thermodynamics of coupled electron and proton transfers. This highlights an curious finding in the experimental results where a low Asp pK(a) is found together with a stabilized reduced heme. The stabilization of a heme oxidation in a model complex by a hydrogen bond to the axial His ligand calculated with continuum electrostatics and with density functional theory were in good agreement.

MCCE2: improving protein pKa calculations with extensive side chain rotamer sampling

Multiconformation continuum electrostatics (MCCE) explores different conformational degrees of freedom in Monte Carlo calculations of protein residue and ligand pK(a)s. Explicit changes in side chain conformations throughout a titration create a position dependent, heterogeneous dielectric response giving a more accurate picture of coupled ionization and position changes. The MCCE2 methods for choosing a group of input heavy atom and proton positions are described. The pK(a)s calculated with different isosteric conformers, heavy atom rotamers and proton positions, with different degrees of optimization are tested against a curated group of 305 experimental pK(a)s in 33 proteins. QUICK calculations, with rotation around Asn and Gln termini, sampling His tautomers and torsion minimum hydroxyls yield an RMSD of 1.34 with 84% of the errors being <1.5 pH units. FULL calculations adding heavy atom rotamers and side chain optimization yield an RMSD of 0.90 with 90% of the errors <1.5 pH unit. Good results are also found for pK(a)s in the membrane protein bacteriorhodopsin. The inclusion of extra side chain positions distorts the dielectric boundary and also biases the calculated pK(a)s by creating more neutral than ionized conformers. Methods for correcting these errors are introduced. Calculations are compared with multiple X-ray and NMR derived structures in 36 soluble proteins. Calculations with X-ray structures give significantly better pK(a)s. Results with the default protein dielectric constant of 4 are as good as those using a value of 8. The MCCE2 program can be downloaded from http://www.sci.ccny.cuny.edu/~mcce.

Analysis of the electrochemistry of hemes with E(m)s spanning 800 mV

The free energy of heme reduction in different proteins is found to vary over more than 18 kcal/mol. It is a challenge to determine how proteins manage to achieve this enormous range of E(m)s with a single type of redox cofactor. Proteins containing 141 unique hemes of a-, b-, and c-type, with bis-His, His-Met, and aquo-His ligation were calculated using Multi-Conformation Continuum Electrostatics (MCCE). The experimental E(m)s range over 800 mV from -350 mV in cytochrome c(3) to 450 mV in cytochrome c peroxidase (vs. SHE). The quantitative analysis of the factors that modulate heme electrochemistry includes the interactions of the heme with its ligands, the solvent, the protein backbone, and sidechains. MCCE calculated E(m)s are in good agreement with measured values. Using no free parameters the slope of the line comparing calculated and experimental E(m)s is 0.73 (R(2) = 0.90), showing the method accounts for 73% of the observed E(m) range. Adding a +160 mV correction to the His-Met c-type hemes yields a slope of 0.97 (R(2) = 0.93). With the correction 65% of the hemes have an absolute error smaller than 60 mV and 92% are within 120 mV. The overview of heme proteins with known structures and E(m)s shows both the lowest and highest potential hemes are c-type, whereas the b-type hemes are found in the middle E(m) range. In solution, bis-His ligation lowers the E(m) by approximately 205 mV relative to hemes with His-Met ligands. The bis-His, aquo-His, and His-Met ligated b-type hemes all cluster about E(m)s which are approximately 200 mV more positive in protein than in water. In contrast, the low potential bis-His c-type hemes are shifted little from in solution, whereas the high potential His-Met c-type hemes are raised by approximately 300 mV from solution. The analysis shows that no single type of interaction can be identified as the most important in setting heme electrochemistry in proteins. For example, the loss of solvation (reaction field) energy, which raises the E(m), has been suggested to be a major factor in tuning in situ E(m)s. However, the calculated solvation energy vs. experimental E(m) shows a slope of 0.2 and R(2) of 0.5 thus correlates weakly with E(m)s. All other individual interactions show even less correlation with E(m). However the sum of these terms does reproduce the range of observed E(m)s. Therefore, different proteins use different aspects of their structures to modulate the in situ heme electrochemistry. This study also shows that the calculated E(m)s are relatively insensitive to different heme partial charges and to the protein dielectric constant used in the simulation.

Stabilization and characterization of a heme-oxy reaction intermediate in inducible nitric-oxide synthase

Nitric-oxide synthases (NOS) are heme-thiolate enzymes that N-hydroxylate L-arginine (L-Arg) to make NO. NOS contain a unique Trp residue whose side chain stacks with the heme and hydrogen bonds with the heme thiolate. To understand its importance we substituted His for Trp188 in the inducible NOS oxygenase domain (iNOSoxy) and characterized enzyme spectral, thermodynamic, structural, kinetic, and catalytic properties. The W188H mutation had relatively small effects on l-Arg binding and on enzyme heme-CO and heme-NO absorbance spectra, but increased the heme midpoint potential by 88 mV relative to wild-type iNOSoxy, indicating it decreased heme-thiolate electronegativity. The protein crystal structure showed that the His188 imidazole still stacked with the heme and was positioned to hydrogen bond with the heme thiolate. Analysis of a single turnover L-Arg hydroxylation reaction revealed that a new heme species formed during the reaction. Its build up coincided kinetically with the disappearance of the enzyme heme-dioxy species and with the formation of a tetrahydrobiopterin (H4B) radical in the enzyme, whereas its subsequent disappearance coincided with the rate of l-Arg hydroxylation and formation of ferric enzyme. We conclude: (i) W188H iNOSoxy stabilizes a heme-oxy species that forms upon reduction of the heme-dioxy species by H4B. (ii) The W188H mutation hinders either the processing or reactivity of the heme-oxy species and makes these steps become rate-limiting for l-Arg hydroxylation. Thus, the conserved Trp residue in NOS may facilitate formation and/or reactivity of the ultimate hydroxylating species by tuning heme-thiolate electronegativity.

Electrostatic basis for the unidirectionality of the primary proton transfer in cytochrome c oxidase

Gaining detailed understanding of the energetics of the proton-pumping process in cytochrome c oxidase (CcO) is one of the challenges of modern biophysics. Despite promising mechanistic proposals, most works have not related the activation barriers of the different assumed steps to the protein structure, and there has not been a physically consistent model that reproduced the barriers needed to create a working pump. This work reevaluates the activation barriers for the primary proton transfer (PT) steps by calculations that reflect all relevant free energy contributions, including the electrostatic energies of the generated charges, the energies of water insertion, and large structural rearrangements of the donor and acceptor. The calculations have reproduced barriers that account for the directionality and sequence of events in the primary PT in CcO. It has also been found that the PT from Glu-286 (E) to the propionate of heme a(3) (Prd) provides a gate for an initial back leakage from the high pH side of the membrane. Interestingly, the rotation of E that brings it closer to Prd appears to provide a way for blocking competing pathways in the primary PT. Our study elucidates and quantifies the nature of the control of the directionality in the primary PT in CcO and provides instructive insight into the role of the water molecules in biological PT, showing that "bridges" of several water molecules in hydrophobic regions present a problem (rather than a solution) that is minimized in the primary PT.

Side-chain protonation and mobility in the sarcoplasmic reticulum Ca2+-ATPase: implications for proton countertransport and Ca2+ release

Protonation of acidic residues in the sarcoplasmic reticulum Ca(2+)-ATPase (SERCA 1a) was studied by multiconformation continuum electrostatic calculations in the Ca(2+)-bound state Ca(2)E1, in the Ca(2+)-free state E2(TG) with bound thapsigargin, and in the E2P (ADP-insensitive phosphoenzyme) analog state with MgF(4)(2-) E2(TG+MgF(4)(2-)). Around physiological pH, all acidic Ca(2+) ligands (Glu(309), Glu(771), Asp(800), and Glu(908)) were unprotonated in Ca(2)E1; in E2(TG) and E2(TG+MgF(4)(2-)) Glu(771), Asp(800), and Glu(908) were protonated. Glu(771) and Glu(908) had calculated pK(a) values larger than 14 in E2(TG) and E2(TG+MgF(4)(2-)), whereas Asp(800) titrated with calculated pK(a) values near 7.5. Glu(309) had very different pK(a) values in the Ca(2+)-free states: 8.4 in E2(TG+MgF(4)(2-)) and 4.7 in E2(TG) because of a different local backbone conformation. This indicates that Glu(309) can switch between a high and a low pK(a) mode, depending on the local backbone conformation. Protonated Glu(309) occupied predominantly two main, very differently orientated side-chain conformations in E2(TG+MgF(4)(2-)): one oriented inward toward the other Ca(2+) ligands and one oriented outward toward a protein channel that seems to be in contact with the cytoplasm. Upon deprotonation, Glu(309) adopted completely the outwardly orientated side-chain conformation. The contact of Glu(309) with the cytoplasm in E2(TG+MgF(4)(2-)) makes this residue unlikely to bind lumenal protons. Instead it might serve as a proton shuttle between Ca(2+)-binding site I and the cytoplasm. Glu(771), Asp(800), and Glu(908) are proposed to take part in proton countertransport.

Protonation and hydrogen bonding of Ca2+ site residues in the E2P phosphoenzyme intermediate of sarcoplasmic reticulum Ca2+-ATPase studied by a combination of infrared spectroscopy and electrostatic calculations

Protonation of the Ca(2+) ligands of the SR Ca(2+)-ATPase (SERCA1a) was studied by a combination of rapid scan FTIR spectroscopy and electrostatic calculations. With FTIR spectroscopy, we investigated the pH dependence of C=O bands of the Ca(2+)-free phosphoenzyme (E2P) and obtained direct experimental evidence for the protonation of carboxyl groups upon Ca(2+) release. At least three of the infrared signals from protonated carboxyl groups of E2P are pH dependent with pK(a) values near 8.3: a band at 1758 cm(-1) characteristic of nonhydrogen-bonded carbonyl groups, a shoulder at 1720 cm(-1), and part of a band at 1710 cm(-1), both characteristic of hydrogen-bonded carbonyl groups. The bands are thus assigned to H(+) binding residues, some of which are involved in H(+) countertransport. At pH 9, bands at 1743 and 1710 cm(-1) remain which we do not attribute to Ca(2+)/H(+) exchange. We also obtained evidence for a pH-dependent conformational change in beta-sheet or turn structures of the ATPase. With MCCE on the E2P analog E2(TG+MgF(4)(2-)), we assigned infrared bands to specific residues and analyzed whether or not the carbonyl groups of the acidic Ca(2+) ligands are hydrogen bonded. The carbonyl groups of Glu(771), Asp(800), and Glu(908) were found to be hydrogen bonded and will thus contribute to the lower wave number bands. The carbonyl group of some side-chain conformations of Asp(800) is left without a hydrogen-bonding partner; they will therefore contribute to the higher wave number band.

XH/pi interactions with the pi system of porphyrin ring in porphyrin-containing proteins

Searching structures of porphyrin-containing proteins from the Protein Data Bank revealed that the pi system of every porphyrin ring is involved in XH/pi interactions, with most of the porphyrins having several interactions. Both five-membered pyrrole rings and six-membered chelate rings are involved in XH/pi interactions; the number of interactions with five-membered rings is larger than the number of interactions with six-membered rings. We found interactions with C-H and N-H groups as hydrogen-atom donors; however, the number of CH/pi interactions is much larger than the number of NH/pi interactions. The amino acids involved in the interactions show a high conservation score. Our results that every porphyrin is involved in XH/pi interactions and that amino acids involved in these interactions are highly conserved demonstrate that XH/pi interactions play an important role in porphyrin-protein stability.

Exploring pathways and barriers for coupled ET/PT in cytochrome c oxidase: a general framework for examining energetics and mechanistic alternatives

Gaining a detailed understanding of the energetics of the proton pumping process in cytochrome c oxidase (CcO) is one of the challenges of modern biophysics. Although there are several current mechanistic proposals, most of these ideas have not been subjected to consistent structure-function considerations. In particular most works have not related the activation barriers for different mechanistic proposals to the protein structure. The present work describes a general approach for exploring the energetics of different feasible models of the action of CcO, using the observed protein structure, established simulation methods and a modified Marcus' formulation. We start by reviewing our methods for evaluation of the energy diagrams for different proton translocation paths and then present a systematic analysis of various constraints that should be imposed on any energy diagram for the pumping process. After the general analysis we turn to the actual computational study, where we construct energy diagrams for forward and backward paths, using the estimated calculated reduction potentials and pK(a) values of all the relevant sites (including internal water molecules). We then explore the relationship between the calculated energy diagrams and key experimental constraints. This comparison allows us to identify some barriers that are not fully consistent with the overall requirement for an efficient pumping. In particular we identify back leakage channels, which are hard to block without stopping the forward channels. This helps to identify open problems that will require further experimental and theoretical studies. We also consider reasonable adjustments of the calculated barriers that may lead to a working pump. Although the present analysis does not establish a unique and workable model for the mechanism of CcO, it presents what is probably the most consistent current analysis of the barriers for different feasible pathways. Perhaps more importantly, the framework developed here should provide a general way for examining any proposal for the action of CcO as well as for the analysis of further experimental findings about the action of this fascinating system.

Calculated proton uptake on anaerobic reduction of cytochrome C oxidase: is the reaction electroneutral?

Cytochrome c oxidase is a transmembrane proton pump that builds an electrochemical gradient using chemical energy from the reduction of O(2). Ionization states of all residues were calculated with Multi-Conformation Continuum Electrostatics (MCCE) in seven anaerobic oxidase redox states ranging from fully oxidized to fully reduced. One long-standing problem is how proton uptake is coupled to the reduction of the active site binuclear center (BNC). The BNC has two cofactors: heme a(3) and Cu(B). If the protein needs to maintain electroneutrality, then 2 protons will be bound when the BNC is reduced by 2 electrons in the reductive half of the reaction cycle. The effective pK(a)s of ionizable residues around the BNC are evaluated in Rhodobacter sphaeroides cytochrome c oxidase. At pH 7, only a hydroxide coordinated to Cu(B) shifts its pK(a) from below 7 to above 7 and so picks up a proton when heme a(3) and Cu(B) are reduced. Glu I-286, Tyr I-288, His I-334, and a second hydroxide on heme a(3) all have pK(a)s above 7 in all redox states, although they have only 1.6-3.5 DeltapK units energy cost for deprotonation. Thus, at equilibrium, they are protonated and cannot serve as proton acceptors. The propionic acids near the BNC are deprotonated with pK(a)s well below 7. They are well stabilized in their anionic state and do not bind a proton upon BNC reduction. This suggests that electroneutrality in the BNC is not maintained during the anaerobic reduction. Proton uptake on reduction of Cu(A), heme a, heme a(3), and Cu(B) shows approximately 2.5 protons bound per 4 electrons, in agreement with prior experiments. One proton is bound by a hydroxyl group in the BNC and the rest to groups far from the BNC. The electrochemical midpoint potential (E(m)) of heme a is calculated in the fully oxidized protein and with 1 or 2 electrons in the BNC. The E(m) of heme a shifts down when the BNC is reduced, which agrees with prior experiments. If the BNC reduction is electroneutral, then the heme a E(m) is independent of the BNC redox state.

Factors influencing the energetics of electron and proton transfers in proteins. What can be learned from calculations

A protein structure should provide the information needed to understand its observed properties. Significant progress has been made in developing accurate calculations of acid/base and oxidation/reduction reactions in proteins. Current methods and their strengths and weaknesses are discussed. The distribution and calculated ionization states in a survey of proteins is described, showing that a significant minority of acidic and basic residues are buried in the protein and that most of these remain ionized. The electrochemistry of heme and quinones are considered. Proton transfers in bacteriorhodopsin and coupled electron and proton transfers in photosynthetic reaction centers, 5-coordinate heme binding proteins and cytochrome c oxidase are highlighted as systems where calculations have provided insight into the reaction mechanism.