key: cord-0733851-3fciye8b authors: Spadaccini, Roberta; Ercole, Carmine; Graziano, Giuseppe; Wechselberger, Rainer; Boelens, Rolf; Picone, Delia title: Mechanism of 3D domain swapping in bovine seminal ribonuclease date: 2013-12-13 journal: FEBS J DOI: 10.1111/febs.12651 sha: d78487218c2b05ca96872e854621088c366ab2e2 doc_id: 733851 cord_uid: 3fciye8b 3D domain swapping (3D‐DS) is a complex protein aggregation process for which no unique mechanism exists. We report an analysis of 3D‐DS in bovine seminal ribonuclease, a homodimeric protein whose subunits are linked by two disulfide bridges, based on NMR and biochemical studies. The presence of the covalent bonds between the subunits stabilizes the unswapped dimer, and allows distinct evaluation of the structural and dynamic effects of the swapping with respect to the dimerization process. In comparison with the monomeric subunit, which, in solution has a compact structure without any propensity for local unfolding, both swapped and unswapped dimers show increased flexibility. NMR analysis, together with urea denaturation and hydrogen–deuterium exchange data, indicates that the two dimers have increased conformational fluctuations. Furthermore, we found that the rate‐limiting step of both the swapping and unswapping pathways is the detachment of the N‐terminal helices from the monomers. These results suggest a new general mechanism in which a dimeric intermediate could facilitate 3D‐DS in globular proteins. STRUCTURED DIGITAL ABSTRACT: http://www.uniprot.org/uniprot/P00669 and http://www.uniprot.org/uniprot/P00669 http://www.ebi.ac.uk/ontology-lookup/?termId=MI:0407 by http://www.ebi.ac.uk/ontology-lookup/?termId=MI:0077 (http://www.ebi.ac.uk/intact/interaction/EBI-8870415). Despite the renewed interest in 3D domain swapping (3D-DS), partly because of its association with the aggregation mechanism of amyloidogenic proteins [1, 2] , the details of the process are not known, and nor can the tendency of a protein to swap be predicted [3] . According to the Eisenberg scheme [4] , 3D-DS is prompted by the partial opening of a monomer structure, which is followed by the association of the partially unfolded proteins: in the swapped dimer, a new contact region between the two subunits is established, defining the so-called open interface, and the two dislocated regions reconstitute the previous inter-face (called closed) . This mechanism appears to be plausible, and only the first step has been revised, because experimental data have suggested two alternative scenarios: local partial unfolding of monomers in several cases [5] , but total monomer unfolding in other cases [6, 7] . Following the partial unfolding scheme, the Gibbs energy balance associated with 3D-DS should consist of: (a) the energy loss necessary to detach the swapping domain and unfold the hinge region, where prolines are often present [8] ; (b) a loss in translational entropy for dimerization; and (c) an energy gain resulting from the new interactions occur-ring at the open interface. On the other hand, in the total unfolding scheme, the energy loss for the first step is significantly larger, and the probability of such an event is considerably lower. However, by comparing the kinetic and thermodynamic parameters of several proteins swapping under different experimental conditions, Liu and Huang recently suggested the general occurrence of total unfolding [9] . They propose that this intermediate step should be necessary not only in the case of proteins requiring chemical or physical stress to swap (such as RNase A and cyanovirin-N), for which global unfolding has already been suggested [6, 7] , but also for the C-terminal domain of severe acute respiratory syndrome coronavirus main protease, for which other authors have found swapping without unfolding [5] . In all of the above-mentioned cases, 3D-DS is always associated with a monomer-to-swapped dimer conversion, which complicates the analysis, because of the need to consider the loss in translational entropy. Therefore, in this study we have investigated the swapping mechanism in an interesting protein representing a special case of intramolecular 3D-DS. Bovine seminal ribonuclease (BS-RNase) is, indeed, biosynthesized as a monomer, but it is converted into a homodimer through two disulfide bonds, between Cys31 and Cys32 of one subunit and Cys32 and Cys31, respectively, of the other subunit. The native protein is isolated as a mixture of unswapped (M=M) and swapped (M9M) isoforms, in a molar equilibrium ratio [M9M]/[M=M] of 2 : 1 [10] . Interestingly, the swapping induces further entangling of the BS-RNase subunits, and elicits additional biological properties of the protein. BS-RNase shows selective cytotoxicity for tumor cells, which is strongly related to the enzymatic activity [11] . Moreover, besides the RNA hydrolytic function, other essential features that endow BS-RNase with these additional properties are the ability to strongly interact with the tumor cell membranes [12, 13] , which results in protein internalization via endocytosis [14] , and the ability to evade the binding of RNase inhibitor [15] , a protein found at high concentrations in the cytosol that strongly binds, and inhibits, most monomeric ribonucleases [16, 17] . It has been proposed that this latter property depends on the swapped isoform, because it retains a dimeric structure even under the reducing cytosol conditions [18, 19] , which hamper the binding by RNase inhibitor. From a structural point of view, the occurrence of the intersubunit disulfide bridges makes BS-RNase a unique case among swapping proteins, because the swapped [Protein Data Bank (PDB) entry: 1BSR] and unswapped (PDB entry: 1R3M) isoforms show very minor structural differences, not only at the closed interface, but also at the open interface. Furthermore, the identical quaternary structure provides the unique opportunity to study the sole effects of the interchange of the swapping elements. In vitro, BS-RNase dislocates its N-terminal helices under very mild temperature and pH conditions, which necessarily involve only local unfolding. In the past, this feature has been attributed to the high flexibility of its 16-22 hinge region [20, 21] , where a proline is present, as in most cases of 3D-DS [8] . However, mutagenesis studies have suggested that the hinge sequence is not responsible for the swapping process [22] . Furthermore, the monomeric derivative of BS-RNase (mBS) in solution behaves like a compact, globular protein, without any detectable propensity for local unfolding [23] . Here, we report an analysis of the interchange of the N-terminal helices in BS-RNase based on a combination of NMR and kinetic data, and we show how the dimerization process can facilitate the 3D-DS process. A comparison of the structural properties of the two dimeric, purified BS-RNase isoforms was performed by means of NMR measurements. Swapped and unswapped conformers showed a superimposable set of resonances (Fig. 1A) , indicative of a very similar structure, in line with X-ray diffraction data. Comparison of backbone chemical shifts with those of mBS revealed that only a few resonances were shifted (Fig. 1A) , most of them belonging to residues located in close proximity to the open interface and the hinge loop. Thus, with simple 1 H/ 15 N-HSQC experiments, it was possible to monitor the swapping process by incubating the samples of the unswapped dimer at 37°C and acquiring spectra at 25°C (to suppress any further swapping) every 2 h over several days. Given the high degree of similarity between the NMR spectra, most signals of the two isoforms overlapped to a high degree. Nevertheless, we were able to monitor the process by following the intensity of a few peaks corresponding to residues located in different regions of the protein surface, i.e. at the open interface (Leu28), and in the b-sheet region (Asp83, Lys98, and Thr99). An overlay of the respective regions in the 1 H/ 15 N-HSQC spectra at different times during the swapping process is shown in Fig. 1B . The time dependence of the signal intensity, over the whole time interval explored, is shown in Fig. 1C . The peak intensities decreased in a similar manner for all of the residues, even in regions far away from the dislocating helices, indicating that the swapping process affects the local environment of the protein residues in a concerted way, but we have not used these data to try to obtain rate constant values. Furthermore, in none of the spectra did we have evidence of populated equilibrium intermediates. In addition, the hinge peptide linking the N-terminal helix to the main protein body was characterized by similar conformational flexibility in both dimers, as indicated by the 1 H/ 15 N heteronuclear NOE signals (Fig. S1 ). On the other hand, the flexibility of this region in both dimers is smaller than in the same region in mBS, as a consequence of the interactions occurring between the side chains of hinge residues and the corresponding partner subunit [22] . This result indicates that the conformational flexibility of the hinge peptide is not sufficient to account for the propensity to dislocate the N-terminal helices, as it has been previously shown that, in aqueous solution, mBS behaves like a stable globular protein, without any tendency to detach the N-terminal helix [23, 24] . As the dimeric forms showed an increased propensity to detach the N-terminal helices, a further analysis of the protein conformational fluctuations in solution was performed by com- paring the early stages of urea-induced denaturation and the hydrogen-deuterium (H-D) exchange rates of both BS-RNase isoforms with those of mBS by means of NMR measurements. Figure 2A shows the effects of the addition of 1 M urea on mBS, by showing the residues whose chemical shifts were displaced more than 0.1 p.p.m. upon addition of 1 M urea, and the residues that underwent a displacement between 0.05 and 0.1 p.p.m. The same data are reported in the form of a bar-graph in Fig. S2 . Figure 2B ,C shows the effect of the addition of 1 M urea on unswapped and swapped BS-RNase. Comparison of the three panels of Fig. 2 indicates that urea accessibility was increased in both dimers with respect to mBS, as, in the presence of 1 M urea, both dimers were affected by urea addition more than mBS. The different behavior between monomer and dimers was more evident at the interface regions, at the so-called deamidation loop (region 65-72), and at the C-terminal region. In addition, the N-terminal helices proved to be more accessible in the unswapped dimer than in the swapped one, in agreement with the higher conformational stability of the latter isoform, which is more abundant in solution. Furthermore, significant differences in regions crucial for swapping between monomer and dimers emerged at higher urea concentration, as shown in Fig. S3 . This shows differences between the two dimers, indicating that the N-terminal helices and the hinge peptides were more affected in the unswapper isoform than in the swapped isoform. Additional evidence for the increased conformational flexibility of the dimeric isoforms emerged from the H-D exchange experiments. On the basis of the protection factors, calculated as described in Experimental procedures, residues were divided into three classes: slow-exchanging, mediumexchanging, and fast-exchanging (Fig. 3) . The comparison of Fig. 3A with Fig. B ,C suggests a significant difference between the monomeric and dimeric isoforms. Most residues in the N-terminal helix passed from the slow-exchange regime in mBS (Fig. 3A) to the medium-exchange regime in the swapped isoform (Fig. 3B) , and even to the fast-exchange regime in the unswapped isoform (Fig. 3C) . Again, among the dimers, the unswapped isoform showed the highest conformational fluctuations in solution. Overall, the experimental data reported in Figs 2 and 3 suggest that dimerization increases the accessibility and the conformational flexibility of the regions playing the fundamental role in the domain swapping of BS-RNase, and that, in particular, the unswapped dimer is the most plastic form in solution. To complete our investigation of swapping in BS-RNase, we investigated the temperature dependence of the rate of conversion of both isoforms by using a well-established biochemical protocol for the quantitative evaluation of isoform amounts, which has already been used to monitor the time course of both the swapping and unswapping processes in native and recombinant BS-RNase [10, 22] , as well as in some mutants [22, 25] . It is worth noting that, as in all other known cases of 3D-DS, the isoform interconversion reactions can be studied solely over a very narrow temperature range, as indicated in Fig. S4A ,B. A standard kinetic analysis allowed the evaluation of the rate constants for the two processes, on the assumption that both swapping and unswapping follow pseudofirst-order kinetics at the different temperatures. These rate constant values were used to construct two Arrhenius plots (Fig. S4C) . Both swapping and unswapping yielded the same activation enthalpy value (99 AE 12 kcal mol À1 and 97 AE 11 kcal mol À1 , respectively). This value corresponds, within the limits of experimental error and with a change of the sign, to twice the enthalpy change associated with the binding of the Speptide to the S-protein: (a) -(43 AE 2) kcal mol À1 at 25°C and pH 5.0, via direct isothermal titration calorimetry measurements [26] ; (b) -(44 AE 3) kcal mol À1 at 39°C and pH 7.0, via a specific analysis of differential scanning calorimetry data for the denaturation of S-peptide-S-protein complexes [27] ; and (c) -(48 AE 3) kcal mol À1 from the kinetics of dissociation of Speptide from S-protein at 30°C and pH 6.7 [28, 29] . This finding indicates that: (a) the rate-limiting step of both the swapping and the unswapping pathways is the detachment of the two N-terminal a-helices from the main body of the monomers; and (b) swapping and unswapping occur via local partial unfolding. This does appear to be plausible, because, for both isoforms, the productive kinetic intermediate should necessarily have the two N-terminal helices detached, and so the magnitude of the activation enthalpy should be similar. However, the kinetics of the swapping process are markedly faster than those of the unswapping process. At 37°C, k (M=M ? M9M) = 0.186 h À1 , whereas k (M9M ?M=M) = 0.041 h À1 . This might be because the two closed interfaces in the swapped isoform are not independent of each other, because they are composed of residues belonging to both monomers, and so the probability of having both the N-terminal helices detached is lower than in the unswapped isoform. It is worth recalling here that the closed interface hosts the catalytic residues of BS-RNase (His12, Lys41, and His119), which, in the swapped isoform, are contributed by the two different monomers. The crosstalk between the two active sites of the swapped isoform is confirmed by the negative cooperativity characterizing its catalytic mechanism; the occupation of one active site hinders the occupation of the other active site by a substrate molecule. This should lead to a lower value of the pre-exponential factor (which is related to the activation entropy) in the Arrhenius equation for the unswapping process (Fig. S4) , giving an increase in the lower kinetic constant values. It should be emphasized that, at equilibrium, the swapped form is favored, having a concentration approximately twice that of the unswapped form (at 37°C), but the standard Gibbs energy difference between the two isoforms amounts to only 0.5 kcal mol À1 at 25°C (i.e. DG = ÀRT ln K). This fact cannot be explained only on a structural basis, because it implies a delicate balance of several factors, among which are the accessibility and dynamics of the solvation shell that, as suggested by the H-D exchange and urea susceptibility data described above, are quite different between the two isoforms. The singularity of BS-RNase, isolated as a mixture of swapped and unswapped homodimeric isoforms that interconvert under mild conditions, allows separate evaluation of the effects of dimerization and swapping processes. This represents an interesting opportunity with respect to all other cases reported in the literature, which are limited to monomer-swapped dimer conversion. On the other hand, the presence of multiple, interconverting, BS-RNase isoforms (monomer, unswapped dimer, swapped dimer, and reducedswapped dimer), together with some possible intermediates [30] , makes a comprehensive study of the features of the different dimeric isoforms in solution very complicated. Our comparative analysis among the structural and dynamic features of the swapped and unswapped dimers with respect to mBS indicates that dimerization promotes swapping, influencing the protein dynamics and conformational fluctuations. This is also in agreement with the increased susceptibility to trypsin cleavage of the two isolated dimeric isoforms in comparison with mBS, reported in a careful biochemical study [31] . According to our data, as summarized in Fig. 4 (pathway A) , we propose that the new interactions established at the open interface should facilitate swapping, not only by preparing the correct scaffold [20] , but also by increasing conformational fluctuations, thus favoring the transient local unfolding of the N-terminal helices and hinge regions reported in Eisenberg's scheme (pathway B), which culminates in the dislocation of N-terminal helices. If this mechanism were also correct for other proteins, the tendency to swap would depend not only on the presence of 'strained regions' in the monomeric structure [32] , but also on the presence of a potentially good open interface. Studies on RNase A derivatives seem to confirm this hypothesis, as mutants engineered with the two interchain disulfide bridges and the BS-RNase hinge yielded unswapped dimers that dislocated N-terminal ends under mild, physiological conditions [33] . The observed increase in RNase A domain-swapped aggregates induced by crowding agents is also in line with this hypothesis [34] . In conclusion, it is not possible to determine whether the BS-RNase pathway of 3D-DS represents a unique case, or whether it can be considered to be a more general scheme, even though it is tempting to speculate that BS-RNase cannot be the only protein that has evolved in this direction, and dimerization could precede and favor swapping in other proteins. This hypothesis emphasizes the role of the two interchain disulfide bridges, because the unswapped isoform may constitute a special 'covalent' case of an unswapped noncovalent dimer that could represent an obligatory intermediate species, even in typical (i.e. monomer-to-dimer) 3D-DS processes. The expression and purification of labeled and unlabeled proteins have been described in detail elsewhere [35] . The amount of swapped protein as a function of temperature was evaluated as described by Ercole et al. [36] . The kinetics of the swapping-unswapping process was followed at four different temperatures: 30°C, 33°C, 35°C, and 37°C. For the kinetic analysis, interconversion rate constants at the various temperatures were obtained by fitting the first data points on the assumption of first-order kinetics. The rate constants were determined at the beginning of the two processes, and so they should not be representative of the situation at thermodynamic equilibrium. At 37°C, the ratio k (M=M ? M9M) = 0.186 h À1 /k (M9M ? M=M) = 0.041 h À1 does not correspond to K = 2. Fig. 4 . Mechanisms of domain swapping. Each monomeric subunit has helical elements involved in the swapping process, respectively indicated by green and pink rectangles. The regions constituting the new ('open') interface in the dimer are also colored. According to our scheme (pathway A), the dimerization establishes intersubunit interactions (represented by faded rectangles) that destabilize the dimer structure and promote the dislocation of exchanging moieties under mild conditions. In this scheme, the 'open interface' is formed before the swapping. For comparison, the 3D-DS mechanism according to the Eisenberg scheme is reported in pathway B. It consists of the partial opening of the monomeric structure, followed by the association of the partially unfolded proteins to form the swapped dimer. In this scheme, both 'open' and 'closed' interfaces are formed in the final step of the process. NMR spectra for structure refinement were acquired at 298 K on Bruker spectrometers operating at 500 MHz, 600 MHz, 700 MHz, and 750 MHz, equipped with tripleresonance gradient probes. Data were processed with NMR-PIPE [37] , and visualization of spectra, peak-picking and analysis were performed in NMRVIEW [38] . Chemical shift assignments were obtained from standard 3D triple-resonance experiments performed on 15 N/ 13 C-labeled samples of unswapped and swapped dimers [39] . To follow the conversion of the swapped dimer in the mixture of unswapped and swapped dimers, a series of 1 H/ 15 N heteronuclear single quantum coherence (HSQC) spectra was recorded at 298 K. Each experiment at 298 K was followed by warming at 310 K for 2 h and subsequent cooling at 298 K to acquire the next HSQC spectrum. The urea-induced unfolding of BS-RNase was performed by adding solid urea (urea d4 ; Sigma Aldrich) directly to the NMR tube. The urea concentration was increased by 1 M at each titration step up to a total final concentration of 7 M. In order to keep the protein concentration and pH constant, we prepared a different protein sample for each point of the titration, and for each of them we adjusted the pH to 5.7. To check for the possible volume increase caused by the addition of solid denaturant, the length of the sample in the NMR tube was measured at each step, and found to be constant (within the limits of experimental error). The unfolding process was monitored with a series of 2D HSQC experiments performed at 298 K. As described by Piccoli et al. [10] , lyophilization affects the stability of the unswapped isoform, for this reason hydrogen exchange experiments were performed by dilution of the 15 N-labeled NMR dimeric protein samples 10-fold with 99.9% D 2 O, and concentration in a stirred ultrafiltration cell (Amicon, Millipore). H-D exchange rates were determined by fitting a first-order exponential equation to the peak heights versus time data: where I(t) represents the cross-peak height at time t, I(0) the original cross-peak height, k obs the observed rate of H-D exchange, and t the time in seconds. The intrinsic H-D exchange rates (k rc ) were calculated with the webserver program SPHERE (http://www.fccc.edu/research/labs/ roder/sphere; [40]). The protection factor (PF) for each detectable NH group was calculated from the following equation: PF = k rc /k obs Accordingly, we divided the exchanging protons into slow-exchanging (PF > 1 9 10 5 ), medium-exchanging (1 9 10 5 < PF < 1 9 10 3 ) and fast-exchanging (PF < 1 9 10 3 ). Additional supporting information may be found in the online version of this article at the publisher's web site: Fig. S1 . 1 H/ 15 N heteronuclear NOEs of swapped (M9M) and unswapped (M=M) BS-RNase isoforms. Fig. S2 . Chemical shift perturbation upon addition of 1 M urea to monomeric, unswapped and swapped BS-RNase. Fig. S3 . Urea denaturation of monomeric, unswapped and swapped BS-RNase. Fig. S4 . Kinetics of swapping and unswapping processes in BS-RNase. Amyloid-like fibrils of ribonuclease A with three-dimensional domain-swapped and native-like structure Deposition diseases and 3D domain swapping Protein acrobatics in pairs -dimerization via domain swapping Domain swapping: entangling alliances between proteins Foldon unfolding mediates the interconversion between M(pro)-C monomer and 3D domain-swapped dimer NMR spectroscopy reveals that RNase A is chiefly denatured in 40% acetic acid, implications for oligomer formation by 3D domain swapping Domain swapping proceeds via complete unfolding: a 19F-and 1H-NMR study of the Cyanovirin-N protein Proline-dependent oligomerization with arm exchange Evidences for the unfolding mechanism of three-dimensional domain swapping The dual-mode quaternary structure of seminal RNase Ribonucleases: Structures and Functions The role of electrostatic interactions in the antitumor activity of dimeric RNases Enforcing the positive charge of N-termini enhances membrane interaction and antitumor activity of bovine seminal ribonuclease Endocytotic internalization as a crucial factor for the cytotoxicity of ribonucleases Evasion of ribonuclease inhibitor as a determinant of ribonuclease cytotoxicity Crystal structure of porcine ribonuclease inhibitor, a protein with leucinerich repeats Mechanism of ribonuclease inhibition by ribonuclease inhibitor protein based on the crystal structure of its complex with ribonuclease A Structure and stability of the non-covalent swapped dimer of bovine seminal ribonuclease: an enzyme tailored to evade ribonuclease protein inhibitor The buried diversity of bovine seminal ribonuclease: shape and cytotoxicity of the swapped noncovalent form of the enzyme Swapping structural determinants of ribonucleases: an energetic analysis of the hinge peptide 16-22 The swapping of terminal arms in ribonucleases: comparison of the solution structure of monomeric bovine seminal and pancreatic ribonucleases The role of the hinge loop in domain swapping. The special case of bovine seminal ribonuclease NMR studies on structure and dynamics of the monomeric derivative of BS-RNase: new insights for 3D domain swapping From ribonuclease A toward bovine seminal ribonuclease: a step by step thermodynamic analysis Role of the hinge peptide and the intersubunit interface in the swapping of N-termini in dimeric bovine seminal RNase Thermodynamics of proteinpeptide interactions in the ribonuclease S system studied by titration calorimetry DSC study of the thermal stability of S-protein and S-peptide/S-protein Kinetic mechanism of a partial folding reaction. 1. Properties of the reaction and effects of denaturants Formation, structure, and dissociation of the ribonuclease S three-dimensional domain-swapped dimer Structurecytotoxicity relationships in bovine seminal ribonuclease: new insights from heat and chemical denaturation studies on variants Trypsin sheds light on the singular case of seminal RNase, a dimer with two quaternary conformations The unfolding story of three-dimensional domain swapping Ribonuclease A can be transformed into a dimeric ribonuclease with antitumor activity Crowding agents and osmolytes provide insight into the formation and dissociation of RNase A oligomers A new mutant of bovine seminal ribonuclease with a reversed swapping propensity Comparison of the structural and functional properties of RNase A and BS-RNase: a stepwise mutagenesis approach NMRPipe: a multidimensional spectral We thank L. Mazzarella for critical reading of the manuscript and helpful discussions. The EU NMR facility of Utrecht is gratefully acknowledged for technical assistance in the framework of the European Union project 'Bio-NMR, a new project for supporting research infrastructures in the provision of transnational access and technological advancements in bio-NMR'.