Effect of cavity-modulating mutations on the stability of Escherichia coli ribonuclease HI

Diverse models of cavity engineering in enzyme modification: Creation, filling, and reshaping

Most enzyme modification strategies focus on designing the active sites or their surrounding structures. Interestingly, a large portion of the enzymes (60%) feature active sites located within spacious cavities. Despite recent discoveries, cavity-mediated enzyme engineering remains crucial for enhancing enzyme properties and unraveling folding-unfolding mechanisms. Cavity engineering influences enzyme stability, catalytic activity, specificity, substrate recognition, and docking. This article provides a comprehensive review of various cavity engineering models for enzyme modification, including cavity creation, filling, and reshaping. Additionally, it also discusses feasible tools for geometric analysis, functional assessment, and modification of cavities, and explores potential future research directions in this field. Furthermore, a promising universal modification strategy for cavity engineering that leverages state-of-the-art technologies and methodologies to tailor cavities according to the specific requirements of industrial production conditions is proposed.

Comparisons of Ribonuclease HI Homologs and Mutants Uncover a Multistate Model for Substrate Recognition

Ribonuclease HI (RNHI) nonspecifically cleaves the RNA strand of RNA:DNA hybrid duplexes in a myriad of biological processes. Several RNHI homologs contain an extended domain, termed the handle region, which is critical to substrate binding. Nuclear magnetic resonance (NMR) spectroscopy and molecular dynamics (MD) simulations have suggested a kinetic model in which the handle region can exist in open (substrate-binding competent) or closed (substrate-binding incompetent) states in homologs containing arginine or lysine at position 88 (using sequence numbering of E. coli RNHI), while the handle region populates states intermediate between the open and closed conformers in homologs with asparagine at residue 88 [Stafford, K. A., et al., PLoS Comput. Biol. 2013, 9, 1-10]. NMR parameters characterizing handle region dynamics are highly correlated with enzymatic activity for RNHI homologs with two-state (open/closed) handle regions [Martin, J. A., et al., Biochemistry 2020, 59, 3201-3205]. The work presented herein shows that homologs containing asparagine 88 display distinct structural features compared with their counterparts containing arginine or lysine 88. Comparisons of RNHI homologs and site-directed mutants with asparagine 88 support a kinetic model for handle region dynamics that includes 12 unique transitions between eight conformations. Overall, these findings present an example of the structure-function relationships of enzymes and spotlight the use of NMR spectroscopy and MD simulations in uncovering fine details of conformational preferences.

Simple models for nonpolar solvation: Parameterization and testing

Implicit solvent models are important for many biomolecular simulations. The polarity of aqueous solvent is essential and qualitatively captured by continuum electrostatics methods like Generalized Born (GB). However, GB does not account for the solvent-induced interactions between exposed hydrophobic sidechains or solute-solvent dispersion interactions. These "nonpolar" effects are often modeled through surface area (SA) energy terms, which lack realism, create mathematical singularities, and have a many-body character. We have explored an alternate, Lazaridis-Karplus (LK) gaussian energy density for nonpolar effects and a dispersion (DI) energy term proposed earlier, associated with GB electrostatics. We parameterized several combinations of GB, SA, LK, and DI energy terms, to reproduce 62 small molecule solvation free energies, 387 protein stability changes due to point mutations, and the structures of 8 protein loops. With optimized parameters, the models all gave similar results, with GBLK and GBDILK giving no performance loss compared to GBSA, and mean errors of 1.7 kcal/mol for the stability changes and 2 Å deviations for the loop conformations. The optimized GBLK model gave poor results in MD of the Trpcage mini-protein, but parameters optimized specifically for MD performed well for Trpcage and three other small proteins. Overall, the LK and DI nonpolar terms are valid alternatives to SA treatments for a range of applications. © 2017 Wiley Periodicals, Inc.

Conformational preferences underlying reduced activity of a thermophilic ribonuclease H

