RNase A oligomerization through 3D domain swapping is favoured by a residue located far from the swapping domains

Bovine pancreatic ribonuclease A forms 3D domain-swapped oligomers by lyophilization from 40% acetic acid solutions or if subjected to various thermally-induced denaturation procedures.Considering that the intrinsic swapping propensity of bovine seminal RNase, the only member of the pancreatic-type RNase super-family that is dimeric in nature, is decreased from 70 to 30% if Arg80 is substituted by Ser (the corresponding residue in native RNase A), we introduced the opposite mutation in position 80 of the pancreatic enzyme. Our aim was to detect if the RNase A tendency to aggregate through domain swapping could increase.Aggregation of the S80R-RNase A mutant was induced either through the ‘classic’ acetic acid lyophilization, or through a thermally-induced method. The results indicate that the S80R mutant aggregates to a higher extent than the native protein, and that the increase occurs especially through N-terminal swapping.Additional investigations on the dimeric and multimeric species formed indicate that the S80R mutation increases their stability against regression to monomer, and does not significantly change their structural and functional features.     

Protein oligomerization through domain swapping: role of inter-molecular interactions and protein concentration

Domain swapping has been shown to be an important mechanism controlling multiprotein assembly and has been suggested recently as a possible mechanism underlying protein aggregation. Understanding oligomerization via domain swapping is therefore of theoretical and practical importance. By using a symmetrized structure-based (Gō) model, we demonstrate that in the free-energy landscape of domain swapping, a large free-energy barrier separates monomeric and domain-swapped dimeric configurations. We investigate the effect of finite monomer concentration, by implementing a new semi-analytical method, which involves computing the second virial coefficient, a thermodynamic indicator of inter-molecular interactions. This method, together with the symmetrized structure-based (Gō) model, minimizes the need for expensive many-protein simulations, providing a convenient framework to investigate concentration effect. Finally, we perform direct simulations of domain-swapped trimer formation, showing that this modeling approach can be used for higher-order oligomers.     

Structures of the two 3D domain-swapped RNase A trimers

When concentrated in mildly acidic solutions, bovine pancreatic ribonuclease (RNase A) forms long-lived oligomers including two types of dimer, two types of trimer, and higher oligomers. In previous crystallographic work, we found that the major dimeric component forms by a swapping of the C-terminal beta-strands between the monomers, and that the minor dimeric component forms by swapping the N-terminal alpha-helices of the monomers. On the basis of these structures, we proposed that a linear RNase A trimer can form from a central molecule that simultaneously swaps its N-terminal helix with a second RNase A molecule and its C-terminal strand with a third molecule. Studies by dissociation are consistent with this model for the major trimeric component: the major trimer dissociates into both the major and the minor dimers, as well as monomers. In contrast, the minor trimer component dissociates into the monomer and the major dimer. This suggests that the minor trimer is cyclic, formed from three monomers that swap their C-terminal beta-strands into identical molecules. These conclusions are supported by cross-linking of lysyl residues, showing that the major trimer swaps its N-terminal helix, and the minor trimer does not. We verified by X-ray crystallography the proposed cyclic structure for the minor trimer, with swapping of the C-terminal beta-strands. This study thus expands the variety of domain-swapped oligomers by revealing the first example of a protein that can form both a linear and a cyclic domain-swapped oligomer. These structures permit interpretation of the enzymatic activities of the RNase A oligomers on double-stranded RNA.     

Binding of a substrate analog to a domain swapping protein: X-ray structure of the complex of bovine seminal ribonuclease with uridylyl(2',5')adenosine

Bovine seminal ribonuclease (BS-RNase) is a unique member of the pancreatic-like ribonuclease superfamily. The native enzyme is a mixture of two dimeric forms with distinct structural features. The most abundant form is characterized by the swapping of N-terminal fragments. In this paper, the crystal structure of the complex between the swapping dimer and uridylyl(2',5')adenosine is reported at 2.06 A resolution. The refined model has a crystallographic R-factor of 0.184 and good stereochemistry. The quality of the electron density maps enables the structure of both the inhibitor and active site residues to be unambiguously determined. The overall architecture of the active site is similar to that of RNase A. The dinucleotide adopts an extended conformation with the pyrimidine and purine base interacting with Thr45 and Asn71, respectively. Several residues (Gln11, His12, Lys41, His119, and Phe120) bind the oxygens of the phosphate group. The structural similarity of the active sites of BS-RNase and RNase A includes some specific water molecules believed to be relevant to catalytic activity. Upon binding of the dinucleotide, small but significant modifications of the tertiary and quaternary structure of the protein are observed. The ensuing correlation of these modifications with the catalytic activity of the enzyme is discussed.     

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.     

