Effect of a pH Gradient on the Protonation States of Cytochrome c Oxidase: A Continuum Electrostatics Study

Protonation-State Dependence of Hydration and Interactions in the Two Proton-Conducting Channels of Cytochrome c Oxidase

Cytochrome c Oxidase (CcO), a membrane protein of the respiratory chain, pumps protons against an electrochemical gradient by using the energy of oxygen reduction to water. The ("chemical") protons required for this reaction and those pumped are taken up via two distinct channels, named D-channel and K-channel, in a step-wise and highly regulated fashion. In the reductive phase of the catalytic cycle, both channels transport protons so that the pumped proton passes the D-channel before the "chemical" proton has crossed the K-channel. By performing molecular dynamics simulations of CcO in the O→E redox state (after the arrival of the first reducing electron) with various combinations of protonation states of the D- and K-channels, we analysed the effect of protonation on the two channels. In agreement with previous work, the amount of water observed in the D-channel was significantly higher when the terminal residue E286 was not (yet) protonated than when the proton arrived at this end of the D-channel and E286 was neutral. Since a sufficient number of water molecules in the channel is necessary for proton transport, this can be understood as E286 facilitating its own protonation. K-channel hydration shows an even higher dependence on the location of the excess proton in the K-channel. Also in agreement with previous work, the K-channel exhibits a very low hydration level that likely hinders proton transfer when the excess proton is located in the lower part of the K-channel, that is, on the N-side of S365. Once the proton has passed S365 (towards the reaction site, the bi-nuclear centre (BNC)), the amount of water in the K-channel provides hydrogen-bond connectivity that renders proton transfer up to Y288 at the BNC feasible. No significant direct effect of the protonation state of one channel on the hydration level, hydrogen-bond connectivity, or interactions between protein residues in the other channel could be observed, rendering proton conductivity in the two channels independent of each other. Regulation of the order of proton uptake and proton passage in the two channels such that the "chemical" proton leaves its channel last must, therefore, be achieved by other means of communication, such as the location of the reducing electron.

Increasing the Realism of in Silico pHLIP Peptide Models with a Novel pH Gradient CpHMD Method

The pH-low insertion peptides (pHLIP) are pH-dependent membrane inserting peptides, whose function depends on the cell microenvironment acidity. Several peptide variants have been designed to improve upon the wt-sequence, particularly the state transition kinetics and the selectivity for tumor pH. The variant 3 (Var3) peptide is a 27 residue long peptide, with a key titrating residue (Asp-13) that, despite showing a modest performance in liposomes (pK^ins ∼ 5.0), excelled in tumor cell experiments. To help rationalize these results, we focused on the pH gradient in the cell membrane, which is one of the crucial properties that are not present in liposomes. We extended our CpHMD-L method and its pH replica-exchange (pHRE) implementation to include a pH gradient and mimic the pHLIP-membrane microenvironment in a cell where the internal pH is fixed (pH 7.2) and the external pH is allowed to change. We showed that, by properly modeling the pH-gradient, we can correctly predict the experimentally observed loss and gain of performance in tumor cells experiments by the wt and Var3 sequences, respectively. In sum, the pH gradient implementation allowed for more accurate and realistic pK_a estimations and was a pivotal step in bridging the in silico data and the in vivo cell experiments.

Determination of proton concentration at cardiolipin-containing membrane interfaces and its relation with the peroxidase activity of cytochrome c

The interface –log[H⁺] defined as pH′ of a mimic inner mitochondrial membrane is ∼3.9 at bulk pH ∼ 6.8, which affects cytochrome c activity.

The activities of biomolecules are affected by the proton concentrations at biological membranes. Here, we succeeded in evaluating the interface proton concentration (–log[H⁺] defined as pH′) of cardiolipin (CL)-enriched membrane models of the inner mitochondrial membrane (IMM) using a spiro-rhodamine-glucose molecule (RHG). According to fluorescence microscopy and ¹H-NMR studies, RHG interacted with the Stern layer of the membrane. The acid/base equilibrium of RHG between its protonated open form (o-RHG) and deprotonated closed spiro-form (c-RHG) at the membrane interface was monitored with UV-vis absorption and fluorescence spectra. The interface pH′ of 25% cardiolipin (CL)-containing large unilamellar vesicles (LUVs), which possess similar lipid properties to those of the IMM, was estimated to be ∼3.9, when the bulk pH was similar to the mitochondrial intermembrane space pH (6.8). However, for the membranes containing mono-anionic lipids, the interface pH′ was estimated to be ∼5.3 at bulk pH 6.8, indicating that the local negative charges of the lipid headgroups in the lipid membranes are responsible for the deviation of the interface pH′ from the bulk pH. The peroxidase activity of cyt c increased 5–7 fold upon lowering the pH to 3.9–4.3 or adding CL-containing (10–25% of total lipids) LUVs compared to that at bulk pH 6.8, indicating that the pH′ decrease at the IMM interface from the bulk pH enhances the peroxidase activity of cyt c. The peroxidase activity of cyt c at the membrane interface of tetraoleoyl CL (TOCL)-enriched (50% of total lipids) LUVs was higher than that estimated from the interface pH′, while the peroxidase activity was similar to that estimated from the interface pH′ for tetramyristoyl CL (TMCL)-enriched LUVs, supporting the hypothesis that when interacting with TOCL (not TMCL), cyt c opens the heme crevice to substrates. The present simple methodology allows us to estimate the interface proton concentrations of complex biological membranes.

Development of constant-pH simulation methods in implicit solvent and applications in biomolecular systems

pH is a critical parameter for biological and technological systems directly related with electrical charges. It can give rise to peculiar electrostatic phenomena, which also makes them more challenging. Due to the quantum nature of the process, involving the forming and breaking of chemical bonds, quantum methods should ideally by employed. Nevertheless, due to the very large number of ionizable sites, different macromolecular conformations, salt conditions, and all other charged species, the CPU time cost simply becomes prohibitive for computer simulations, making this a quite complex problem. Simplified methods based on Monte Carlo sampling have been devised and will be reviewed here, highlighting the updated state-of-the-art of this field, advantages, and limitations of different theoretical protocols for biomolecular systems (proteins and nucleic acids). Following a historical perspective, the discussion will be associated with the applications to protein interactions with other proteins, polyelectrolytes, and nanoparticles.