The conformational basis for reduced activity of the thermophilic ribonuclease HI enzyme from Thermus thermophilus, compared to its mesophilic homolog from Escherichia coli, is elucidated using a combination of NMR spectroscopy and molecular dynamics (MD) simulations. Explicit-solvent all-atom MD simulations of the two wild-type proteins and an E. coli mutant in which a glycine residue is inserted after position 80 to mimic the T. thermophilus protein reproduce the differences in conformational dynamics determined from (15)N spin-relaxation NMR spectroscopy of three loop regions that surround the active site and contain functionally important residues: the glycine-rich region, the handle region, and the β5/αE loop. Examination of the MD trajectories indicates that the thermophilic protein samples conformations productive for substrate binding and activity less frequently than the mesophilic enzyme, although these differences may manifest as either increased or decreased relative flexibility of the different regions. Additional MD simulations indicate that mutations increasing activity of the T. thermophilus enzyme at mesophilic temperatures do so by reconfiguring the local environments of the mutated sites to more closely resemble active conformations. Taken together, the results show that both locally increased and decreased flexibility contribute to an overall reduction in activity of T. thermophilus ribonuclease H compared to its mesophilic E. coli homolog.

Inferring stabilizing mutations from protein phylogenies: application to influenza hemagglutinin

One selection pressure shaping sequence evolution is the requirement that a protein fold with sufficient stability to perform its biological functions. We present a conceptual framework that explains how this requirement causes the probability that a particular amino acid mutation is fixed during evolution to depend on its effect on protein stability. We mathematically formalize this framework to develop a Bayesian approach for inferring the stability effects of individual mutations from homologous protein sequences of known phylogeny. This approach is able to predict published experimentally measured mutational stability effects (ΔΔG values) with an accuracy that exceeds both a state-of-the-art physicochemical modeling program and the sequence-based consensus approach. As a further test, we use our phylogenetic inference approach to predict stabilizing mutations to influenza hemagglutinin. We introduce these mutations into a temperature-sensitive influenza virus with a defect in its hemagglutinin gene and experimentally demonstrate that some of the mutations allow the virus to grow at higher temperatures. Our work therefore describes a powerful new approach for predicting stabilizing mutations that can be successfully applied even to large, complex proteins such as hemagglutinin. This approach also makes a mathematical link between phylogenetics and experimentally measurable protein properties, potentially paving the way for more accurate analyses of molecular evolution.

Mutating a protein frequently causes a change in its stability. As scientists, we often care about these changes because we would like to engineer a protein's stability or understand how its stability is impacted by a naturally occurring mutation. Evolution also cares about mutational stability changes, because a basic evolutionary requirement is that proteins remain sufficiently stable to perform their biological functions. Our work is based on the idea that it should be possible to use the fact that evolution selects for stability to infer from related proteins the effects of specific mutations. We show that we can indeed use protein evolutionary histories to computationally predict previously measured mutational stability changes more accurately than methods based on either of the two main existing strategies. We then test whether we can predict mutations that increase the stability of hemagglutinin, an influenza protein whose rapid evolution is partly responsible for the ability of this virus to cause yearly epidemics. We experimentally create viruses carrying predicted stabilizing mutations and find that several do in fact improve the virus's ability to grow at higher temperatures. Our computational approach may therefore be of use in understanding the evolution of this medically important virus.

SIMPLE estimate of the free energy change due to aliphatic mutations: Superior predictions based on first principles

The bioinformatics revolution of the last decade has been instrumental in the development of empirical potentials to quantitatively estimate protein interactions for modeling and design. Although computationally efficient, these potentials hide most of the relevant thermodynamics in 5-to-40 parameters that are fitted against a large experimental database. Here, we revisit this longstanding problem and show that a careful consideration of the change in hydrophobicity, electrostatics, and configurational entropy between the folded and unfolded state of aliphatic point mutations predicts 20-30% less false positives and yields more accurate predictions than any published empirical energy function. This significant improvement is achieved with essentially no free parameters, validating past theoretical and experimental efforts to understand the thermodynamics of protein folding. Our first principle analysis strongly suggests that both the solute-solute van der Waals interactions in the folded state and the electrostatics free energy change of exposed aliphatic mutations are almost completely compensated by similar interactions operating in the unfolded ensemble. Not surprisingly, the problem of properly accounting for the solvent contribution to the free energy of polar and charged group mutations, as well as of mutations that disrupt the protein backbone remains open.

Identification of RNase HII from psychrotrophic bacterium, Shewanella sp. SIB1 as a high-activity type RNase H