Propensity for C-terminal domain swapping correlates with increased regional flexibility in the C-terminus of RNase A

Domain swapping is a type of oligomerization in which monomeric proteins exchange a structural element, resulting in oligomers whose subunits recapitulate the native, monomeric fold. It has been implicated as a potential mechanism for protein aggregation, which provides a strong impetus to understand the structural determinants and folding mechanisms that trigger domain swapping. Bovine pancreatic ribonuclease A (RNase A) is a well-studied protein known to domain swap under extreme conditions, such as lyophilization from acetic acid. The major domain-swapped dimer form of RNase A exchanges a β-strand at its C-terminus to form a C-terminal domain-swapped dimer. To study the mechanism by which C-terminal swapping occurs, we used a variant of RNase A containing a P114G mutation that readily domain swaps under physiological conditions. Using NMR and hydrogen-deuterium exchange, we find that the P114G variant has decreased protection from hydrogen exchange compared to the wild-type protein near the C-terminal hinge region. Our results suggest that domain swapping occurs via a local high-energy fluctuation at the C-terminus.     

Chain termini cross-talk in the swapping process of bovine pancreatic ribonuclease

3D domain swapping is the process by which two or more protein molecules exchange part of their structure to form intertwined dimers or higher oligomers. Bovine pancreatic ribonuclease (RNase A) is able to swap the N-terminal α-helix (residues 1-13) and/or the C-terminal β-strand (residues 116-124), thus forming a variety of oligomers, including two different dimers. Cis-trans isomerization of the Asn113-Pro114 peptide group was observed when the protein formed the C-terminal swapped dimer. To study the effect of the substitution of Pro114 on the swapping process of RNase A, we have prepared and characterized the P114A monomeric and dimeric variants of the enzyme. In contrast with previous reports, the crystal structure and NMR data on the monomer reveals a mixed cis-trans conformation for the Asn113-Ala114 peptide group, whereas the X-ray structure of the C-terminal swapped dimer of the variant is very close to that of the corresponding dimer of RNase A. The mutation at the C-terminus affects the capability of the N-terminal α-helix to swap and the stability of both dimeric forms. The present results underscore the importance of the hydration shell in determining the cross-talk between the chain termini in the swapping process of RNase A.     

The unswapped chain of bovine seminal ribonuclease: Crystal structure of the free and liganded monomeric derivative

Bovine seminal ribonuclease, a homodimeric enzyme joined covalently by two interchain disulphide bonds, is an equilibrium mixture of two conformational isomers, MxM and M=M. The major form, MxM, whose crystal structure has been previously determined at 1.9 A resolution, presents the swapping of the N-terminal segments (residues 1-15) and composite active sites formed by residues of different chains. The three-dimensional domain swapping does not occur in the M=M form. The different fold of each N-terminal tail is directed by the hinge loop (residue 16-22) connecting the swapping domain to the body of the protein. Reduction and alkylation of interchain disulphide bridges produce a monomeric derivative and a noncovalent swapped dimer, which are both active. The free and nucleotide-bound forms of the monomer have been crystallized at an alkaline pH and refined at 1.45 and 1.65 A resolution, respectively. In both cases, the N-terminal fragment is folded on the main body of the protein to produce an intact active site and a chain architecture very similar to that of bovine pancreatic ribonuclease. In this new fold of the seminal chain, the hinge loop is disordered. Despite the difference between the tertiary structure of the monomer and that of the chains in the MxM form, the active sites of the two enzymes are virtually indistinguishable. Furthermore, the structure of the liganded enzyme represents the first example of a ribonuclease complex studied at an alkaline pH and provides new information on the binding of a nucleotide when the catalytic histidines are deprotonated.     

Structural features for the mechanism of antitumor action of a dimeric human pancreatic ribonuclease variant

