Segment swapping aided the evolution of enzyme function: The case of uroporphyrinogen III synthase

In an earlier study, we showed that two‐domain segment‐swapped proteins can evolve by domain swapping and fusion, resulting in a protein with two linkers connecting its domains. We proposed that a potential evolutionary advantage of this topology may be the restriction of interdomain motions, which may facilitate domain closure by a hinge‐like movement, crucial for the function of many enzymes. Here, we test this hypothesis computationally on uroporphyrinogen III synthase, a two‐domain segment‐swapped enzyme essential in porphyrin metabolism. To compare the interdomain flexibility between the wild‐type, segment‐swapped enzyme (having two interdomain linkers) and circular permutants of the same enzyme having only one interdomain linker, we performed geometric and molecular dynamics simulations for these species in their ligand‐free and ligand‐bound forms. We find that in the ligand‐free form, interdomain motions in the wild‐type enzyme are significantly more restricted than they would be with only one interdomain linker, while the flexibility difference is negligible in the ligand‐bound form. We also estimated the entropy costs of ligand binding associated with the interdomain motions, and find that the change in domain connectivity due to segment swapping results in a reduction of this entropy cost, corresponding to ～20% of the total ligand binding free energy. In addition, the restriction of interdomain motions may also help the functional domain‐closure motion required for catalysis. This suggests that the evolution of the segment‐swapped topology facilitated the evolution of enzyme function for this protein by influencing its dynamic properties. Proteins 2016; 85:46–53. © 2016 Wiley Periodicals, Inc.     

Functional differences between human and yeast protein disulfide isomerase family proteins

Previously, it has been reported that a mammalian protein disulfide isomerase (PDI), when expressed on a single copy number plasmid, can rescue growth of a PDI1-disrupted yeast. However, here, for the first time we demonstrated by tetrad analysis that human PDI (hPDI) is unable to replace yeast PDI (yPDI) when hPDI cDNA is integrated into the yeast chromosome. This observation indicates that hPDI is not functionally equivalent to yPDI. Estimation of the actual copy number of the plasmid, as well as comparison of isomerase and chaperone activities between human and yeast PDI homologues, indicates that one copy of hPDI cDNA is not sufficient to rescue the PDI1-disrupted strain. Notably, the isomerase activities of yPDI family proteins, Mpd1p, Mpd2p, and Eug1p, were extremely low, although yPDI itself exhibited twice as much isomerase activity as hPDI in vitro. Moreover, with the exception of Mpd1p, all hPDI and yPDI family proteins had chaperone activity, this being particularly strong in the case of yPDI and Mpd2p. These observations indicate that the growth of Saccharomyces cerevisiae is completely dependent on the isomerase activity of yPDI.     

Functional properties of the two redox-active sites in yeast protein disulphide isomerase in vitro and in vivo

