Protonation State Research Articles

Perdeuterated GbpA Enables Neutron Scattering Experiments of a Lytic Polysaccharide Monooxygenase.

Lytic polysaccharide monooxygenases (LPMOs) are surface-active redox enzymes that catalyze the degradation of recalcitrant polysaccharides, making them important tools for energy production from renewable sources. In addition, LPMOs are important virulence factors for fungi, bacteria, and viruses. However, many knowledge gaps still exist regarding their catalytic mechanism and interaction with their insoluble, crystalline substrates. Moreover, conventional structural biology techniques, such as X-ray crystallography, usually do not reveal the protonation state of catalytically important residues. In contrast, neutron crystallography is highly suited to obtain this information, albeit with significant sample volume requirements and challenges associated with hydrogen's large incoherent scattering signal. We set out to demonstrate the feasibility of neutron-based techniques for LPMOs using N-acetylglucosamine-binding protein A (GbpA) from Vibrio cholerae as a target. GbpA is a multifunctional protein that is secreted by the bacteria to colonize and degrade chitin. We developed an efficient deuteration protocol, which yields >10 mg of pure 97% deuterated protein per liter expression media, which was scaled up further at international facilities. The deuterated protein retains its catalytic activity and structure, as demonstrated by small-angle X-ray and neutron scattering studies of full-length GbpA and X-ray crystal structures of its LPMO domain (to 1.1 Å resolution), setting the stage for neutron scattering experiments with its substrate chitin.

Insights into the protonation state and spin structure for the g = 2 multiline electron paramagnetic resonance signal of the oxygen-evolving complex.

In photosystem II (PSII), one-electron oxidation of the most stable oxidation state of the Mn4CaO5 cluster (S1) leads to formation of two distinct states, the open-cubane S2 conformation [Mn1(III)Mn2(IV)Mn3(IV)Mn4(IV)] with low spin and the closed-cubane S2 conformation [Mn1(IV)Mn2(IV)Mn3(IV)Mn4(III)] with high spin. In electron paramagnetic resonance (EPR) spectroscopy, the open-cubane S2 conformation exhibits a g = 2 multiline signal. However, its protonation state remains unclear. Here, we investigated the protonation state of the open-cubane S2 conformation by calculating exchange couplings in the presence of the PSII protein environment and simulating the pulsed electron-electron double resonance (PELDOR). When a ligand water molecule, which forms an H-bond with D1-Asp61 (W1), is deprotonated at dangling Mn4(IV), the first-exited energy (34 cm-1) in manifold spin excited states aligns with the observed value in temperature-dependent pulsed EPR analyses, and the PELDOR signal is best reproduced. Consequently, the g = 2 multiline signal observed in EPR corresponds to the open-cubane S2 conformation with the deprotonated W1 (OH-).

Fentanyl and the Fluorinated Fentanyl Derivative NFEPP Elicit Distinct Hydrogen-Bond Dynamics of the Opioid Receptor.

The development of safe therapeutics to manage pain is of central interest for biomedical applications. The fluorinated fentanyl derivative N-(3-fluoro-1-phenethylpiperidin-4-yl)-N-phenylpropionamide (NFEPP) is potentially a safer alternative to fentanyl because unlike fentanyl─which binds to the μ-opioid receptor (MOR) at both physiological and acidic pH─NFEPP might bind to the MOR only at acidic pH typical of inflamed tissue. Knowledge of the protonation-coupled dynamics of the receptor-drug interactions is thus required to understand the molecular mechanism by which receptor activation initiates cell signaling to silence pain. To this end, here we have carried out extensive atomistic simulations of the MOR in different protonation states, in the absence of opioid drugs, and in the presence of fentanyl vs NFEPP. We used graph-based analyses to characterize internal hydrogen-bond networks that could contribute to the activation of the MOR. We find that fentanyl and NFEPP prefer distinct binding poses and that, in their binding poses, fentanyl and NFEPP partake in distinct internal hydrogen-bond networks, leading to the cytoplasmic G-protein-binding region. Moreover, the protonation state of functionally important aspartic and histidine side chains impacts hydrogen-bond networks that extend throughout the receptor, such that the ligand-bound MOR presents at its cytoplasmic G-protein-binding side, a hydrogen-bonding environment where dynamics depend on whether fentanyl or NFEPP is bound, and on the protonation state of specific MOR groups. The exquisite sensitivity of the internal protein-water hydrogen-bond network to the protonation state and to details of the drug binding could enable the MOR to elicit distinct pH- and opioid-dependent responses at its cytoplasmic G-protein-binding site.