The gene encoding RNase HII from the psychrotrophic bacterium, Shewanella sp. SIB1 was cloned, overexpressed in Escherichia coli, and the recombinant protein was purified and biochemically characterized. SIB1 RNase HII is a monomeric protein with 212 amino acid residues and shows an amino acid sequence identity of 64% to E. coli RNase HII. The enzymatic properties of SIB1 RNase HII, such as metal ion preference, pH optimum, and cleavage mode of substrate, were similar to those of E. coli RNase HII. SIB1 RNase HII was less stable than E. coli RNase HII, but the difference was marginal. The half-lives of SIB1 and E. coli RNases HII at 30 degrees C were approximately 30 and 45 min, respectively. The midpoint of the urea denaturation curve and optimum temperature of SIB1 RNase HII were lower than those of E. coli RNase HII by approximately 0.2 M and approximately 5 degrees C, respectively. However, SIB1 RNase HII was much more active than E. coli RNase HII at all temperatures studied. The specific activity of SIB1 RNase HII at 30 degrees C was 20 times that of E. coli RNase HII. Because SIB1 RNase HII was also much more active than SIB1 RNase HI, RNases HI and HII represent low- and high-activity type RNases H, respectively, in SIB1. In contrast, RNases HI and HII represent high- and low-activity type RNases H, respectively, in E. coli. We propose that bacterial cells usually contain low- and high-activity type RNases H, but these types are not correlated with RNase H families.

Stabilization of E. coli Ribonuclease HI by the ‘stability profile of mutant protein’ (SPMP)-inspired random and non-random mutagenesis

The change in the structural stability of Escherichia coli ribonuclease HI (RNase HI) due to single amino acid substitutions has been estimated computationally by the stability profile of mutant protein (SPMP) [Ota, M., Kanaya, S. Nishikawa, K., 1995. Desk-top analysis of the structural stability of various point mutations introduced into ribonuclease H. J. Mol. Biol. 248, 733-738]. As well, an effective strategy using random mutagenesis and genetic selection has been developed to obtain E. coli RNase HI mutants with enhanced thermostability [Haruki, M., Noguchi, E., Akasako, A., Oobatake, M., Itaya, M., Kanaya, S., 1994. A novel strategy for stabilization of Escherichia coli ribonuclease HI involving a screen for an intragenic suppressor of carboxyl-terminal deletions. J. Biol. Chem. 269, 26904-26911]. In this study, both methods were combined: random mutations were individually introduced to Lys99-Val101 on the N-terminus of the alpha-helix IV and the preceding beta-turn, where substitutions of other amino acid residues were expected to significantly increase the stability from SPMP, and then followed by genetic selection. Val101 to Ala, Gln, and Arg mutations were selected by genetic selection. The Val101-->Ala mutation increased the thermal stability of E. coli RNase HI by 2.0 degrees C in Tm at pH 5.5, whereas the Val101-->Gln and Val101-->Arg mutations decreased the thermostability. Separately, the Lys99-->Pro and Asn100-->Gly mutations were also introduced directly. The Lys99-->Pro mutation increased the thermostability of E. coli RNase HI by 1.8 degrees C in Tm at pH 5.5, whereas the Asn100-->Gly mutation decreased the thermostability by 17 degrees C. In addition, the Lys99-->Pro mutation altered the dependence of the enzymatic activity on divalent metal ions.

Multiple time scale backbone dynamics of homologous thermophilic and mesophilic ribonuclease HI enzymes

Backbone conformational fluctuations on multiple time scales in a cysteine-free Thermus thermophilus ribonuclease HI mutant (ttRNH(*)) are quantified using (15)N nuclear magnetic spin relaxation. Laboratory-frame relaxation data acquired at 310 K and at static magnetic field strengths of 11.7, 14.1 and 18.8 T are analysed using reduced spectral density mapping and model-free approaches. Chemical exchange line broadening is characterized using Hahn-echo transverse and multiple quantum relaxation data acquired over a temperature range of 290-320 K and at a static magnetic field strength of 14.1 T. Results for ttRNH(*) are compared to previously published data for a mesophilic homologue, Escherichia coli ribonuclease HI (ecRNH). Intramolecular conformational fluctuations on the picosecond-to-nanosecond time scale generally are similar for ttRNH(*) and ecRNH. beta-Strands 3 and 5 and the glycine-rich region are more rigid while the substrate-binding handle region and C-terminal tail are more flexible in ttRNH(*) than in ecRNH. Rigidity in the two beta-strands and the glycine-rich region, located along the periphery of the central beta-sheet, may be associated with the increased thermodynamic stability of the thermophilic enzyme. Chemical exchange line broadening, reflecting microsecond-to-millisecond time scale conformational changes, is more pronounced in ttRNH(*) than in ecRNH, particularly for residues in the handle and surrounding the catalytic site. The temperature dependence of chemical exchange show an increase of approximately 15 kJ/mol in the apparent activation energies for ttRNH(*) residues in the handle compared to ecRNH. Increased activation barriers, coupled with motion between alpha-helices B and C not present in ecRNH, may be associated with the reduced catalytic activity of the thermophilic enzyme at 310 K.

