Domain-domain Interface Research Articles

Generation of a high confidence set of domain-domain interface types to guide protein complex structure predictions by AlphaFold.

While the release of AlphaFold (AF) represented a breakthrough for the prediction of protein complex structures, its sensitivity, especially when using full length protein sequences, still remains limited. Modeling success rates might increase if AF predictions were guided by likely interacting protein fragments. This approach requires available sets of highly confident protein-protein interface types. Computational resources, such as 3did, infer interacting globular domain types from observed contacts in protein structures. Assessing the accuracy of these predicted interface types is difficult because we lack hand-curated reference sets of verified domain-domain interface (DDI) types. To improve protein complex modeling of DDIs by AF, we manually inspected 80 randomly selected DDI types from the 3did resource to generate a first reference set of DDI types. Identified cases of DDI type nonapproval (40%) primarily resulted from inaccurate Pfam domain matches, crystal contacts, and synthetic protein constructs. Using logistic regression, we predicted a subset of 2411 out of 5724 considered DDI types in 3did to be of high confidence, which we subsequently applied to 53000 human-protein interactions to predict DDIs followed by AF modeling. We obtained highly confident AF models for 604 out of 1129 predicted DDIs. Of note, for 47% of them no confident AF structural model could be obtained using full length protein sequences. Code is available at https://github.com/KatjaLuckLab/DDI_manuscript.

ProBID-Net: a deep learning model for protein-protein binding interface design.

Protein-protein interactions are pivotal in numerous biological processes. The computational design of these interactions facilitates the creation of novel binding proteins, crucial for advancing biopharmaceutical products. With the evolution of artificial intelligence (AI), protein design tools have swiftly transitioned from scoring-function-based to AI-based models. However, many AI models for protein design are constrained by assuming complete unfamiliarity with the amino acid sequence of the input protein, a feature most suited for de novo design but posing challenges in designing protein-protein interactions when the receptor sequence is known. To bridge this gap in computational protein design, we introduce ProBID-Net. Trained using natural protein-protein complex structures and protein domain-domain interface structures, ProBID-Net can discern features from known target protein structures to design specific binding proteins based on their binding sites. In independent tests, ProBID-Net achieved interface sequence recovery rates of 52.7%, 43.9%, and 37.6%, surpassing or being on par with ProteinMPNN in binding protein design. Validated using AlphaFold-Multimer, the sequences designed by ProBID-Net demonstrated a close correspondence between the design target and the predicted structure. Moreover, the model's output can predict changes in binding affinity upon mutations in protein complexes, even in scenarios where no data on such mutations were provided during training (zero-shot prediction). In summary, the ProBID-Net model is poised to significantly advance the design of protein-protein interactions.

Peptidase Activation by a Leader Peptide-Bound RiPP Recognition Element.

The RiPP precursor recognition element (RRE) is a conserved domain found in many prokaryotic ribosomally synthesized and post-translationally modified peptide (RiPP) biosynthetic gene clusters (BGCs). RREs bind with high specificity and affinity to a recognition sequence within the N-terminal leader region of RiPP precursor peptides. Lasso peptide biosynthesis involves an RRE-dependent leader peptidase, which is discretely encoded or fused to the RRE as a di-domain protein. Here we leveraged thousands of predicted BGCs to define the RRE:leader peptidase interaction through evolutionary covariance analysis. Each interacting domain contributes a three-stranded β-sheet to form a hydrophobic β-sandwich-like interface. The bioinformatics-guided predictions were experimentally confirmed using proteins from discrete and fused lasso peptide BGC architectures. Support for the domain-domain interface derived from chemical shift perturbation, paramagnetic relaxation enhancement experiments, and rapid variant activity screening using cell-free biosynthesis. Further validation of selected variants was performed with purified proteins. We developed a p-nitroanilide-based leader peptidase assay to illuminate the role of RRE domains. Our data show that RRE domains play a dual function. RRE domains deliver the precursor peptide to the leader peptidase, and the rate is saturable as expected for a substrate. RRE domains also partially compose the elusive S2 proteolytic pocket that binds the penultimate threonine of lasso leader peptides. Because the RRE domain is required to form the active site, leader peptidase activity is greatly diminished when the RRE domain is supplied at substoichiometric levels. Full proteolytic activation requires RRE engagement with the recognition sequence-containing portion of the leader peptide. Together, our observations define a new mechanism for protease activity regulation.

