The role of the hinge loop in domain swapping. The special case of bovine seminal ribonuclease

Human FoxP Transcription Factors as Tractable Models of the Evolution and Functional Outcomes of Three-Dimensional Domain Swapping

The association of two or more proteins to adopt a quaternary complex is one of the most widespread mechanisms by which protein function is modulated. In this scenario, three-dimensional domain swapping (3D-DS) constitutes one plausible pathway for the evolution of protein oligomerization that exploits readily available intramolecular contacts to be established in an intermolecular fashion. However, analysis of the oligomerization kinetics and thermodynamics of most extant 3D-DS proteins shows its dependence on protein unfolding, obscuring the elucidation of the emergence of 3D-DS during evolution, its occurrence under physiological conditions, and its biological relevance. Here, we describe the human FoxP subfamily of transcription factors as a feasible model to study the evolution of 3D-DS, due to their significantly faster dissociation and dimerization kinetics and lower dissociation constants in comparison to most 3D-DS models. Through the biophysical and functional characterization of FoxP proteins, relevant structural aspects highlighting the evolutionary adaptations of these proteins to enable efficient 3D-DS have been ascertained. Most biophysical studies on FoxP suggest that the dynamics of the polypeptide chain are crucial to decrease the energy barrier of 3D-DS, enabling its fast oligomerization under physiological conditions. Moreover, comparison of biophysical parameters between human FoxP proteins in the context of their minute sequence differences suggests differential evolutionary strategies to favor homoassociation and presages the possibility of heteroassociations, with direct impacts in their gene regulation function.

Monothiol and dithiol glutaredoxin-1 from Clostridium oremlandii: identification of domain-swapped structures by NMR, X-ray crystallography and HDX mass spectrometry

Protein dimerization or oligomerization resulting from swapping part of the protein between neighboring polypeptide chains is known to play a key role in the regulation of protein function and in the formation of protein aggregates. Glutaredoxin-1 from Clostridium oremlandii (cGrx1) was used as a model to explore the formation of multiple domain-swapped conformations, which were made possible by modulating several hinge-loop residues that can form a pivot for domain swapping. Specifically, two alternative domain-swapped structures were generated and analyzed using nuclear magnetic resonance (NMR), X-ray crystallography, circular-dichroism spectroscopy and hydrogen/deuterium-exchange (HDX) mass spectrometry. The first domain-swapped structure (β3-swap) was formed by the hexameric cGrx1-cMsrA complex. The second domain-swapped structure (β1-swap) was formed by monothiol cGrx1 (C16S) alone. In summary, the first domain-swapped structure of an oxidoreductase in a hetero-oligomeric complex is presented. In particular, a single point mutation of a key cysteine residue to serine led to the formation of an intramolecular disulfide bond, as opposed to an intermolecular disulfide bond, and resulted in modulation of the underlying free-energy landscape of protein oligomerization.

Structure, stability and aggregation propensity of a Ribonuclease A-Onconase chimera

Structural roles of loop regions are frequently overlooked in proteins. Nevertheless, they may be key players in the definition of protein topology and in the self-assembly processes occurring through domain swapping. We here investigate the effects on structure and stability of replacing the loop connecting the last two β-strands of RNase A with the corresponding region of the more thermostable Onconase. The crystal structure of this chimeric variant (RNaseA-ONC) shows that its terminal loop size better adheres to the topological rules for the design of stabilized proteins, proposed by Baker and coworkers [43]. Indeed, RNaseA-ONC displays a thermal stability close to that of RNase A, despite the lack of Pro at position 114, which, due to its propensity to favor a cis peptide bond, has been identified as an important stabilizing factor of the native protein. Accordingly, RNaseA-ONC is significantly more stable than RNase A variants lacking Pro114; RNaseA-ONC also displays a higher propensity to form oligomers in native conditions when compared to either RNase A or Onconase. This finding demonstrates that modifications of terminal loops should to be carefully controlled in terms of size and sequence to avoid unwanted and/or potentially harmful aggregation processes.

