Allosteric Coupling of Drug Binding and Intracellular Signaling in the A2A Adenosine Receptor.

Signaling across cellular membranes, the 826 human G protein-coupled receptors (GPCRs) govern a wide range of vital physiological processes, making GPCRs prominent drug targets. X-ray crystallography provided GPCR molecular architectures, which also revealed the need for additional structural dynamics data to support drug development. Here, nuclear magnetic resonance (NMR) spectroscopy with the wild-type-like A2A adenosine receptor (A2AAR) in solution provides a comprehensive characterization of signaling-related structural dynamics. All six tryptophan indole and eight glycine backbone 15N-1H NMR signals in A2AAR were individually assigned. These NMR probes provided insight into the role of Asp522.50 as an allosteric link between the orthosteric drug binding site and the intracellular signaling surface, revealing strong interactions with the toggle switch Trp 2466.48, and delineated the structural response to variable efficacy of bound drugs across A2AAR. The present data support GPCR signaling based on dynamic interactions between two semi-independent subdomains connected by an allosteric switch at Asp522.50.

Structural Studies of Amyloid Proteins at the Molecular Level.

Dozens of proteins are known to convert to the aggregated amyloid state. These include fibrils associated with systemic and neurodegenerative diseases and cancer, functional amyloid fibrils in microorganisms and animals, and many denatured proteins. Amyloid fibrils can be much more stable than other protein assemblies. In contrast to globular proteins, a single protein sequence can aggregate into several distinctly different amyloid structures, termed polymorphs, and a given polymorph can reproduce itself by seeding. Amyloid polymorphs may be the molecular basis of prion strains. Whereas the Protein Data Bank contains some 100,000 globular protein and 3,000 membrane protein structures, only a few dozen amyloid protein structures have been determined, and most of these are short segments of full amyloid-forming proteins. Regardless, these amyloid structures illuminate the architecture of the amyloid state, including its stability and its capacity for formation of polymorphs.

Structure of FUS Protein Fibrils and Its Relevance to Self-Assembly and Phase Separation of Low-Complexity Domains.

Polymerization and phase separation of proteins containing low-complexity (LC) domains are important factors in gene expression, mRNA processing and trafficking, and localization of translation. We have used solid-state nuclear magnetic resonance methods to characterize the molecular structure of self-assembling fibrils formed by the LC domain of the fused in sarcoma (FUS) RNA-binding protein. From the 214-residue LC domain of FUS (FUS-LC), a segment of only 57 residues forms the fibril core, while other segments remain dynamically disordered. Unlike pathogenic amyloid fibrils, FUS-LC fibrils lack hydrophobic interactions within the core and are not polymorphic at the molecular structural level. Phosphorylation of core-forming residues by DNA-dependent protein kinase blocks binding of soluble FUS-LC to FUS-LC hydrogels and dissolves phase-separated, liquid-like FUS-LC droplets. These studies offer a structural basis for understanding LC domain self-assembly, phase separation, and regulation by post-translational modification.

Recruitment of trimeric eIF2 by phosphatase non-catalytic subunit PPP1R15B.

Regulated protein phosphorylation controls most cellular processes. The protein phosphatase PP1 is the catalytic subunit of many holoenzymes that dephosphorylate serine/threonine residues. How these enzymes recruit their substrates is largely unknown. Here, we integrated diverse approaches to elucidate how the PP1 non-catalytic subunit PPP1R15B (R15B) captures its full trimeric eIF2 substrate. We found that the substrate-recruitment module of R15B is largely disordered with three short helical elements, H1, H2, and H3. H1 and H2 form a clamp that grasps the substrate in a region remote from the phosphorylated residue. A homozygous N423D variant, adjacent to H1, reducing substrate binding and dephosphorylation was discovered in a rare syndrome with microcephaly, developmental delay, and intellectual disability. These findings explain how R15B captures its 125 kDa substrate by binding the far end of the complex relative to the phosphosite to present it for dephosphorylation by PP1, a paradigm of broad relevance.

A universal allosteric mechanism for G protein activation.

