Molecular dissection of the pseudokinase ZED1 expands effector recognition to the tomato immune receptor ZAR1.

The highly conserved angiosperm immune receptor HOPZ-ACTIVATED RESISTANCE 1 (ZAR1) is a bacterial pathogen recognition hub that mediates resistance by guarding host kinases for modification by pathogen effectors. The pseudokinase HOPZ-ETI DEFICIENT 1 (ZED1) is the only known ZAR1-guarded protein that interacts directly with a pathogen effector, HopZ1a, from the bacterial pathogen Pseudomonas syringae, making it a promising system for rational design of effector recognition for plant immunity. Here, we conducted an in-depth molecular analysis of ZED1. We generated a library of 164 random ZED1 mutants and identified 50 mutants that could not recognize the effector HopZ1a when transiently expressed in Nicotiana benthamiana. Based on our random mutants, we generated a library of 27 point mutants and found evidence of minor functional divergence between Arabidopsis (Arabidopsis thaliana) and N. benthamiana ZAR1 orthologs. We leveraged our point mutant library to identify regions in ZED1 critical for ZAR1 and HopZ1a interactions and identified two likely ZED1-HopZ1a binding conformations. We explored ZED1 nucleotide and cation binding activity and showed that ZED1 is a catalytically dead pseudokinase, functioning solely as an allosteric regulator upon effector recognition. We used our library of ZED1 point mutants to identify the ZED1 activation loop regions as the most likely cause of interspecies ZAR1-ZED1 incompatibility. Finally, we identified a mutation that abolished ZAR1-ZED1 interspecies incompatibility while retaining the ability to mediate HopZ1a recognition, which enabled recognition of HopZ1a through tomato (Solanum lycopersicum) ZAR1. This provides an example of expanded effector recognition through a ZAR1 ortholog from a non-model species.

Structural mapping of PEAK pseudokinase interactions identifies 14-3-3 as a molecular switch for PEAK3 signaling.

PEAK pseudokinases regulate cell migration, invasion and proliferation by recruiting key signaling proteins to the cytoskeleton. Despite lacking catalytic activity, alteration in their expression level is associated with several aggressive cancers. Here, we elucidate the molecular details of key PEAK signaling interactions with the adapter proteins CrkII and Grb2 and the scaffold protein 14-3-3. Our findings rationalize why the dimerization of PEAK proteins has a crucial function in signal transduction and provide biophysical and structural data to unravel binding specificity within the PEAK interactome. We identify a conserved high affinity 14-3-3 motif on PEAK3 and demonstrate its role as a molecular switch to regulate CrkII binding and signaling via Grb2. Together, our studies provide a detailed structural snapshot of PEAK interaction networks and further elucidate how PEAK proteins, especially PEAK3, act as dynamic scaffolds that exploit adapter proteins to control signal transduction in cell growth/motility and cancer.

Production and purification of the PEAK pseudokinases for structural and functional studies.

The PEAK family of pseudokinases, which comprises PEAK1, PEAK2 and PEAK3, are newly identified scaffolds that dynamically assemble oncogenic signaling pathways known to contribute to the development of several aggressive cancers. A striking feature of this unique family of pseudokinase scaffolds is their large multi-domain structure, which allows them to achieve protein complex assemblies through their structural plasticity and functional versatility. Recent structural advances have begun to reveal the critical regulatory elements that control their function. Specifically, the dimer-dependent scaffolding activity of PEAK pseudokinases is emerging as a critical mechanism for their signaling function, in addition to their ability to hetero-associate to form higher-order regulatory networks to diversify and amplify their signaling output. Here, we present a suite of techniques that enable the efficient expression and purification of PEAK proteins for functional characterization.

Crystal structure of the putative cell-wall lipoglycan biosynthesis protein LmcA from Mycobacterium smegmatis.