Understanding protein domain-swapping using structure-based models of protein folding

In domain-swapping, two or more identical protein monomers exchange structural elements and fold into dimers or multimers whose units are structurally similar to the original monomer. Domain-swapping is of biotechnological interest because inhibiting domain-swapping can reduce disease-causing fibrillar protein aggregation. To achieve such inhibition, it is important to understand both the energetics that stabilize the domain-swapped structure and the protein dynamics that enable the swapping. Structure-based models (SBMs) encode the folded structure of the protein in their potential energy functions. SBMs have been successfully used to understand diverse aspects of monomer folding. Symmetrized SBMs model interactions between two identical protein chains using only intra-monomer interactions. Molecular dynamics simulations of such symmetrized SBMs have been used to correctly predict the domain-swapped structure and to understand the mechanism of domain-swapping. Here, we review such models and illustrate that monomer topology determines key aspects of domain-swapping. However, in some proteins, specifics of local energetic interactions modulate domain-swapping and these need to be added to the symmetrized SBMs. We then summarize some general principles of the mechanism of domain-swapping that emerge from the symmetrized SBM simulations. Finally, using our own results, we explore how symmetrized SBMs could be used to design domain-swapping in proteins.

Intermolecular β-strand networks avoid hub residues and favor low interconnectedness: a potential protection mechanism against chain dissociation upon mutation

Altogether few protein oligomers undergo a conformational transition to a state that impairs their function and leads to diseases. But when it happens, the consequences are not harmless and the so-called conformational diseases pose serious public health problems. Notorious examples are the Alzheimer's disease and some cancers associated with a conformational change of the amyloid precursor protein (APP) and of the p53 tumor suppressor, respectively. The transition is linked with the propensity of β-strands to aggregate into amyloid fibers. Nevertheless, a huge number of protein oligomers associate chains via β-strand interactions (intermolecular β-strand interface) without ever evolving into fibers. We analyzed the layout of 1048 intermolecular β-strand interfaces looking for features that could provide the β-strands resistance to conformational transitions. The interfaces were reconstructed as networks with the residues as the nodes and the interactions between residues as the links. The networks followed an exponential decay degree distribution, implying an absence of hubs and nodes with few links. Such layout provides robustness to changes. Few links per nodes do not restrict the choices of amino acids capable of making an interface and maintain high sequence plasticity. Few links reduce the “bonding” cost of making an interface. Finally, few links moderate the vulnerability to amino acid mutation because it entails limited communication between the nodes. This confines the effects of a mutation to few residues instead of propagating them to many residues via hubs. We propose that intermolecular β-strand interfaces are organized in networks that tolerate amino acid mutation to avoid chain dissociation, the first step towards fiber formation. This is tested by looking at the intermolecular β-strand network of the p53 tetramer.

Mechanism of 3D domain swapping in bovine seminal ribonuclease

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)

Bovine seminal ribonuclease triggers Beclin1-mediated autophagic cell death in pancreatic cancer cells

Among the large number of variants belonging to the pancreatic-type secretory ribonuclease (RNase) superfamily, bovine pancreatic ribonuclease (RNase A) is the proto-type and bovine seminal RNase (BS-RNase) represents the unique natively dimeric member. In the present manuscript, we evaluate the anti-tumoral property of these RNases in pancreatic adenocarcinoma cell lines and in nontumorigenic cells as normal control. We demonstrate that BS-RNase stimulates a strong anti-proliferative and pro-apoptotic effect in cancer cells, while RNase A is largely ineffective. Notably, we reveal for the first time that BS-RNase triggers Beclin1-mediated autophagic cancer cell death, providing evidences that high proliferation rate of cancer cells may render them more susceptible to autophagy by BS-RNase treatment. Notably, to improve the autophagic response of cancer cells to BS-RNase we used two different strategies: the more basic (as compared to WT enzyme) G38K mutant of BS-RNase, known to interact more strongly than wt with the acidic membrane of cancer cells, or BS-RNase oligomerization (tetramerization or formation of larger oligomers). Both mutant BS-RNase and BS-RNase oligomers potentiated autophagic cell death as compared to WT native dimer of BS-RNase, while the various RNase A oligomers remained completely ineffective. Altogether, our results shed more light on the mechanisms lying at the basis of BS-RNase antiproliferative effect in cancer cells, and support its potential use to develop new anti-cancer strategies.

