Simulating biological charge transfer: Continuum dielectric theory or molecular dynamics?

Protein charge transfer far from equilibrium: a theoretical perspective

Potential differences for protein-assisted electron transfer across lipid bilayers or in bio-nano setups can amount to several 100 mV; they lie far outside the range of linear response theory. We describe these situations by Pauli-master equations that are based on Marcus theory of charge transfer between self-trapped electrons and that obey Kirchhoff's current law. In addition, we take on-site blockade effects and a full non-linear response of the local potentials into account. We present analytical and numerical current-potential curves and electron populations for multi-site model systems and biological electron transfer chains. Based on these, we provide empirical rules for electron populations and chemical potentials along the chain. The Pauli-master mean-field results are validated by kinetic Monte Carlo simulations. We briefly discuss the biochemical and evolutionary aspects of our findings.

Spectrochemistry of Firefly Bioluminescence

The chemical reactions underlying the emission of light in fireflies and other bioluminescent beetles are some of the most thoroughly studied processes by scientists worldwide. Despite these remarkable efforts, fierce academic arguments continue around even some of the most fundamental aspects of the reaction mechanism behind the beetle bioluminescence. In an attempt to reach a consensus, we made an exhaustive search of the available literature and compiled the key discoveries on the fluorescence and chemiluminescence spectrochemistry of the emitting molecule, the firefly oxyluciferin, and its chemical analogues reported over the past 50+ years. The factors that affect the light emission, including intermolecular interactions, solvent polarity, and electronic effects, were analyzed in the context of both the reaction mechanism and the different colors of light emitted by different luciferases. The collective data points toward a combined emission of multiple coexistent forms of oxyluciferin as the most probable explanation for the variation in color of the emitted light. We also highlight realistic research directions to eventually address some of the remaining questions related to firefly bioluminescence. It is our hope that this extensive compilation of data and detailed analysis will not only consolidate the existing body of knowledge on this important phenomenon but will also aid in reaching a wider consensus on some of the mechanistic details of firefly bioluminescence.

Monte Carlo simulation and thermodynamic integration applied to protein charge transfer

We introduce a combination of Monte Carlo simulation and thermodynamic integration methods to address a model problem in free energy computations, electron transfer in proteins. The feasibility of this approach is tested using the ferredoxin protein from Clostridium acidurici. The results are compared to numerical solutions of the Poisson-Boltzmann equation and data from recent molecular dynamics simulations on charge transfer in a protein complex, the NrfHA nitrite reductase of Desulfovibrio vulgaris. Despite the conceptual and computational simplicity of the Monte Carlo approach, the data agree well with those obtained by other methods. A link to experiments is established via the cytochrome subunit of the bacterial photosynthetic reaction center of Rhodopseudomonas viridis.

Long-range electron-electron interaction and charge transfer in protein complexes: a numerical approach

With application to the nitrite reductase hexameric protein complex of Desulfovibrio vulgaris, NrfH2A4, we suggest a strategy to compute the energy landscape of electron transfer in large systems of biochemical interest. For small complexes, the energy of all electronic configurations can be scanned completely on the level of a numerical solution of the Poisson-Boltzmann equation. In contrast, larger systems have to be treated using a pair approximation, which is verified here. Effective Coulomb interactions between neighbouring sites of excess electron localization may become as large as 200 meV, and they depend in a nontrivial manner on the intersite distance. We discuss the implications of strong Coulomb interactions on the thermodynamics and kinetics of charging and decharging a protein complex. Finally, we turn to the effect of embedding the system into a biomembrane.