Lysine 362 in cytochrome c oxidase regulates opening of the K-channel via changes in pKA and conformation.

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.

Computing pK(A) values of hexa-aqua transition metal complexes.

Aqueous pKA values for 15 hexa-aqua transition metal complexes were computed using a combination of quantum chemical and electrostatic methods. Two different structure models were considered optimizing the isolated complexes in vacuum or in presence of explicit solvent using a QM/MM approach. They yield very good agreement with experimentally measured pKA values with an overall root mean square deviation of about 1 pH unit, excluding a single but different outlier for each of the two structure models. These outliers are hexa-aqua Cr(III) for the vacuum and hexa-aqua Mn(III) for the QM/MM structure model. Reasons leading to the deviations of the outlier complexes are partially explained. Compared to previous approaches from the same lab the precision of the method was systematically improved as discussed in this study. The refined methods to obtain the appropriate geometries of the complexes, developed in this work, may allow also the computation of accurate pKA values for multicore transition metal complexes in different oxidation states.

Proton transfer in the K-channel analog of B-type Cytochrome c oxidase from Thermus thermophilus.

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.

Two-photon-induced singlet fission in rubrene single crystal.

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.

Computations of 36 tautomer/isomer equilibria of different lactams.

Gas-phase energies of 36 tautomer/isomer pairs of 18 six-membered N-heterocyclic compounds were computed quantum chemically. Among the considered B3LYP, BH&HLYP, BH&HLYP(G), and PW6B95 DFT functionals, the latter two provide accurate tautomer/isomer pair energies with root-mean-square deviations (rmsd) relative to experiments of 0.2 and 0.3 kcal/mol, respectively. Since only few (namely five) experimental data are available, 15 tautomer/isomer pair energies were computed with the very precise QCISD(T)(quadruple-ζ) method serving as reference. Relative to this reference the PW6B95 DFT functional is slightly superior to the BH&HLYP(G) functional, yielding an rmsd of 0.7 and 0.8 kcal/mol, respectively. In contrast to BH&HLYP(G), the PW6B95 DFT functional yields also accurate tautomer/isomer pair energies if zwitterionic structures are involved. The tautomer/isomer pair states possess different amounts of aromaticity. This is characterized by nucleus-independent chemical shift (NICS) values. The tautomer/isomer pair reference energies, from which the energies computed with PW6B95 are subtracted, correlate linearly with the corresponding differences in the NICS values. This correlation is used to construct a correction term for the pair energies computed with PW6B95, yielding tautomer/isomer pair energies with rmsd of 0.3 kcal/mol with respect to the more CPU time demanding QCISD(T)(quadruple-ζ) method.

Understanding selectin counter-receptor binding from electrostatic energy computations and experimental binding studies.

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.

Exploring the possible role of Glu286 in CcO by electrostatic energy computations combined with molecular dynamics.

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.