Unambiguous Assignment of Reduction Potentials in Diheme Cytochromes

Excitation landscape of the CP43 photosynthetic antenna complex from multiscale simulations

Photosystem II (PSII), the principal enzyme of oxygenic photosynthesis, contains two integral light harvesting proteins (CP43 and CP47) that bind chlorophylls and carotenoids. The two intrinsic antennae play crucial roles in excitation energy transfer and photoprotection. CP43 interacts most closely with the reaction center of PSII, specifically with the branch of the reaction center (D1) that is responsible for primary charge separation and electron transfer. Deciphering the function of CP43 requires detailed atomic-level insights into the properties of the embedded pigments. To advance this goal, we employ a range of multiscale computational approaches to determine the site energies and excitonic profile of CP43 chlorophylls, using large all-atom models of a membrane-bound PSII monomer. In addition to time-dependent density functional theory (TD-DFT) used in the context of a quantum-mechanics/molecular-mechanics setup (QM/MM), we present a thorough analysis using the perturbed matrix method (PMM), which enables us to utilize information from long-timescale molecular dynamics simulations of native PSII-complexed CP43. The excited state energetics and excitonic couplings have both similarities and differences compared with previous experimental fits and theoretical calculations. Both static TD-DFT and dynamic PMM results indicate a layered distribution of site energies and reveal specific groups of chlorophylls that have shared contributions to low-energy excitations. Importantly, the contribution to the lowest energy exciton does not arise from the same chlorophylls at each system configuration, but rather changes as a function of conformational dynamics. An unexpected finding is the identification of a low-energy charge-transfer excited state within CP43 that involves a lumenal (C2) and the central (C10) chlorophyll of the complex. The results provide a refined basis for structure-based interpretation of spectroscopic observations and for further deciphering excitation energy transfer in oxygenic photosynthesis.

Evidence of a Distinctive Enantioselective Binding Mode for the Photoinduced Radical Cyclization of α-Chloroamides in Ene-Reductases

We demonstrate here through molecular simulations and mutational studies the origin of the enantioselectivity in the photoinduced radical cyclization of α-chloroacetamides catalyzed by ene-reductases, in particular the Gluconobacter oxidans ene-reductase and the Old Yellow Enzyme 1, which show opposite enantioselectivity. Our results reveal that neither the π-facial selectivity model nor a protein-induced selective stabilization of the transition states is able to explain the enantioselectivity of the radical cyclization in the studied flavoenzymes. We propose a new enantioinduction scenario according to which enantioselectivity is indeed controlled by transition-state stability; however, the relative stability of the prochiral transition states is not determined by direct interaction with the protein but is rather dependent on an inherent degree of freedom within the substrate itself. This intrinsic degree of freedom, distinct from the traditional π-facial exposure mode, can be controlled by the substrate conformational selection upon binding to the enzyme.

Diheme cytochromes: Effect of mixed-axial ligation on the electronic structure and electrochemical properties with cobalt porphyrin dimer

A series of six-coordinate diCo(III) porphyrin dimers, as synthetic analogues of diheme cytochromes, have been reported here having bis(imidazole), bis(pyridine) and mixed thiophenolate-pyridine/imidazole axial ligands. In the X-ray structures of bis(imidazole) and bis(pyridine) complexes, the axial ligands are in perpendicular orientation while they are parallelly oriented in their monomeric analog. The porphyrin rings are also highly ruffle-distorted in dimer but planar in monomer which reflect the effect of intramolecular interaction between two Co(porphyrin) units in dimers. In the X-ray structure of diCo(III) thiophenolate-pyridine mixed-ligated complex, the axial Co-S and Co-N(py) distances are 2.256(1) and 2.063(2) Å, respectively. The Co-N(py) distance of 2.063(2) Å is much longer than the distances of 1.961(3) and 1.972(3) Å observed in bis(pyridine) complex and the Co-S distance is larger than Co-N(py) in the mixed ligated complex which results in a displacement of Co by 0.15 Å towards the pyridine ligand from the mean porphyrin plane. Indeed, this is the first X-ray structure of a metalloporphyrin with mixed thiophenolate-pyridine axial ligands. The effect of mixed-axial ligation is demonstrated by a blue-shift of the Soret band in the UV-visible spectroscopy and also a positive shift of the Co(III)/Co(II) redox couple as compared to their bis(pyridine) analogue. The redox potentials are shifted to a large negative value just upon replacing the metal from iron to cobalt. The present investigation emphasizes the role of axial ligation, metal ions, and also the effect of heme-heme interaction in controlling the spectral and electrochemical properties.

Alternative Fast and Slow Primary Charge-Separation Pathways in Photosystem II