The bacterial genus Mycobacterium includes important pathogens, most notably M. tuberculosis, which infects one-quarter of the entire human population, resulting in around 1.4 million deaths from tuberculosis each year. Mycobacteria, and the closely related corynebacteria, synthesize a class of abundant glycolipids, the phosphatidyl-myo-inositol mannosides (PIMs). PIMs serve as membrane anchors for hyperglycosylated species, lipomannan (LM) and lipoarabinomannan (LAM), which are surface-exposed and modulate the host immune response. Previously, in studies using the model species Corynebacterium glutamicum, NCgl2760 was identified as a novel membrane protein that is required for the synthesis of full-length LM and LAM. Here, the first crystal structure of its ortholog in Mycobacterium smegmatis, MSMEG_0317, is reported at 1.8 Šresolution. The structure revealed an elongated ß-barrel fold enclosing two distinct cavities and one α-helix extending away from the ß-barrel core, resembling a `cone with a flake' arrangement. Through xenon derivatization and structural comparison with AlphaFold2-derived predictions of the M. tuberculosis homolog Rv0227c, structural elements were identified that may undergo conformational changes to switch from `closed' to `open' conformations, allowing cavity access. An AlphaFold2-derived NCgl2760 model predicted a smaller ß-barrel core with an enclosed central cavity, suggesting that all three proteins, which were collectively termed LmcA, may have a common mechanism of ligand binding through these cavities. These findings provide new structural insights into the biosynthetic pathway for a family of surface lipoglycans with important roles in mycobacterial pathogenesis.

Structural basis for small molecule targeting of Doublecortin Like Kinase 1 with DCLK1-IN-1.

Doublecortin-like kinase 1 (DCLK1) is an understudied bi-functional kinase with a proven role in tumour growth and development. However, the presence of tissue-specific spliced DCLK1 isoforms with distinct biological functions have challenged the development of effective strategies to understand the role of DCLK1 in oncogenesis. Recently, DCLK1-IN-1 was reported as a highly selective DCLK1 inhibitor, a powerful tool to dissect DCLK1 biological functions. Here, we report the crystal structures of DCLK1 kinase domain in complex with DCLK1-IN-1 and its precursors. Combined, our data rationalises the structure-activity relationship that informed the development of DCLK1-IN-1 and provides the basis for the high selectivity of DCLK1-IN-1, with DCLK1-IN-1 inducing a drastic conformational change of the ATP binding site. We demonstrate that DCLK1-IN-1 binds DCLK1 long isoforms but does not prevent DCLK1's Microtubule-Associated Protein (MAP) function. Together, our work provides an invaluable structural platform to further the design of isoform-specific DCLK1 modulators for therapeutic intervention.

Development and application of a high-throughput screening assay for identification of small molecule inhibitors of the P. falciparum reticulocyte binding-like homologue 5 protein.

The P. falciparum parasite, responsible for the disease in humans known as malaria, must invade erythrocytes to provide an environment for self-replication and survival. For invasion to occur, the parasite must engage several ligands on the host erythrocyte surface to enable adhesion, tight junction formation and entry. Critical interactions include binding of erythrocyte binding-like ligands and reticulocyte binding-like homologues (Rhs) to the surface of the host erythrocyte. The reticulocyte binding-like homologue 5 (Rh5) is the only member of this family that is essential for invasion and it binds to the basigin host receptor. The essential nature of Rh5 makes it an important vaccine target, however to date, Rh5 has not been targeted by small molecule intervention. Here, we describe the development of a high-throughput screening assay to identify small molecules which interfere with the Rh5-basigin interaction. To validate the utility of this assay we screened a known drug library and the Medicines for Malaria Box and demonstrated the reproducibility and robustness of the assay for high-throughput screening purposes. The screen of the known drug library identified the known leukotriene antagonist, pranlukast. We used pranlukast as a model inhibitor in a post screening evaluation cascade. We procured and synthesised analogues of pranlukast to assist in the hit confirmation process and show which structural moieties of pranlukast attenuate the Rh5 - basigin interaction. Evaluation of pranlukast analogues against P. falciparum in a viability assay and a schizont rupture assay show the parasite activity was not consistent with the biochemical inhibition of Rh5, questioning the developability of pranlukast as an antimalarial. The high-throughput assay developed from this work has the capacity to screen large collections of small molecules to discover inhibitors of P. falciparum Rh5 for future development of invasion inhibitory antimalarials.

Mapping and functional analysis of heterochromatin protein 1 phosphorylation in the malaria parasite Plasmodium falciparum.

