Model for Proton Transport Coupled to Protein Conformational Change: Application to Proton Pumping in the Bacteriorhodopsin Photocycle

Quantifying Bacteriorhodopsin Activity as a Function of its Local Environment with a Raman-Based Assay

Bacteriorhodopsin (bR) is a transmembrane protein that functions as a light-driven proton pump in halophilic archaea. The bR photocycle has been well-characterized; however, these measurements almost exclusively measured purified bR, outside of its native membrane. To investigate what effect the cellular environment has on the bR photocycle, we have developed a Raman-based assay that can monitor the activity of the bR in a variety of conditions, including in its native membrane. The assay uses two continuous-wave lasers, one to initiate photochemistry and one to monitor bR activity. The excitation leads to the steady-state depletion of ground-state bR, which directly relates to the population of photocycle intermediate states. We have used this assay to monitor bR activity both in vitro and in vivo. Our in vitro measurements confirm that our assay is sensitive to bulk environmental changes reported in the literature. Our in vivo measurements show a decrease in bR activity with increasing extracellular pH for bR in its native membrane. The difference in activity with increasing pH indicates that the native membrane environment affects the function of bR. This assay opens the door to future measurements into understanding how the local environment of this transmembrane protein affects function.

Concerted Motions and Molecular Function: What Physical Chemistry We Can Learn from Light-Driven Ion-Pumping Rhodopsins

Transmembrane ion gradients are generated and maintained by ion-pumping proteins in cells. Light-driven ion-pumping rhodopsins are retinal-containing proteins found in archaea, bacteria, and eukarya. Photoisomerization of the retinal chromophore induces structural changes in the protein, allowing the transport of ions in a particular direction. Understanding unidirectional ion transport by ion-pumping rhodopsins is an exciting challenge for biophysical chemistry. Concerted changes in ion-binding affinities of the ion-binding sites in proteins are key to unidirectional ion transport, as is the coupling between the chromophore and the protein moiety to drive the concerted motions regulating ion-binding affinities. The commonality of ion-pumping rhodopsin protein structures and the diversity of their ion-pumping functions suggest universal principles governing ion transport, which would be widely applicable to molecular systems. In this Perspective, I review the insights obtained from previous studies on rhodopsins and discuss future perspectives.

Thermodynamics of transport through the ammonium transporter Amt-1 investigated with free energy calculations

Amt-1 from Archaeoglobus fulgidus (AfAmt-1) belongs to the Amt/Rh family of ammonium/ammonia transporting membrane proteins. The transport mode and the precise microscopic permeation mechanism utilized by these proteins are intensely debated. Open questions concern the identity of the transported substrate (ammonia and/or ammonium) and whether the transport is passive or active. To address these questions, we studied the overall thermodynamics of the different transport modes as a function of the environmental conditions. Then, we investigated the thermodynamics of the underlying microscopic transport mechanisms with free energy calculations within a continuum electrostatics model. The formalism developed for this purpose is of general utility in the calculation of binding free energies for ligands with multiple protonation forms or other binding forms. The results of our calculations are compared to the available experimental and theoretical data on Amt/Rh proteins and discussed in light of the current knowledge on the physiological conditions experienced by microorganisms and plants. We found that microscopic models of electroneutral and electrogenic transport modes are in principle thermodynamically viable. However, only the electrogenic variants have a net thermodynamic driving force under the physiological conditions experienced by microorganisms and plants. Thus, the transport mechanism of AfAmt-1 is most likely electrogenic.

GMCT : a Monte Carlo simulation package for macromolecular receptors

Generalized Monte Carlo titration (GMCT) is a versatile suite of computer programs for the efficient simulation of complex macromolecular receptor systems as for example proteins. The computational model of the system is based on a microstate description of the receptor and an average description of its surroundings in terms of chemical potentials. The receptor can be modeled in great detail including conformational flexibility and many binding sites with multiple different forms that can bind different ligand types. Membrane embedded systems can be modeled including electrochemical potential gradients. Overall properties of the receptor as well as properties of individual sites can be studied with a variety of different Monte Carlo (MC) simulation methods. Metropolis MC, Wang-Landau MC and efficient free energy calculation methods are included. GMCT is distributed as free open source software at www.bisb.uni-bayreuth.de under the terms of the GNU Affero General Public License.

Artificial photosynthetic reaction centers coupled to light-harvesting antennas

We analyze a theoretical model for energy and electron transfer in an artificial photosynthetic system. The photosystem consists of a molecular triad (i.e., with a donor, a photosensitive unit, and an acceptor) coupled to four accessory light-harvesting-antenna pigments. The resonant energy transfer from the antennas to the artificial reaction center (the molecular triad) is described here by the Förster mechanism. We consider two different kinds of arrangements of the accessory light-harvesting pigments around the reaction center. The first arrangement allows direct excitation transfer to the reaction center from all the surrounding pigments. The second configuration transmits energy via a cascade mechanism along a chain of light-harvesting chromophores, where only one chromophore is connected to the reaction center. We show that the artificial photosynthetic system using the cascade energy transfer absorbs photons in a broader wavelength range and converts their energy into electricity with a higher efficiency than the system based on direct couplings between all the antenna chromophores and the reaction center.