A specialized class of RNases shows a high cytotoxicity toward tumor cell lines, which is critically dependent on their ability to reach the cytosol and to evade the action of the ribonuclease inhibitor (RI). The cytotoxicity and antitumor activity of bovine seminal ribonuclease (BSRNase), which exists in the native state as an equilibrium mixture of a swapped and an unswapped dimer, are peculiar properties of the swapped form. A dimeric variant (HHP2‐RNase) of human pancreatic RNase, in which the enzyme has been engineered to reproduce the sequence of BSRNase helix‐II (Gln28→Leu, Arg31→Cys, Arg32→Cys, and Asn34→Lys) and to eliminate a negative charge on the surface (Glu111→Gly), is also extremely cytotoxic. Surprisingly, this activity is associated also to the unswapped form of the protein. The crystal structure reveals that on this molecule the hinge regions, which are highly disordered in the unswapped form of BSRNase, adopt a very well‐defined conformation in both subunits. The results suggest that the two hinge peptides and the two Leu28 side chains may provide an anchorage to a transient noncovalent dimer, which maintains Cys31 and Cys32 of the two subunits in proximity, thus stabilizing a quaternary structure, similar to that found for the noncovalent swapped dimer of BSRNase, that allows the molecule to escape RI and/or to enhance the formation of the interchain disulfides.     

Stability strengths and weaknesses in protein structures detected by statistical potentials: Application to bovine seminal ribonuclease

We present an in silico method to estimate the contribution of each residue in a protein to its overall stability using three database‐derived statistical potentials that are based on inter‐residue distances, backbone torsion angles and solvent accessibility, respectively. Residues that contribute very unfavorably to the folding free energy are defined as stability weaknesses, whereas residues that show a highly stabilizing contribution are called stability strengths. Strengths and/or weaknesses on residues that are in spatial contact are clustered into 3‐dimensional (3D) stability patches. The identification and analysis of strength‐ and weakness‐containing regions in a protein may reveal structural or functional characteristics, and/or interesting spots to introduce mutations. To illustrate the power of our method, we apply it to bovine seminal ribonuclease. This enzyme catalyzes the degradation of RNA strands, and has the peculiarity of undergoing 3D domain swapping in physiological conditions. The weaknesses and strengths were compared among the monomeric, dimeric and swapped dimeric forms. We identified weaknesses among the catalytic residues and a mixture of weaknesses and strengths among the substrate‐binding residues in the three forms. In the regions involved in 3D swapping, we observed an accumulation of weaknesses in the monomer, which disappear in the dimer and especially in the swapped dimer. Moreover, monomeric homologous proteins were found to exhibit less weaknesses in these regions, whereas mutants known to favor unswapped dimerization appear stabilized in this form. Our method has several perspectives for functional annotation, rational prediction of targeted mutations, and mapping of stability changes upon conformational rearrangements. Proteins 2016; 84:143–158. © 2015 Wiley Periodicals, Inc.     

Crystal structure of the UBR‐box from UBR6/FBXO11 reveals domain swapping mediated by zinc binding

The UBR‐box is a 70‐residue zinc finger domain present in the UBR family of E3 ubiquitin ligases that directly binds N‐terminal degradation signals in substrate proteins. UBR6, also called FBXO11, is an UBR‐box containing E3 ubiquitin ligase that does not bind N‐terminal signals. Here, we present the crystal structure of the UBR‐box domain from human UBR6. The dimeric crystal structure reveals a unique form of domain swapping mediated by zinc coordination, where three independent protein chains come together to regenerate the topology of the monomeric UBR‐box fold. Analysis of the structure suggests that the absence of N‐terminal residue binding arises from the lack of an amino acid binding pocket.     

Population shift vs induced fit: the case of bovine seminal ribonuclease swapping dimer

Bovine seminal ribonuclease (BS-RNase) is a unique member of the pancreatic-like ribonuclease superfamily. This enzyme exists as two conformational isomers with distinctive biological properties. The structure of the major isomer is characterized by the swapping of the N-terminal segment (MxM BS-RNase). In this article, the crystal structures of the ligand-free MxM BS-RNase and its complex with 2'-deoxycitidylyl(3',5')-2'-deoxyadenosine derived from isomorphous crystals have been refined. Interestingly, the comparison between this novel ligand-free form and the previously published sulfate-bound structure reveals significant differences. In particular, the ligand-free MxM BS-RNase is closer to the structure of MxM BS-RNase productive complexes than to the sulfate-bound form. These results reveal that MxM BS-RNase presents a remarkable flexibility, despite the structural constraints of the interchain disulfide bridges and the swapping of the N-terminal helices. These findings have important implications to the ligand binding mechanism of MxM BS-RNase. Indeed, a population shift rather than a substrate-induced conformational transition may occur in the MxM BS-RNase ligand binding process.     

Three-dimensional domain swapping as a mechanism to lock the active conformation in a super-active octamer of SARS-CoV main protease