The Structure of Saccharomyces cerevisiae Arginyltransferase 1 (ATE1)

Eukaryotic post-translational arginylation, mediated by the family of enzymes known as the arginyltransferases (ATE1s), is an important post-translational modification that can alter protein function and even dictate cellular protein half-life. Multiple major biological pathways are linked to the fidelity of this process, including neural and cardiovascular developments, cell division, and even the stress response. Despite this significance, the structural, mechanistic, and regulatory mechanisms that govern ATE1 function remain enigmatic. To that end, we have used X-ray crystallography to solve the crystal structure of ATE1 from the model organism Saccharomyces cerevisiae ATE1 (ScATE1) in the apo form. The three-dimensional structure of ScATE1 reveals a bilobed protein containing a GCN5-related N-acetyltransferase (GNAT) fold, and this crystalline behavior is faithfully recapitulated in solution based on size-exclusion chromatography-coupled small angle X-ray scattering (SEC-SAXS) analyses and cryo-EM 2D class averaging. Structural superpositions and electrostatic analyses point to this domain and its domain-domain interface as the location of catalytic activity and tRNA binding, and these comparisons strongly suggest a mechanism for post-translational arginylation. Additionally, our structure reveals that the N-terminal domain, which we have previously shown to bind a regulatory [Fe-S] cluster, is dynamic and disordered in the absence of metal bound in this location, hinting at the regulatory influence of this region. When taken together, these insights bring us closer to answering pressing questions regarding the molecular-level mechanism of eukaryotic post-translational arginylation.

Computational Design of an Allosteric Antibody Switch by Deletion and Rescue of a Complex Structural Constellation.

Therapeutic monoclonal antibodies have transformed medicine, especially with regards to treating cancers and disorders of the immune system. More than 50 antibody-derived drugs have already reached the clinic, the majority of which target cytokines or cell-surface receptors. Unfortunately, many of these targets have pleiotropic functions: they serve multiple different roles, and often not all of these roles are disease-related. This can be problematic because antibodies act throughout the body, and systemic neutralization of such targets can lead to safety concerns. To address this, we have developed a strategy whereby an antibody’s ability to recognize its antigen is modulated by a second layer of control, relying on addition of an exogenous small molecule. In previous studies, we began to explore this idea by introducing a deactivating tryptophan-to-glycine mutation in the domain–domain interface of a single-chain variable fragment (scFv), and then restoring activity by adding back indole to fit the designed cavity. Here, we now describe a novel computational strategy for enumerating larger cavities that can be formed by simultaneously introducing multiple adjacent large-to-small mutations; we then carry out a complementary virtual screen to identify druglike compounds to match each candidate cavity. We first demonstrate the utility of this strategy in a fluorescein-binding single-chain variable fragment (scFv) and experimentally characterize a triple mutant with reduced antigen-binding (Rip-3) that can be rescued using a complementary ligand (Stitch-3). Because our design is built upon conserved residues in the antibody framework, we then show that the same mutation/ligand pair can also be used to modulate antigen-binding in an scFv build from a completely unrelated framework. This set of residues is present in many therapeutic antibodies as well, suggesting that this mutation/ligand pair may serve as a general starting point for introducing ligand-dependence into many clinically relevant antibodies.

Structure and function of a novel periplasmic chitooligosaccharide-binding protein from marine Vibrio bacteria

