3D domain swapping: as domains continue to swap

Three-dimensional (3D) domain swapping creates a bond between two or more protein molecules as they exchange their identical domains. Since the term '3D domain swapping' was first used to describe the dimeric structure of diphtheria toxin, the database of domain-swapped proteins has greatly expanded. Analyses of the now about 40 structurally characterized cases of domain-swapped proteins reveal that most swapped domains are at either the N or C terminus and that the swapped domains are diverse in their primary and secondary structures. In addition to tabulating domain-swapped proteins, we describe in detail several examples of 3D domain swapping which show the swapping of more than one domain in a protein, the structural evidence for 3D domain swapping in amyloid proteins, and the flexibility of hinge loops. We also discuss the physiological relevance of 3D domain swapping and a possible mechanism for 3D domain swapping. The present state of knowledge leads us to suggest that 3D domain swapping can occur under appropriate conditions in any protein with an unconstrained terminus. As domains continue to swap, this review attempts not only a summary of the known domain-swapped proteins, but also a framework for understanding future findings of 3D domain swapping.     

Replacement of the Human Topoisomerase Linker Domain with the Plasmodial Counterpart Renders the Enzyme Camptothecin Resistant

A human/plasmodial hybrid enzyme, generated by swapping the human topoisomerase IB linker domain with the corresponding domain of the Plasmodium falciparum enzyme, has been produced and characterized. The hybrid enzyme displays a relaxation activity comparable to the human enzyme, but it is characterized by a much faster religation rate. The hybrid enzyme is also camptothecin resistant. A 3D structure of the hybrid enzyme has been built and its structural-dynamical properties have been analyzed by molecular dynamics simulation. The analysis indicates that the swapped plasmodial linker samples a conformational space much larger than the corresponding domain in the human enzyme. The large linker conformational variability is then linked to important functional properties such as an increased religation rate and a low drug reactivity, demonstrating that the linker domain has a crucial role in the modulation of the topoisomerase IB activity.     

Structural and functional modularity of proteins in the de novo purine biosynthetic pathway

It is generally accepted that naturally existing functional domains can serve as building blocks for complex protein structures, and that novel functions can arise from assembly of different combinations of these functional domains. To inform our understanding of protein evolution and explore the modular nature of protein structure, two model enzymes were chosen for study, purT‐encoded glycinamide ribonucleotide formyltransferase (PurT) and purK‐encoded N⁵‐carboxylaminoimidazole ribonucleotide synthetase (PurK). Both enzymes are found in the de novo purine biosynthetic pathway of Escherichia coli. In spite of their low sequence identity, PurT and PurK share significant similarity in terms of tertiary structure, active site organization, and reaction mechanism. Their characteristic three domain structures categorize both PurT and PurK as members of the ATP‐grasp protein superfamily. In this study, we investigate the exchangeability of individual protein domains between these two enzymes and the in vivo and in vitro functional properties of the resulting hybrids. Six domain‐swapped hybrids were unable to catalyze full wild‐type reactions, but each hybrid protein could catalyze partial reactions. Notably, an additional loop replacement in one of the domain‐swapped hybrid proteins was able to restore near wild‐type PurK activity. Therefore, in this model system, domain‐swapped proteins retained the ability to catalyze partial reactions, but further modifications were required to efficiently couple the reaction intermediates and achieve catalysis of the full reaction. Implications for understanding the role of domain swapping in protein evolution are discussed.     

Open interface and large quaternary structure movements in 3D domain swapped proteins: insights from molecular dynamics simulations of the C-terminal swapped dimer of ribonuclease A

