In-silico assessment of protein-protein electron transfer. a case study: cytochrome c peroxidase--cytochrome c

Why Do Most Aromatics Fail to Support Hole Hopping in the Cytochrome c Peroxidase-Cytochrome c Complex?

Electron transport through aromatic species (especially tryptophan and tyrosine) plays a central role in water splitting, redox signaling, oxidative damage protection, and bioenergetics. The cytochrome c peroxidase (CcP)-cytochrome c (Cc) complex (CcP:Cc) is used widely to study interprotein electron transfer (ET) mechanisms. Tryptophan 191 (Trp191) of CcP supports hole hopping charge recombination in the CcP:Cc complex. Experimental studies find that when Trp191 is substituted by tyrosine, phenylalanine, or redox-active aniline derivatives bound in the W191G cavity, enzymatic activity and charge recombination rates both decrease. Theoretical analysis of these CcP:Cc complexes finds that the ET kinetics depend strongly on the chemistry of the modified Trp site. The computed electronic couplings in the W191F and W191G species are orders of magnitude smaller than in the native protein, due largely to the absence of a hopping intermediate and the large tunneling distance. Small molecules bound in the W191G cavity are weakly coupled electronically to the Cc heme, and the structural disorder of the guest molecule in the binding pocket may contribute further to the lack of enzymatic activity. The couplings in W191Y are not substantially weakened compared to the native species, but the redox potential difference for tyrosine vs tryptophan oxidation accounts for the slower rate in the Tyr mutant. Thus, theoretical analysis explains why only the native Trp supports rapid hole hopping in the CcP:Cc complex. Favorable free energies and electronic couplings are essential for establishing an efficient hole hopping relay in this protein-protein complex.

Structural-Based Modeling in Protein Engineering. A Must Do

Biotechnological solutions will be a key aspect in our immediate future society, where optimized enzymatic processes through enzyme engineering might be an important solution for waste transformation, clean energy production, biodegradable materials, and green chemistry, for example. Here we advocate the importance of structural-based bioinformatics and molecular modeling tools in such developments. We summarize our recent experiences indicating a great prediction/success ratio, and we suggest that an early in silico phase should be performed in enzyme engineering studies. Moreover, we demonstrate the potential of a new technique combining Rosetta and PELE, which could provide a faster and more automated procedure, an essential aspect for a broader use.

Characterization of a Dye-Decolorizing Peroxidase from Irpex lacteus Expressed in Escherichia coli: An Enzyme with Wide Substrate Specificity Able to Transform Lignosulfonates

A dye-decolorizing peroxidase (DyP) from Irpex lacteus was cloned and heterologously expressed as inclusion bodies in Escherichia coli. The protein was purified in one chromatographic step after its in vitro activation. It was active on ABTS, 2,6-dimethoxyphenol (DMP), and anthraquinoid and azo dyes as reported for other fungal DyPs, but it was also able to oxidize Mn2+ (as manganese peroxidases and versatile peroxidases) and veratryl alcohol (VA) (as lignin peroxidases and versatile peroxidases). This corroborated that I. lacteus DyPs are the only enzymes able to oxidize high redox potential dyes, VA and Mn+2. Phylogenetic analysis grouped this enzyme with other type D-DyPs from basidiomycetes. In addition to its interest for dye decolorization, the results of the transformation of softwood and hardwood lignosulfonates suggest a putative biological role of this enzyme in the degradation of phenolic lignin.

Tuning Radical Relay Residues by Proton Management Rescues Protein Electron Hopping

Transient tyrosine and tryptophan radicals play key roles in the electron transfer (ET) reactions of photosystem (PS) II, ribonucleotide reductase (RNR), photolyase, and many other proteins. However, Tyr and Trp are not functionally interchangeable, and the factors controlling their reactivity are often unclear. Cytochrome c peroxidase (CcP) employs a Trp191^•+ radical to oxidize reduced cytochrome c (Cc). Although a Tyr191 replacement also forms a stable radical, it does not support rapid ET from Cc. Here we probe the redox properties of CcP Y191 by non-natural amino acid substitution, altering the ET driving force and manipulating the protic environment of Y191. Higher potential fluorotyrosine residues increase ET rates marginally, but only addition of a hydrogen bond donor to Tyr191^• (via Leu232His or Glu) substantially alters activity by increasing the ET rate by nearly 30-fold. ESR and ESEEM spectroscopies, crystallography, and pH-dependent ET kinetics provide strong evidence for hydrogen bond formation to Y191^• by His232/Glu232. Rate measurements and rapid freeze quench ESR spectroscopy further reveal differences in radical propagation and Cc oxidation that support an increased Y191^• formal potential of ∼200 mV in the presence of E232. Hence, Y191 inactivity results from a potential drop owing to Y191^•+ deprotonation. Incorporation of a well-positioned base to accept and donate back a hydrogen bond upshifts the Tyr^• potential into a range where it can effectively oxidize Cc. These findings have implications for the Y_Z/Y_D radicals of PS II, hole-hopping in RNR and cryptochrome, and engineering proteins for long-range ET reactions.