Previous studies in model eukaryotes have demonstrated that phosphorylation of heterochromatin protein 1 (HP1) is important for dynamically regulating its various functions. However, in the malaria parasite Plasmodium falciparum both the function of HP1 phosphorylation and the identity of the protein kinases targeting HP1 are still elusive. In order to functionally analyze phosphorylation of P. falciparum HP1 (PfHP1), we first mapped PfHP1 phosphorylation sites by liquid chromatography tandem mass spectrometry (LC-MS/MS) analysis of native PfHP1, which identified motifs from which potential kinases could be predicted; in particular, several phosphorylated residues were embedded in motifs rich in acidic residues, reminiscent of targets for P. falciparum casein kinase 2 (PfCK2). Secondly, we tested recombinant PfCK2 and a number of additional protein kinases for their ability to phosphorylate PfHP1 in in vitro kinase assays. These experiments validated our prediction that PfHP1 acts as a substrate for PfCK2. Furthermore, LC-MS/MS analysis showed that PfCK2 phosphorylates three clustered serine residues in an acidic motif within the central hinge region of PfHP1. To study the role of PfHP1 phosphorylation in live parasites we used CRISPR/Cas9-mediated genome editing to generate a number of conditional PfHP1 phosphomutants based on the DiCre/LoxP system. Our studies revealed that neither PfCK2-dependent phosphorylation of PfHP1, nor phosphorylation of the hinge domain in general, affect PfHP1's ability to localize to heterochromatin, and that PfHP1 phosphorylation in this region is dispensable for the proliferation of P. falciparum blood stage parasites.

Conformational switching of the pseudokinase domain promotes human MLKL tetramerization and cell death by necroptosis.

Necroptotic cell death is mediated by the most terminal known effector of the pathway, MLKL. Precisely how phosphorylation of the MLKL pseudokinase domain activation loop by the upstream kinase, RIPK3, induces unmasking of the N-terminal executioner four-helix bundle (4HB) domain of MLKL, higher-order assemblies, and permeabilization of plasma membranes remains poorly understood. Here, we reveal the existence of a basal monomeric MLKL conformer present in human cells prior to exposure to a necroptotic stimulus. Following activation, toggling within the MLKL pseudokinase domain promotes 4HB domain disengagement from the pseudokinase domain αC helix and pseudocatalytic loop, to enable formation of a necroptosis-inducing tetramer. In contrast to mouse MLKL, substitution of RIPK3 substrate sites in the human MLKL pseudokinase domain completely abrogated necroptotic signaling. Therefore, while the pseudokinase domains of mouse and human MLKL function as molecular switches to control MLKL activation, the underlying mechanism differs between species.

FSHD2- and BAMS-associated mutations confer opposing effects on SMCHD1 function.

Structural maintenance of chromosomes flexible hinge domain-containing 1 (Smchd1) plays important roles in epigenetic silencing and normal mammalian development. Recently, heterozygous mutations in SMCHD1 have been reported in two disparate disorders: facioscapulohumeral muscular dystrophy type 2 (FSHD2) and Bosma arhinia microphthalmia syndrome (BAMS). FSHD2-associated mutations lead to loss of function; however, whether BAMS is associated with loss- or gain-of-function mutations in SMCHD1 is unclear. Here, we have assessed the effect of SMCHD1 missense mutations from FSHD2 and BAMS patients on ATP hydrolysis activity and protein conformation and the effect of BAMS mutations on craniofacial development in a Xenopus model. These data demonstrated that FSHD2 mutations only result in decreased ATP hydrolysis, whereas many BAMS mutations can result in elevated ATPase activity and decreased eye size in Xenopus Interestingly, a mutation reported in both an FSHD2 patient and a BAMS patient results in increased ATPase activity and a smaller Xenopus eye size. Mutations in the extended ATPase domain increased catalytic activity, suggesting critical regulatory intramolecular interactions and the possibility of targeting this region therapeutically to boost SMCHD1's activity to counter FSHD.

Probing the folding pathway of a consensus serpin using single tryptophan mutants.

Conserpin is an engineered protein that represents the consensus of a sequence alignment of eukaryotic serpins: protease inhibitors typified by a metastable native state and a structurally well-conserved scaffold. Previously, this protein has been found to adopt a native inhibitory conformation, possess an atypical reversible folding pathway and exhibit pronounced resistance to inactivation. Here we have designed a version of conserpin, cAT, with the inhibitory specificity of α1-antitrypsin, and generated single-tryptophan variants to probe its folding pathway in more detail. cAT exhibited similar thermal stability to the parental protein, an inactivation associated with oligomerisation rather a transition to the latent conformation, and a native state with pronounced kinetic stability. The tryptophan variants reveal the unfolding intermediate ensemble to consist of an intact helix H, a distorted helix F and 'breach' region structurally similar to that of a mesophilic serpin intermediate. A combination of intrinsic fluorescence, circular dichroism, and analytical gel filtration provide insight into a highly cooperative folding pathway with concerted changes in secondary and tertiary structure, which minimises the accumulation of two directly-observed aggregation-prone intermediate species. This functional conserpin variant represents a basis for further studies of the relationship between structure and stability in the serpin superfamily.