Bovine pancreatic ribonuclease (RNase A) forms two three-dimensional (3D) domain swapped dimers. Crystallographic investigations have revealed that these dimers display completely different quaternary structures: one dimer (N-dimer), which presents the swapping of the N-terminal helix, is characterized by a compact structure, whereas the other (C-dimer), which is stabilized by the exchange of the C-terminal end, shows a rather loose assembly of the two subunits. The dynamic properties of monomeric RNase A and of the N-dimer have been extensively characterized. Here, we report a molecular dynamics investigation carried out on the C-dimer. This computational experiment indicates that the quaternary structure of the C-dimer undergoes large fluctuations. These motions do not perturb the proper folding of the two subunits, which retain the dynamic properties of RNase A and the N-dimer. Indeed, the individual subunits of the C-dimer display the breathing motion of the beta-sheet structure, which is important for the enzymatic activity of pancreatic-like ribonucleases. In contrast to what has been observed for the N-dimer, the breathing motion of the two subunits of the C-dimer is not coupled. This finding suggests that the intersubunit communications in a 3D domain swapped dimer strongly rely on the extent of the interchain interface. Furthermore, the observation that the C-dimer is endowed with a high intrinsic flexibility holds interesting implications for the specific properties of 3D domain swapped dimers. Indeed, a survey of the quaternary structures of the other 3D domain swapped dimers shows that large variations are often observed when the structural determinations are conducted in different experimental conditions. The 3D domain swapping phenomenon coupled with the high flexibility of the quaternary structure may be relevant for protein-protein recognition, and in particular for the pathological aggregations.     

Design,selection, and characterization of a split chorismate mutase

Split proteins are versatile tools for detecting protein–protein interactions and studying protein folding. Here, we report a new, particularly small split enzyme, engineered from a thermostable chorismate mutase (CM). Upon dissecting the helical‐bundle CM from Methanococcus jannaschii into a short N‐terminal helix and a 3‐helix segment and attaching an antiparallel leucine zipper dimerization domain to the individual fragments, we obtained a weakly active heterodimeric mutase. Using combinatorial mutagenesis and in vivo selection, we optimized the short linker sequences connecting the leucine zipper to the enzyme domain. One of the selected CMs was characterized in detail. It spontaneously assembles from the separately inactive fragments and exhibits wild‐type like CM activity. Owing to the availability of a well characterized selection system, the simple 4‐helix bundle topology, and the small size of the N‐terminal helix, the heterodimeric CM could be a valuable scaffold for enzyme engineering efforts and as a split sensor for specifically oriented protein–protein interactions.     

A new method for modeling large-scale rearrangements of protein domains

A method for modeling large-scale rearrangements of protein domains connected by a single- or a double-stranded linker is proposed. Multidomain proteins may undergo substantial domain displacements, while their intradomain structure remains essentially unchanged. The method allows automatic identification of an interdomain linker and builds an all-atom model of a protein structure in internal coordinates. Torsion angles belonging to the interdomain linkers and side chains potentially able to form domain interfaces are set free while all remaining torsions, bond lengths, and bond angles are fixed. Large-scale sampling of the reduced torsion conformational subspace is effected with the “biased probability Monte Carlo-minimization” method [Abagyan, R.A., Totrov, M.M. (1994): J. Mol. Biol. 235, 983–1002]. Solvation and side-chain entropic contributions are added to the energy function. A special procedure has been developed to generate concerted deformations of a double-stranded interdomain linker in such a way that the polypeptide chain continuity is preserved. The method was tested on Bence-Jones protein with a single-stranded linker and lysine/arginine/ornithine-binding (LAO) protein with a double-stranded linker. For each protein, structurally diverse low-energy conformations with ideal covalent geometry were generated, and an overlap between two sets of conformations generated starting from the crystallographically determined “closed” and “open” forms was found. One of the low-energy conformations generated in a run starting from the LAO “closed” form was only 2.2 Å away from the structure of the “open” form. The method can be useful in predicting the scope of possible domain rearrangements of a multidomain protein. Proteins 27:410–424, 1997. © 1997 Wiley-Liss, Inc.     

Effect of backbone cyclization on protein folding stability: chain entropies of both the unfolded and the folded states are restricted