Computational Design of Protein Function

The computational design of enzymes is a tremendous challenge for both chemistry and biochemistry. The ability to design stable and functional biocatalysts that could operate under different conditions to perform chemical reactions without precedent in nature, allowing the large-scale production of chemicals à la carte, would revolutionise both synthetic, pharmacologic and materials chemistry. Despite the great advances achieved, this highly multidisciplinary area of research is still in its infancy. This chapter describes the ‘inside-out’ protocol for computational enzyme design and both the achievements and limitations of the current technology are highlighted. Furthermore, molecular dynamics simulations have proved to be invaluable in the enzyme design process, constituting an important tool for discovering elusive catalytically relevant conformations of the engineered or designed enzyme. As a complement to the ‘inside-out’ design protocol, different examples where hybrid QM/MM approaches have been directly applied to discover beneficial mutations in rational computational enzyme design are described.

On the role of residue phosphorylation in 14-3-3 partners: AANAT as a case study

Twenty years ago, a novel concept in protein structural biology was discovered: the intrinsically disordered regions (IDRs). These regions remain largely unstructured under native conditions and the more are studied, more properties are attributed to them. Possibly, one of the most important is their ability to conform a new type of protein-protein interaction. Besides the classical domain-to-domain interactions, IDRs follow a 'fly-casting' model including 'induced folding'. Unfortunately, it is only possible to experimentally explore initial and final states. However, the complete movie of conformational changes of protein regions and their characterization can be addressed by in silico experiments. Here, we simulate the binding of two proteins to describe how the phosphorylation of a single residue modulates the entire process. 14-3-3 protein family is considered a master regulator of phosphorylated proteins and from a modern point-of-view, protein phosphorylation is a three component system, with writers (kinases), erasers (phosphatases) and readers. This later biological role is attributed to the 14-3-3 protein family. Our molecular dynamics results show that phosphorylation of the key residue Thr31 in a partner of 14-3-3, the aralkylamine N-acetyltransferase, releases the fly-casting mechanism during binding. On the other hand, the non-phosphorylation of the same residue traps the proteins, systematically and repeatedly driving the simulations into wrong protein-protein conformations.

Constraints on the Radical Cation Center of Cytochrome c Peroxidase for Electron Transfer from Cytochrome c

The tryptophan 191 cation radical of cytochrome c peroxidase (CcP) compound I (Cpd I) mediates long-range electron transfer (ET) to cytochrome c (Cc). Here we test the effects of chemical substitution at position 191. CcP W191Y forms a stable tyrosyl radical upon reaction with peroxide and produces spectral properties similar to those of Cpd I but has low reactivity toward reduced Cc. CcP W191G and W191F variants also have low activity, as do redox ligands that bind within the W191G cavity. Crystal structures of complexes between Cc and CcP W191X (X = Y, F, or G), as well as W191G with four bound ligands reveal similar 1:1 association modes and heme pocket conformations. The ligands display structural disorder in the pocket and do not hydrogen bond to Asp235, as does Trp191. Well-ordered Tyr191 directs its hydroxyl group toward the porphyrin ring, with no basic residue in the range of interaction. CcP W191X (X = Y, F, or G) variants substituted with zinc-porphyrin (ZnP) undergo photoinduced ET with Cc(III). Their slow charge recombination kinetics that result from loss of the radical center allow resolution of difference spectra for the charge-separated state [ZnP(+), Cc(II)]. The change from a phenyl moiety at position 191 in W191F to a water-filled cavity in W191G produces effects on ET rates much weaker than the effects of the change from Trp to Phe. Low net reactivity of W191Y toward Cc(II) derives either from the inability of ZnP(+) or the Fe-CcP ferryl to oxidize Tyr or from the low potential of the resulting neutral Tyr radical.

A theoretical multiscale treatment of protein-protein electron transfer: The ferredoxin/ferredoxin-NADP(+) reductase and flavodoxin/ferredoxin-NADP(+) reductase systems