Diffusion-controlled generation of a proton-motive force across a biomembrane

Respiration in bacteria involves a sequence of energetically coupled electron and proton transfers creating an electrochemical gradient of protons (a proton-motive force) across the inner bacterial membrane. With a simple kinetic model, we analyze a redox loop mechanism of proton-motive force generation mediated by a molecular shuttle diffusing inside the membrane. This model, which includes six electron-binding and two proton-binding sites, reflects the main features of nitrate respiration in E. coli bacteria. We describe the time evolution of the proton translocation process. We find that the electron-proton electrostatic coupling on the shuttle plays a significant role in the process of energy conversion between electron and proton components. We determine the conditions where the redox loop mechanism is able to translocate protons against the transmembrane voltage gradient above 200 mV with a thermodynamic efficiency of about 37%, in the physiologically important range of temperatures from 250 to 350 K.

Simulating the proton transfer in gramicidin A by a sequential dynamical Monte Carlo method

The large interest in long-range proton transfer in biomolecules is triggered by its importance for many biochemical processes such as biological energy transduction and drug detoxification. Since long-range proton transfer occurs on a microsecond time scale, simulating this process on a molecular level is still a challenging task and not possible with standard simulation methods. In general, the dynamics of a reactive system can be described by a master equation. A natural way to describe long-range charge transfer in biomolecules is to decompose the process into elementary steps which are transitions between microstates. Each microstate has a defined protonation pattern. Although such a master equation can in principle be solved analytically, it is often too demanding to solve this equation because of the large number of microstates. In this paper, we describe a new method which solves the master equation by a sequential dynamical Monte Carlo algorithm. Starting from one microstate, the evolution of the system is simulated as a stochastic process. The energetic parameters required for these simulations are determined by continuum electrostatic calculations. We apply this method to simulate the proton transfer through gramicidin A, a transmembrane proton channel, in dependence on the applied membrane potential and the pH value of the solution. As elementary steps in our reaction, we consider proton uptake and release, proton transfer along a hydrogen bond, and rotations of water molecules that constitute a proton wire through the channel. A simulation of 8 mus length took about 5 min on an Intel Pentium 4 CPU with 3.2 GHz. We obtained good agreement with experimental data for the proton flux through gramicidin A over a wide range of pH values and membrane potentials. We find that proton desolvation as well as water rotations are equally important for the proton transfer through gramicidin A at physiological membrane potentials. Our method allows to simulate long-range charge transfer in biological systems at time scales, which are not accessible by other methods.

A fast and accurate computational approach to protein ionization

We report a very fast and accurate physics-based method to calculate pH-dependent electrostatic effects in protein molecules and to predict the pK values of individual sites of titration. In addition, a CHARMm-based algorithm is included to construct and refine the spatial coordinates of all hydrogen atoms at a given pH. The present method combines electrostatic energy calculations based on the Generalized Born approximation with an iterative mobile clustering approach to calculate the equilibria of proton binding to multiple titration sites in protein molecules. The use of the GBIM (Generalized Born with Implicit Membrane) CHARMm module makes it possible to model not only water-soluble proteins but membrane proteins as well. The method includes a novel algorithm for preliminary refinement of hydrogen coordinates. Another difference from existing approaches is that, instead of monopeptides, a set of relaxed pentapeptide structures are used as model compounds. Tests on a set of 24 proteins demonstrate the high accuracy of the method. On average, the RMSD between predicted and experimental pK values is close to 0.5 pK units on this data set, and the accuracy is achieved at very low computational cost. The pH-dependent assignment of hydrogen atoms also shows very good agreement with protonation states and hydrogen-bond network observed in neutron-diffraction structures. The method is implemented as a computational protocol in Accelrys Discovery Studio and provides a fast and easy way to study the effect of pH on many important mechanisms such as enzyme catalysis, ligand binding, protein-protein interactions, and protein stability.

Dynamics of voltage profile in enzymatic ion transporters, demonstrated in electrokinetics of proton pumping rhodopsin

H(+)-pumping rhodopsins mediate a primordial conversion of light to metabolic energy. Bacteriorhodopsin from Halobacterium salinarium is the first identified and (biochemically) best-studied H(+)-pumping rhodopsin. The electrical properties of H(+)-pumping rhodopsins, however, are known in more detail for the homolog Acetabularia rhodopsin, isolated from the eukaryotic green alga Acetabularia acetabulum. Based on data from Acetabularia rhodopsin we present a general reaction kinetic model of H(+)-pumping rhodopsins with only seven independent parameters, which fits the kinetic properties of photocurrents as functions of light, transmembrane voltage, internal and external pH, and time. The model describes fast photoisomerization of retinal with simultaneous H(+) transfer to an H(+) acceptor, reprotonation of retinal from the intracellular face via an H(+) donor, and proton release to the extracellular space via an H(+) release complex. The voltage sensitivities of the individual reaction steps and their temporal changes are treated here by a novel approach, whereby--as in an Ohmic voltage divider--the effective portions of the total transmembrane voltage decrease with the relative velocities of the individual reaction steps. This analysis quantitatively infers dynamic changes of the voltage profile and of the pK values of the H(+)-binding sites involved.