The multiple forms of bovine seminal ribonuclease: structure and stability of a C-terminal swapped dimer

Bovine seminal ribonuclease (BS-RNase) acquires an interesting anti-tumor activity associated with the swapping on the N-terminal. The first direct experimental evidence on the formation of a C-terminal swapped dimer (C-dimer) obtained from the monomeric derivative of BS-RNase, although under non-native conditions, is here reported. The X-ray model of this dimer reveals a quaternary structure different from that of the C-dimer of RNase A, due to the presence of three mutations in the hinge peptide 111-116. The mutations increase the hinge peptide flexibility and decrease the stability of the C-dimer against dissociation. The biological implications of the structural data are also discussed.

Double domain swapping in bovine seminal RNase: formation of distinct N- and C-swapped tetramers and multimers with increasing biological activities

Bovine seminal (BS) RNase, the unique natively dimeric member of the RNase super-family, represents a special case not only for its additional biological actions but also for the singular features of 3D domain swapping. The native enzyme is indeed a mixture of two isoforms: M = M, a dimer held together by two inter-subunit disulfide bonds, and MxM, 70% of the total, which, besides the two mentioned disulfides, is additionally stabilized by the swapping of its N-termini.

When lyophilized from 40% acetic acid, BS-RNase oligomerizes as the super-family proto-type RNase A does. In this paper, we induced BS-RNase self-association and analyzed the multimers by size-exclusion chromatography, cross-linking, electrophoresis, mutagenesis, dynamic light scattering, molecular modelling. Finally, we evaluated their enzymatic and cytotoxic activities.

Several BS-RNase domain-swapped oligomers were detected, including two tetramers, one exchanging only the N-termini, the other being either N- or C-swapped. The C-swapping event, confirmed by results on a BS-K113N mutant, has been firstly seen in BS-RNase here, and probably stabilizes also multimers larger than tetramers.

Interestingly, all BS-RNase oligomers are more enzymatically active than the native dimer and, above all, they display a cytotoxic activity that definitely increases with the molecular weight of the multimers. This latter feature, to date unknown for BS-RNase, suggests again that the self-association of RNases strongly modulates their biological and potentially therapeutic properties.

Functional repertoire, molecular pathways and diseases associated with 3D domain swapping in the human proteome

Background

3D domain swapping is a novel structural phenomenon observed in diverse set of protein structures in oligomeric conformations. A distinct structural feature, where structural segments in a protein dimer or higher oligomer were shared between two or more chains of a protein structure, characterizes 3D domain swapping. 3D domain swapping was observed as a key mediator of numerous functional mechanisms and play pathogenic role in various diseases including conformational diseases like amyloidosis, Alzheimer's disease, Parkinson's disease and prion diseases. We report the first study with a focus on identifying functional classes, pathways and diseases mediated by 3D domain swapping in the human proteome.

Methods

We used a panel of four enrichment tools with two different ontologies and two annotations database to derive biological and clinical relevant information associated with 3D domain swapping. Protein domain enrichment analysis followed by Gene Ontology (GO) term enrichment analysis revealed the functional repertoire of proteins involved in swapping. Pathway analysis using KEGG annotations revealed diverse pathway associations of human proteins involved in 3D domain swapping. Disease Ontology was used to find statistically significant associations with proteins in swapped conformation and various disease categories (P-value < 0.05).

Results

We report meta-analysis results of a literature-curated dataset of human gene products involved in 3D domain swapping and discuss new insights about the functional repertoire, pathway associations and disease implications of proteins involved in 3D domain swapping.

Conclusions

