Theoretical and computational studies of vectorial processes in biomolecular systems

Microscopic basis for kinetic gating in Cytochrome c oxidase: insights from QM/MM analysis

Understanding the mechanism of vectorial proton pumping in biomolecules requires establishing the microscopic basis for the regulation of both thermodynamic and kinetic features of the relevant proton transfer steps.

Understanding the mechanism of vectorial proton pumping in biomolecules requires establishing the microscopic basis for the regulation of both thermodynamic and kinetic features of the relevant proton transfer steps. For the proton pump cytochrome c oxidase, while the regulation of thermodynamic driving force for key proton transfers has been discussed in great detail, the microscopic basis for the control of proton transfer kinetics has been poorly understood. Here we carry out extensive QM/MM free energy simulations to probe the kinetics of relevant proton transfer steps and analyze the effects of local structure and hydration level. We show that protonation of the proton loading site (PLS, taken to be a propionate of heme a₃) requires a concerted process in which a key glutamic acid (Glu286H) delivers the proton to the PLS while being reprotonated by an excess proton coming from the D-channel. The concerted nature of the mechanism is a crucial feature that enables the loading of the PLS before the cavity containing Glu286 is better hydrated to lower its pK_a to experimentally measured range; the charged rather than dipolar nature of the process also ensures a tight coupling with heme a reduction, as emphasized by Siegbahn and Blomberg. In addition, we find that rotational flexibility of the PLS allows its protonation before that of the binuclear center (the site where oxygen gets reduced to water). Together with our recent study (P. Goyal, et al., Proc. Natl. Acad. Sci. U. S. A., 2013, 110, 18886–18891) that focused on the modulation of Glu286 pK_a, the current work suggests a mechanism that builds in a natural sequence for the protonation of the PLS prior to that of the binuclear center. This provides microscopic support to the kinetic constraints revealed by kinetic network analysis as essential elements that ensure an efficient vectorial proton transport in cytochrome c oxidase.

Changing hydration level in an internal cavity modulates the proton affinity of a key glutamate in cytochrome c oxidase

Cytochrome c oxidase contributes to the transmembrane proton gradient by removing two protons from the high-pH side of the membrane each time the binuclear center active site is reduced. One proton goes to the binuclear center, whereas the other is pumped to the low-pH periplasmic space. Glutamate 286 (Glu286) has been proposed to serve as a transiently deprotonated proton donor. Using unrestrained atomistic molecular dynamics simulations, we show that the size of and water distribution in the hydrophobic cavity that holds Glu286 is controlled by the protonation state of the propionic acid of heme a3, a group on the proton outlet pathway. Protonation of the propionate disrupts hydrogen bonding to two side chains, allowing a loop to swing open. Continuum electrostatics and atomistic free-energy perturbation calculations show that the resultant changes in hydration and electrostatic interactions lower the Glu proton affinity by at least 5 kcal/mol. These changes in the internal hydration level occur in the absence of major conformational transitions and serve to stabilize needed transient intermediates in proton transport. The trigger is not the protonation of the Glu of interest, but rather the protonation of a residue ∼10 Å away. Thus, unlike local water penetration to stabilize a new charge, this finding represents a specific role for water molecules in the protein interior, mediating proton transfers and facilitating ion transport.

The mechanism for proton pumping in cytochrome c oxidase from an electrostatic and quantum chemical perspective

The mechanism for proton pumping in cytochrome c oxidase in the respiratory chain, has for decades been one of the main unsolved problems in biochemistry. However, even though several different suggested mechanisms exist, many of the steps in these mechanisms are quite similar and constitute a general consensus framework for discussing proton pumping. When these steps are analyzed, at least three critical gating situations are found, and these points are where the suggested mechanisms in general differ. The requirements for gating are reviewed and analyzed in detail, and a mechanism is suggested, where solutions for all the gating situations are formulated. This mechanism is based on an electrostatic analysis of a kinetic experiment fior the O to E transition. The key component of the mechanism is a positively charged transition state. An electron on heme a opens the gate for proton transfer from the N-side to a pump loading site (PLS). When the negative charge of the electron is compensated by a chemical proton, the positive transition state prevents backflow from the PLS to the N-side at the most critical stage of the pumping process. The mechanism has now been tested by large model DFT calculations, and these calculations give strong support for the suggested mechanism.

Microscopic pKa analysis of Glu286 in cytochrome c oxidase (Rhodobacter sphaeroides): toward a calibrated molecular model