Circular versions of a large number of proteins have been designed by connecting the N and C termini via peptide linkers. A motivation for these designs is the assumed enhancement in folding stability, because backbone cyclization reduces the chain entropy of the unfolded state. Here, it is recognized that backbone cyclization also reduces the chain entropy of a flexible peptide linker in the folded state. Specifically, the end-to-end distance of the linker is restricted to fluctuations around the average displacement between the N and C termini of the folded protein. The balance of the chain-entropy reductions in the folded and unfolded states is used to predict the change in the unfolding free energy, deltadeltaG(cycl), by backbone cyclization. Predicted values of deltadeltaG(cycl) are in quantitative agreement with results of a careful study on cyclizing the 34 residue PIN1 WW domain by linkers with two to seen residues. The experimental results of an optimal linker length l=4 and a maximum stabilization of 1.7 kcal/mol are reproduced. Calculations of deltadeltaG(cycl) for a broad selection of circular proteins suggest that the stabilizing effect of backbone cyclization is modest, reflecting entropy reductions in both the unfolded and the folded states.     

Comprehensive analysis of motions in molecular dynamics trajectories of the actin capping protein and its inhibitor complexes

The actin capping protein (CP) binds to actin filaments to block further elongation. The capping activity is inhibited by proteins V‐1 and CARMIL interacting with CP via steric and allosteric mechanisms, respectively. The crystal structures of free CP, CP/V‐1, and CP/CARMIL complexes suggest that the binding of CARMIL alters the flexibility of CP rather than the overall structure of CP, and this is an allosteric inhibition mechanism. Here, we performed molecular dynamics (MD) simulations of CP in the free form, and in complex with CARMIL or V‐1. The resulting trajectories were analyzed exhaustively using Motion Tree, which identifies various rigid‐body motions ranging from small local motions to large domain motions. After enumerating all the motions, CP flexibilities with different ligands were characterized by a list of frequencies for 20 dominant rigid‐body motions, some of which were not identified in previous studies. The comparative analysis highlights the influence of the binding of the CARMIL peptide to CP flexibility. In free CP and the CP/V‐1 complex, domain motions around a large crevice between the N‐stalk and the CP‐S domain occur frequently. The CARMIL peptide binds the crevice and suppresses the motions effectively. In addition, the binding of the CARMIL peptide enhances and alters local motions around the pocket that participates in V‐1 binding. These newly identified motions are likely to suppress the binding of V‐1 to CP. The observed changes in CP motion provide insights that describe the mechanism of allosteric regulation by CARMIL through modulating CP flexibility. Proteins 2016; 84:948–956. © 2016 Wiley Periodicals, Inc.     

A proposed structural model for amyloid fibril elongation: domain swapping forms an interdigitating beta-structure polymer

We propose a model illustrating how proteins, which differ in their overall sequences and structures, can form the propagating, twisted beta-sheet conformations, characteristic of amyloids. Some cases of amyloid formation can be explained through a "domain swapping" event, where the swapped segment is either a beta-hairpin or an unstable conformation which can partially unfold and assume a beta-hairpin structure. As in domain swapping, here the swapped beta-hairpin is at the edge of the structure, has few (if any) salt bridges and hydrogen bonds connecting it to the remainder of the structure and variable extents of buried non-polar surface areas. Additionally, in both cases the swapped piece constitutes a transient "building block" of the structure, with a high population time. Whereas in domain swapping the swapped fragment has been shown to be an alpha-helix, loop, strand or an entire domain, but so far not a beta-hairpin, despite the large number of cases in which it was already detected, here swapping may involve such a structural motif. We show how the swapping of beta-hairpins would form an interdigitated, twisted beta-sheet conformation, explaining the remarkable high stability of the protofibril in vitro. Such a swapping mechanism is attractive as it involves a universal mechanism in proteins, critical for their function, namely hinge-bending motions. Our proposal is consistent with structural superpositioning of mutational variants. While the overall r.m.s.d.s of the wild-type and mutants are small, the proposed hinge-bending region consistently shows larger deviations. These larger deviations illustrate that this region is more prone to respond to the mutational changes, regardless of their location in the sequence or in the structure. Nevertheless, above all, we stress that this proposition is hypothetical, since it is based on assumptions lacking definitive experimental support.     

Formation and carbon monoxide‐dependent dissociation of Allochromatium vinosum cytochrome c′ oligomers using domain‐swapped dimers