Periplasmic solute-binding proteins in bacteria are involved in the active transport of nutrients into the cytoplasm. In marine bacteria of the genus Vibrio, a chitooligosaccharide-binding protein (CBP) is thought to be the major solute-binding protein controlling the rate of chitin uptake in these bacteria. However, the molecular mechanism of the CBP involvement in chitin metabolism has not been elucidated. Here, we report the structure and function of a recombinant chitooligosaccharide-binding protein from Vibrio harveyi, namely VhCBP, expressed in Escherichia coli Isothermal titration calorimetry revealed that VhCBP strongly binds shorter chitooligosaccharides ((GlcNAc) n , where n = 2, 3, and 4) with affinities that are considerably greater than those for glycoside hydrolase family 18 and 19 chitinases but does not bind longer ones, including insoluble chitin polysaccharides. We also found that VhCBP comprises two domains with flexible linkers and that the domain-domain interface forms the sugar-binding cleft, which is not long extended but forms a small cavity. (GlcNAc)2 bound to this cavity, apparently triggering a closed conformation of VhCBP. Trp-363 and Trp-513, which stack against the two individual GlcNAc rings, likely make a major contribution to the high affinity of VhCBP for (GlcNAc)2 The strong chitobiose binding, followed by the conformational change of VhCBP, may facilitate its interaction with an active-transport system in the inner membrane of Vibrio species.

An Alternatively Packed Dry Molten Globule-like Intermediate in the Native State Ensemble of a Multidomain Protein.

It has been difficult to quantify the degree of side-chain conformational heterogeneity in the native (N) state ensemble of proteins and the relative energetic contributions of the side-chain packing and the hydrophobic effect in protein stability. Here, we show using multiple site-specific spectroscopic probes and tools of thermodynamics that the N state ensemble of a multidomain protein contains an equilibrium intermediate (I) whose interdomain region resembles a dry molten globule. In the I state, a tryptophan residue in the interdomain region is alternatively packed, but its secondary structure and intradomain packing are N-like. The I state also has a larger interdomain distance, but the domain-domain interface is dry and molten. Our results indicate that hydrophobic desolvation and side-chain packing are decoupled during protein folding and that interdomain packing interactions have an important energetic contribution in protein stability. Dynamic interconversion between alternatively packed N-like states could be important for multiple allosteric and ligand binding functions of this protein.

Structure of Chorismate Mutase-like Domain of DAHPS from Bacillus subtilis Complexed with Novel Inhibitor Reveals Conformational Plasticity of Active Site

3-deoxy-D-arabino-heptulosonate-7-phosphate-synthase (DAHPS) is the first enzyme of the shikimate pathway and is responsible for the synthesis of aromatic amino acids in microorganisms. This pathway is an attractive target for antimicrobial drugs. In Bacillus subtilis, the N-terminal domain of the bifunctional DAHPS enzyme belongs to an AroQ class of chorismate mutase and is functionally homologous to the downstream AroH class chorismate mutase. This is the first structure of chorismate mutase, AroQ (BsCM_2) enzyme from Bacillus subtilis in complex with citrate and chlorogenic acid at 1.9 Å and 1.8 Å resolution, respectively. This work provides the structural basis of ligand binding into the active site of AroQ class of chorismate mutase, while accompanied by the conformational flexibility of active site loop. Molecular dynamics results showed that helix H2′ undergoes uncoiling at the first turn and increases the mobility of loop L1′. The side chains of Arg45, Phe46, Arg52 and Lys76 undergo conformational changes, which may play an important role in DAHPS regulation by the formation of the domain-domain interface. Additionally, binding studies showed that the chlorogenic acid binds to BsCM_2 with a higher affinity than chorismate. These biochemical and structural findings could lead to the development of novel antimicrobial drugs.

Full-Length Anion Exchanger 1 Structure and Interactions with Ankyrin-1 Determined by Zero Length Crosslinking of Erythrocyte Membranes