Our integrated bioinformatics pipeline comprising of four different enrichment tools, two ontologies and two annotations revealed new insights into the functional and disease correlations with 3D domain swapping. GO term enrichment were used to infer terms associated with three different GO categories. Protein domain enrichment was used to identify conserved domains enriched in swapped proteins. Pathway enrichment analysis using KEGG annotations revealed that proteins with swapped conformations are present in all six classes of KEGG BRITE hierarchy and significantly enriched KEGG pathways were observed in five classes. Five major classes of disease were found to be associated with 3D domain swapping using functional disease ontology based enrichment analysis. Five classes of human diseases: cancer, diseases of the respiratory or pulmonary system, degenerative diseases of the central nervous system, vascular disease and encephalitis were found to be significant. In conclusion, our study shows that bioinformatics based analytical approaches using curated data can enhance the understanding of functional and disease implications of 3D domain swapping.

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.

NMR studies on structure and dynamics of the monomeric derivative of BS-RNase: new insights for 3D domain swapping

Three-dimensional domain swapping is a common phenomenon in pancreatic-like ribonucleases. In the aggregated state, these proteins acquire new biological functions, including selective cytotoxicity against tumour cells. RNase A is able to dislocate both N- and C-termini, but usually this process requires denaturing conditions. In contrast, bovine seminal ribonuclease (BS-RNase), which is a homo-dimeric protein sharing 80% of sequence identity with RNase A, occurs natively as a mixture of swapped and unswapped isoforms. The presence of two disulfides bridging the subunits, indeed, ensures a dimeric structure also to the unswapped molecule. In vitro, the two BS-RNase isoforms interconvert under physiological conditions. Since the tendency to swap is often related to the instability of the monomeric proteins, in these paper we have analysed in detail the stability in solution of the monomeric derivative of BS-RNase (mBS) by a combination of NMR studies and Molecular Dynamics Simulations. The refinement of NMR structure and relaxation data indicate a close similarity with RNase A, without any evidence of aggregation or partial opening. The high compactness of mBS structure is confirmed also by H/D exchange, urea denaturation, and TEMPOL mapping of the protein surface. The present extensive structural and dynamic investigation of (monomeric) mBS did not show any experimental evidence that could explain the known differences in swapping between BS-RNase and RNase A. Hence, we conclude that the swapping in BS-RNase must be influenced by the distinct features of the dimers, suggesting a prominent role for the interchain disulfide bridges.

Understanding the molecular mechanism of enzyme dynamics of ribonuclease A through protonation/deprotonation of HIS48

Molecular dynamics simulation is carried out to investigate the enzyme dynamics of RNase A with the HIS48 in three different states (HIP48 (protonated), HID48 (deprotonated), and H48A mutant). Insights derived from the current theoretical study, combined with the available experimental observations, enabled us to provide a microscopic picture for the efficient enzyme dynamics. Specifically, in the "closed" state or HIP48, the N-terminal hinge loop is intact and the enzyme remains in a relatively stable conformation which is preferred for catalytic reaction. Deprotonation of HIS48 induces the denaturing of this hinge-loop into a 3(10)-helix, causing it to break the original interaction network around the loop-1 and drive the partial unfolding of the N-terminal. The enhanced dynamic motion of the N-terminal helix facilitates the release of the catalytic product (the rate limiting step) and speeds up the overall catalytic process. The current study established that HIS49 acts as a modulator for the transformation of conformational states through the perturbing of hydrogen bond networks across loop-1, the N-terminal helix, and other residues nearby. Our study suggests that HIS48 may also serve to transport loop-1's kinetic energy to the reaction center.

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.

Toward an antitumor form of bovine pancreatic ribonuclease: the crystal structure of three noncovalent dimeric mutants