The number of artificial protein supramolecules has been increasing; however, control of protein oligomer formation remains challenging. Cytochrome c′ from Allochromatium vinosum (AVCP) is a homodimeric protein in its native form, where its protomer exhibits a four‐helix bundle structure containing a covalently bound five‐coordinate heme as a gas binding site. AVCP exhibits a unique reversible dimer–monomer transition according to the absence and presence of CO. Herein, domain‐swapped dimeric AVCP was constructed and utilized to form a tetramer and high‐order oligomers. The X‐ray crystal structure of oxidized tetrameric AVCP consisted of two monomer subunits and one domain‐swapped dimer subunit, which exchanged the region containing helices αA and αB between protomers. The active site structures of the domain‐swapped dimer subunit and monomer subunits in the tetramer were similar to those of the monomer subunits in the native dimer. The subunit–subunit interactions at the interfaces of the domain‐swapped dimer and monomer subunits in the tetramer were also similar to the subunit–subunit interaction in the native dimer. Reduced tetrameric AVCP dissociated to a domain‐swapped dimer and two monomers upon CO binding. Without monomers, the domain‐swapped dimers formed tetramers, hexamers, and higher‐order oligomers in the absence of CO, whereas the oligomers dissociated to domain‐swapped dimers in the presence of CO, demonstrating that the domain‐swapped dimer maintains the CO‐induced subunit dissociation behavior of native ACVP. These results suggest that protein oligomer formation may be controlled by utilizing domain swapping for a dimer–monomer transition protein.     

Probing the primary structural determinants of streptokinase inter-domain linkers by site-specific substitution and deletion mutagenesis

The bacterial protein streptokinase (SK) contains three independently folded domains (α, β and γ), interconnected by two flexible linkers with noticeable sequence homology. To investigate their primary structure requirements, the linkers were swapped amongst themselves i.e. linker 1 (between α and β domains) was swapped with linker 2 (between β and γ domains) and vice versa. The resultant construct exhibited very low activity essentially due to an enhanced proteolytic susceptibility. However, a SK mutant with two linker 1 sequences, which was proteolytically as stable as WT-rSK retained about 10% of the plasminogen activator activity of rSK When the native sequence of each linker was substituted with 9 consecutive glycine sequences, in case of the linker 1 substitution mutant substantial activity was seen to survive, whereas the linker 2 mutant lost nearly all its activity. The optimal length of linkers was then studied through deletion mutagenesis experiments, which showed that deletion beyond three residues in either of the linkers resulted in virtually complete loss of activator activity. The effect of length of the linkers was then also examined by insertion of extraneous pentapeptide sequences having a propensity for adopting either an extended conformation or a relatively rigid conformation. The insertion of poly-Pro sequences into native linker 2 sequence caused up to 10-fold reduction in activity, whereas its effect in linker 1 was relatively minor. Interestingly, most of the linker mutants could form stable 1:1 complexes with human plasminogen. Taken together, these observations suggest that (i) the functioning of the inter-domain linkers of SK requires a critical minimal length, (ii) linker 1 is relatively more tolerant to insertions and sequence alterations, and appears to function primarily as a covalent connector between the α and β domains, and (iii) the native linker 2 sequence is virtually indispensable for the activity of SK probably because of structural and/or flexibility requirements in SK action during catalysis.     

Effect of modification of the length and flexibility of the acyl carrier protein-thioesterase interdomain linker on functionality of the animal fatty acid synthase

A natural linker of approximately 20 residues connects the acyl carrier protein with the carboxy-terminal thioesterase domain of the animal fatty acid synthase. This study examines the effects of changes in the length and amino acid composition of this linker on catalytic activity, product composition, and segmental motion of the thioesterase domain. Deletion of 10 residues, almost half of the interdomain linker, had no effect on either mobility of the thioesterase domain, estimated from fluorescence polarization of a pyrenebutyl methylphosphono moiety bound covalently to the active site serine residue, or functionality of the fatty acid synthase; further shortening of the linker limited mobility of the thioesterase domain and resulted in reduced fatty acid synthase activity and an increase in product chain length from 16 to 18 and 20 carbon atoms. Surprisingly, however, even when the entire linker region was deleted, the fatty acid synthase retained 28% activity. Lengthening of the linker, by insertion of an unusually long acyl carrier protein-thioesterase linker from a modular polyketide synthase, increased mobility of the thioesterase domain without having any significant effect on catalytic properties of the complex. Interdomain linkers could also be used to tether, to the acyl carrier protein domain of the fatty acid synthase, a thioesterase active toward shorter chain length acyl thioesters generating novel short-chain fatty acid synthases. These studies reveal that although truncation of the interdomain linker partially impacts the ability of the thioesterase domain to terminate growth of the acyl chain, the overall integrity of the fatty acid synthase is quite tolerant to moderate changes in linker length and flexibility. The retention of fatty acid synthesizing activity on deletion of the entire linker region implies that the inherent flexibility of the phosphopantetheine "swinging arm" also contributes significantly to the successful docking of the long-chain acyl moiety in the thioesterase active site.     