G proteins play a central role in signal transduction and pharmacology. Signaling is initiated by cell-surface receptors, which promote guanosine triphosphate (GTP) binding and dissociation of Gα from the Gßγ subunits. Structural studies have revealed the molecular basis of subunit association with receptors, RGS proteins, and downstream effectors. In contrast, the mechanism of subunit dissociation is poorly understood. We use cell signaling assays, molecular dynamics (MD) simulations, and biochemistry and structural analyses to identify a conserved network of amino acids that dictates subunit release. In the presence of the terminal phosphate of GTP, a glycine forms a polar network with an arginine and glutamate, putting torsional strain on the subunit binding interface. This "G-R-E motif" secures GTP and, through an allosteric link, discharges the Gßγ dimer. Replacement of network residues prevents subunit dissociation regardless of agonist or GTP binding. These findings reveal the molecular basis of the final committed step of G protein activation.

The Solution Structure of FUS Bound to RNA Reveals a Bipartite Mode of RNA Recognition with Both Sequence and Shape Specificity.

Fused in sarcoma (FUS) is an RNA binding protein involved in regulating many aspects of RNA processing and linked to several neurodegenerative diseases. Transcriptomics studies indicate that FUS binds a large variety of RNA motifs, suggesting that FUS RNA binding might be quite complex. Here, we present solution structures of FUS zinc finger (ZnF) and RNA recognition motif (RRM) domains bound to RNA. These structures show a bipartite binding mode of FUS comprising of sequence-specific recognition of a NGGU motif via the ZnF and an unusual shape recognition of a stem-loop RNA via the RRM. In addition, sequence-independent interactions via the RGG repeats significantly increase binding affinity and promote destabilization of structured RNA conformation, enabling additional binding. We further show that disruption of the RRM and ZnF domains abolishes FUS function in splicing. Altogether, our results rationalize why deciphering the RNA binding mode of FUS has been so challenging.

An Insulin-Responsive Sensor in the SIRT1 Disordered Region Binds DBC1 and PACS-2 to Control Enzyme Activity.

Current models of SIRT1 enzymatic regulation primarily consider the effects of fluctuating levels of its co-substrate NAD+, which binds to the stably folded catalytic domain. By contrast, the roles of the sizeable disordered N- and C-terminal regions of SIRT1 are largely unexplored. Here we identify an insulin-responsive sensor in the SIRT1 N-terminal region (NTR), comprising an acidic cluster (AC) and a 3-helix bundle (3HB), controlling deacetylase activity. The allosteric assistor DBC1 removes a distal N-terminal shield from the 3-helix bundle, permitting PACS-2 to engage the acidic cluster and the transiently exposed helix 3 of the 3-helix bundle, disrupting its structure and inhibiting catalysis. The SIRT1 activator (STAC) SRT1720 binds and stabilizes the 3-helix bundle, protecting SIRT1 from inhibition by PACS-2. Identification of the SIRT1 insulin-responsive sensor and its engagement by the DBC1 and PACS-2 regulatory hub provides important insight into the roles of disordered regions in enzyme regulation and the mode by which STACs promote metabolic fitness.

Interaction networks within disease-associated GαS variants characterized by an integrative biophysical approach.

Activation of G proteins through nucleotide exchange initiates intracellular signaling cascades essential for life processes. Under normal conditions, nucleotide exchange is regulated by the formation of G protein-G protein-coupled receptor complexes. Single point mutations in the Gα subunit of G proteins bypass this interaction, leading to loss of function or constitutive gain of function, which is closely linked with the onset of multiple diseases. Despite the recognized significance of Gα mutations in disease pathology, structural information for most variants is lacking, potentially due to inherent protein dynamics that pose challenges for crystallography. To address this, we leveraged an integrative spectroscopic and computational approach to structurally characterize seven of the most frequently observed and clinically relevant mutations in the stimulatory Gα subunit, GαS. A previously proposed allosteric model of Gα activation linked structural changes in the nucleotide-binding pocket with functionally important changes in interactions between switch regions. We investigated this allosteric connection in GαS by integrating data from variable temperature CD spectroscopy, which measured changes in global protein structure and stability, and molecular dynamics simulations, which observed changes in interaction networks between GαS switch regions. Additionally, saturation-transfer difference NMR spectroscopy was applied to observe changes in nucleotide interactions with residues within the nucleotide binding site. These data have enabled testing of predictions regarding how mutations in GαS result in loss or gain of function and evaluation of proposed structural mechanisms. The integration of experimental and computational data allowed us to propose a more nuanced classification of mechanisms underlying GαS gain-of-function and loss-of-function mutations.