Proteolytic processing of viral polyproteins is indispensible for the lifecycle of coronaviruses. The main protease (M^pro) of SARS-CoV is an attractive target for anti-SARS drug development as it is essential for the polyprotein processing. M^pro is initially produced as part of viral polyproteins and it is matured by autocleavage. Here, we report that, with the addition of an N-terminal extension peptide, M^pro can form a domain-swapped dimer. After complete removal of the extension peptide from the dimer, the mature M^pro self-assembles into a novel super-active octamer (AO-M^pro). The crystal structure of AO-M^pro adopts a novel fold with four domain-swapped dimers packing into four active units with nearly identical conformation to that of the previously reported M^pro active dimer, and 3D domain swapping serves as a mechanism to lock the active conformation due to entanglement of polypeptide chains. Compared with the previously well characterized form of M^pro, in equilibrium between inactive monomer and active dimer, the stable AO-M^pro exhibits much higher proteolytic activity at low concentration. As all eight active sites are bound with inhibitors, the polyvalent nature of the interaction between AO-M^pro and its polyprotein substrates with multiple cleavage sites, would make AO-M^pro functionally much more superior than the M^pro active dimer for polyprotein processing. Thus, during the initial period of SARS-CoV infection, this novel active form AO-M^pro should play a major role in cleaving polyproteins as the protein level is extremely low. The discovery of AO-M^pro provides new insights about the functional mechanism of M^pro and its maturation process.     

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.     

A single mutation induces amyloid aggregation in the alpha-spectrin SH3 domain: analysis of the early stages of fibril formation

The Src-homology region 3 domain of chicken alpha-spectrin (Spc-SH3) is a small two-state folding protein, which has never been described to form amyloid fibrils under any condition investigated so far. We show here that the mutation of asparagine 47 to alanine at the distal loop, which destabilises similarly the native and folding transition states of the domain, induces the formation of amyloid fibrils under mild acid conditions. Amyloid aggregation of the mutant is enhanced by the increase in temperature, protein concentration and NaCl concentration. The early stages of amyloid formation have been monitored as a function of time and temperature using a variety of biophysical methods. Differential scanning calorimetry experiments under conditions of amyloid formation have allowed the identification of different thermal transitions corresponding to conformational and aggregation processes as well as to the high-temperature disaggregation and unfolding of the amyloid fibrils. Aggregation is preceded by a rapid conformational change in the monomeric domain involving about 40% of the global unfolding enthalpy, considerable change in secondary structure, large loss of tertiary structure and exposure of hydrophobic patches to the solvent. The conformational change is followed by formation of a majority of oligomeric species with apparent hydrodynamic radius between 2.5 nm and 10nm, depending on temperature, together with the appearance and progressive growth of protofibrillar aggregates. After these early aggregation stages, long and curved fibrils of up to several micrometers start to develop by elongation of the protofibrils. The calorimetric data indicate that the specific enthalpy of fibril disaggregation and unfolding is relatively low, suggesting a low density of interactions within the fibril structure as compared to the native protein and a main entropy contribution to the stability of the amyloid fibrils.     

A new crystal form of the NK1 splice variant of HGF/SF demonstrates extensive hinge movement and suggests that the NK1 dimer originates by domain swapping

NK1 is a splice variant of the polypeptide growth factor HGF/SF that consists of the N terminal (N) and first kringle (K) domains and retains receptor binding and signalling. While NK1 behaves as a monomer in solution, two independent crystallographic structures have previously shown an identical, tightly packed dimer. Here we report a novel orthorhombic crystal form of NK1 at 2.5 A resolution in which four NK1 protomers are packed in two distinct dimers in the asymmetric unit. Although the basic architecture of the new NK1 dimers is similar to the two described earlier, the new crystal form demonstrates extensive hinge movement between the N and K domain that leads to re-orientation of the receptor-binding sites. The hinge bending is evidence of the paucity of strong interactions between domains within the protomer, in contrast to the extensive interactions between protomers in the dimer. These observations are consistent with domain swapping in the dimer, such that the interdomain interactions of the monomer are replaced by equivalent interprotomer interactions in the dimer and offer a route for protein engineering of NK1 variants which may act as receptor antagonists.     

Interferon-γ activates the cleavage of double-stranded RNA by bovine seminal ribonuclease