In silico evidence of Remdesivir action in blood coagulation cascade modulation in COVID-19 treatment.

Coronavirus disease 2019 (COVID-19) has demonstrated several clinical manifestations which include not only respiratory system issues but also liver, kidney, and other organ injuries. One of these abnormalities is coagulopathies, including thrombosis and disseminated intravascular coagulation. Because of this, the administration of low molecular weight heparin is required for patients that need to be hospitalized. In addition, Remdesivir is an antiviral that was used against Middle East Acute Respiratory Syndrome, Ebola, Acute Respiratory Syndrome, and other diseases, showing satisfactory results on recovery. Besides, there is evidence suggesting that this medication can provide a better prognosis for patients with COVID-19. To investigate in silico the interaction between Remdesivir and clotting factors, pursuing a possibility of using it as medicine. In this in silico study, the 3D structures of angiotensin-converting enzyme 2 (ACE2), Factor I (fibrinogen), Factor II (prothrombin), Factor III (thromboplastin), Factor V (proaccelerin), Factor VII (proconvertin), Factor VIII (antihemophilic factor A), Factor IX (antihemophilic factor B), Factor X (Stuart-Prower factor), and Factor XI (precursor of thromboplastin (these structures are technically called receptors) were selected from the Protein Data Bank. The structures of the antivirals Remdesivir and Osetalmivir (these structures are called ligands) were selected from the PubChem database, while the structure of Atazanavir was selected from the ZINC database. The software AutoDock Tools (ADT) was used to prepare the receptors for molecular docking. Ions, peptides, water molecules, and other ones were removed from each ligand, and then, hydrogen atoms were added to the structures. The grid box was delimited and calculated using the same software ADT. A physiological environment with pH 7.4 is needed to make the ligands interact with the receptors, and still the software Marvin sketch® (ChemAxon®) was used to forecast the protonation state. To perform molecular docking, ADT and Vina software was connected. Using PyMol® software and Discovery studio® software from BIOVIA, it was possible to analyze the amino acid residues from receptors that were involved in the interactions with the ligands. Ligand tortions, atoms that participated in the interactions, and the type, strength, and duration of the interactions were also analyzed using those software. Molecular docking analysis showed that Remdesivir and ACE2 had an affinity energy of -8.8 kcal/moL, forming a complex with eight hydrogen bonds involving seven atoms of Remdesivir and five amino acid residues of ACE2. Remdesivir and prothrombin had an interaction with six hydrogen bonds involving atoms of the drug and five amino acid residues of the clotting factor. Similar to that, Remdesivir and thromboplastin presented interactions via seven hydrogen bonds involving five atoms of the drug and four residues of the clotting factor. While Remdesivir and Factor V established a complex with seven hydrogen bonds between six antiviral atoms and six amino acid residues from the factor, and Factor VII connected with the drug by four hydrogen bonds, which involved three atoms of the drug and three residues of amino acids of the factor. The complex between Remdesivir and Factor IX formed an interaction via 11 hydrophilic bonds with seven atoms of the drug and seven residues of the clotting factor, plus one electrostatic bond and three hydrophobic interactions. Factor X and Remdesivir had an affinity energy of -9.6 kcal/moL, and the complex presented 10 hydrogen bonds and 14 different hydrophobic interactions which involved nine atoms of the drug and 16 amino acid residues of the clotting factor. The interaction between Remdesivir and Factor XI formed five hydrogen bonds involving five amino acid residues of the clotting factor and five of the antiviral atoms. Because of the in silico significant affinity, Remdesivir possibly could act in the severe acute respiratory syndrome coronavirus 2 infection blockade by interacting with ACE2 and concomitantly act in the modulation of the coagulation cascade preventing the hypercoagulable state.

Search for toverline{t}H/Ato toverline{t}toverline{t} production in the multilepton final state in proton–proton collisions at sqrt{s} = 13 TeV with the ATLAS detector

A search for a new heavy scalar or pseudo-scalar Higgs boson (H/A) produced in association with a pair of top quarks, with the Higgs boson decaying into a pair of top quarks (H/A → toverline{t} ) is reported. The search targets a final state with exactly two leptons with same-sign electric charges or at least three leptons. The analysed dataset corresponds to an integrated luminosity of 139 fb−1 of proton–proton collisions collected at a centre-of-mass energy of 13 TeV with the ATLAS detector at the LHC. Two multivariate classifiers are used to separate the signal from the background. No significant excess of events over the Standard Model expectation is observed. The results are interpreted in the context of a type-II two-Higgs-doublet model. The observed (expected) upper limits at 95% confidence level on the toverline{t}H/A production cross-section times the branching ratio of H/A → toverline{t} range between 14 (10) fb and 6 (5) fb for a heavy Higgs boson with mass between 400 GeV and 1000 GeV, respectively. Assuming that only one particle, either the scalar H or the pseudo-scalar A, contributes to the toverline{t}toverline{t} final state, values of tan β below 1.2 or 0.5 are excluded for a mass of 400 GeV or 1000 GeV, respectively. These exclusion ranges increase to tan β below 1.6 or 0.6 when both particles are considered.

Structural and dynamic properties of urinary human serum albumin fragments: a molecular dynamics study

A microalbuminuria level acts as a good index to screen and monitor diabetes and renal failure. However, the urinary albumin loss after sample preservation and storage is the major bottleneck to obtain the accurate microalbuminuria test. Such loss is due to the rapid albumin fragmentation by urinary proteases. Some fragments were suggested to be bioactive biomarkers of diabetes and renal disease, but no structural and dynamical properties of albumin fragments are available. Thus, in this work, the structural and dynamical properties of reported albumin fragments are revealed using molecular dynamics simulations. The properties of nine fragments (F1–F9) discovered recently were studied at the real pH conditions of urine samples (pH 4.5, 7 and 8). The complete loss of secondary structure is found in short fragments (F1–F6), while large-sized polypeptides (F7–F9) can somehow maintain their folds. Especially, F8 (subdomain IIIB) is the most stable fragment. The difference in histidine protonation states has no impact on the structural stability of albumin fragments. The ability of F8 (subdomain IIIB) to maintain its stability and folds suggests it as an alternative albumin biomarker in urine. An insight obtained here will become the fundamental importance for understanding clinical assays for albumin detection, sample stability and peptidomics analysis of urine. Communicated by Ramaswamy H. Sarma

PH-Dependent Binding and Releasing Mechanism of Acetate in the Inner Water Cavity of Heliorhodopsin.

The high-resolution structure of heliorhodopsin crystallized at low pH reveals the presence of a planar triangle molecule, acetate, in the inner water cavity. Here, we investigate how the acetate molecule is stabilized at the counterion Glu107 moiety, using molecular dynamics (MD) simulations and a quantum mechanical/molecular mechanical (QM/MM) approach. QM/MM calculations indicate that the density is best described as acetate among triangle acids, including nitric acid and bicarbonate. The calculated protonation state indicates that protonated acetate donates an H-bond to deprotonated Glu107 in the low-pH crystal structure. The observed red-shift of ∼30 nm in the absorption wavelength with pKa ≈ 4 is likely due to the His23/His80 protonation, rather than the Glu107 protonation. MD simulations also show that acetate can exist at the Glu107 moiety only when it is protonated. When ionized, acetate is released from the Glu107 moiety via Asn101 at the channel bottleneck and Arg91 on the intracellular protein surface. These observations could explain how acetate binds at low pH and releases at high pH.

Search for resonant WZ production in the fully leptonic final state in proton–proton collisions at mathbf {sqrt{s} = 13} TeV with the ATLAS detector

A search for a WZ resonance, in the fully leptonic final state (electrons or muons), is performed using 139 fb^{-1} of data collected at a centre-of-mass energy of 13 TeV by the ATLAS detector at the Large Hadron Collider. The results are interpreted in terms of a singly charged Higgs boson of the Georgi–Machacek model, produced by WZ fusion, and of a Heavy Vector Triplet, with the resonance produced by WZ fusion or the Drell–Yan process. No significant excess over the Standard Model prediction is observed and limits are set on the production cross-section times branching ratio as a function of the resonance mass for these processes.

Molecular Structure Related to an S = 5/2 High-Spin S2 State Manganese Cluster of Photosystem II Investigated by Q-Band Pulse EPR Spectroscopy.

The high-spin S2 state of the photosynthetic oxygen-evolving cluster Mn4CaO5, corresponding to the g = 4.1 signal for X-band electron paramagnetic resonance (EPR), was investigated using Q-band pulsed EPR, which detected a main peak at g = 3.10 and satellite peaks at 5.25, 4.55, and 2.80. We evaluated the spin state as the zero-field splitting of D = 0.465 cm-1 and E/D = 0.245 with S = 5/2. The temperature dependence of the T1 relaxation time revealed that the excited-state energy was 28.7 cm-1 higher than that of the S = 5/2 ground state. By comparing present quantum mechanical (QM) calculation models, a closed-cubane structure with the protonation state of two oxygens, W1 (= OH-) and W2 (= H2O), was the most probable structure for the S = 5/2 state. The three-pulse electron spin-echo envelope modulation (ESEEM) detected the nuclear signal, which was assigned to nitrogen as His332 ligated to the Mn1 ion. The obtained hyperfine constant for the nitrogen signal was significantly reduced from that in the S = 1/2 low-spin state. These results indicate that the S = 5/2 spin state arises from the closed-cubane structure.

Automated Parameterization of Coarse-Grained Polyethylenimine under a Martini Framework.

As a versatile polymer in many applications, synthesized polyethylenimine (PEI) is polydisperse with diverse branched structures that attain pH-dependent protonation states. Understanding the structure-function relationship of PEI is necessary for enhancing its efficacy in various applications. Coarse-grained (CG) simulations can be performed at length and time scales directly comparable with experimental data while maintaining the molecular perspective. However, manually developing CG forcefields for complex PEI structures is time-consuming and prone to human errors. This article presents a fully automated algorithm that can coarse-grain any branched architecture of PEI from its all-atom (AA) simulation trajectories and topology. The algorithm is demonstrated by coarse-graining a branched 2 kDa PEI, which can replicate the AA diffusion coefficient, radius of gyration, and end-to-end distance of the longest linear chain. Commercially available 25 and 2 kDa Millipore-Sigma PEIs are used for experimental validation. Specifically, branched PEI architectures are proposed, coarse-grained using the automated algorithm, and then simulated at different mass concentrations. The CG PEIs can reproduce existing experimental data on PEI's diffusion coefficient and Stokes-Einstein radius at infinite dilution as well as its intrinsic viscosity. This suggests a strategy where probable chemical structures of synthetic PEIs can be inferred computationally using the developed algorithm. The coarse-graining methodology presented here can also be extended to other polymers.

Visualization of the human manganese superoxide dismutase product inhibition mechanism and protonation states

Visualization of the human manganese superoxide dismutase product inhibition mechanism and protonation states

Dual Piperidine-Based Histamine H3 and Sigma-1 Receptor Ligands in the Treatment of Nociceptive and Neuropathic Pain.

In search of new dual-acting histamine H3/sigma-1 receptor ligands, we designed a series of compounds structurally based on highly active in vivo ligands previously studied and described by our team. However, we kept in mind that within the previous series, a pair of closely related compounds, KSK67 and KSK68, differing only in the piperazine/piperidine moiety in the structural core showed a significantly different affinity at sigma-1 receptors (σ1Rs). Therefore, we first focused on an in-depth analysis of the protonation states of piperazine and piperidine derivatives in the studied compounds. In a series of 16 new ligands, mainly based on the piperidine core, we selected three lead structures (3, 7, and 12) for further biological evaluation. Compound 12 showed a broad spectrum of analgesic activity in both nociceptive and neuropathic pain models based on the novel molecular mechanism.

Visualizing protonation states in serine hydroxymethyltransferase with neutron crystallography

Visualizing protonation states in serine hydroxymethyltransferase with neutron crystallography

Molecular details of ovalbumin solvation by an aqueous solution of xylitol in different pH environment:Ultrasonic and molecular simulation studies

Artificial sweetener-protein interactions are playing a major role in enhancing the functional properties of protein-based food and pharmaceutical products. This article discusses the molecular details of ovalbumin (OA) solvation in an aqueous solution of xylitol (XY) in different pH environments. This task can be accomplished by both Ultrasonic and Molecular Simulation studies. In ultrasonic studies, density, ultrasound velocity, viscosity, and surface tension have been measured for OA in various pH of phosphate buffer with and without XY. Further, some of the calculated parameters such as adiabatic compressibility, intermolecular free length, free volume, etc have been worked out. From ultrasonic results, pH dependent preferential hydration of OA was noticed in the presence of XY in an aqueous environment. In Molecular Dynamics Simulation (MDS) studies, water and XY accumulation to OA surface for different pH models were analyzed. The different states of protonation of the OA surface promote two distinct effects on the solvent and cosolvent orientation.

Monitoring Conformation and Protonation States of Glutathione by Raman Optical Activity and Molecular Dynamics.

Glutathione (GSH) is a common antioxidant and its biological activity depends on the conformation and protonation state. We used molecular dynamics, Raman and Raman optical activity (ROA) spectroscopies to investigate GSH structural changes in a broad pH range. Factor analysis of the spectra provided protonation constants (2.05, 3.45, 8.62, 9.41) in good agreement with previously published values. Following the analysis, spectra of differently protonated forms were obtained by extrapolation. The complete deprotonation of the thiol group above pH 11 was clearly visible in the spectra; however, many spectral features did not change much with pH. Experimental spectra at various pH values were decomposed into the simulated ones, which allowed us to study the conformer populations and quality of molecular dynamics (MD). According to this combined ROA/MD analysis conformation of the GSH backbone is affected by the pH changes only in a limited way. The combination of ROA with the computations thus has the potential to improve the MD force field and obtain more accurate populations of the conformer species. The methodology can be used for any molecule, but for a more detailed insight better computational techniques are needed in the future.

Molecular Dynamics Simulation Study of the Protonation State Dependence of Glutamic Acid Transport through a Cyclic Peptide Nanotube.

The effect of the protonation state of glutamic acid on its translocation through cyclic peptide nanotubes (CPNs) was assessed by using molecular dynamics (MD) simulations. Anionic (GLU-), neutral zwitterionic (GLU0), and cationic (GLU+) forms of glutamic acid were selected as three different protonation states for an analysis of energetics and diffusivity for acid transport across a cyclic decapeptide nanotube. Based on the solubility-diffusion model, permeability coefficients for the three protonation states of the acid were calculated and compared with experimental results for CPN-mediated glutamate transport through CPNs. Potential of mean force (PMF) calculations reveal that, due to the cation-selective nature of the lumen of CPNs, GLU-, so-called glutamate, shows significantly high free energy barriers, while GLU+ displays deep energy wells and GLU0 has mild free energy barriers and wells inside the CPN. The considerable energy barriers for GLU- inside CPNs are mainly attributed to unfavorable interactions with DMPC bilayers and CPNs and are reduced by favorable interactions with channel water molecules through attractive electrostatic interactions and hydrogen bonding. Unlike the distinct PMF curves, position-dependent diffusion coefficient profiles exhibit comparable frictional behaviors regardless of the charge status of three protonation states due to similar confined environments imposed by the lumen of the CPN. The calculated permeability coefficients for the three protonation states clearly demonstrate that glutamic acid has a strong protonation state dependence for its transport through CPNs, as determined by the energetics rather than the diffusivity of the protonation state. In addition, the permeability coefficients also imply that GLU- is unlikely to pass through a CPN due to the high energy barriers inside the CPN, which is in disagreement with experimental measurements, where a considerable amount of glutamate permeating through the CPN was detected. To resolve the discrepancy between this work and the experimental observations, several possibilities are proposed, including a large concentration gradient of glutamate between the inside and outside of lipid vesicles and bilayers in the experiments, the glutamate activity difference between our MD simulations and experiments, an overestimation of energy barriers due to the artifacts imposed in MD simulations, and/or finally a transformation of the protonation state from GLU- to GLU0 to reduce the energy barriers. Overall, our study demonstrates that the protonation state of glutamic acid has a strong effect on the transport of the acid and suggests a possible protonation state change for glutamate permeating through CPNs.

Effect of Polypeptide Complex Coacervate Microenvironment on Protonation of a Guest Molecule.

Complex coacervate droplets formed by the liquid-liquid phase separation of polyelectrolyte solutions capture several important features of membraneless organelles including their ability to accumulate guest molecules and to provide distinct microenvironments. Here, we examine how polyions in complex coacervates can influence localized guest molecules, leading to a shifted protonation state of the guest molecule in response to its electrostatic environment. A fluorescent ratiometric pH indicator dye was used as a model guest molecule able to report its protonation state in the coacervate phase. Experimentally observed differences in dye-reported local apparent pH inside versus outside of coacervate droplets were largest for polyion pairs having lower salt stabilities and/or larger polyion length mismatch, which we attribute to the relative concentration of open sites on polyions within the coacervates based on theoretical calculations. Using the transfer matrix method, we confirmed that theoretical phase diagrams and critical salt stabilities generated for each polyion pair were consistent with experimental turbidity measurements and estimated the amount of available binding sites on polyions for guest molecules. We conclude that dye molecules likely experience an effective pKa shift due to interactions with coacervate polyions rather than reporting directly on local proton concentrations. Such a local pKa shift can also be anticipated for other guest molecules having protonatable groups, including, for example, many metabolites, ligands, and/or drug molecules that partition into coacervates or membraneless organelles based on ion pairing interactions.

Constructing the Mechanism of Dinoflagellate Luciferin Bioluminescence Using Computation.

Dinoflagellate luciferin bioluminescence is unique since it does not rely on decarboxylation but is poorly understood compared to that of firefly, bacteria, and coelenterata luciferins. Here we computationally investigate possible protonation states, stereoisomers, a chemical mechanism, and the dynamics of the bioluminescence intermediate that is responsible for chemiexcitation. Using semiempirical dynamics, time-dependent density functional theory static calculations, and a correlation diagram, we find that the intermediate's functional group that is likely responsible for chemiexcitation is a 4-member ring, a dioxetanol, that undergoes [2π + 2π] cycloreversion and the biolumiphore is the cleaved structure. The simulated emission spectra and luciferase-dependent absorbance spectra agree with the experimental data, giving support to our proposed mechanism and biolumiphore. We also compute circular dichroism spectra of the intermediate's four stereoisomers to guide future experiments in differentiating them.

A Remarkably Unsymmetric Hexairon Core Embraced by Two High-Symmetry Tripodal Oligo-α-pyridylamido Ligands.

Oligo-α-pyridylamides offer an appealing route to polyiron complexes with short Fe-Fe separations and large room-temperature magnetic moments. A derivative of tris(2-aminoethyl)amine (H6tren) containing three oligo-α-pyridylamine branches and 13 nitrogen donors (H6L) reacts with [Fe2(Mes)4] to yield an organic nanocage built up by two tripodal ligands with interdigitated branches (HMes = mesitylene). The nanocage has crystallographic D3 symmetry but hosts a remarkably unsymmetric hexairon-oxo core, with a central Fe5(μ5-O) square pyramid, two oxygen donors bridging basal sites, and an additional Fe center residing in one of the two tren-like pockets. Bond valence sum (BVS) analysis, density functional theory (DFT) calculations, and electrochemical data were then used to establish the protonation state of oxygen atoms and the formal oxidation states of the metals. For this purpose, a specialized set of BVS parameters was devised for Fe2+-N3- bonds with nitrogen donors of oligo-α-pyridylamides. This allowed us to formulate the compound as [Fe6O2(OH)(H3L)L], with nominally four FeII ions and two FeIII ions. Mössbauer spectra indicate that the compound contains two unique FeII sites, identified as a pair of closely spaced hydroxo-bridged metal ions in the central Fe5(μ5-O) pyramid, and a substantially valence-delocalized FeII2FeIII2 unit. Broken-symmetry DFT calculations predict strong ferromagnetic coupling between the two iron(II) ions, leading to a local S = 4 state that persists to room temperature and explaining the large magnetic moment measured at 300 K. The compound behaves as a single-molecule magnet, with magnetization dynamics detectable in zero static field and dominated by an Orbach-like mechanism with activation parameters Ueff/kB = 49(2) K and τ0 = 4(2) × 10-10 s.

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.