The enzymatic oxygen sensor cysteamine dioxygenase binds its protein substrates through their N-termini.

The non-heme iron-dependent dioxygenase 2-aminoethanethiol dioxygenase (ADO) has recently been identified as an enzymatic oxygen sensor that coordinates cellular changes to hypoxia by regulating the stability of proteins bearing an N-terminal cysteine (Nt-cys) through the N-degron pathway. It catalyses Nt-cys sulfinylation, which promotes O2-dependent proteasomal degradation of the target. Only a few ADO substrates have been verified, including regulators of G-protein signalling (RGS) 4 and 5, and the pro-inflammatory cytokine interleukin-32 (IL32), all of which exhibit cell and/or tissue specific expression patterns. ADO, in contrast, is ubiquitously expressed, suggesting it can regulate the stability of additional Nt-cys proteins in an O2-dependent manner. Furthermore, the role of individual chemical groups, active site metal, amino acid composition and globular structure on protein substrate association remains elusive. To help identify new targets and examine the underlying biochemistry of the system, we conducted a series of biophysical experiments to investigate the binding requirements of established ADO substrates RGS5 and IL32. We demonstrate, using surface plasmon response (SPR) and enzyme assays, that a free, unmodified Nt-thiol and Nt-amine are vital for substrate engagement through active site metal coordination, with residues next to Nt-cys moderately impacting association and catalytic efficiency. Additionally, we show, through 1H-15N heteronuclear single quantum coherence (15N-HSQC) nuclear magnetic resonance (NMR) titrations, that the globular portion of RGS5 has limited impact on ADO association, with interactions restricted to the N-terminus. This work establishes key features involved in ADO substrate binding, which will help identify new protein targets and, subsequently, elucidate its role in hypoxic adaptation.

Identification of small-molecule ligand-binding sites on and in the ARNT PAS-B domain.

Transcription factors are challenging to target with small-molecule inhibitors due to their structural plasticity and lack of catalytic sites. Notable exceptions include naturally ligand-regulated transcription factors, including our prior work with the hypoxia-inducible factor (HIF)-2 transcription factor, showing that small-molecule binding within an internal pocket of the HIF-2α Per-Aryl hydrocarbon Receptor Nuclear Translocator (ARNT)-Sim (PAS)-B domain can disrupt its interactions with its dimerization partner, ARNT. Here, we explore the feasibility of targeting small molecules to the analogous ARNT PAS-B domain itself, potentially opening a promising route to modulate several ARNT-mediated signaling pathways. Using solution NMR fragment screening, we previously identified several compounds that bind ARNT PAS-B and, in certain cases, antagonize ARNT association with the transforming acidic coiled-coil containing protein 3 transcriptional coactivator. However, these ligands have only modest binding affinities, complicating characterization of their binding sites. We address this challenge by combining NMR, molecular dynamics simulations, and ensemble docking to identify ligand-binding "hotspots" on and within the ARNT PAS-B domain. Our data indicate that the two ARNT/transforming acidic coiled-coil containing protein 3 inhibitors, KG-548 and KG-655, bind to a ß-sheet surface implicated in both HIF-2 dimerization and coactivator recruitment. Furthermore, while KG-548 binds exclusively to the ß-sheet surface, KG-655 can additionally bind within a water-accessible internal cavity in ARNT PAS-B. Finally, KG-279, while not a coactivator inhibitor, exemplifies ligands that preferentially bind only to the internal cavity. All three ligands promoted ARNT PAS-B homodimerization, albeit to varying degrees. Taken together, our findings provide a comprehensive overview of ARNT PAS-B ligand-binding sites and may guide the development of more potent coactivator inhibitors for cellular and functional studies.

Structural basis of binding the unique N-terminal domain of microtubule-associated protein 2c to proteins regulating kinases of signaling pathways.