The cytotoxic action of bovine seminal ribonuclease (BS-RNase) depends on its noncovalent swapped dimeric form (NCD-BS), which presents a compact structure that allows the molecule to escape ribonuclease inhibitor (RI). A key role in the acquisition of this structure has been attributed to the concomitant presence of a proline in position 19 and a leucine in position 28. The introduction of Leu28, Cys31, and Cys32 and, in addition, of Pro19 in the sequence of bovine pancreatic ribonuclease (RNase A) has produced two dimeric variants LCC and PLCC, which do exhibit a cytotoxic activity, though at a much lower level than BS-RNase. The crystal structure analysis of the noncovalent swapped form (NCD) of LCC and PLCC, complexed with the substrate analogue 2 '-deoxycytidylyl(3 ',5 ')-2 '-deoxyguanosine, has revealed that, differently from NCD-BS, the dimers adopt an opened quaternary structure, with the two Leu residues fully exposed to the solvent, that does not hinder the binding of RI. Similar results have been obtained for a third mutant of the pancreatic enzyme, engineered with the hinge peptide sequence of the seminal enzyme (residues 16-22) and the two cysteines in position 31 and 32, but lacking the hydrophobic Leu residue in position 28. The comparison of these three structures with those previously reported for other ribonuclease swapped dimers strongly suggests that, in addition to Pro19 and Leu28, the presence of a glycine at the N-terminal end of the hinge peptide is also important to push the swapped form of RNase A dimer into the compact quaternary organization observed for NCD-BS.

Mechanism and energy landscape of domain swapping in the B1 domain of protein G

Three-dimensional domain swapping has emerged as a ubiquitous process for homo-oligomer formation in many unrelated proteins, but the molecular mechanism of this process is still poorly understood. Here we present a mechanism for the swapping reaction in the B1 domain of the immunoglobulin G binding protein from group G of Streptococcus (GB1). This is a particularly attractive system for investigating the swapping process, as the swapped dimer formed by the quadruple mutant (L5V/F30V/Y33F/A34F) of GB1 was recently shown to exist in equilibrium with a monomer-like conformation over time scales of minutes. According to our mechanism, swapping in GB1 starts from the C-terminus of the polypeptide chain and progresses by exchanging an increasing portion of the chains until a stable conformational state is reached. This exchange process does not involve unfolding. Rather, the conformational changes of individual monomers and their association are tightly coupled to minimize solvent exposure and maximize the total number of native contacts at all times, thereby closely approximating the minimum energy path of the reaction. Using detailed atomic descriptions, we compute the complete free-energy profiles of the exchange reaction for the GB1 quadruple mutant that forms swapped dimers and for the wild-type protein, which is monomeric. In both GB1 forms, intermediates sample a surprisingly wide range of nearly isoenergetic association modes and hinge conformations, indicating that the exchange reaction is a non-specific process akin to encounter complex formation where the amino acid sequence plays a marginal role. The main role of the mutations in the swapping process is to destabilize the GB1 monomer state, while stabilizing the swapped dimer conformation, with non-native intersubunit interactions, fostered by mutant side chains, contributing significantly to this stabilization. Our findings are rationalized in terms of a generic swapping mechanism that involves the association of activated molecular species, and it is argued that a similar mechanism may apply to swapping in other protein systems.

Overcoming residual frustration in domain-swapping: the roles of disulfide bonds in dimerization and aggregation

The prevalence of domain-swapping in nature is a manifestation of the principle of minimal frustration in that the interactions designed by evolution to stabilize the protein are also involved in this mode of binding. We previously demonstrated that the Symmetrized-Go potential accurately predicts the experimentally observed domain-swapped structure of Eps8 based solely on the structure of the monomer. There can be, however, multiple modes of domain-swapping, reflecting a higher level of frustration, which is a consequence of symmetry. The human prion and cyanovirin-N are too frustrated to form unique domain-swapped structures on the basis of the Symmetrized-Go potential. However, supplementing the completely symmetric model with intermolecular and intramolecular disulfide bonds in the prion and cyanovirin-N proteins, respectively, yielded unique domain-swapped structures with a remarkable similarity to the experimentally observed ones. These results suggest that the disulfide bonds may sometimes be critical in overcoming the intrinsic frustration of the symmetrized energy landscapes for domain-swapping. We also discuss the implications of intermolecular disulfide bonds in the formation of mammalian prion aggregates.