Photosystem-II (PSII) is a multi-subunit protein complex that harvests sunlight to perform oxygenic photosynthesis. Initial light-activated charge separation takes place at a reaction centre consisting of four chlorophylls and two pheophytins. Understanding the processes following light excitation remains elusive due to spectral congestion, the ultrafast nature, and multi-component behaviour of the charge-separation process. Here, using advanced computational multiscale approaches which take into account the large-scale configurational flexibility of the system, we identify two possible primary pathways to radical-pair formation that differ by three orders of magnitude in their kinetics. The fast (short-range) pathway is dominant, but the existence of an alternative slow (long-range) charge-separation pathway hints at the evolution of redundancy that may serve other purposes, adaptive or protective, related to formation of the unique oxidative species that drives water oxidation in PSII.

Evidence for a high pKa of an aspartic acid residue in the active site of CALB by a fully atomistic multiscale approach

Candida antarctica Lipase B (CALB) is a paradigm for the family of lipases. At pH 7, the optimal pH for catalysis, the protonation state of an aspartic acid of the active site (Asp134) could not be conclusively assigned. In fact, the pK_a estimate provided by a widely used computational tool, namely PropKa, that predicts pKa values of ionizable groups in proteins based on the crystallographic structure, is only slightly above 7 (pKa = 7.25). This, along with the lack of an experimental evaluation, makes the assignment of its protonation state at neutral pH challenging. Here, we calculate the pK_a of Asp134 by means of a fully atomistic multiscale computational approach based on classical molecular dynamics (MD) simulation and the perturbed matrix method (PMM), namely the MD-PMM approach. MD-PMM is able to take into account the dynamics of the system and, at the same time, to treat the deprotonation step at the quantum level. The calculations provide a pK_a value of 8.9 ± 1.1, hence suggesting that Asp134 in CALB should be protonated at neutral, and even at slightly basic, pH.

Theoretical Modeling of Redox Potentials of Biomolecules

The estimation of the redox potentials of biologically relevant systems by means of theoretical-computational approaches still represents a challenge. In fact, the size of these systems typically does not allow a full quantum-mechanical treatment needed to describe electron loss/gain in such a complex environment, where the redox process takes place. Therefore, a number of different theoretical strategies have been developed so far to make the calculation of the redox free energy feasible with current computational resources. In this review, we provide a survey of such theoretical-computational approaches used in this context, highlighting their physical principles and discussing their advantages and limitations. Several examples of these approaches applied to the estimation of the redox potentials of both proteins and nucleic acids are described and critically discussed. Finally, general considerations on the most promising strategies are reported.

Heme-Heme Interactions in Diheme Cytochromes: Effect of Mixed-Axial Ligation on the Electronic Structure and Electrochemical Properties

Diheme cytochromes, the simplest members in the multiheme family, play substantial biochemical roles in enzymatic catalysis as well as in electron transfer. A series of diiron(III) porphyrin dimers have been synthesized as active site analogues of diheme cytochromes. The complexes contain six-coordinated iron(III) having thiophenol and imidazole at the fifth and sixth coordination sites, respectively. The iron centers in the complexes have been found to be in a low-spin state, as confirmed through solid-state Mössbauer and electron paramagnetic resonance (EPR) spectroscopic investigations. Mössbauer quadrupole splitting of complexes having mixed ligands is substantially larger than that observed when both axial ligands are the same. Rhombic types of EPR spectra with narrow separation between g_x, g_y, and g_z clearly distinguish heme thiolate coordination compared to bis(imidazole)-ligated low-spin heme centers. The redox potential in diheme cytochromes has been found to be tuned by interheme interactions along with the nature of axial ligands. The effect of mixed-axial ligation within the diiron(III) porphyrin dimers is demonstrated by a positive shift in the Fe(III)/Fe(II) redox couple upon thiophenolate coordination compared to their bis(imidazole) analogues. The pK_a of the imidazole also decides the extent of the shift for the Fe(III)/Fe(II) couple, while the potential of the mixed-ligated diiron(III) porphyrin dimer is more positive compared to their monomeric analogue. A variation of around 1.1 V for the Fe(III)/Fe(II) redox potential in the diiron(III) porphyrin dimer can be achieved with the combined effect of axial ligation and a metal spin state, while such a large variation in the redox potential, compared to their monomeric analogues, is attributed to the heme-heme interactions observed in dihemes. Moreover, theoretical calculations also support the experimental shifts in the redox potential values.

Exploring Coupled Redox and pH Processes with a Force-Field-Based Approach: Applications to Five Different Systems

