Undulating Free Energy Landscapes Buffer Redox Chains from Environmental Fluctuations

Roller-coaster or undulating free energy landscapes, with alternating high and low potential cofactors, occur frequently in biological redox chains. Yet, there is little understanding of the possible advantages created by these landscapes. We examined the tetraheme subunit associated with Blastochloris viridis reaction centers, comparing the dynamics of the native protein and of hypothetical (in silico) mutants. We computed the variation in the total number of electrons in wild type (WT) and mutant tetrahemes connected to an electron reservoir in the presence of a time-varying potential, as a model for a fluctuating redox environment. We found that roller-coaster free energy landscapes buffer the redox cofactor populations from these fluctuations. The WT roller-coaster landscape slows forward and backward electron transfer in the face of fluctuations, and may offer the advantage of sustaining the reduction of essential cofactors, such as the chlorophyll special pair in photosynthesis, even though an undulating landscape introduces thermodynamically uphill steps in multistep redox chains.

Virtual Model Compound Approach for Calculating Redox Potentials of [Fe2S2]-Cys4 Centers in Proteins - Structure Quality Matters

The midpoint potential of the [Fe2S2]-Cys4-cluster in proteins is known to vary between -200 and -450 mV. This variation is caused by the different electrostatic environment of the cluster in the respective proteins. Continuum electrostatics can quantify the impact of the protein environment on the redox potential. Thus, if the redox potential of a [Fe2S2]-Cys4-cluster model compound in aqueous solution would be known, then redox potentials in various protein complexes could be calculated. However, [Fe2S2]-Cys4-cluster models are not water-soluble, and thus, their redox potential can not be measured in aqueous solution. To overcome this problem, we introduce a method that we call Virtual Model Compound Approach (VMCA) to extrapolate the model redox potential from known redox potentials of proteins. We carefully selected high-resolution structures for our analysis and divide them into a fit set, for fitting the model redox potential, and an independent test set, to check the validity of the model redox potential. However, from our analysis, we realized that the some structures can not be used as downloaded from the PDB but had to be re-refined in order to calculate reliable redox potentials. Because of the re-refinement, we were able to significantly reduce the standard deviation of our derived model redox potential for the [Fe2S2]-Cys4-cluster from 31 mV to 10 mV. As the model redox potential, we obtained -184 mV. This model redox potential can be used to analyze the redox behavior of [Fe2S2]-Cys4-clusters in larger protein complexes.

Dissecting the Ultrafast Stepwise Bidirectional Proton Relay in a Blue-Light Photoreceptor

Proton relays through H-bond networks are essential in realizing the functionality of protein machines such as in photosynthesis and photoreceptors. It has been challenging to dissect the rates and energetics of individual proton-transfer steps during the proton relay. Here, we have designed a proton rocking blue light using a flavin (BLUF) domain with the flavin mononucleotide (FMN)-glutamic acid (E)-tryptophan (W) triad and have resolved the four individual proton-transfer steps kinetically using ultrafast spectroscopy. We have found that after the photo-induced charge separation forming FMN^·-/E-COOH/WH^·+, the proton first rapidly jumps from the bridging E-COOH to FMN^- (τ_fPT2 = 3.8 ps; KIE = 1.0), followed by a second proton transfer from WH^·+ to E-COO^- (τ_fPT1 = 336 ps; KIE = 2.6) which immediately rocks back to W^· (τ_rPT1 = 85 ps; KIE = 6.7), followed by a proton return from FMNH^· to E-COO^- (τ_rPT2 = 34 ps; KIE = 3.3) with the final charge recombination between FMN^·- and WH^·+ to close the reaction cycle. Our results revisited the Grotthuss mechanism on the ultrafast timescale using the BLUF domain as a paradigm protein.

Nonlinear Migration Dynamics of Excess Electrons along Linear Oligopeptides Controlled by an Applied Electric Field