Bovine seminal ribonuclease (BS-RNase), a dimeric homologue of RNase A, cleaves both single- and double-stranded RNA and inhibits the growth of tumor cells. Its catalytic activity against double-stranded RNA, either homopolymeric ([³H]polyA/polyU) or mixed sequence, is enhanced by bovine or human recombinant interferon-γ (IFN-γ). Activation is seen with as little as 4–10 interferon units per assay. Enhancing the degradation of double-stranded RNA, an intermediate in the growth cycle of many viruses, could contribute to IFN-γ's ability to control cell growth and induce an antiviral state.     

3D structure of bovine pancreatic ribonuclease A in aqueous solution: An approach to tertiary structure determination from a small basis of1H NMR NOE correlations

Summary A method is proposed to generate initial structures in cases where the distance geometry method may fail, such as when the set of¹H NMR NOE-based distance constraints is small in relation to the size of the protein. The method introduces an initial correlation between the and backbone angles (based on empirical observations) which is relaxed in later stages of the calculation. The obtained initial structures are refined by well-established methods of energy minimization and restrained molecular dynamics. The method is applied to determine the solution structure of Ribonuclease A (124 residues) from a NOE basis consisting of 467 NOE cross-correlations (97 intra-residue, 206 sequential, 23 medium-range and 141 long-range) obtained at 360 MHz. The global shape and backbone overall fold of the eight final refined structures are close to those shown by the crystal structure. A meaningful difference in the positioning of the catalytically important His¹¹⁹ side chain in the solution and crystal structures has been detected.     

The folding mechanism of a two-domain protein: folding kinetics and domain docking of the gene-3 protein of phage fd

The gene-3 protein (G3P) of filamentous phages is essential for the infection of Escherichia coli. The carboxy-terminal domain anchors this protein in the phage coat, whereas the two amino-terminal domains N1 and N2 protrude from the phage surface. We analyzed the folding mechanism of the two-domain fragment N1-N2 of G3P (G3P(*)) and the interplay between folding and domain assembly. For this analysis, a variant of G3P(*) was used that contained four stabilizing mutations (IIHY-G3P(*)). The observed refolding kinetics extend from 10 ms to several hours. Domain N1 refolds very rapidly (with a time constant of 9.4 ms at 0.5 M guanidinium chloride, 25 degrees C) both as a part of IIHY-G3P(*) and as an isolated protein fragment. The refolding of domain N2 is slower and involves two reactions with time constants of seven seconds and 42 seconds. These folding reactions of the individual domains are followed by a very slow, spectroscopically silent docking process, which shows a time constant of 6200 seconds. This reaction was detected by a kinetic unfolding assay for native molecules. Before docking, N1 and N2 unfold fast and independently, after docking they unfold slowly in a correlated fashion. A high energy barrier is thus created by domain docking, which protects G3P kinetically against unfolding. The slow domain docking is possibly important for the infection of E.coli by the phage. Upon binding to the F pilus, the N2 domain separates from N1 and the binding site for TolA on domain N1 is exposed. Since domain reassembly is so slow, this binding site remains accessible until pilus retraction has brought N1 close to TolA on the bacterial surface.     

Small-angle X-ray scattering of reduced ribonuclease A: effects of solution conditions and comparisons with a computational model of unfolded proteins

The disulfide-reduced form of bovine ribonuclease A, with the Cys thiols irreversibly blocked, was characterized by small-angle x-ray scattering. To help resolve the conflicting results and interpretations from previous studies of this model unfolded protein, we measured scattering profiles using a range of solution conditions and compared them with the profiles predicted by a computational model for a random-coil polypeptide. Analysis of the simulated and experimental profiles reveals that scattering intensities at intermediate angles, corresponding to interatomic distances in the range of 5-20 Å, are particularly sensitive to changes in solvation and can be used to assess the internal scaling behavior of the polypeptide chain, expressed as a mass fractal dimension, D_m. This region of the scattering curve is also much less sensitive to experimental artifacts than is the very small angle regime (the Guinier region) that has been more typically used to characterize unfolded proteins. The experimental small-angle x-ray scattering profiles closely matched those predicted by the computational model assuming relatively small solvation energies. The scaling behavior of the polypeptide approaches that of a well-solvated polymer under conditions where it has a large net charge and at high urea concentrations. At lower urea concentrations and neutral pH, the behavior of the chain approaches that expected for θ-conditions, where the effects of slightly unfavorable interactions with solvent balance those of excluded volume, leading to scaling behavior comparable to that of an idealized random walk chain. Though detectable, the shift toward more compact conformations at lower urea concentrations does not correspond to a transition to a globule state and is associated with little or no reduction in conformational entropy. This type of collapse, therefore, is unlikely to greatly reduce the conformational search for the native state.