Anion exchanger 1 (AE1) is a critical transporter and the primary structural scaffold for large macromolecular complexes responsible for erythrocyte membrane flexibility and integrity. We used zero-length crosslinking and mass spectrometry to probe AE1 structures and interactions in intact erythrocyte membranes. An experimentally verified full-length model of AE1 dimers was developed by combining crosslink-defined distance constraints with homology modeling. Previously unresolved cytoplasmic loops in the AE1 C-terminal domain are packed at the domain-domain interface on the cytoplasmic face of the membrane where they anchor the N-terminal domain's location and prevent it from occluding the ion channel. Crosslinks between AE1 dimers and ankyrin-1 indicate the likely topology for AE1 tetramers and suggest that ankyrin-1 wraps around AE1 tetramers, which may stabilize this oligomer state. This interaction and interactions of AE1 with other major erythrocyte membrane proteins show that protein-protein contacts are often substantially more extensive than previously reported.

Structure and Function of a C-C Bond Cleaving Oxygenase in Atypical Angucycline Biosynthesis.

C-C bond ring cleaving oxygenases represent a unique family of enzymes involved in the B ring cleavage reaction only observed in atypical angucycline biosynthesis. B ring cleavage is the key reaction leading to dramatic divergence in the final structures of atypical angucyclines. Here, we present the crystal structure of AlpJ, the first structure of this family of enzymes. AlpJ has been verified as the enzyme catalyzing C-C bond cleavage in kinamycin biosynthesis. The crystal structure of the AlpJ monomer resembles the dimeric structure of ferredoxin-like proteins. The N- and C-terminal halves of AlpJ are homologous, and both contain a putative hydrophobic substrate binding pocket in the "closed" and "open" conformations, respectively. Structural comparison of AlpJ with ActVA-Orf6 and protein-ligand docking analysis suggest that the residues including Asn60, Trp64, and Trp181 are possibly involved in substrate recognition. Site-directed mutagenesis results supported our hypothesis, as mutation of these residues led to nearly a complete loss of the activity of AlpJ. Structural analysis also revealed that AlpJ possesses an intramolecular domain-domain interface, where the residues His50 and Tyr178 form a hydrogen bond that probably stabilizes the three-dimensional structure of AlpJ. Site-directed mutagenesis showed that the two residues, His50 and Tyr178, were vital for the activity of AlpJ. Our findings shed light on the structure and catalytic mechanism of the AlpJ family of oxygenases, which presumably involves two active sites that might function in a cooperative manner.

An Unusual Cation-Binding Site and Distinct Domain-Domain Interactions Distinguish Class II Enolpyruvylshikimate-3-phosphate Synthases.

Enolpyruvylshikimate-3-phosphate synthase (EPSPS) catalyzes a critical step in the biosynthesis of a number of aromatic metabolites. An essential prokaryotic enzyme and the molecular target of the herbicide glyphosate, EPSPSs are the subject of both pharmaceutical and commercial interest. Two EPSPS classes that exhibit low sequence homology, differing substrate/glyphosate affinities, and distinct cation activation properties have previously been described. Here, we report structural studies of the monovalent cation-binding class II Coxiella burnetii EPSPS (cbEPSPS). Three cbEPSPS crystal structures reveal that the enzyme undergoes substantial conformational changes that alter the electrostatic potential of the active site. A complex with shikimate-3-phosphate, inorganic phosphate (Pi), and K(+) reveals that ligand induced domain closure produces an unusual cation-binding site bordered on three sides by the N-terminal domain, C-terminal domain, and the product Pi. A crystal structure of the class I Vibrio cholerae EPSPS (vcEPSPS) clarifies the basis of differential class I and class II cation responsiveness, showing that in class I EPSPSs a lysine side chain occupies the would-be cation-binding site. Finally, we identify distinct patterns of sequence conservation at the domain-domain interface and propose that the two EPSPS classes have evolved to differently optimize domain opening-closing dynamics.

Chemical cross-linking/mass spectrometry maps the amyloid β peptide binding region on both apolipoprotein E domains.