Investigating the mechanisms of photosynthetic proteins using continuum electrostatics

Computational methods based on continuum electrostatics are widely used in theoretical biochemistry to analyze the function of proteins. Continuum electrostatic methods in combination with quantum chemical and molecular mechanical methods can help to analyze even very complex biochemical systems. In this article, applications of these methods to proteins involved in photosynthesis are reviewed. After giving a short introduction to the basic concepts of the continuum electrostatic model based on the Poisson-Boltzmann equation, we describe the application of this approach to the docking of electron transfer proteins, to the comparison of isofunctional proteins, to the tuning of absorption spectra, to the analysis of the coupling of electron and proton transfer, to the analysis of the effect of membrane potentials on the energetics of membrane proteins, and to the kinetics of charge transfer reactions. Simulations as those reviewed in this article help to analyze molecular mechanisms on the basis of the structure of the protein, guide new experiments, and provide a better and deeper understanding of protein functions.

The energetics of the primary proton transfer in bacteriorhodopsin revisited: it is a sequential light-induced charge separation after all

The light-induced proton transport in bacteriorhodopsin has been considered as a model for other light-induced proton pumps. However, the exact nature of this process is still unclear. For example, it is not entirely clear what the driving force of the initial proton transfer is and, in particular, whether it reflects electrostatic forces or other effects. The present work simulates the primary proton transfer (PT) by a specialized combination of the EVB and the QCFF/PI methods. This combination allows us to obtain sufficient sampling and a quantitative free energy profile for the PT at different protein configurations. The calculated profiles provide new insight about energetics of the primary PT and its coupling to the protein conformational changes. Our finding confirms the tentative analysis of an earlier work (A. Warshel, Conversion of light energy to electrostatic energy in the proton pump of Halobacterium halobium, Photochem. Photobiol. 30 (1979) 285-290) and determines that the overall PT process is driven by the energetics of the charge separation between the Schiff base and its counterion Asp85. Apparently, the light-induced relaxation of the steric energy of the chromophore leads to an increase in the ion-pair distance, and this drives the PT process. Our use of the linear response approximation allows us to estimate the change in the protein conformational energy and provides the first computational description of the coupling between the protein structural changes and the PT process. It is also found that the PT is not driven by twist-modulated changes of the Schiff base's pKa, changes in the hydrogen bond directionality, or other non-electrostatic effects. Overall, based on a consistent use of structural information as the starting point for converging free energy calculations, we conclude that the primary event should be described as a light-induced formation of an unstable ground state, whose relaxation leads to charge separation and to the destabilization of the ion-pair state. This provides the driving force for the subsequent PT steps.

Proton pathways in a [NiFe]-hydrogenase: A theoretical study

We present here a theoretical study to investigate possible proton pathways in the [NiFe]-hydrogenase from Desulfovibrio gigas. The approach used in this study consists of a combination of Poisson-Boltzmann and Monte Carlo simulations together with a distance-based network analysis to find possible groups involved in the proton transfer. Results obtained at different pH values show a reasonable number of proton active residues distributed by the protein interior and surface, with a concentration around the metal centres. The electrostatic interactions in this protein are strong, as shown by the unusual shape of the titration curves of several sites. Some residue pairs show strongly correlated protonations, indicating the sharing and probably exchange of a proton between them. The conjugation of the PB and MC simulations with the distance-based analysis allows a detailed characterization of the possible proton pathways. We discuss previous suggestions and propose a new complete pathway for the proton transfer between the active site and the surface. This pathway is mainly composed of histidines and glutamic acid residues.

Simulation of the Electron Transfer between the Tetraheme Subunit and the Special Pair of the Photosynthetic Reaction Center Using a Microstate Description

Charge transfer through biological macromolecules is essential for many biological processes such as, for instance, photosynthesis and respiration. Protons or electrons are transferred between titratable residues or redox-active cofactors, respectively. Transfer rates between these sites depend on the current charge configuration of neighboring sites. Here, we formulate the kinetics of charge-transfer systems in a microstate formalism. A unique transfer rate constant can be assigned to the interconversion of microstates. Mutual interactions between sites participating in the transfer reactions are naturally taken into account. The formalism is applied to the kinetics of electron transfer in the tetraheme subunit and the special pair of the reaction center of Blastochloris viridis. It is shown that continuum electrostatic calculations can be used in combination with an existing empirical rate law to obtain electron-transfer rate constants. The re-reduction kinetics of the photo-oxidized special pair simulated in a microstate formalism is shown to be in good agreement with experimental data. A flux analysis is used to follow the individual electron-transfer steps.