In the photosynthetic electron transfer (ET) chain, two electrons transfer from photosystem I to the flavin-dependent ferredoxin-NADP(+) reductase (FNR) via two sequential independent ferredoxin (Fd) electron carriers. In some algae and cyanobacteria (as Anabaena), under low iron conditions, flavodoxin (Fld) replaces Fd as single electron carrier. Extensive mutational studies have characterized the protein-protein interaction in FNR/Fd and FNR/Fld complexes. Interestingly, even though Fd and Fld share the interaction site on FNR, individual residues on FNR do not participate to the same extent in the interaction with each of the protein partners, pointing to different electron transfer mechanisms. Despite of extensive mutational studies, only FNR/Fd X-ray structures from Anabaena and maize have been solved; structural data for FNR/Fld remains elusive. Here, we present a multiscale modelling approach including coarse-grained and all-atom protein-protein docking, the QM/MM e-Pathway analysis and electronic coupling calculations, allowing for a molecular and electronic comprehensive analysis of the ET process in both complexes. Our results, consistent with experimental mutational data, reveal the ET in FNR/Fd proceeding through a bridge-mediated mechanism in a dominant protein-protein complex, where transfer of the electron is facilitated by Fd loop-residues 40-49. In FNR/Fld, however, we observe a direct transfer between redox cofactors and less complex specificity than in Fd; more than one orientation in the encounter complex can be efficient in ET.

Catalytic surface radical in dye-decolorizing peroxidase: a computational, spectroscopic and site-directed mutagenesis study

Dye-decolorizing peroxidase (DyP) of Auricularia auricula-judae has been expressed in Escherichia coli as a representative of a new DyP family, and subjected to mutagenic, spectroscopic, crystallographic and computational studies. The crystal structure of DyP shows a buried haem cofactor, and surface tryptophan and tyrosine residues potentially involved in long-range electron transfer from bulky dyes. Simulations using PELE (Protein Energy Landscape Exploration) software provided several binding-energy optima for the anthraquinone-type RB19 (Reactive Blue 19) near the above aromatic residues and the haem access-channel. Subsequent QM/MM (quantum mechanics/molecular mechanics) calculations showed a higher tendency of Trp-377 than other exposed haem-neighbouring residues to harbour a catalytic protein radical, and identified the electron-transfer pathway. The existence of such a radical in H₂O₂-activated DyP was shown by low-temperature EPR, being identified as a mixed tryptophanyl/tyrosyl radical in multifrequency experiments. The signal was dominated by the Trp-377 neutral radical contribution, which disappeared in the W377S variant, and included a tyrosyl contribution assigned to Tyr-337 after analysing the W377S spectra. Kinetics of substrate oxidation by DyP suggests the existence of high- and low-turnover sites. The high-turnover site for oxidation of RB19 (k(cat) > 200 s⁻¹) and other DyP substrates was assigned to Trp-377 since it was absent from the W377S variant. The low-turnover site/s (RB19 k(cat) ~20 s⁻¹) could correspond to the haem access-channel, since activity was decreased when the haem channel was occluded by the G169L mutation. If a tyrosine residue is also involved, it will be different from Tyr-337 since all activities are largely unaffected in the Y337S variant.

Insights into Laccase Engineering from Molecular Simulations: Toward a Binding-Focused Strategy

Understanding the molecular determinants of enzyme performance is of primary importance for the rational design of ad hoc mutants. A novel approach, which combines efficient conformational sampling and quick reactivity scoring, is used here to shed light on how substrate oxidation was improved during the directed evolution experiment of a fungal laccase (from Pycnoporus cinnabarinus), an industrially relevant class of oxidoreductases. It is found that the enhanced activity of the evolved enzyme is mainly the result of substrate arrangement in the active site, with no important change in the redox potential of the T1 copper. Mutations at the active site shift the binding mode into a more buried substrate position and provide a more favorable electrostatic environment for substrate oxidation. As a consequence, engineering the binding event seems to be a viable way to in silico evolution of oxidoreductases.

Could tyrosine and tryptophan serve multiple roles in biological redox processes?

Single-step electron tunnelling reactions can transport charges over distances of 15-20 Åin proteins. Longer-range transfer requires multi-step tunnelling processes along redox chains, often referred to as hopping. Long-range hopping via oxidized radicals of tryptophan and tyrosine, which has been identified in several natural enzymes, has been demonstrated in artificial constructs of the blue copper protein azurin. Tryptophan and tyrosine serve as hopping way stations in high-potential charge transport processes. It may be no coincidence that these two residues occur with greater-than-average frequency in O(2)- and H(2)O(2)-reactive enzymes. We suggest that appropriately placed tyrosine and/or tryptophan residues prevent damage from high-potential reactive intermediates by reduction followed by transfer of the oxidizing equivalent to less harmful sites or out of the protein altogether.

Charge transfer in dynamical biosystems, or the treachery of (static) images

The image is not the thing. Just as a pipe rendered in an oil painting cannot be smoked, quantum mechanical coupling pathways rendered on LCDs do not convey electrons. The aim of this Account is to examine some of our recent discoveries regarding biological electron transfer (ET) and transport mechanisms that emerge when one moves beyond treacherous static views to dynamical frameworks.