Apolipoprotein E (apoE) binds the amyloid β peptide (Aβ), one of the major culprits in Alzheimer's disease development. The formation of apoE:Aβ complexes is implicated in both Aβ clearance and fibrillization. However, the binding interface between apoE and Aβ is poorly defined despite substantial previous research efforts, and the exact role of apoE in the pathology of Alzheimer's disease remains largely elusive. Here, we compared the three main isoforms of apoE (E2, E3, and E4) for their interaction with Aβ1-42 in an early stage of aggregation and at near physiological conditions. Using electron microscopy and Western blots, we showed that all three isoforms are able to prevent Aβ fibrillization and form a noncovalent complex, with one molecule of Aβ bound per apoE. Using chemical cross-linking coupled to mass spectrometry, we further examined the interface of interaction between apoE2/3/4 and Aβ. Multiple high-confidence intermolecular apoE2/3/4:Aβ cross-links confirmed that Lys16 is located in the region of Aβ binding to apoE2/3/4. Further, we demonstrated that both N- and C-terminal domains of apoE2/3/4 are interacting with Aβ. The cross-linked sites were mapped onto and evaluated in light of a recent structure of apoE. Our results support binding of the hydrophobic Aβ at the apoE domain-domain interaction interface, which would explain how apoE is able to stabilize Aβ and thereby prevent its subsequent aggregation.

Structural Insights into the Nature of the Unique Anion Binding Site within the Cardiac Ryanodine Receptor N-Terminal Region and Associated Disease Mutations

Ryanodine receptors (RyRs) are large intracellular calcium-release channels of the endo/sarcoplasmic reticulum that play a critical role in mediating excitation-contraction coupling. Mutations within the skeletal muscle isoform (RyR1) are associated with the pharmacogenetic disorder known as Malignant Hyperthermia (MH) as well as congenital myopathies. Concurrently, mutations within the cardiac (RyR2) isoform have been linked to life-threatening arrhythmias, such as catecholaminergic polymorphic ventricular tachycardia (CPVT), a potential precursor of sudden cardiac death. While recent high-resolution structures have provided some mechanistic insights as to how mutations may impact receptor function, this subject requires more work. The RyR1 amino-terminal region contains an elaborate ionic interaction network embedded between its three domains which is key to the protein's stability. Despite all the ionic interaction partners being conserved within RyR2, the replacement of two histidine residues with Tyr125 and Arg420 within RyR2 is enough to abolish the ionic network and introduce a unique central anion binding site within RyR2. Here, we present structural insights into the nature of the RyR2 anion binding site. A 2.6A crystal structure has been solved of the amino-terminal RyR2 region (residues 1-547) into which two histidines, Y125H and R420H, were introduced which abolished chloride binding but which failed to restore the ionic network of its counterpart in RyR1. The amino-terminal intrasubunit domain-domain interfaces of both RyR1 and RyR2 are each targeted by over twenty disease-associated mutations. Here, we provide the biochemical and structural characterization of a novel MH RyR1 mutation, H113Q, and a potential CPVT-causing mutant, RyR2 R420W, both of which reside at a domain-domain interface.

Examining variable domain orientations in antigen receptors gives insight into TCR-like antibody design.

The variable domains of antibodies and T-Cell receptors (TCRs) share similar structures. Both molecules act as sensors for the immune system but recognise their respective antigens in different ways. Antibodies bind to a diverse set of antigenic shapes whilst TCRs only recognise linear peptides presented by a major histocompatibility complex (MHC). The antigen specificity and affinity of both receptors is determined primarily by the sequence and structure of their complementarity determining regions (CDRs). In antibodies the binding site is also known to be affected by the relative orientation of the variable domains, VH and VL. Here, the corresponding property for TCRs, the Vβ-Vα orientation, is investigated and compared with that of antibodies. We find that TCR and antibody orientations are distinct. General antibody orientations are found to be incompatible with binding to the MHC in a canonical TCR-like mode. Finally, factors that cause the orientation of TCRs and antibodies to be different are investigated. Packing of the long Vα CDR3 in the domain-domain interface is found to be influential. In antibodies, a similar packing affect can be achieved using a bulky residue at IMGT position 50 on the VH domain. Along with IMGT VH 50, other positions are identified that may help to promote a TCR-like orientation in antibodies. These positions should provide useful considerations in the engineering of therapeutic TCR-like antibodies.