Isoforms of microtubule-associated protein 2 (MAP2) differ from their homolog Tau in the sequence and interactions of the N-terminal region. Binding of the N-terminal region of MAP2c (N-MAP2c) to the dimerization/docking domains of the regulatory subunit RIIα of cAMP-dependent protein kinase (RIIDD2) and to the Src-homology domain 2 (SH2) of growth factor receptor-bound protein 2 (Grb2) have been described long time ago. However, the structural features of the complexes remained unknown due to the disordered nature of MAP2. Here, we provide structural description of the complexes. We have solved solution structure of N-MAP2c in complex with RIIDD2, confirming formation of an amphiphilic α-helix of MAP2c upon binding, defining orientation of the α-helix in the complex and showing that its binding register differs from previous predictions. Using chemical shift mapping, we characterized the binding interface of SH2-Grb2 and rat MAP2c phosphorylated by the tyrosine kinase Fyn in their complex and proposed a model explaining differences between SH2-Grb2 complexes with rat MAP2c and phosphopeptides with a Grb2-specific sequence. The results provide the structural basis of a potential role of MAP2 in regulating cAMP-dependent phosphorylation cascade via interactions with RIIDD2 and Ras signaling pathway via interactions with SH2-Grb2.

A selection and optimization strategy for single-domain antibodies targeting the PHF6 linear peptide within the tau intrinsically disordered protein.

The use of variable domain of the heavy-chain of the heavy-chain-only antibodies (VHHs) as disease-modifying biomolecules in neurodegenerative disorders holds promises, including targeting of aggregation-sensitive proteins. Exploitation of their clinical values depends however on the capacity to deliver VHHs with optimal physico-chemical properties for their specific context of use. We described previously a VHH with high therapeutic potential in a family of neurodegenerative diseases called tauopathies. The activity of this promising parent VHH named Z70 relies on its binding within the central region of the tau protein. Accordingly, we carried out random mutagenesis followed by yeast two-hybrid screening to obtain optimized variants. The VHHs selected from this initial screen targeted the same epitope as VHH Z70 as shown using NMR spectroscopy and had indeed improved binding affinities according to dissociation constant values obtained by surface plasmon resonance spectroscopy. The improved affinities can be partially rationalized based on three-dimensional structures and NMR data of three complexes consisting of an optimized VHH and a peptide containing the tau epitope. Interestingly, the ability of the VHH variants to inhibit tau aggregation and seeding could not be predicted from their affinity alone. We indeed showed that the in vitro and in cellulo VHH stabilities are other limiting key factors to their efficacy. Our results demonstrate that only a complete pipeline of experiments, here described, permits a rational selection of optimized VHH variants, resulting in the selection of VHH variants with higher affinities and/or acting against tau seeding in cell models.

Arid5a uses disordered extensions of its core ARID domain for distinct DNA- and RNA-recognition and gene regulation.

AT-rich interacting domain (ARID)-containing proteins, Arids, are a heterogeneous DNA-binding protein family involved in transcription regulation and chromatin processing. For the member Arid5a, no exact DNA-binding preference has been experimentally defined so far. Additionally, the protein binds to mRNA motifs for transcript stabilization, supposedly through the DNA-binding ARID domain. To date, however, no unbiased RNA motif definition and clear dissection of nucleic acid-binding through the ARID domain have been undertaken. Using NMR-centered biochemistry, we here define the Arid5a DNA preference. Further, high-throughput in vitro binding reveals a consensus RNA-binding motif engaged by the core ARID domain. Finally, transcriptome-wide binding (iCLIP2) reveals that Arid5a has a weak preference for (A)U-rich regions in pre-mRNA transcripts of factors related to RNA processing. We find that the intrinsically disordered regions flanking the ARID domain modulate the specificity and affinity of DNA binding, while they appear crucial for RNA interactions. Ultimately, our data suggest that Arid5a uses its extended ARID domain for bifunctional gene regulation and that the involvement of IDR extensions is a more general feature of Arids in interacting with different nucleic acids at the chromatin-mRNA interface.

Identification of potential C1-binding sites in the immunoglobulin CL domains.