Migration of an excess electron along linear oligopeptides governed by the external electric field (E_ex ) which is against the inner dipole electric field is theoretically investigated, including the effects of E_ex on the structural and electronic properties of electron migration. Two structural properties including electron-binding ability and the dipole moment of linear oligopeptides are sensitive to the E_ex values and can be largely modulated by E_ex due to the competition of E_ex and the inner electric field and electron transfer caused by E_ex . In the case of low E_ex values, two structural properties decrease slightly_, while for high E_ex values, the electron-binding ability continually increases strongly, with dipole moments firstly increasing significantly and then increasing more slowly at higher E_ex . Additionally, linear oligopeptides of different chain lengths influence the modulation extent of E_ex and the longer the chain length is, the more sensitive modulation of E_ex is. In addition, electronic properties represented by electron spin densities and singly occupied molecular orbital distributions vary with E_ex intensities, leading to an unusual electron migration behavior. As E_ex increases, an excess electron transfers from the N-terminus to the C-terminus and jumps over a neighboring dipole unit of two termini to other units, respectively, instead of transferring by means of a one-by-one dipole unit hopping mechanism. These findings not only promote a deeper understanding of the connection between E_ex and structural and electronic properties of electron transfer behavior in peptides, but also provide a new insight into the modulation of electron migration along the oligopeptides.

Computational Biochemistry-Enzyme Mechanisms Explored

Understanding enzyme mechanisms is a major task to achieve in order to comprehend how living cells work. Recent advances in biomolecular research provide huge amount of data on enzyme kinetics and structure. The analysis of diverse experimental results and their combination into an overall picture is, however, often challenging. Microscopic details of the enzymatic processes are often anticipated based on several hints from macroscopic experimental data. Computational biochemistry aims at creation of a computational model of an enzyme in order to explain microscopic details of the catalytic process and reproduce or predict macroscopic experimental findings. Results of such computations are in part complementary to experimental data and provide an explanation of a biochemical process at the microscopic level. In order to evaluate the mechanism of an enzyme, a structural model is constructed which can be analyzed by several theoretical approaches. Several simulation methods can and should be combined to get a reliable picture of the process of interest. Furthermore, abstract models of biological systems can be constructed combining computational and experimental data. In this review, we discuss structural computational models of enzymatic systems. We first discuss various models to simulate enzyme catalysis. Furthermore, we review various approaches how to characterize the enzyme mechanism both qualitatively and quantitatively using different modeling approaches.

The open, the closed, and the empty: time-resolved fluorescence spectroscopy and computational analysis of RC-LH1 complexes from Rhodopseudomonas palustris

We studied the time-resolved fluorescence of isolated RC-LH1 complexes from Rhodopseudomonas palustris as a function of the photon fluence and the repetition rate of the excitation laser. Both parameters were varied systematically over 3 orders of magnitude. On the basis of a microstate description we developed a quantitative model for RC-LH1 and obtained very good agreement between experiments and elaborate simulations based on a global master equation approach. The model allows us to predict the relative population of RC-LH1 complexes with the special pair in the neutral state or in the oxidized state P(+) and those complexes that lack a reaction center.

Storage, transport, release: heme versatility in nitrite reductase electron transfer studied by molecular dynamics simulations

Using molecular dynamics simulations of the thermodynamic integration type, we study the energetics and kinetics of electron transfer through the nitrite reductase enzyme of Sulfurospirillum deleyianum, Wolinella succinogenes and Campylobacter jejuni.

Influence of protein interactions on oxidation/reduction midpoint potentials of cofactors in natural and de novo metalloproteins

As discussed throughout this special issue, oxidation and reduction reactions play critical roles in the function of many organisms. In photosynthetic organisms, the conversion of light energy drives oxidation and reduction reactions through the transfer of electrons and protons in order to create energy-rich compounds. These reactions occur in proteins such as cytochrome c, a heme-containing water-soluble protein, the bacteriochlorophyll-containing reaction center, and photosystem II where water is oxidized at the manganese cluster. A critical measure describing the ability of cofactors in proteins to participate in such reactions is the oxidation/reduction midpoint potential. In this review, the basic concepts of oxidation/reduction reactions are reviewed with a summary of the experimental approaches used to measure the midpoint potential of metal cofactors. For cofactors in proteins, the midpoint potential not only depends upon the specific chemical characteristics of cofactors but also upon interactions with the surrounding protein, such as the nature of the coordinating ligands and protein environment. These interactions can be tailored to optimize an oxidation/reduction reaction carried out by the protein. As examples, the midpoint potentials of hemes in cytochromes, bacteriochlorophylls in reaction centers, and the manganese cluster of photosystem II are discussed with an emphasis on the influence that protein interactions have on these potentials. This article is part of a Special Issue entitled: Metals in Bioenergetics and Biomimetics Systems.

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.

Cucurbit[8]uril mediated donor-acceptor ternary complexes: a model system for studying charge-transfer interactions