Conformational stability as a design target to control protein aggregation.

Non-native protein aggregation is a prevalent problem occurring in many biotechnological manufacturing processes and can compromise the biological activity of the target molecule or induce an undesired immune response. Additionally, some non-native aggregation mechanisms lead to amyloid fibril formation, which can be associated with debilitating diseases. For natively folded proteins, partial or complete unfolding is often required to populate aggregation-prone conformational states, and therefore one proposed strategy to mitigate aggregation is to increase the free energy for unfolding (ΔGunf) prior to aggregation. A computational design approach was tested using human γD crystallin (γD-crys) as a model multi-domain protein. Two mutational strategies were tested for their ability to reduce/increase aggregation rates by increasing/decreasing ΔGunf: stabilizing the less stable domain and stabilizing the domain-domain interface. The computational protein design algorithm, RosettaDesign, was implemented to identify point variants. The results showed that although the predicted free energies were only weakly correlated with the experimental ΔGunf values, increased/decreased aggregation rates for γD-crys correlated reasonably well with decreases/increases in experimental ΔGunf, illustrating improved conformational stability as a possible design target to mitigate aggregation. However, the results also illustrate that conformational stability is not the sole design factor controlling aggregation rates of natively folded proteins.

Computational study on the inhibitor binding mode and allosteric regulation mechanism in hepatitis C virus NS3/4A protein.

HCV NS3/4A protein is an attractive therapeutic target responsible for harboring serine protease and RNA helicase activities during the viral replication. Small molecules binding at the interface between the protease and helicase domains can stabilize the closed conformation of the protein and thus block the catalytic function of HCV NS3/4A protein via an allosteric regulation mechanism. But the detailed mechanism remains elusive. Here, we aimed to provide some insight into the inhibitor binding mode and allosteric regulation mechanism of HCV NS3/4A protein by using computational methods. Four simulation systems were investigated. They include: apo state of HCV NS3/4A protein, HCV NS3/4A protein in complex with an allosteric inhibitor and the truncated form of the above two systems. The molecular dynamics simulation results indicate HCV NS3/4A protein in complex with the allosteric inhibitor 4VA adopts a closed conformation (inactive state), while the truncated apo protein adopts an open conformation (active state). Further residue interaction network analysis suggests the communication of the domain-domain interface play an important role in the transition from closed to open conformation of HCV NS3/4A protein. However, the inhibitor stabilizes the closed conformation through interaction with several key residues from both the protease and helicase domains, including His57, Asp79, Asp81, Asp168, Met485, Cys525 and Asp527, which blocks the information communication between the functional domains interface. Finally, a dynamic model about the allosteric regulation and conformational changes of HCV NS3/4A protein was proposed and could provide fundamental insights into the allosteric mechanism of HCV NS3/4A protein function regulation and design of new potent inhibitors.

The interdomain interface in bifunctional enzyme protein 3/4A (NS3/4A) regulates protease and helicase activities

