Destabilizing mutations alter the hydrogen exchange mechanism in ribonuclease A

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

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

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.

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.

Structural and kinetic mapping of side-chain exposure onto the protein energy landscape

Identification and characterization of structural fluctuations that occur under native conditions is crucial for understanding protein folding and function, but such fluctuations are often rare and transient, making them difficult to study. Native-state hydrogen exchange (NSHX) has been a powerful tool for identifying such rarely populated conformations, but it generally reveals no information about the placement of these species along the folding reaction coordinate or the barriers separating them from the folded state and provides little insight into side-chain packing. To complement such studies, we have performed native-state alkyl-proton exchange, a method analogous to NSHX that monitors cysteine modification rather than backbone amide exchange, to examine the folding landscape of Escherichia coli ribonuclease H, a protein well characterized by hydrogen exchange. We have chosen experimental conditions such that the rate-limiting barrier acts as a kinetic partition: residues that become exposed only upon crossing the unfolding barrier are modified in the EX1 regime (alkylation rates report on the rate of unfolding), while those exposed on the native side of the barrier are modified predominantly in the EX2 regime (alkylation rates report on equilibrium populations). This kinetic partitioning allows for identification and placement of partially unfolded forms along the reaction coordinate. Using this approach we detect previously unidentified, rarely populated conformations residing on the native side of the barrier and identify side chains that are modified only upon crossing the unfolding barrier. Thus, in a single experiment under native conditions, both sides of the rate-limiting barrier are investigated.

X-ray crystallographic studies of RNase A variants engineered at the most destabilizing positions of the main hydrophobic core: further insight into protein stability

To investigate the structural origin of decreased pressure and temperature stability, the crystal structure of bovine pancreatic ribonuclease A variants V47A, V54A, V57A, I81A, I106A, and V108A was solved at 1.4-2.0 A resolution and compared with the structure of wild-type protein. The introduced mutations had only minor influence on the global structure of ribonuclease A. The structural changes had individual character that depends on the localization of mutated residue, however, they seemed to expand from mutation site to the rest of the structure. Several different parameters have been evaluated to find correlation with decrease of free energy of unfolding DeltaDeltaG(T), and the most significant correlation was found for main cavity volume change. Analysis of the difference distance matrices revealed that the ribonuclease A molecule is organized into five relatively rigid subdomains with individual response to mutation. This behavior consistent with results of unfolding experiments is an intrinsic feature of ribonuclease A that might be surviving remnants of folding intermediates and reflects the dynamic nature of the molecule.

Mapping the stability clusters in bovine pancreatic ribonuclease A

In the present work, we have thermodynamically characterized the thermally induced unfolding of 20 variants of bovine pancreatic ribonuclease A (RNase A) to experimentally describe the residues and the regions that are critical for the stability of the enzyme. The achieved results, complemented with previous studies by our group, allowed us to define the significance of the two hydrophobic nuclei present in the RNase A structure, as well as the contribution of the participating residues within each nucleus, to the global enzyme stability. We propose a structural model for the major and the minor hydrophobic nuclei of RNase A. The major nucleus is composite and located in the cavity delimited by alpha-helices 1 and 3, and the beta-sheet that is formed by strands 2, 3, 5, and 6. It consists of a central tight packed part constituted by residues Phe8, Met13, Val54, Val57, Ile106, Val108, and Phe120. This central part is surrounded by a layer formed by residues Val63, Tyr73, Met79, Ile107, Val116, and Val118. The minor nucleus, although less complex, is also constituted by a tight packing that involves the side chains of residues Tyr25, Met29, Met30, Leu35, Phe46, and Tyr97, which fill the cavity that originates the beta-sheet formed by beta-strands 1, 4, and 5 together with alpha-helix2.

The (1)H, (13)C, (15)N resonance assignment, solution structure, and residue level stability of eosinophil cationic protein/RNase 3 determined by NMR spectroscopy

Eosinophil cationic protein (ECP)/human RNase 3, a member of the RNase A family, is a remarkably cytotoxic protein implicated in asthma and allergies. These activities are probably due to ECP's ability to interact with and disrupt membranes and depend on two Trp, 19 Arg, and possibly an extremely high conformational stability. Here, we have used NMR spectroscopy to assign essentially all (1)H, (15)N, and backbone (13)C resonances, to solve the 3D structure in aqueous solution and to quantify its residue-level stability. The NMR solution structure was determined on the basis of 2316 distance constraints and is well-defined (backbone RMSD = 0.81 A). The N-terminus and the loop composed of residues 114-123 are relatively well-ordered; in contrast, conformational diversity is observed for the loop segments 17-22, 65-68, and 92-95 and most exposed sidechains. The side chain NH groups of the two Trp and 19 Arg showed no significant protection against hydrogen/deuterium exchange. The most protected NH groups belong to the first and last two beta-strands, and curiously, the first alpha-helix. Analysis of their exchange rates reveals a strikingly high global stability of 11.8 kcal/mol. This value and other stability measurements are used to better quantify ECP's unfolding thermodynamics.