Contribution of engineered electrostatic interactions to the stability of cytosolic malate dehydrogenase

Protein engineering is a promising tool to obtain stable proteins. Comparison between homologous thermophilic and mesophilic enzymes from a given structural family can reveal structural features responsible for the enhanced stability of thermophilic proteins. Structures from pig heart cytosolic and Thermus flavus malate dehydrogenases (cMDH, Tf MDH), two proteins showing a 55% sequence homology, were compared with the aim of increasing cMDH stability using features from the Thermus flavus enzyme. Three potential salt bridges from Tf MDH were selected on the basis of their location in the protein (surface R176-D200, inter-subunit E57-K168 and intrasubunit R149-E275) and implemented on cMDH using site-directed mutagenesis. Mutants containing E275 were not produced in any detectable amount, which shows that the energy penalty of introducing a charge imbalance in a region that was not exposed to solvent was too unfavourable to allow proper folding of the protein. The salt bridge R149-E275, if formed, would not enhance stability enough to overcome this effect. The remaining mutants were expressed and active and no differences from wild-type other than stability were found. Of the mutants assayed, Q57E/L168K led to a stability increase of 0.4 kcal/mol, as determined by either guanidinium chloride denaturalization or thermal inactivation experiments. This results in a 15 degrees C shift in the optimal temperature, thus confirming that the inter-subunit salt bridge initially present in the T.flavus enzyme was formed in the cMDH structure and that the extra energy obtained is transformed into an increase in protein stability. These results indicate that the use of structural features of thermophilic enzymes, revealed by a detailed comparison of three-dimensional structures, is a valid strategy to improve the stability of mesophilic malate dehydrogenases.

Pressure-denatured state of Escherichia coli ribonuclease HI as monitored by Fourier transform infrared and NMR spectroscopy

Pressure denaturation of Escherichia coli ribonuclease HI (RNase HI) was studied by Fourier transform infrared (FTIR) and two-dimensional NMR spectroscopy at pD* 3.0 and 25 degrees C. A reversible transition in the pressure range of 0.1-1090 MPa was observed with second-derivative FTIR experiments. A cooperative and gradual denaturation, involving both the secondary and tertiary structures, was observed between 240 and 450 MPa. The two peaks at 1629 and 1652 cm(-1), due to beta-strands and alpha-helices, respectively, did not fully disappear after the denaturation, and are different from the spectra of the random coil peptides. The hydrogen-deuterium exchange rates of the individual backbone amide protons were determined by heteronuclear NMR combined with the pressure-jump technique at 500, 650, and 850 MPa. Although most of the amides protected in the native structure are also highly protected in the pressure-denatured state, the rate constants (0.048 +/- 0.007 min(-1)) for the amide protons at 500 MPa are similar regardless of their locations, which is an indication of the EX1 mechanism of hydrogen-deuterium exchange. The pressure-denatured state of RNase HI at 500 MPa represents a novel denatured state, which is different from a typical molten globule state at atmospheric pressure (0.1 MPa), from the viewpoint of the homogeneous rate constants. The observations at 650 MPa are essentially the same as those at 500 MPa. However, at 850 MPa, the amide exchange rates for the highly hydrophobic C-terminal half of alpha-helix I are significantly slower than those for the other part of the protein, which can be interpreted as a hydrophobic collapse centered at the C-terminal half of alpha-helix I.

