The metabolism of aerobic life uses the conversion of molecular oxygen to water as an energy source. This reaction is catalyzed by cytochrome e oxidase (CeO) consuming four electrons and four protons, which move along specific routes. While all four electrons are transferred via the same cofactors to the binuclear reaction center (BNC), the protons take two different routes in the A-type CeO, i.e., two of the four chemical protons consumed in the reaction arrive via the D-channel in the oxidative first half starting after oxygen binding. The other two chemical protons enter via the K-channel in the reductive second half of the reaction cycle. To date, the mechanism behind these separate proton transport pathways has not been understood. In this study, we propose a model that can explain the reaction-step specific opening and closing of the K-channel by conformational and pKA changes of its central lysine 362. Molecular dynamics simulations reveal an upward movement of Lys362 towards the BNC, which had already been supposed by several experimental studies. Redox state-dependent pKA calculations provide evidence that Lys362 may protonate transiently, thereby opening the K-channel only in the reductive second half of the reaction cycle. From our results, we develop a model that assigns a key role to Lys362 in the proton gating between the two proton input channels of the A-type CeO.

The two-photon-induced singlet fission was observed in rubrene single crystal and studied by use of femtosecond pump-probe spectroscopy. The location of two-photon excited states was obtained from the nondegenerate two-photon absorption (TPA) spectrum. Time evolution of the two-photon-induced transient absorption spectra reveals the direct singlet fission from the two-photon excited states. The TPA absorption coefficient of rubrene single crystal is 52 cm∕GW at 740 nm, as obtained from Z-scan measurements. Quantum chemical calculations based on time-dependent density functional theory support our experimental data.

Cytochrome c oxidase (CcO) is a central enzyme in aerobic life catalyzing the conversion of molecular oxygen to water and utilizing the chemical energy to pump protons and establish an electrochemical gradient. Despite intense research, it is not understood how CcO achieves unidirectional proton transport and avoids short circuiting the proton pump. Within this work, we analyzed the potential role of Glu286 as a proton valve. We performed unconstrained MD simulations of CcO with an explicit membrane for up to 80 ns. Those MD simulations revealed that deprotonated Glu286 (Glu286-) is repelled by the negatively charged propionic acid PRD of heme a3. Thus, it destabilizes a potential linear chain of waters in the hydrophobic cavity connecting Glu286 with PRD and the binuclear center (BNC). Conversely, protonated Glu286 (Glu286H) may remain in an upward position (oriented toward PRD) and can stabilize the connecting linear water chain in the hydrophobic cavity. We calculated the pKa of Glu286 under physiological conditions to be above 12, but this value decreases to about 9 under increased water accessibility of Glu286. The latter value is in accordance with experimental measurements. In the time course of MD simulation, we also observed conformations where Glu286 bridges between water molecules located on both sides (the D channel being connected to the N side and the hydrophobic cavity), which might lead to proton backflow.

A key enzyme in aerobic metabolism is cytochrome c oxidase (CcO), which catalyzes the reduction of molecular oxygen to water in the mitochondrial and bacterial membranes. Substrate electrons and protons are taken up from different sides of the membrane and protons are pumped across the membrane, thereby generating an electrochemical gradient. The well-studied A-type CcO uses two different entry channels for protons: the D-channel for all pumped and two consumed protons, and the K-channel for the other two consumed protons. In contrast, the B-type CcO uses only a single proton input channel for all consumed and pumped protons. It has the same location as the A-type K-channel (and thus is named the K-channel analog) without sharing any significant sequence homology. In this study, we performed molecular-dynamics simulations and electrostatic calculations to characterize the K-channel analog in terms of its energetic requirements and functionalities. The function of Glu-15B as a proton sink at the channel entrance is demonstrated by its rotational movement out of the channel when it is deprotonated and by its high pKA value when it points inside the channel. Tyr-244 in the middle of the channel is identified as the valve that ensures unidirectional proton transfer, as it moves inside the hydrogen-bond gap of the K-channel analog only while being deprotonated. The electrostatic energy landscape was calculated for all proton-transfer steps in the K-channel analog, which functions via proton-hole transfer. Overall, the K-channel analog has a very stable geometry without large energy barriers.

Higher organisms defend themselves against invading micro-organisms and harmful substances with their immune system. Key players of the immune system are the white blood cells (WBC), which in case of infection move in an extravasation process from blood vessels toward infected tissue promoting inflammation. This process starts with the attachment of the WBC to the blood vessel wall, mediated by protein pair interactions of selectins and counter-receptors (C-R). Individual selectin C-R binding is weak and varies only moderately between the three selectin types. Multivalency enhances such small differences, rendering selectin-binding type specific. In this work, we study selectin C-R binding, the initial step of extravasation. We performed electrostatic energy computations based on the crystal structure of one selectin type co-crystallized with the ligating part of the C-R. The agreement with measured free energies of binding is satisfactory. Additionally, we modeled selectin mutant structures in order to explain differences in binding of the different selectin types. To verify our modeling procedures, surface plasmon resonance data were measured for several mutants and compared with computed binding affinities. Binding affinities computed with soaked rather than co-crystallized selectin C-R structures do not agree with measured data. Hence, these structures are inappropriate to describe the binding mode. The analysis of selectin/C-R binding unravels the role played by individual molecular components in the binding event. This opens new avenues to prevent immune system malfunction, designing drugs that can control inflammatory processes by moderating selectin C-R binding.