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

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.

Interactions crucial for three-dimensional domain swapping in the HP-RNase variant PM8

The structural determinants that are responsible for the formation of higher order associations of folded proteins remain unknown. We have investigated the role on the dimerization process of different residues of a domain-swapped dimer human pancreatic ribonuclease variant. This variant is a good model to study the dimerization and swapping processes because dimer and monomer forms interconvert, are easily isolated, and only one dimeric species is produced. Thus, simple models for the swapping process can be proposed. The dimerization (dissociation constant) and swapping propensity have been studied using different variants with changes in residues that belong to different putative molecular determinants of dimerization. Using NMR spectroscopy, we show that these mutations do not substantially alter the overall conformation and flexibility, but affect the residue level stability. Overall, the most critical residues for the swapping process are those of one subunit that interact with the hinge loop of another one-subunit residue, stabilizing it in a conformation that favors the interchange. Tyr(25), Gln(101), and Pro(19), with Asn(17), Ser(21), and Ser(23), are found to be the most significant; notably, Glu(103) and Arg(104), which were postulated to form salt bridges that would stabilize the dimer, are not critical for dimerization.

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.

Functional and structural analyses of N-acylsulfonamide-linked dinucleoside inhibitors of RNase A

Molecular probes are useful for both studying and controlling the functions of enzymes and other proteins. The most useful probes have high affinity for their target, along with small size and resistance to degradation. Here, we report on new surrogates for nucleic acids that fulfill these criteria. Isosteres in which phosphoryl [R–O–P(O₂⁻)–O–R′] groups are replaced with N-acylsulfonamidyl [R–C(O)–N⁻–S(O₂)–R′] or sulfonimidyl [R–S(O₂)–N⁻–S(O₂)–R′] groups increase the number of nonbridging oxygens from two (phosphoryl) to three (N-acylsulfonamidyl) or four (sulfonimidyl). Six such isosteres were found to be more potent inhibitors of catalysis by bovine pancreatic RNase A than are parent compounds containing phosphoryl groups. The atomic structures of two RNase A·N-acylsulfonamide complexes were determined at high resolution by X-ray crystallography. The N-acylsulfonamidyl groups were observed to form more hydrogen bonds with active site residues than did the phosphoryl groups in analogous complexes. These data encourage the further development and use of N-acylsulfonamides and sulfonimides as antagonists of nucleic acid-binding proteins.

A new RNase sheds light on the RNase/angiogenin subfamily from zebrafish

Recently, extracellular RNases of the RNase A superfamily, with the characteristic CKxxNTF sequence signature, have been identified in fish. This has led to the recognition that these RNases are present in the whole vertebrate subphylum. In fact, they comprise the only enzyme family unique to vertebrates. Four RNases from zebrafish (Danio rerio) have been previously reported and have a very low RNase activity; some of these are endowed, like human angiogenin, with powerful angiogenic and bactericidal activities. In the present paper, we report the three-dimensional structure, the thermodynamic behaviour and the biological properties of a novel zebrafish RNase, ZF-RNase-5. The investigation of its structural and functional properties, extended to all other subfamily members, provides an inclusive description of the whole zebrafish RNase subfamily.

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.

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.

Conformational changes below the Tm: molecular dynamics studies of the thermal pretransition of ribonuclease A

Recent work suggests that some native conformations of proteins can vary with temperature. To obtain an atomic-level description of this structural and conformational variation, we have performed all-atom, explicit-solvent molecular dynamics simulations of bovine pancreatic ribonuclease A (RNase A) up to its melting temperature (Tm approximately 337 K). RNase A has a thermal pretransition near 320 K [Stelea, S. D., Pancoska, P., Benight, A. S., and Keiderling, T. A. (2001) Protein Sci. 10, 970-978]. Our simulations identify a conformational change that coincides with this pretransition. Between 310 and 320 K, there is a small but significant decrease in the number of native contacts, beta-sheet hydrogen bonding, and deviation of backbone conformation from the starting structure, and an increase in the number of nonnative contacts. Native contacts are lost in beta-sheet regions and in alpha1, partially due to movement of alpha1 away from the beta-sheet core. At 330 and 340 K, a nonnative helical segment of residues 15-20 forms, corresponding to a helix observed in the N-terminal domain-swapped dimer [Liu, Y. S., Hart, P. J., Schulnegger, M. P., and Eisenberg, D. (1998) Proc. Natl. Acad. Sci. U.S.A. 95, 3437-3432]. The conformations observed at the higher temperatures possess nativelike topology and overall conformation, with many native contacts, but they have a disrupted active site. We propose that these conformations may represent the native state at elevated temperature, or the N' state. These simulations show that subtle, functionally important changes in protein conformation can occur below the Tm.

The evolution of seminal ribonuclease: pseudogene reactivation or multiple gene inactivation events?

Two approaches, one novel, are applied to analyze the divergent evolution of ruminant seminal ribonucleases (RNases), paralogs of the well-known pancreatic RNases of mammals. Here, the goal was to identify periods of divergence of seminal RNase under functional constraints, periods of divergence as a pseudogene, and periods of divergence driven by positive selection pressures. The classical approach involves the analysis of nonsynonymous to synonymous replacements ratios (omega) for the branches of the seminal RNase evolutionary tree. The novel approach coupled these analyses with the mapping of substitutions on the folded structure of the protein. These analyses suggest that seminal RNase diverged during much of its history after divergence from pancreatic RNase as a functioning protein, followed by homoplastic inactivations to create pseudogenes in multiple ruminant lineages. Further, they are consistent with adaptive evolution only in the most recent episode leading to the gene in modern oxen. These conclusions contrast sharply with the view, cited widely in the literature, that seminal RNase decayed after its formation by gene duplication into an inactive pseudogene, whose lesions were repaired in a reactivation event. Further, the 2 approaches, omega estimation and mapping of replacements on the protein structure, were compared by examining their utility for establishing the functional status of the seminal RNase genes in 2 deer species. Hog and roe deer share common lesions, which strongly suggests that the gene was inactive in their last common ancestor. In this specific example, the crystallographic approach made the correct implication more strongly than the omega approach. Studies of this type should contribute to an integrated framework of tools to assign functional and nonfunctional episodes to recently created gene duplicates and to understand more broadly how gene duplication leads to the emergence of proteins with novel functions.