Protein folding catalysed by protein disulphide isomerase (PDI) has been studied both in vivo and in vitro using different assays. PDI contains a CGHC active site in each of its two catalytic domains (a and a'). The relative importance of each active site in PDI from Saccharomyces cerevisiae (yPDI) has been analysed by exchanging the active-site cysteine residues for serine residues. The activity of the mutant forms of yPDI was determined quantitatively by following the refolding of bovine pancreatic trypsin inhibitor in vitro. In this assay the activity of the wild-type yPDI is quite similar to that of human PDI, both in rearrangement and oxidation reactions. However, while the a domain active site of the human enzyme is more active than the a'-site, the reverse is the case for yPDI. This prompted us to set up an assay to investigate whether the situation would be different with a native yeast substrate, procarboxypeptidase Y. In this assay, however, the a' domain active site also appeared to be much more potent than the a-site. These results were unexpected, not only because of the difference with human PDI, but also because analysis of folding of procarboxypeptidase Y in vivo had shown the a-site to be most important. We furthermore show that the apparent difference between in vivo and in vitro activities is not due to catalytic contributions from the other PDI homologues found in yeast.     

Reduction of protein disulfide isomerase results in open conformations and stimulates dynamic exchange between structural ensembles

Human protein disulfide isomerase (PDI) is an essential redox-regulated enzyme required for oxidative protein folding. It comprises four thioredoxin domains, two catalytically active (a, a’) and two inactive (b, b’), organized to form a flexible abb’a’ U-shape. Snapshots of unbound oxidized and reduced PDI have been obtained by X-ray crystallography. Yet, how PDI’s structure changes in response to the redox environment and inhibitor binding remains controversial. Here, we used multiparameter confocal single-molecule FRET to track the movements of the two catalytic domains with high temporal resolution. We found that at equilibrium, PDI visits three structurally distinct conformational ensembles, two “open” (O₁ and O₂) and one “closed” (C). We show that the redox environment dictates the time spent in each ensemble and the rate at which they exchange. While oxidized PDI samples O₁, O₂, and C more evenly and in a slower fashion, reduced PDI predominantly populates O₁ and O₂ and exchanges between them more rapidly, on the submillisecond timescale. These findings were not expected based on crystallographic data. Using mutational analyses, we further demonstrate that the R300-W396 cation-π interaction and active site cysteines dictate, in unexpected ways, how the catalytic domains relocate. Finally, we show that irreversible inhibitors targeting the active sites of reduced PDI did not abolish these protein dynamics but rather shifted the equilibrium toward the closed ensemble. This work introduces a new structural framework that challenges current views of PDI dynamics, helps rationalize its multifaceted role in biology, and should be considered when designing PDI-targeted therapeutics.     

Protein disulfide–isomerase,a folding catalyst and a redox-regulated chaperone

Protein disulfide–isomerase (PDI) was the first protein-folding catalyst to be characterized, half a century ago. It plays critical roles in a variety of physiological events by displaying oxidoreductase and redox-regulated chaperone activities. This review provides a brief history of the identification of PDI as both an enzyme and a molecular chaperone and of the recent advances in studies on the structure and dynamics of PDI, the substrate binding and release, and the cooperation with its partners to catalyze oxidative protein folding and maintain ER redox homeostasis. In this review, we highlight the structural features of PDI, including the high interdomain flexibility, the multiple binding sites, the two synergic active sites, and the redox-dependent conformational changes.     

KCTD5 is endowed with large,functionally relevant,interdomain motions

The KCTD family is an emerging class of proteins that are involved in important biological processes whose biochemical and structural properties are rather poorly characterized or even completely undefined. We here used KCTD5, the only member of the family with a known three-dimensional structure, to gain insights into the intrinsic structural stability of the C-terminal domain (CTD) and into the mutual dynamic interplay between the two domains of the protein. Molecular dynamics (MD) simulations indicate that in the simulation timescale (120 ns), the pentameric assembly of the CTD is endowed with a significant intrinsic stability. Moreover, MD analyses also led to the identification of exposed β-strand residues. Being these regions intrinsically sticky, they could be involved in the substrate recognition. More importantly, simulations conducted on the full-length protein provide interesting information of the relative motions between the BTB domain and the CTD of the protein. Indeed, the dissection of the overall motion of the protein is indicative of a large interdomain twisting associated with limited bending movements. Notably, MD data indicate that the entire interdomain motion is pivoted by a single residue (Ser150) of the hinge region that connects the domains. The functional relevance of these motions was evaluated in the context of the functional macromolecular machinery in which KCTD5 is involved. This analysis indicates that the interdomain twisting motion here characterized may be important for the correct positioning of the substrate to be ubiquitinated with respect to the other factors of the ubiquitination machinery.     

Activation of nanoscale allosteric protein domain motion revealed by neutron spin echo spectroscopy

NHERF1 is a multidomain scaffolding protein that assembles signaling complexes, and regulates the cell surface expression and endocytic recycling of a variety of membrane proteins. The ability of the two PDZ domains in NHERF1 to assemble protein complexes is allosterically modulated by the membrane-cytoskeleton linker protein ezrin, whose binding site is located as far as 110 Ångstroms away from the PDZ domains. Here, using neutron spin echo (NSE) spectroscopy, selective deuterium labeling, and theoretical analyses, we reveal the activation of interdomain motion in NHERF1 on nanometer length-scales and on submicrosecond timescales upon forming a complex with ezrin. We show that a much-simplified coarse-grained model suffices to describe interdomain motion of a multidomain protein or protein complex. We expect that future NSE experiments will benefit by exploiting our approach of selective deuteration to resolve the specific domain motions of interest from a plethora of global translational and rotational motions. Our results demonstrate that the dynamic propagation of allosteric signals to distal sites involves changes in long-range coupled domain motions on submicrosecond timescales, and that these coupled motions can be distinguished and characterized by NSE.     

Estimation of interdomain flexibility of N-terminus of factor H using residual dipolar couplings

Characterization of segmental flexibility is needed to understand the biological mechanisms of the very large category of functionally diverse proteins, exemplified by the regulators of complement activation, that consist of numerous compact modules or domains linked by short, potentially flexible, sequences of amino acid residues. The use of NMR-derived residual dipolar couplings (RDCs), in magnetically aligned media, to evaluate interdomain motion is established but only for two-domain proteins. We focused on the three N-terminal domains (called CCPs or SCRs) of the important complement regulator, human factor H (i.e., FH1-3). These domains cooperate to facilitate cleavage of the key complement activation-specific protein fragment, C3b, forming iC3b that no longer participates in the complement cascade. We refined a three-dimensional solution structure of recombinant FH1-3 based on nuclear Overhauser effects and RDCs. We then employed a rudimentary series of RDC data sets, collected in media containing magnetically aligned bicelles (disklike particles formed from phospholipids) under three different conditions, to estimate interdomain motions. This circumvents a requirement of previous approaches for technically difficult collection of five independent RDC data sets. More than 80% of conformers of this predominantly extended three-domain molecule exhibit flexions of <40°. Such segmental flexibility (together with the local dynamics of the hypervariable loop within domain 3) could facilitate recognition of C3b via initial anchoring and eventual reorganization of modules to the conformation captured in the previously solved crystal structure of a C3b:FH1-4 complex.     

Quality control in the endoplasmic reticulum: PDI mediates the ER retention of unassembled procollagen C-propeptides

Quality control within the endoplasmic reticulum (ER) is thought to be mediated by the interaction of a folding protein with one or several resident ER proteins [1]. Protein disulphide isomerase (PDI) is one such ER resident protein that has been previously shown to interact with proteins during their folding and assembly pathways [2, 3]. It has been assumed that, as a consequence of this interaction, unassembled proteins are retained within the ER. Here, we experimentally show that this is indeed the case. We have taken advantage of our previous finding that PDI interacts with procollagen chains early on in their assembly pathway [2] to address the role of this protein in directly retaining unassembled chains within the ER. Our experimental approach involved expressing individual C-propeptide domains from different procollagen chains in mammalian cells and determining the ability of these domains to interact with PDI and to be secreted. The C-propeptide from the proalpha2(I) chain was retained within the cell, where it formed a complex with PDI. Conversely, the C-propeptide from the proalpha1(III) chain did not form a complex with PDI and was secreted. Both domains were secreted, however, from a stable cell line expressing a secreted form of PDI lacking its ER retrieval signal. Hence, we have demonstrated directly that the intracellular retention of one substrate for ER quality control is due to an interaction with PDI.     

Model-free methods of analyzing domain motions in proteins from simulation: A comparison of normal mode analysis and molecular dynamics simulation of lysozyme

Model-free methods are introduced to determine quantities pertaining to protein domain motions from normal mode analyses and molecular dynamics simulations. For the normal mode analysis, the methods are based on the assumption that in low frequency modes, domain motions can be well approximated by modes of motion external to the domains. To analyze the molecular dynamics trajectory, a principal component analysis tailored specifically to analyze interdomain motions is applied. A method based on the curl of the atomic displacements is described, which yields a sharp discrimination of domains, and which defines a unique interdomain screw-axis. Hinge axes are defined and classified as twist or closure axes depending on their direction. The methods have been tested on lysozyme. A remarkable correspondence was found between the first normal mode axis and the first principal mode axis, with both axes passing within 3 Å of the alpha-carbon atoms of residues 2, 39, and 56 of human lysozyme, and near the interdomain helix. The axes of the first modes are overwhelmingly closure axes. A lesser degree of correspondence is found for the second modes, but in both cases they are more twist axes than closure axes. Both analyses reveal that the interdomain connections allow only these two degrees of freedom, one more than provided by a pure mechanical hinge. Proteins 27:425–437, 1997. © 1997 Wiley-Liss, Inc.     

Interpretation of NMR relaxation properties of Pin1, a two-domain protein, based on Brownian dynamic simulations

Many important proteins contain multiple domains connected by flexible linkers. Inter-domain motion is suggested to play a key role in many processes involving molecular recognition. Heteronuclear NMR relaxation is sensitive to motions in the relevant time scales and could provide valuable information on the dynamics of multi-domain proteins. However, the standard analysis based on the separation of global tumbling and fast local motions is no longer valid for multi-domain proteins undergoing internal motions involving complete domains and that take place on the same time scale than the overall motion.The complexity of the motions experienced even for the simplest two-domain proteins are difficult to capture with simple extensions of the classical Lipari-Szabo approach. Hydrodynamic effects are expected to dominate the motion of the individual globular domains, as well as that of the complete protein. Using Pin1 as a test case, we have simulated its motion at the microsecond time scale, at a reasonable computational expense, using Brownian Dynamic simulations on simplified models. The resulting trajectories provide insight on the interplay between global and inter-domain motion and can be analyzed using the recently published method of isotropic Reorientational Mode Dynamics which offer a way of calculating their contribution to heteronuclear relaxation rates. The analysis of trajectories computed with Pin1 models of different flexibility provides a general framework to understand the dynamics of multi-domain proteins and explains some of the observed features in the relaxation rate profile of free Pin1.     

Interdomain flexibility and interfacial integrity of β-lactamase inhibitory protein (BLIP) modulate its binding to class A β-lactamases

β-Lactamase inhibitory protein (BLIP) consists of a tandem repeat of αβ domains conjugated by an interdomain loop and can effectively bind and inactivate class A β-lactamases, which are responsible for resistance of bacteria to β-lactam antibiotics. The varied ability of BLIP to bind different β-lactamases and the structural determinants for significant enhancement of BLIP variants with a point mutation are poorly understood. Here, we investigated the conformational dynamics of BLIP upon binding to three clinically prevalent class A β-lactamases (TEM1, SHV1, and PC1) with dissociation constants between subnanomolar and micromolar. Hydrogen deuterium exchange mass spectrometry revealed that the flexibility of the interdomain region was significantly suppressed upon strong binding to TEM1, but was not significantly changed upon weak binding to SHV1 or PC1. E73M and K74G mutations in the interdomain region improved binding affinity toward SHV1 and PC1, respectively, showing significantly increased flexibility of the interdomain region compared to the wild-type and favorable conformational changes upon binding. In contrast, more rigidity of the interfacial loop 135–145 was observed in these BLIP mutants in both free and bound states. Consistently, molecular dynamics simulations of BLIP exhibited drastic changes in the flexibility of the loop 135–145 in all complexes. Our results indicated for the first time that higher flexibility of the interdomain linker, as well as more rigidity of the interfacial loop 135–145, could be desirable determinants for enhancing inhibition of BLIP to class A β-lactamases. Together, these findings provide unique insights into the design of enhanced inhibitors.     

Interaction interfaces of protein domains are not topologically equivalent across families within superfamilies: Implications for metabolic and signaling pathways

Using a data set of aligned protein domain superfamilies of known three-dimensional structure, we compared the location of interdomain interfaces on the tertiary folds between members of distantly related protein domain superfamilies. The data set analyzed is comprised of interdomain interfaces, with domains occurring within a polypeptide chain and those between two polypeptide chains. We observe that, in general, the interfaces between protein domains are formed entirely in different locations on the tertiary folds in such pairs. This variation in the location of interface happens in protein domains involved in a wide range of functions, such as enzymes, adapters, and domains that bind protein ligands, or cofactors. While basic biochemical functionality is preserved at the domain superfamily level, the effect of biochemical function on protein assemblies is different in these protein domains related by superfamily. The divergence between proteins, in most cases, is coupled with domain recruitment, with different modes of interaction with the recruited domain. This is in complete contrast to the observation that in closely related homologous protein domains, almost always the interaction interfaces are topologically equivalent. In a small subset of interacting domains within proteins related by remote homology, we observe that the relative positioning of domains with respect to one another is preserved. Based on the analysis of multidomain proteins of known or unknown structure, we suggest that variation in protein-protein interactions in members within a superfamily could serve as diverging points in otherwise parallel metabolic or signaling pathways. We discuss a few representative cases of diverging pathways involving domains in a superfamily.     

On the Role of Protein Disulfide Isomerase in the Retrograde Cell Transport of Secreted Phospholipases A2

Following the finding that ammodytoxin (Atx), a neurotoxic secreted phospholipase A₂ (sPLA₂) in snake venom, binds specifically to protein disulfide isomerase (PDI) in vitro we show that these proteins also interact in living rat PC12 cells that are able to internalize this group IIA (GIIA) sPLA₂. Atx and PDI co-localize in both differentiated and non-differentiated PC12 cells, as shown by fluorescence microscopy. Based on a model of the complex between Atx and yeast PDI (yPDI), a three-dimensional model of the complex between Atx and human PDI (hPDI) was constructed. The Atx binding site on hPDI is situated between domains b and b’. Atx interacts hPDI with an extensive area on its interfacial binding surface. The mammalian GIB, GIIA, GV and GX sPLA₂s have the same fold as Atx. The first three sPLA₂s have been detected intracellularly but not the last one. The models of their complexes with hPDI were constructed by replacement of Atx with the respective mammalian sPLA₂ in the Atx—hPDI complex and molecular docking of the structures. According to the generated models, mammalian GIB, GIIA and GV sPLA₂s form complexes with hPDI very similar to that with Atx. The contact area between GX sPLA₂ and hPDI is however different from that of the other sPLA₂s. Heterologous competition of Atx binding to hPDI with GV and GX sPLA₂s confirmed the model-based expectation that GV sPLA₂ was a more effective inhibitor than GX sPLA₂, thus validating our model. The results suggest a role of hPDI in the (patho)physiology of some snake venom and mammalian sPLA₂s by assisting the retrograde transport of these molecules from the cell surface. The sPLA₂–hPDI model constitutes a valuable tool to facilitate further insights into this process and into the (patho)physiology of sPLA₂s in relation to their action intracellularly.     

Substrate-mediated electron transfer in peptidylglycine alpha-hydroxylating monooxygenase.

Peptide amidation is a ubiquitous posttranslational modification of bioactive peptides. Peptidylglycine alpha-hydroxylating monooxygenase (PHM; EC 1.14.17.3), the enzyme that catalyzes the first step of this reaction, is composed of two domains, each of which binds one copper atom. The coppers are held 11 A apart on either side of a solvent-filled interdomain cleft, and the PHM reaction requires electron transfer between these sites. A plausible mechanism for electron transfer might involve interdomain motion to decrease the distance between the copper atoms. Our experiments show that PHM catalytic core (PHMcc) is enzymatically active in the crystal phase, where interdomain motion is not possible. Instead, structures of two states relevant to catalysis indicate that water, substrate and active site residues may provide an electron transfer pathway that exists only during the PHM catalytic cycle.     

Structure and functional analysis of the siderophore periplasmic binding protein from the fuscachelin gene cluster of Thermobifida fusca

Iron acquisition is a complex, multicomponent process critical for most organisms' survival and virulence. Small iron chelating molecules, siderophores, mediate transport as key components of common pathways for iron assimilation in many microorganisms. The chemistry and biology of the extraordinary tight and specific metal binding siderophores is of general interest in terms of host/guest chemistry and is a potential target toward the development of therapeutic treatments for microbial virulence. The siderophore pathway of the moderate thermophile, Thermobifida fusca, is an excellent model system to study the process in Gram‐positive bacteria. Here we describe the structure and characterization of the siderophore periplasmic binding protein, FscJ from the fuscachelin gene cluster of T. fusca. The structure shows a di‐domain arrangement connected with a long α‐helix hinge. Several X‐ray structures detail ligand‐free conformational changes at different pH values, illustrating complex interdomain flexibility of the siderophore receptors. We demonstrated that FscJ has a unique recognition mechanism and details the binding interaction with ferric‐fuscachelin A through ITC and docking analysis. The presented work provides a structural basis for the complex molecular mechanisms of siderophore recognition and transportation. Proteins 2016; 84:118–128. © 2015 Wiley Periodicals, Inc.