Studies over the last two decades introduced both atomistic detail and macromolecule dynamics to the description of biological ET. The first model to move beyond the structureless square-barrier tunneling description is the Pathway model, which predicts how protein secondary motifs and folding-induced through-bond and through-space tunneling gaps influence kinetics. Explicit electronic structure theory is applied routinely now to elucidate ET mechanisms, to capture pathway interferences, and to treat redox cofactor electronic structure effects. Importantly, structural sampling of proteins provides an understanding of how dynamics may change the mechanisms of biological ET, as ET rates are exponentially sensitive to structure. Does protein motion average out tunneling pathways? Do conformational fluctuations gate biological ET? Are transient multistate resonances produced by energy gap fluctuations? These questions are becoming accessible as the static view of biological ET recedes and dynamical viewpoints take center stage.

This Account introduces ET reactions at the core of bioenergetics, summarizes our team’s progress toward arriving at an atomistic-level description, examines how thermal fluctuations influence ET, presents metrics that characterize dynamical effects on ET, and discusses applications in very long (micrometer scale) bacterial nanowires. The persistence of structural effects on the ET rates in the face of thermal fluctuations is considered. Finally, the flickering resonance (FR) view of charge transfer is presented to examine how fluctuations control low-barrier transport among multiple groups in van der Waals contact. FR produces exponential distance dependence in the absence of tunneling; the exponential character emerges from the probability of matching multiple vibronically broadened electronic energies within a tolerance defined by the rms coupling among interacting groups. FR thus produces band like coherent transport on the nanometer length scale, enabled by conformational fluctuations. Taken as a whole, the emerging context for ET in dynamical biomolecules provides a robust framework to design and interpret the inner workings of bioenergetics from the molecular to the cellular scale and beyond, with applications in biomedicine, biocatalysis, and energy science.

Basidiomycete DyPs: Genomic diversity, structural-functional aspects, reaction mechanism and environmental significance

The first enzyme with dye-decolorizing peroxidase (DyP) activity was described in 1999 from an arthroconidial culture of the fungus Bjerkandera adusta. However, the first DyP sequence had been deposited three years before, as a peroxidase gene from a culture of an unidentified fungus of the family Polyporaceae (probably Irpex lacteus). Since the first description, fewer than ten basidiomycete DyPs have been purified and characterized, but a large number of sequences are available from genomes. DyPs share a general fold and heme location with chlorite dismutases and other DyP-type related proteins (such as Escherichia coli EfeB), forming the CDE superfamily. Taking into account the lack of an evolutionary relationship with the catalase-peroxidase superfamily, the observed heme pocket similarities must be considered as a convergent type of evolution to provide similar reactivity to the enzyme cofactor. Studies on the Auricularia auricula-judae DyP showed that high-turnover oxidation of anthraquinone type and other DyP substrates occurs via long-range electron transfer from an exposed tryptophan (Trp377, conserved in most basidiomycete DyPs), whose catalytic radical was identified in the H2O2-activated enzyme. The existence of accessory oxidation sites in DyP is suggested by the residual activity observed after site-directed mutagenesis of the above tryptophan. DyP degradation of substituted anthraquinone dyes (such as Reactive Blue 5) most probably proceeds via typical one-electron peroxidase oxidations and product breakdown without a DyP-catalyzed hydrolase reaction. Although various DyPs are able to break down phenolic lignin model dimers, and basidiomycete DyPs also present marginal activity on nonphenolic dimers, a significant contribution to lignin degradation is unlikely because of the low activity on high redox-potential substrates.

Structural implications of the C-terminal tail in the catalytic and stability properties of manganese peroxidases from ligninolytic fungi

The genome of Ceriporiopsis subvermispora includes 13 manganese peroxidase (MnP) genes representative of the three subfamilies described in ligninolytic fungi, which share an Mn(2+)-oxidation site and have varying lengths of the C-terminal tail. Short, long and extralong MnPs were heterologously expressed and biochemically characterized, and the first structure of an extralong MnP was solved. Its C-terminal tail surrounds the haem-propionate access channel, contributing to Mn(2+) oxidation by the internal propionate, but prevents the oxidation of 2,2'-azino-bis(3-ethylbenzothiazoline-6-sulfonate) (ABTS), which is only oxidized by short MnPs and by shortened-tail variants from site-directed mutagenesis. The tail, which is anchored by numerous contacts, not only affects the catalytic properties of long/extralong MnPs but is also associated with their high acidic stability. Cd(2+) binds at the Mn(2+)-oxidation site and competitively inhibits oxidation of both Mn(2+) and ABTS. Moreover, mutations blocking the haem-propionate channel prevent substrate oxidation. This agrees with molecular simulations that position ABTS at an electron-transfer distance from the haem propionates of an in silico shortened-tail form, while it cannot reach this position in the extralong MnP crystal structure. Only small differences exist between the long and the extralong MnPs, which do not justify their classification as two different subfamilies, but they significantly differ from the short MnPs, with the presence/absence of the C-terminal tail extension being implicated in these differences.