Kinetic and stoichiometric analysis for the binding of Escherichia coli ribonuclease HI to RNA-DNA hybrids using surface plasmon resonance

To understand how ribonucleases H recognize RNA-DNA hybrid substrates, we analyzed kinetic parameters of binding of Escherichia coli RNase HI to RNA-DNA hybrids ranging in length from 18 to 36 base pairs (bp) using surface plasmon resonance (BIAcoreTM). The kon and koff values for the binding of the enzyme to the 36-bp substrate were 1.5 x 10(6) M-1 s-1 and 3.2 x 10(-2) s-1, respectively. Similar values were obtained with the shorter substrates. Using uncleavable 2'-O-methylated RNA-DNA substrates, values for kon and koff were 2.1 x 10(5) M-1 s-1 and 1.3 x 10(-1) s-1 in the absence of Mg2+ that were further reduced in the presence of Mg2+ to 7.4 x 10(3) M-1 s-1 and 2.6 x 10(-2) s-1. Kinetic parameters similar to the wild-type enzyme were obtained using an active-site mutant enzyme, Asp134 replaced by Ala, whereas a greatly reduced on-rate was observed for another inactive mutant enzyme, in which the basic protrusion is eliminated, thereby distinguishing between poor catalysis and inability to bind to the substrate. Stoichiometric analyses of RNase HI binding to substrates of 18, 24, 30, and 36 bp are consistent with previous reports suggesting that RNase HI binds to 9-10 bp of RNA-DNA hybrid.

Conformational stabilities of Escherichia coli RNase HI variants with a series of amino acid substitutions at a cavity within the hydrophobic core

Escherichia coli ribonuclease HI has a cavity within the hydrophobic core. Two core residues, Ala52 and Val74, resided at both ends of this cavity. We have constructed a series of single mutant proteins at Ala52, and double mutant proteins, in which Ala52 was replaced by Gly, Val, Ile, Leu, or Phe, and Val74 was replaced by Ala or Leu. All of these mutant proteins, except for A52W, A52R, and A52G/V74A, were overproduced and purified. Measurement of the thermal denaturations of the proteins at pH 3.2 by CD suggests that the cavity is large enough to accommodate three methyl or methylene groups without creating serious strains. A correlation was observed between the protein stability and the hydrophobicity of the substituted residue. As a result, a number of the mutant proteins were more stable than the wild-type protein. The stabilities of the mutant proteins with charged or extremely bulky residues at the cavity were lower than those expected from the hydrophobicities of the substituted residues, suggesting that considerable strains are created at the mutation sites in these mutant proteins. However, examination of the far- and near-UV CD spectra and the enzymatic activities suggest that all of the mutant proteins have structures similar to that of the wild-type protein. These results suggest that the cavity in the hydrophobic core of E. coli RNase HI is conformationally fairly stable. This may be the reason why the cavity-filling mutations effectively increase the thermal stability of this protein.

Engineering proteins for environmental applications

Recently, significant new insight has been obtained into the structure and catalytic mechanism of enzymes that convert environmental pollutants. Recent advances in protein engineering make it possible to use this information for improving the catalytic performance of such enzymes to achieve increased stability and expanded substrate range.

Investigating the role of conserved residue Asp134 in Escherichia coli ribonuclease HI by site-directed random mutagenesis

The role of the conserved Asp134 residue in Escherichia coli ribonuclease HI, which is located at the center of the alpha V helix and lies close to the active site, was analyzed by means of site-directed random mutagenesis. Mutant rnhA genes encoding proteins with ribonuclease H activities were screened by their ability to suppress the ribonuclease-H-dependent, temperature-sensitive growth phenotype of E. coli strain MIC3001. Based on the DNA sequences, nine mutant proteins were predicted to have ribonuclease H activity in vivo. All of these mutant proteins were purified to homogeneity and examined for enzymic activity and protein stability. Among them, only the mutant proteins [D134H]RNase H and [D134N]RNase H were shown to have considerable ribonuclease H activities. Determination of the kinetic parameters revealed that replacement of Asp134 by amino acid residues other than asparagine and histidine dramatically decreased the enzymic activity without seriously affecting the substrate binding. Determination of the CD spectra indicated that none of the mutations seriously affected secondary and tertiary structure. The protein stability was determined from the thermal denaturation curves. All mutant proteins were more stable than the wild-type protein. Such stabilization effects would be a result of a reduction in the negative charge repulsion between Asp134 and the active-site residues, and/or an enhancement of the stability of the alpha V helix. These results strongly suggest that Asp134 does not contribute to the maintenance of the molecular architecture but the carboxyl oxygen at its delta 1 position impacts catalysis.