Immunoglobulin G (IgG) molecules that bind antigens on the membrane of target cells spontaneously form hexameric rings, thus recruiting C1 to initiate the complement pathway. However, our previous report indicated that a mouse IgG mutant lacking the Cγ1 domain activates the pathway independently of antigen presence through its monomeric interaction with C1q via the CL domain, as well as Fc. In this study, we investigated the potential interaction between C1q and human CL isoforms. Quantitative single-molecule observations using high-speed atomic force microscopy revealed that human Cκ exhibited comparable C1q binding capabilities with its mouse counterpart, surpassing the Cλ types, which have a higher isoelectric point than the Cκ domains. Nuclear magnetic resonance and mutation experiments indicated that the human and mouse Cκ domains share a common primary binding site for C1q, centred on Glu194, a residue conserved in the Cκ domains but absent in the Cλ domains. Additionally, the Cγ1 domain, with its high isoelectric point, can cause electrostatic repulsion to the C1q head and impede the C1q-interaction adjustability of the Cκ domain in Fab. The removal of the Cγ1 domain is considered to eliminate these factors and thus promote Cκ interaction with C1q with the potential risk of uncontrolled activation of the complement pathway in vivo in the absence of antigen. However, this research underscores the presence of potential subsites in Fab for C1q binding, offering promising targets for antibody engineering to refine therapeutic antibody design.

Raf Kinase Inhibitory Protein regulates the cAMP-dependent protein kinase signaling pathway through a positive feedback loop.

Raf Kinase Inhibitory Protein (RKIP) maintains cellular robustness and prevents the progression of diseases such as cancer and heart disease by regulating key kinase cascades including MAP kinase and protein kinase A (PKA). Phosphorylation of RKIP at S153 by Protein Kinase C (PKC) triggers a switch from inhibition of Raf to inhibition of the G protein coupled receptor kinase 2 (GRK2), enhancing signaling by the ß-adrenergic receptor (ß-AR) that activates PKA. Here we report that PKA-phosphorylated RKIP promotes ß-AR-activated PKA signaling. Using biochemical, genetic, and biophysical approaches, we show that PKA phosphorylates RKIP at S51, increasing S153 phosphorylation by PKC and thereby triggering feedback activation of PKA. The S51V mutation blocks the ability of RKIP to activate PKA in prostate cancer cells and to induce contraction in primary cardiac myocytes in response to the ß-AR activator isoproterenol, illustrating the functional importance of this positive feedback circuit. As previously shown for other kinases, phosphorylation of RKIP at S51 by PKA is enhanced upon RKIP destabilization by the P74L mutation. These results suggest that PKA phosphorylation at S51 may lead to allosteric changes associated with a higher-energy RKIP state that potentiates phosphorylation of RKIP at other key sites. This allosteric regulatory mechanism may have therapeutic potential for regulating PKA signaling in disease states.

Precision engineering for localization, validation, and modification of allergenic epitopes.

The allergen-IgE interaction is essential for the genesis of allergic responses, yet investigation of the molecular basis of these interactions is in its infancy. Precision engineering has unveiled the molecular features of allergen-antibody interactions at the atomic level. High-resolution technologies, including x-ray crystallography, nuclear magnetic resonance spectroscopy, and cryo-electron microscopy, determine allergen-antibody structures. X-ray crystallography of an allergen-antibody complex localizes in detail amino acid residues and interactions that define the epitope-paratope interface. Multiple structures involving murine IgG mAbs have recently been resolved. The number of amino acids forming the epitope broadly correlates with the epitope area. The production of human IgE mAbs from B cells of allergic subjects is an exciting recent development that has for the first time enabled an actual IgE epitope to be defined. The biologic activity of defined IgE epitopes can be validated in vivo in animal models or by measuring mediator release from engineered basophilic cell lines. Finally, gene-editing approaches using the Clustered Regularly Interspaced Short Palindromic Repeats technology to either remove allergen genes or make targeted epitope engineering at the source are on the horizon. This review presents an overview of the identification and validation of allergenic epitopes by precision engineering.

Four-Order Power Reduction in Nanoscale Electron-Nuclear Double Resonance with a Nitrogen-Vacancy Center in Diamonds.

Detecting nuclear spins using single nitrogen-vacancy (NV) centers is of particular importance in nanoscale science and engineering but often suffers from the heating effect of microwave fields for spin manipulation, especially under high magnetic fields. Here, we realize an energy-efficient nanoscale nuclear-spin detection using a phase-modulation electron-nuclear double resonance scheme. The microwave field can be reduced to 1/250 of the previous requirements, and the corresponding power is over four orders lower. Meanwhile, the microwave-induced broadening to the line-width of the spectroscopy is significantly canceled, and we achieve a nuclear-spin spectrum with a resolution down to 2.1 kHz under a magnetic field at 1840 Gs. The spectral resolution can be further improved by upgrading the experimental control precision. This scheme can also be used in sensing microwave fields and can be extended to a wide range of applications in the future.