Protein flexibility using constraints from molecular dynamics simulations

Proteins are held together in the native state by hydrophobic interactions, hydrogen bonds and interactions with the surrounding water, whose strength as well as spatial and temporal distribution affects protein flexibility and hence function. We study these effects using 10 ns molecular dynamics simulations of pure water and of two proteins, the glutamate receptor ligand binding domain and barnase. We find that most of the noncovalent interactions flicker on and off over typically nanoseconds, and so we can obtain good statistics from the molecular dynamics simulations. Based on this information, a topological network of rigid bonds corresponding to a protein structure with covalent and noncovalent bonds is constructed, with account being taken of the influence of the flickering hydrogen bonds. We define the duty cycle for the noncovalent interactions as the percentage of time a given interaction is present, which we use as an input to investigate flexibility/rigidity patterns, in the algorithm FIRST which constructs and analyses topological networks.

The Importance of Dynamic Effects on the Enzyme Activity

Onconase (ONC), a member of the RNase A superfamily extracted from oocytes of Rana pipiens, is an effective cancer killer. It is currently used in treatment of various forms of cancer. ONC antitumor properties depend on its ribonucleolytic activity that is low in comparison with other members of the superfamily. The most damaging side effect from Onconase treatment is renal toxicity, which seems to be caused by the unusual stability of the enzyme. Therefore, mutants with reduced thermal stability and/or increased catalytic activity may have significant implications for human cancer chemotherapy. In this context, we have determined the crystal structures of two Onconase mutants (M23L-ONC and C87S,des103-104-ONC) and performed molecular dynamic simulations of ONC and C87S,des103-104-ONC with the aim of explaining on structural grounds the modifications of the activity and thermal stability of the mutants. The results also provide the molecular bases to explain the lower catalytic activity of Onconase compared with RNase A and the unusually high thermal stability of the amphibian enzyme.

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

Bovine seminal ribonuclease (BS-RNase) is a covalent homodimeric enzyme homologous to pancreatic ribonuclease (RNase A), endowed with a number of special biological functions. It is isolated as an equilibrium mixture of swapped (MxM) and unswapped (M=M) dimers. The interchanged N termini are hinged on the main bodies through the peptide 16-22, which changes conformation in the two isomers. At variance with other proteins, domain swapping in BS-RNase involves two dimers having a similar and highly constrained quaternary association, mainly dictated by two interchain disulfide bonds. This provides the opportunity to study the intrinsic ability to swap as a function of the hinge sequence, without additional effects arising from dissociation or quaternary structure modifications. Two variants, having Pro19 or the whole sequence of the hinge replaced by the corresponding residues of RNase A, show equilibrium and kinetic parameters of the swapping similar to those of the parent protein. In comparison, the x-ray structures of MxM indicate, within a substantial constancy of the quaternary association, a greater mobility of the hinge residues. The relative insensitivity of the swapping tendency to the substitutions in the hinge region, and in particular to the replacement of Pro19 by Ala, contrasts with the results obtained for other swapped proteins and can be rationalized in terms of the unique features of the seminal enzyme. Moreover, the results indirectly lend credit to the hypothesis that the major role of Pro19 resides in directing the assembly of the non-covalent dimer, the species produced by selective reduction of the interchain disulfides and considered responsible for the special biological functions of BS-RNase.

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.

Structure and Stability of the Non-covalent Swapped Dimer of Bovine Seminal Ribonuclease

A growing number of pancreatic-type ribonucleases (RNases) present cytotoxic activity against malignant cells. The cytoxicity of these enzymes is related to their resistance to the ribonuclease protein inhibitor (RI). In particular, bovine seminal ribonuclease (BS-RNase) is toxic to tumor cells both in vitro and in vivo. BS-RNase is a covalent dimer with two intersubunit disulfide bridges between Cys(31) of one chain and Cys(32) of the second and vice versa. The native enzyme is an equilibrium mixture of two isomers, MxM and M=M. In the former the two subunits swap their N-terminal helices. The cytotoxic action is a peculiar property of MxM. In the reducing environment of cytosol, M=M dissociates into monomers, which are strongly inhibited by RI, whereas MxM remains as a non-covalent dimer (NCD), which evades RI. We have solved the crystal structure of NCD, carboxyamidomethylated at residues Cys(31) and Cys(32) (NCD-CAM), in a complex with 2'-deoxycitidylyl(3'-5')-2'-deoxyadenosine. The molecule reveals a quaternary structural organization much closer to MxM than to other N-terminal-swapped non-covalent dimeric forms of RNases. Model building of the complexes between these non-covalent dimers and RI reveals that NCD-CAM is the only dimer equipped with a quaternary organization capable of interfering seriously with the binding of the inhibitor. Moreover, a detailed comparative structural analysis of the dimers has highlighted the residues, which are mostly important in driving the quaternary structure toward that found in NCD-CAM.