Coupled redox and pH-driven processes are at the core of many important biological mechanisms. As the distribution of protonation and redox states in a system is associated with the pH and redox potential of the solution, having efficient computational tools that can simulate under these conditions becomes very important. Such tools have the potential to provide information that complement and drive experiments. In previous publications we have presented the implementation of the constant pH and redox potential molecular dynamics (C(pH,E)MD) method in AMBER and we have shown how multidimensional replica exchange can be used to significantly enhance the convergence efficiency of our simulations. In the current work, after an improvement in our C(pH,E)MD approach that allows a given residue to be simultaneously pH- and redox-active, we have employed our methodologies to study five different systems of interest in the literature. We present results for capped tyrosine dipeptide, two maquette systems containing one pH- and redox-active tyrosine (α₃Y and peptide A), and two proteins that contain multiple heme groups (diheme cytochrome c from Rhodobacter sphaeroides and Desulfovibrio vulgaris Hildenborough cytochrome c₃). We show that our results can provide new insights into previous theoretical and experimental findings by using a fully force-field-based and GPU-accelerated approach, which allows the simulations to be executed with high computational performance.

Evidence of a Thermodynamic Ramp for Hole Hopping to Protect a Redox Enzyme from Oxidative Damage

Redox proteins and enzymes are at risk of irreversible oxidative damage from highly oxidizing intermediates generated in the active site in the case of unsuccessful functional reaction. Chains of tyrosine and/or tryptophan residues have been recently proposed to provide protection to the active site and the whole protein by delivering oxidizing equivalents (holes) out of the protein via a multistep hopping mechanism. In the present work we use a hybrid quantum/classical theoretical-computational methodology based on the perturbed matrix method and on molecular dynamics simulations to calculate the oxidation potential difference along a chain of tyrosine and tryptophan residues in a human redox enzyme of major importance, a superoxide dismutase, which acts as antioxidant defense. We show that the hole hopping is thermodynamically favored along such a chain and that the hopping propensity is strongly affected by the protein environment and in particular by the active site and its second coordination sphere.

Multiheme proteins: effect of heme–heme interactions

This Frontier illustrates a brief personal account on the effect of heme–heme interactions in dihemes which thereby discloses some of the evolutionary design principles involved in multiheme proteins for their diverse structures and functions.

Extending the essential dynamics analysis to investigate molecular properties: application to the redox potential of proteins

Here, a methodology is proposed to investigate the collective fluctuation modes of an arbitrary set of observables, maximally contributing to the fluctuation of another functionally relevant observable. The methodology, based on the analysis of fully classical molecular dynamics (MD) simulations, exploits the essential dynamics (ED) method, originally developed to analyse the collective motions in proteins. We apply this methodology to identify the residues that are more relevant for determining the reduction potential (E(0)) of a redox-active protein. To this aim, the fluctuation modes of the single-residue electrostatic potentials mostly contributing to the fluctuations of the total electrostatic potential (the main determinant of E(0)) are investigated for wild-type azurin and two of its mutants with a higher E(0). By comparing the results here obtained with a previous study on the same systems [Zanetti-Polzi et al., Org. Biomol. Chem., 2015, 13, 11003] we show that the proposed methodology is able to identify the key sites that determine E(0). This information can be used for a general deeper understanding of the molecular mechanisms on the basis of the redox properties of the proteins under investigation, as well as for the rational design of mutants with a higher or lower E(0). From the results of the present analysis we propose a new azurin mutant that, according to our calculations, shows a further increase of E(0).

The dynamics of the β-propeller domain in Kelch protein KLHL40 changes upon nemaline myopathy-associated mutation

The nemaline myopathy-associated E528K mutation in the KLHL40 alters the communication between the Kelch propeller blades.

A few key residues determine the high redox potential shift in azurin mutants

The wide range of variability of the reduction potential (E(0)) of blue-copper proteins has been the subject of a large number of studies in the past several years. In particular, a series of azurin mutants have been recently rationally designed tuning E(0) over a very broad range (700 mV) without significantly altering the redox-active site [Marshall et al., Nature, 2009, 462, 113]. This clearly suggests that interactions outside the primary coordination sphere are relevant to determine E(0) in cupredoxins. However, the molecular determinants of the redox potential variability are still undisclosed. Here, by means of atomistic molecular dynamics simulations and hybrid quantum/classical calculations, the mechanisms that determine the E(0) shift of two azurin mutants with high potential shifts are unravelled. The reduction potentials of native azurin and of the mutants are calculated obtaining results in good agreement with the experiments. The analysis of the simulations reveals that only a small number of residues (including non-mutated ones) are relevant in determining the experimentally observed E(0) variation via site-specific, but diverse, mechanisms. These findings open the path to the rational design of new azurin mutants with different E(0).