A supramolecular self-assembly approach is described which allows for the convenient preparation of a wide range of charge-transfer (CT) donor-acceptor complexes in aqueous solutions. When one equiv of the macrocyclic host cucurbit[8]uril (CB[8]) is added to an aqueous donor and acceptor solution, a heteroternary complex forms inside the host's cavity with a well-defined face-to-face π-π-stacking geometry of the donor and acceptor. This heteroternary, CB[8]-mediated complex offers the opportunity to study the CT phenomena at low concentrations and free from complications arising from any donor-donor and acceptor-acceptor interactions as a result of the large binding affinities and the very high selectivity over the formation of these homoternary complexes. Thus, this supramolocular self-assembly strategy is a practical donor-acceptor mix-and-match approach with synthetic advantages over much more cumbersome tethering schemes. While the characteristic UV/vis features of a few CB[8] ternary systems had been described as a CT band, we present for the first time systematic evidence for the existence of CT interactions between several donor-acceptor pairs that are mediated by the host CB[8]. Correlation of the experimentally obtained CT λ(max) to computed HOMO-LUMO energies demonstrated that the CT process in the host's cavity can be described by the Mulliken model. Furthermore, the literature claim of a "CT driving force" for the formation of CB[8] ternary complexes was scrutinized and evaluated by calorimetric (ITC) and ESI-MS measurements. The findings indicated that neither in the aqueous medium nor in the "gas-phase" is CT of energetic relevance to the Gibbs free binding energy. In contrast, electrostatic considerations combined with solvation effects are much better suited to rationalize the observed trends in binding affinities. Additionally, the CT λ(max) was found to be much more red-shifted (≥75 nm) inside the CB[8] cavity than in any polar organic solvents or water, indicating a significant stabilization of the CT excited state within the host cavity, further demonstrating the unique electrostatic, polar properties of the host cavity.

Distinguishing two groups of flavin reductases by analyzing the protonation state of an active site carboxylic acid

Flavin-containing reductases are involved in a wide variety of physiological reactions such as photosynthesis, nitric oxide synthesis, and detoxification of foreign compounds, including therapeutic drugs. Ferredoxin-NADP(H)-reductase (FNR) is the prototypical enzyme of this family. The fold of this protein is highly conserved and occurs as one domain of several multidomain enzymes such as the members of the diflavin reductase family. The enzymes of this family have emerged as fusion of a FNR and a flavodoxin. Although the active sites of these enzymes are very similar, different enzymes function in opposite directions, that is, some reduce oxidized nicotinamide adenine dinucleotide phosphate (NADP(+)) and some oxidize reduced nicotinamide adenine dinucleotide phosphate (NADPH). In this work, we analyze the protonation behavior of titratable residues of these enzymes through electrostatic calculations. We find that a highly conserved carboxylic acid in the active site shows a different titration behavior in different flavin reductases. This residue is deprotonated in flavin reductases present in plastids, but protonated in bacterial counterparts and in diflavin reductases. The protonation state of the carboxylic acid may also influence substrate binding. The physiological substrate for plastidic enzymes is NADP(+), but it is NADPH for the other mentioned reductases. In this article, we discuss the relevance of the environment of this residue for its protonation and its importance in catalysis. Our results allow to reinterpret and explain experimental data.

Engineering enzymes

Fundamental research into bioinorganic catalysis of the kind presented at this Faraday Discussion has the potential to turn inspiration drawn from impressive natural energy and chemical transformations into artificial catalyst constructions useful to mankind. Creating bio-inspired artificial constructions requires a level of understanding well beyond simple description of structures and mechanisms of natural enzymes. To be useful, such description must be augmented by a practical sense of structural and energetic engineering tolerances of the mechanism. Significant barriers to achieving an engineering understanding of enzyme mechanisms arise from natural protein complexity. In certain cases we can surmount these barriers to understanding, such as natural electron tunneling, coupling of electron tunneling to light capture and proton exchange as well as simpler bond breaking redox catalysis. Hope for similar solutions of more complex bioinorganic enzymes is indicated in several papers presented in this Discussion. Armed with an engineering understanding of mechanism, the current serious frustrations to successful creation of functional artificial proteins that are rooted in protein complexity can fall away. Here we discuss the genetic and biological roots of protein complexity and show how to dodge and minimize the effects of complexity. In the best-understood cases, artificial enzymes can be designed from scratch using the simplest of protein scaffolds.