Interdomain orientation of cardiac Troponin C characterized by paramagnetic relaxation enhancement NMR reveals a compact state

Cardiac troponin C (cTnC) is the calcium binding subunit of the troponin complex that triggers the thin filament response to calcium influx into the sarcomere. cTnC consists of two globular EF-hand domains (termed the N- and C-domains) connected by a flexible linker. While the conformation of each domain of cTnC has been thoroughly characterized through NMR studies involving either the isolated N-domain (N-cTnC) or C-domain (C-cTnC), little attention has been paid to the range of interdomain orientations possible in full-length cTnC that arises as a consequence of the flexibility of the domain linker. Flexibility in the domain linker of cTnC is essential for effective regulatory function of troponin. We have therefore utilized paramagnetic relaxation enhancement (PRE) NMR to assess the interdomain orientation of cTnC. Ensemble fitting of our interdomain PRE measurements reveals that isolated cTnC has considerable interdomain flexibility and preferentially adopts a bent conformation in solution, with a defined range of relative domain orientations.     

Amino‐acid composition after loop deletion drives domain swapping

Rational engineering of a protein to enable domain swapping requires an understanding of the sequence, structural and energetic factors that favor the domain‐swapped oligomer over the monomer. While it is known that the deletion of loops between β‐strands can promote domain swapping, the spliced sequence at the position of the loop deletion is thought to have a minimal role to play in such domain swapping. Here, two loop‐deletion mutants of the non‐domain‐swapping protein monellin, frame‐shifted by a single residue, were designed. Although the spliced sequence in the two mutants differed by only one residue at the site of the deletion, only one of them (YEIKG) promoted domain swapping. The mutant containing the spliced sequence YENKG was entirely monomeric. This new understanding that the domain swapping propensity after loop deletion may depend critically on the chemical composition of the shortened loop will facilitate the rational design of domain swapping.     

ATPase Domain and Interdomain Linker Play a Key Role in Aggregation of Mitochondrial Hsp70 Chaperone Ssc1

The co-chaperone Hep1 is required to prevent the aggregation of mitochondrial Hsp70 proteins. We have analyzed the interaction of Hep1 with mitochondrial Hsp70 (Ssc1) and the determinants in Ssc1 that make it prone to aggregation. The ATPase and peptide binding domain (PBD) of Hsp70 proteins are connected by a linker segment that mediates interdomain communication between the domains. We show here that the minimal Hep1 binding entity of Ssc1 consists of the ATPase domain and the interdomain linker. In the absence of Hep1, the ATPase domain with the interdomain linker had the tendency to aggregate, in contrast to the ATPase domain with the mutated linker segment or without linker, and in contrast to the PBD. The closest homolog of Ssc1, bacterial DnaK, and a Ssc1 chimera, in which a segment of the ATPase domain of Ssc1 was replaced by the corresponding segment from DnaK, did not aggregate in Δhep1 mitochondria. The propensity to aggregate appears to be a specific property of the mitochondrial Hsp70 proteins. The ATPase domain in combination with the interdomain linker is crucial for aggregation of Ssc1. In conclusion, our results suggest that interdomain communication makes Ssc1 prone to aggregation. Hep1 counteracts aggregation by binding to this aggregation-prone conformer.     

Probing conservation of HAMP linker structure and signal transduction mechanism through analysis of hybrid sensor kinases