Structure of SgK223 pseudokinase reveals novel mechanisms of homotypic and heterotypic association.

The mammalian pseudokinase SgK223, and its structurally related homologue SgK269, are oncogenic scaffolds that nucleate the assembly of specific signalling complexes and regulate tyrosine kinase signalling. Both SgK223 and SgK269 form homo- and hetero-oligomers, a mechanism that underpins a diversity of signalling outputs. However, mechanistic insights into SgK223 and SgK269 homo- and heterotypic association are lacking. Here we present the crystal structure of SgK223 pseudokinase domain and its adjacent N- and C-terminal helices. The structure reveals how the N- and C-regulatory helices engage in a novel fold to mediate the assembly of a high-affinity dimer. In addition, we identified regulatory interfaces on the pseudokinase domain required for the self-assembly of large open-ended oligomers. This study highlights the diversity in how the kinase fold mediates non-catalytic functions and provides mechanistic insights into how the assembly of these two oncogenic scaffolds is achieved in order to regulate signalling output.

Smoothing a rugged protein folding landscape by sequence-based redesign.

The rugged folding landscapes of functional proteins puts them at risk of misfolding and aggregation. Serine protease inhibitors, or serpins, are paradigms for this delicate balance between function and misfolding. Serpins exist in a metastable state that undergoes a major conformational change in order to inhibit proteases. However, conformational labiality of the native serpin fold renders them susceptible to misfolding, which underlies misfolding diseases such as α1-antitrypsin deficiency. To investigate how serpins balance function and folding, we used consensus design to create conserpin, a synthetic serpin that folds reversibly, is functional, thermostable, and polymerization resistant. Characterization of its structure, folding and dynamics suggest that consensus design has remodeled the folding landscape to reconcile competing requirements for stability and function. This approach may offer general benefits for engineering functional proteins that have risky folding landscapes, including the removal of aggregation-prone intermediates, and modifying scaffolds for use as protein therapeutics.

Homo- and Heterotypic Association Regulates Signaling by the SgK269/PEAK1 and SgK223 Pseudokinases.

SgK269/PEAK1 is a pseudokinase and scaffolding protein that plays a critical role in regulating growth factor receptor signal output and is implicated in the progression of several cancers, including those of the breast, colon, and pancreas. SgK269 is structurally related to SgK223, a human pseudokinase that also functions as a scaffold but recruits a distinct repertoire of signaling proteins compared with SgK269. Structural similarities between SgK269 and SgK223 include a predicted α-helical region (designated CH) immediately preceding the conserved C-terminal pseudokinase (PK) domain. Structure-function analyses of SgK269 in MCF-10A mammary epithelial cells demonstrated a critical role for the CH and PK regions in promoting cell migration and Stat3 activation. Characterization of the SgK269 "interactome" by mass spectrometry-based proteomics identified SgK223 as a novel binding partner, and association of SgK269 with SgK223 in cells was dependent on the presence of the CH and PK domains of both pseudokinases. Homotypic association of SgK269 and SgK223 was also demonstrated and exhibited the same structural requirements. Further analysis using pulldowns and size-exclusion chromatography underscored the critical role of the CH region in SgK269/SgK223 association. Importantly, although SgK269 bridged SgK223 to Grb2, it was unable to activate Stat3 or efficiently enhance migration in SgK223 knock-out cells generated by CRISPR/Cas9. These results reveal previously unrecognized interplay between two oncogenic scaffolds and demonstrate a novel signaling mechanism for pseudokinases whereby homotypic and heterotypic association is used to assemble scaffolding complexes with distinct binding properties and hence qualitatively regulate signal output.

Biochemical and Structural Insights into Doublecortin-like Kinase Domain 1.

Doublecortin-like kinase 1 (DCLK1) is a serine/threonine kinase that belongs to the family of microtubule-associated proteins. Originally identified for its role in neurogenesis, DCLK1 has recently been shown to regulate biological processes outside of the CNS. DCLK1 is among the 15 most common putative driver genes for gastric cancers and is highly mutated across various other human cancers. However, our present understanding of how DCLK1 dysfunction leads to tumorigenesis is limited. Here, we provide evidence that DCLK1 kinase activity negatively regulates microtubule polymerization. We present the crystal structure of the DCLK1 kinase domain at 1.7 Å resolution, providing detailed insight into the ATP-binding site that will serve as a framework for future drug design. This structure also allowed for the mapping of cancer-causing mutations within the kinase domain, suggesting that a loss of kinase function may contribute to tumorigenesis.