Hepatitis C (HCV) protein 3/4A (NS3/4A) is a bifunctional enzyme comprising two separate domains with protease and helicase activities, which are essential for viral propagation. Both domains are stable and have enzymatic activity separately, and the relevance and implications of having protease and helicase together as a single protein remains to be explored. Altered in vitro activities of isolated domains compared with the full-length NS3/4A protein suggest the existence of interdomain communication. The molecular mechanism and extent of this communication was investigated by probing the domain-domain interface observed in HCV NS3/4A crystal structures. We found in molecular dynamics simulations that the two domains of NS3/4A are dynamically coupled through the interface. Interestingly, mutations designed to disrupt this interface did not hinder the catalytic activities of either domain. In contrast, substrate cleavage and DNA unwinding by these mutants were mostly enhanced compared with the wild-type protein. Disrupting the interface did not significantly alter RNA unwinding activity; however, the full-length protein was more efficient in RNA unwinding than the isolated protease domain, suggesting a more direct role in RNA processing independent of the interface. Our findings suggest that HCV NS3/4A adopts an "extended" catalytically active conformation, and interface formation acts as a switch to regulate activity. We propose a unifying model connecting HCV NS3/4A conformational states and protease and helicase function, where interface formation and the dynamic interplay between the two enzymatic domains of HCV NS3/4A potentially modulate the protease and helicase activities in vivo.

Mechanical Stability of the Antibody Domain CH3 Homodimer in Different Oxidation States

The CH3 homodimer at the C-terminal end of the antibody heavy chain is the key noncovalent interaction stabilizing antibody proteins. Here, we use single-molecule force spectroscopy to investigate the dissociation mechanics of CH3 as a proxy for antibody mechanical stability. We find the CH3 homodimer to be a highly stable complex, and its dissociation force of >150 pN at a loading rate of ≈5500 pN/s exceeds the stability of most protein-protein interactions studied to date. Separated CH3 monomers, on the other hand, are mechanically labile and only short-lived. Each CH3 monomer contains a conserved buried disulfide bridge, and we find that the successive reduction of one or both disulfide bridges in the dimer results in a stepwise decrease of the dissociation force. This suggests a structural role of the disulfide bridges helping to mold the high-affinity domain-domain interface, even though they are neither required for nor directly involved in dimerization. Taken together, our results set a limit on how much force a single antibody can bear and reveal the CH3 homodimer as a mechanical fastener that prevents antibody dissociation.

Protein-Peptide Complex Prediction through Fragment Interaction Patterns

The number of protein-peptide interactions in a cell isso large that experimental determination of all these complex structures would be a daunting task. Although homology modeling and refinement protocols have vastly improved the number and quality of predicted structural models, ab initio methods are still challenged by both the large number of possible docking sites and the conformational space accessible to flexible peptides. We present a method that addresses these challenges by sampling the entire accessible surface of a protein with a reduced conformational space of interacting backbone fragment pairs from unrelated structures. We demonstrate its potential by predicting ab initio the bound structure for a variety of protein-peptide complexes. In addition, we show the potential of our method forthe discovery of domain interaction sites and domain-domain docking.

Structural and mechanistic insights into the electron flow through protein for cytochrome c-tethering copper nitrite reductase

Copper-containing nitrite reductases (CuNiRs), which catalyse the reversible one-electron reduction of nitrite to nitric oxide, are members of a large family of multi-copper enzymes that require an interprotein electron transfer (ET) reaction with redox partner proteins. Here, we show that the naturally fused type of CuNiR tethering a cytochrome c (Cyt c) at the C-terminus folds as a unique trimeric domain-swapped structure and has a self-sufficient electron flow system. The C-terminal Cyt c domain is located at the surface of the type 1 copper (T1Cu) site in the N-terminal CuNiR domain from the adjacent subunit, the heme-to-Cu distance (10.6 Å) of which is comparable to the transient ET complex of normal CuNiR with Cyt c. The structural aspects for the domain-domain interface and the ET kinetics indicate that the Cyt c-CuNiR domain interaction should be highly transient. The further electrochemical analysis of the interprotein ET reaction with a cognate redox partner protein suggested that an electron is directly transferred from the partner to the T1Cu. Structural and mechanistic comparisons of Cyt c-CuNiR with another cupredoxin-tethering CuNiR highlight the behaviours of extra domains on the fusion types of CuNiRs required for ET through proteins.