The HAMP linker, a predicted structural element observed in many sensor kinases and methyl-accepting chemotaxis proteins, transmits signals between sensory input modules and output modules. HAMP linkers are located immediately inside the cytoplasmic membrane and are predicted to form two short amphipathic alpha-helices (AS-1 and AS-2) joined by an unstructured connector. HAMP linkers are found in the Escherichia coli nitrate- and nitrite-responsive sensor kinases NarX and NarQ (which respond to ligand by increasing kinase activity) and the sensor kinase CpxA (which responds to ligand by decreasing kinase activity). We constructed a series of hybrids with fusion points throughout the HAMP linker, in which the sensory modules of NarX or NarQ are fused to the transmitter modules of NarX, NarQ, or CpxA. A hybrid of the NarX sensor module and the CpxA HAMP linker and transmitter module (NarX-CpxA-1) responded to nitrate by decreasing kinase activity, whereas a hybrid in which the HAMP linker of NarX was replaced by that of CpxA (NarX-CpxA-NarX-1) responded to nitrate by increasing kinase activity. However, sequence variations between HAMP linkers do not allow free exchange of HAMP linkers or their components. Certain deletions in the NarX HAMP linker resulted in characteristic abnormal responses to ligand; similar deletions in the NarQ and NarX-CpxA-1 HAMP linkers resulted in responses to ligand generally similar to those seen in NarX. We conclude that the structure and action of the HAMP linker are conserved and that the HAMP linker transmits a signal to the output domain that ligand is bound.     

Solution structure of calcium-saturated cardiac troponin C bound to cardiac troponin I

Cardiac troponin C (TnC) is composed of two globular domains connected by a flexible linker. In solution, linker flexibility results in an ill defined orientation of the two globular domains relative to one another. We have previously shown a decrease in linker flexibility in response to cardiac troponin I (cTnI) binding. To investigate the relative orientation of calcium-saturated TnC domains when bound to cTnI, (1)H-(15)N residual dipolar couplings were measured in two different alignment media. Similarity in alignment tensor orientation for the two TnC domains supports restriction of domain motion in the presence of cTnI. The relative spatial orientation of TnC domains bound to TnI was calculated from measured residual dipolar couplings and long-range distance restraints utilizing a rigid body molecular dynamics protocol. The relative domain orientation is such that hydrophobic pockets face each other, forming a latch to constrain separate helical segments of TnI. We have utilized this structure to successfully explain the observed functional consequences of linker region deletion mutants. Together, these studies suggest that, although linker plasticity is important, the ability of TnC to function in muscle contraction can be correlated with a preferred domain orientation and interdomain distance.     

Structural mobility of the extracellular ligand-binding core of an ionotropic glutamate receptor. Analysis of NMR relaxation dynamics

Ionotropic glutamate receptors play important roles in a variety of neuronal processes and have been implicated in multiple neurodegenerative diseases. The extracellular ligand-binding (S1S2) core of the GluR2 subtype can be expressed in bacteria as a soluble, monomeric protein with binding properties essentially identical to those of the intact receptor. The crystal structure of this protein has been determined in the presence and absence of various agonists and antagonists [Armstrong, N., Sun, Y., Chen, G. Q., and Gouaux, E. (1998) Nature 395, 913-917; Armstrong, N., and Gouaux, E. (2000) Neuron 28, 165-181]. The protein consists of two lobes, with the S1 segment composing the majority of lobe 1 and the S2 segment composing most of lobe 2. A domain closure upon ligand binding has been postulated, but details of intradomain motions have not been investigated. In this paper, the backbone motions of the ligand-binding core of GluR2 bound to glutamate were studied using (15)N longitudinal (T1) and transverse (T2) relaxation measurements as well as [1H]-15N nuclear Overhauser effects at 500 and 600 MHz. Residues in the agonist-binding pocket exhibited two main classes of motion. Those contacting the alpha-substituents of the ligand glutamate exhibited minimal internal motion, while those contacting the gamma-constituents exhibited exchange dynamics, indicating two dynamically distinct portions of the binding pocket. Also, two residues in transdomain linkers between lobes 1 and 2 show exchange, lending new insight into the previously proposed domain closure hypothesis. Finally, concerted motion of helix F suggests a pathway for ligand dissociation without the necessity of domain reopening.