Molecular Mechanism of CCAAT-Enhancer Binding Protein Recruitment by the TRIB1 Pseudokinase.

CCAAT-enhancer binding proteins (C/EBPs) are transcription factors that play a central role in the differentiation of myeloid cells and adipocytes. Tribbles pseudokinases govern levels of C/EBPs by recruiting them to the COP1 ubiquitin ligase for ubiquitination. Here, we present the first crystal structure of a Tribbles protein, which reveals a catalytically inactive TRIB1 pseudokinase domain with a unique adaptation in the αC helix. A second crystal structure and biophysical studies of TRIB1 with its C-terminal extension, which includes the COP1-binding motif, show that the C-terminal extension is sequestered at a site formed by the modified TRIB1 αC helix. In addition, we have identified and characterized the TRIB1 substrate-recognition sequence within C/EBPα, which is evolutionarily conserved in C/EBP transcription factors. Binding studies indicate that C/EBPα recruitment is weaker in the presence of the C-terminal COP1-binding motif, but the magnitude of this effect suggests that the two bind distinct rather directly overlapping binding sites.

The roles of helix I and strand 5A in the folding, function and misfolding of α1-antitrypsin.

α(1)-Antitrypsin, the archetypal member of the serpin superfamily, is a metastable protein prone to polymerization when exposed to stressors such as elevated temperature, low denaturant concentrations or through the presence of deleterious mutations which, in a physiological context, are often associated with disease. Experimental evidence suggests that α(1)-Antitrypsin can polymerize via several alternative mechanisms in vitro. In these polymerization mechanisms different parts of the molecule are proposed to undergo conformational change. Both strand 5 and helix I are proposed to adopt different conformations when forming the various polymers, and possess a number of highly conserved residues however their role in the folding and misfolding of α(1)-Antitrypsin has never been examined. We have therefore created a range of α(1)Antitypsin variants in order to explore the role of these conserved residues in serpin folding, misfolding, stability and function. Our data suggest that key residues in helix I mediate efficient folding from the folding intermediate and residues in strand 5A ensure native state stability in order to prevent misfolding. Additionally, our data indicate that helix I is involved in the inhibitory process and that both structural elements undergo differing conformational rearrangements during unfolding and misfolding. These findings suggest that the ability of α(1)-Antitrypsin to adopt different types of polymers under different denaturing conditions may be due to subtle conformational differences in the transiently populated structures adopted prior to the I and M* states.

Expression, purification and characterization of recombinant Z alpha(1)-antitrypsin--the most common cause of alpha(1)-antitrypsin deficiency.

Alpha(1)-antitrypsin (alpha(1)AT), the most abundant proteinase inhibitor circulating in the blood, protects extracellular matrix proteins of the lung against proteolytic destruction by neutrophil elastase. alpha(1)AT deficiency predisposes patients to emphysema, juvenile cirrhosis and hepatocellular carcinoma. Over 90% of clinical cases of severe alpha(1)AT deficiency are caused by the Z variant (E342K) of alpha(1)AT. The presence of the Z mutation results in misfolding and polymerization of alpha(1)AT. Due to its inherent propensity to polymerize there are no reported cases of recombinant Z alpha(1)AT production. This has created a major impediment to studying the effect of the Z mutation on alpha(1)AT. Here we report our attempts to produce recombinant Z alpha(1)AT using both Escherichia coli and Pichia pastoris as host systems. Using a range of expression vectors in E. coli we were unable to produce soluble active Z alpha(1)AT. Cytosolic expression of the Z alpha(1)AT gene in P. pastoris was successful. Monomeric and active recombinant Z alpha(1)AT was purified from the yeast cytosol using affinity chromatography and anion exchange chromatography. Biochemical analyses demonstrated that the recombinant Z alpha(1)AT has identical properties to its native counterpart purified from plasma of patients homozygous for the Z allele. A recombinant source of pathological Z alpha(1)AT will increase the chances of elucidating the mechanism of its polymerization and thus the development of therapeutic strategies.

A family of E. coli expression vectors for laboratory scale and high throughput soluble protein production.