Nonsevere Burn Induces a Prolonged Systemic Metabolic Phenotype Indicative of a Persistent Inflammatory Response Postinjury.

Globally, burns are a significant cause of injury that can cause substantial acute trauma as well as lead to increased incidence of chronic comorbidity and disease. To date, research has primarily focused on the systemic response to severe injury, with little in the literature reported on the impact of nonsevere injuries (<15% total burn surface area; TBSA). To elucidate the metabolic consequences of a nonsevere burn injury, longitudinal plasma was collected from adults (n = 35) who presented at hospital with a nonsevere burn injury at admission, and at 6 week follow up. A cross-sectional baseline sample was also collected from nonburn control participants (n = 14). Samples underwent multiplatform metabolic phenotyping using 1H nuclear magnetic resonance spectroscopy and liquid chromatography-mass spectrometry to quantify 112 lipoprotein and glycoprotein signatures and 852 lipid species from across 20 subclasses. Multivariate data modeling (orthogonal projections to latent structures-discriminate analysis; OPLS-DA) revealed alterations in lipoprotein and lipid metabolism when comparing the baseline control to hospital admission samples, with the phenotypic signature found to be sustained at follow up. Univariate (Mann-Whitney U) testing and OPLS-DA indicated specific increases in GlycB (p-value < 1.0e-4), low density lipoprotein-2 subfractions (variable importance in projection score; VIP > 6.83e-1) and monoacyglyceride (20:4) (p-value < 1.0e-4) and decreases in circulating anti-inflammatory high-density lipoprotein-4 subfractions (VIP > 7.75e-1), phosphatidylcholines, phosphatidylglycerols, phosphatidylinositols, and phosphatidylserines. The results indicate a persistent systemic metabolic phenotype that occurs even in cases of a nonsevere burn injury. The phenotype is indicative of an acute inflammatory profile that continues to be sustained postinjury, suggesting an impact on systems health beyond the site of injury. The phenotypes contained metabolic signatures consistent with chronic inflammatory states reported to have an elevated incidence postburn injury. Such phenotypic signatures may provide patient stratification opportunities, to identify individual responses to injury, personalize intervention strategies, and improve acute care, reducing the risk of chronic comorbidity.

NMR characterization of an assembling RHIM (RIP homotypic interaction motif) amyloid reveals a cryptic region for self-recognition.

The RIP homotypic interaction motif (RHIM) is an essential protein motif in inflammatory signaling and certain cell death pathways. RHIM signaling occurs following the assembly of functional amyloids, and while the structural biology of such higher-order RHIM complexes has started to emerge, the conformations and dynamics of nonassembled RHIMs remain unknown. Here, using solution NMR spectroscopy, we report the characterization of the monomeric form of the RHIM in receptor-interacting protein kinase 3 (RIPK3), a fundamental protein in human immunity. Our results establish that the RHIM of RIPK3 is an intrinsically disordered protein motif, contrary to prediction, and that exchange dynamics between free monomers and amyloid-bound RIPK3 monomers involve a 20-residue stretch outside the RHIM that is not incorporated within the structured cores of the RIPK3 assemblies determined by cryo-EM or solid-state NMR. Thus, our findings expand on the structural characterization of RHIM-containing proteins, specifically highlighting conformational dynamics involved in assembly processes.

Hemes on a string: insights on the functional mechanisms of PgcA from Geobacter sulfurreducens.

Microbial extracellular reduction of insoluble compounds requires soluble electron shuttles that diffuse in the environment, freely diffusing cytochromes, or direct contact with cellular conductive appendages that release or harvest electrons to assure a continuous balance between cellular requirements and environmental conditions. In this work, we produced and characterized the three cytochrome domains of PgcA, an extracellular triheme cytochrome that contributes to Fe(III) and Mn(IV) oxides reduction in Geobacter sulfurreducens. The three monoheme domains are structurally homologous, but their heme groups show variable axial coordination and reduction potential values. Electron transfer experiments monitored by NMR and visible spectroscopy show the variable extent to which the domains promiscuously exchange electrons while reducing different electron acceptors. The results suggest that PgcA is part of a new class of cytochromes - microbial heme-tethered redox strings - that use low-complexity protein stretches to bind metals and promote intra- and intermolecular electron transfer events through its cytochrome domains.