As stringent tests for the molecular model and computational protocol, microscopic pK(a) calculations are performed for the key residue, Glu286, in cytochrome c oxidase (CcO) using a combined quantum mechanical/molecular mechanical (QM/MM) potential and a thermodynamic integration protocol. The impact of the number of water molecules in the hydrophobic cavity and protonation state of several key residues (e.g., His334, Cu(B)-bound water, and PRD(a3)) on the computed microscopic pK(a) values of Glu286 has been systematically examined. To help evaluate the systematic errors in the QM/MM-based protocol, microscopic pK(a) calculations have also been carried out for sites in a soluble protein (Asp70 in T4 lysozyme) and a better-characterized membrane protein (Asp85 in bacteriorhodopsin). Overall, the results show a significant degree of internal consistency and reproducibility that support the effectiveness of the computational framework. Although the number of water molecules in the hydrophobic cavity does not greatly influence the computed pK(a) of Glu286, the protonation states of several residues, some of which are rather far away, have more significant impacts. Adopting the standard protonation state for all titratable residues leaves a large net charge on the system and a significantly elevated pK(a) for Glu286, highlighting that any attempt to address the energetics of proton transfers in CcO at a microscopic level should carefully select the protonation state of residues, even those not in the immediate neighborhood of the active site. The calculations indirectly argue against the deprotonation of His334 for the proton pumping process, although further studies that explicitly compute its pK(a) are required for a more conclusive statement. Finally, the deprotonated Glu286 is found to be in a stable water-mediated connection with PRD(a3) for at least several nanoseconds when this presumed pumping site is protonated. This does not support the proposed role of Glu286 as a robust gating valve that prevents proton leakage, although a conclusive statement awaits a more elaborate characterization of the Glu286-PRD(a3) connectivity with free energy simulations and a protonated PRD(a3). The large sets of microscopic simulations performed here have provided useful guidance to the establishment of a meaningful molecular model and effective computational protocol for explicitly analyzing the proton transfer kinetics in CcO, which is required for answering key questions regarding the pumping function of this fascinating and complex system.

Extension of the Self-Consistent-Charge Density-Functional Tight-Binding Method: Third-Order Expansion of the Density Functional Theory Total Energy and Introduction of a Modified Effective Coulomb Interaction

The standard self-consistent-charge density-functional-tight-binding (SCC-DFTB) method (Phys. Rev. B 1998, 58, 7260) is derived by a second-order expansion of the density functional theory total energy expression, followed by an approximation of the charge density fluctuations by charge monopoles and an effective damped Coulomb interaction between the atomic net charges. The central assumptions behind this effective charge-charge interaction are the inverse relation of atomic size and chemical hardness and the use of a fixed chemical hardness parameter independent of the atomic charge state. While these approximations seem to be unproblematic for many covalently bound systems, they are quantitatively insufficient for hydrogen-bonding interactions and (anionic) molecules with localized net charges. Here, we present an extension of the SCC-DFTB method to incorporate third-order terms in the charge density fluctuations, leading to chemical hardness parameters that are dependent on the atomic charge state and a modification of the Coulomb scaling to improve the electrostatic treatment within the second-order terms. These modifications lead to a significant improvement in the description of hydrogen-bonding interactions and proton affinities of biologically relevant molecules.

pKa analysis for the zinc-bound water in human carbonic anhydrase II: Benchmark for "multiscale" QM/MM simulations and mechanistic implications

To quantitatively explore the applicability of the generalized solvent boundary potential (GSBP) based QM/MM approach as a "multiscale" framework for studying chemical reactions in biomolecules, the structural and energetic properties of the Human Carbonic Anhydrase II (CAII) are analyzed and compared to those from periodic boundary condition (PBC) simulations and available experimental data. Although the atomic fluctuations from GSBP based simulations are consistently lower compared to those from PBC simulations or crystallographic data, the fluctuations and internal coordinate distributions for residues in the proximity of the active site as well as diffusion constants of active-site water molecules are fairly well described by GSBP simulations. The pKa of the zinc-bound water, calculated with a SCC-DFTB/MM-GSBP based thermodynamic integration approach, agrees well with experiments for the wild type CAII. For the E106Q mutant, however, a 9 pKa unit downward shift relative to the wild type is found in contrast with previous experiments that found little change. This dramatic discrepancy signals a possible change in the mechanism for the interconversion between CO2/HCO3- in the E106Q mutant, which may be similar to the bicarbonate mediated mechanism proposed for the Co2+ substituted CAII (J. Am. Chem. Soc. 2001, 123, 5861).1 The study highlights pKa analyses as a valuable approach for quantitatively validating the computational model for complex biomolecules as well as for revealing energetic properties intimately related to the chemical process of interest.

Mechanochemical coupling in the myosin motor domain. I. Insights from equilibrium active-site simulations

Although the major structural transitions in molecular motors are often argued to couple to the binding of Adenosine triphosphate (ATP), the recovery stroke in the conventional myosin has been shown to be dependent on the hydrolysis of ATP. To obtain a clearer mechanistic picture for such "mechanochemical coupling" in myosin, equilibrium active-site simulations with explicit solvent have been carried out to probe the behavior of the motor domain as functions of the nucleotide chemical state and conformation of the converter/relay helix. In conjunction with previous studies of ATP hydrolysis with different active-site conformations and normal mode analysis of structural flexibility, the results help establish an energetics-based framework for understanding the mechanochemical coupling. It is proposed that the activation of hydrolysis does not require the rotation of the lever arm per se, but the two processes are tightly coordinated because both strongly couple to the open/close transition of the active site. The underlying picture involves shifts in the dominant population of different structural motifs as a consequence of changes elsewhere in the motor domain. The contribution of this work and the accompanying paper [] is to propose the actual mechanism behind these "population shifts" and residues that play important roles in the process. It is suggested that structural flexibilities at both the small and large scales inherent to the motor domain make it possible to implement tight couplings between different structural motifs while maintaining small free-energy drops for processes that occur in the detached states, which is likely a feature shared among many molecular motors. The significantly different flexibility of the active site in different X-ray structures with variable level arm orientations supports the notation that external force sensed by the lever arm may transmit into the active site and influence the chemical steps (nucleotide hydrolysis and/or binding).