Stabilization of creatinase from Pseudomonas putida by random mutagenesis

Creatinase (creatine amidinohydrolase, EC 3.5.3.3) from Pseudomonas putida is a homodimer of 45 kDa subunit molecular mass, the three-dimensional structure of which is known at 1.9 A resolution. Three point mutants, A109V, V355M, and V182I, as well as one double mutant combining A109V and V355M, and the triple mutant with all three replacements, were compared with wild-type creatinase regarding their physical and enzymological properties. High-resolution crystal data for wild-type creatinase and the first two mutants suggest isomorphism at least for these three proteins (R. Huber, pers. comm.). Physicochemical measurements confirm this prediction, showing that the mutations have no effect either on the quaternary structure and gross conformation or the catalytic properties as compared to wild-type creatinase. The replacement of V182 (at the solvent-exposed end of the first helix of the C-terminal domain) does not cause significant differences in comparison with the wild-type enzyme. The other point mutations stabilize the first step in the biphasic denaturation transition without affecting the second one. In sum, the enhanced stability seems to reflect slight improvements in the local packing without creating new well-defined bonds. The increase in hydrophobicity generated by the introduction of additional methyl groups (A109V, V182I) must be compensated by minor readjustments of the global structure. Secondary or quaternary interactions are not affected. In going from single to double and triple mutants, to a first approximation, the increments of stabilization are additive.

Prediction and analysis of structure, stability and unfolding of thermolysin-like proteases

Bacillus neutral proteases (NPs) form a group of well-characterized homologous enzymes, that exhibit large differences in thermostability. The three-dimensional (3D) structures of several of these enzymes have been modelled on the basis of the crystal structures of the NPs of B. thermoproteolyticus (thermolysin) and B. cereus. Several new techniques have been developed to improve the model-building procedures. Also a 'model-building by mutagenesis' strategy was used, in which mutants were designed just to shed light on parts of the structures that were particularly hard to model. The NP models have been used for the prediction of site-directed mutations aimed at improving the thermostability of the enzymes. Predictions were made using several novel computational techniques, such as position-specific rotamer searching, packing quality analysis and property-profile database searches. Many stabilizing mutations were predicted and produced: improvement of hydrogen bonding, exclusion of buried water molecules, capping helices, improvement of hydrophobic interactions and entropic stabilization have been applied successfully. At elevated temperatures NPs are irreversibly inactivated as a result of autolysis. It has been shown that this denaturation process is independent of the protease activity and concentration and that the inactivation follows first-order kinetics. From this it has been conjectured that local unfolding of (surface) loops, which renders the protein susceptible to autolysis, is the rate-limiting step. Despite the particular nature of the thermal denaturation process, normal rules for protein stability can be applied to NPs. However, rather than stabilizing the whole protein against global unfolding, only a small region has to be protected against local unfolding. In contrast to proteins in general, mutational effects in proteases are not additive and their magnitude is strongly dependent on the location of the mutation. Mutations that alter the stability of the NP by a large amount are located in a relatively weak region (or more precisely, they affect a local unfolding pathway with a relatively low free energy of activation). One weak region, that is supposedly important in the early steps of NP unfolding, has been determined in the NP of B. stearothermophilus. After eliminating this weakest link a drastic increase in thermostability was observed and the search for the second-weakest link, or the second-lowest energy local unfolding pathway is now in progress. Hopefully, this approach can be used to unravel the entire early phase of unfolding.

Mechanism-based strategies for protein thermostabilization

Strategies for stabilizing enzymes can be derived from a two-step model of irreversible inactivation that involves preliminary reversible unfolding, followed by an irreversible step. Reversible unfolding is best prevented by covalent immobilization, whereas methods such as covalent modification of amino acid residues or 'medium engineering' (by the addition of low-molecular-weight compounds) are effective against irreversible 'incorrect' refolding. Genetic modification of the protein sequence is the most effective approach for preventing chemical deterioration.