BACKGROUND: In the past few years, both automated and manual high-throughput protein expression and purification has become an accessible means to rapidly screen and produce soluble proteins for structural and functional studies. However, many of the commercial vectors encoding different solubility tags require different cloning and purification steps for each vector, considerably slowing down expression screening. We have developed a set of E. coli expression vectors with different solubility tags that allow for parallel cloning from a single PCR product and can be purified using the same protocol. RESULTS: The set of E. coli expression vectors, encode for either a hexa-histidine tag or the three most commonly used solubility tags (GST, MBP, NusA) and all with an N-terminal hexa-histidine sequence. The result is two-fold: the His-tag facilitates purification by immobilised metal affinity chromatography, whilst the fusion domains act primarily as solubility aids during expression, in addition to providing an optional purification step. We have also incorporated a TEV recognition sequence following the solubility tag domain, which allows for highly specific cleavage (using TEV protease) of the fusion protein to yield native protein. These vectors are also designed for ligation-independent cloning and they possess a high-level expressing T7 promoter, which is suitable for auto-induction. To validate our vector system, we have cloned four different genes and also one gene into all four vectors and used small-scale expression and purification techniques. We demonstrate that the vectors are capable of high levels of expression and that efficient screening of new proteins can be readily achieved at the laboratory level. CONCLUSION: The result is a set of four rationally designed vectors, which can be used for streamlined cloning, expression and purification of target proteins in the laboratory and have the potential for being adaptable to a high-throughput screening.

Different conformational changes within the F-helix occur during serpin folding, polymerization, and proteinase inhibition.

The intrinsic metastability of the serpin native state is the thermodynamic driving force for both proteinase inhibition and the formation of inactive polymers. A number of mechanisms has been proposed to explain how both these conformational changes are achieved. However, one aspect that has received little attention is the movement of the F-helix, which physically impedes both these events. We have applied a protein engineering approach to investigate the conformational changes of this helix during proteinase inhibition, serpin folding, and polymerization. We systematically mutated two highly conserved hydrophobic residues on the F-helix, V161 and I157, and in addition, removed a hydrogen bond between D149 and the first turn of the helix. Our data demonstrate that while all three interactions are important for the stability and folding of the molecule, their contribution during inhibition and polymerization differ. The presence of I157 is crucial to all conformational changes as its loss results in inactivation of the serpin and rapid polymerization. The replacement of D149 does not affect activity but significantly increases the polymerization rate. The interactions formed by V161 play an important role only in maintaining the native conformation. Taken together, these data suggest that the F-helix undergoes a reversible conformational change in both its N- and C-termini during proteinase inhibition only the C-terminus undergoes changes during polymerization, but there is a global change required for folding.

X-ray structure of methanol dehydrogenase from Paracoccus denitrificans and molecular modeling of its interactions with cytochrome c-551i.

The X-ray structure of methanol dehydrogenase (MEDH) from Paracoccus denitrificans (MEDH-PD) was determined at 2.5 A resolution using molecular replacement based on the structure of MEDH from Methylophilus methylotrophus W3A1 (MEDH-WA). The overall structures from the two bacteria are similar to each other except that the former has a longer C-terminal tail in each subunit and shows local differences in several insertion regions. The "X-ray sequence" of the segment alphaGly444-alphaLeu452 was established, including one insertion and seven replacements compared with the reported sequence. The primary electron acceptor of MEDH-PD is cytochrome c-551i (Cyt c551i). Based on the crystal structure of MEDH-PD and of the published structure of Cyt c551i, their interactions were investigated by molecular modeling. As a guide and starting point, the covalently attached cytochrome and PQQ domains of the alcohol dehydrogenase from Pseudomonas putida HK5 (ADH2B) were used. In the modeling, two molecules of Cyt c551i could be accommodated in their interaction with the MEDH heterotetramer in accordance with the two-fold molecular symmetry of the latter. Two models are proposed, in both of which electrostatic and hydrogen bonding interactions make major contributions to inter-protein binding. One of these models involves salt bridges from alphaArg99 of MEDH to the heme propionic acids of Cyt c551i and the other involves salt bridges from alphaArg426 of MEDH to Glu112 of Cyt c551i. Both involve salt bridges from alphaLys93 of MEDH to Asp75 of Cyt c551i. The size and nature of the cytochrome/quinoprotein heterodimer interfaces and calculations of electronic coupling and electron transfer rates favor one of these models over the other.