Utilization of low-stability variants in protein evolutionary engineering.

Evolutionary engineering involves repeated mutations and screening and is widely used to modify protein functions. However, it is important to diversify evolutionary pathways to eliminate the bias and limitations of the variants by using traditionally unselected variants. In this study, we focused on low-stability variants that are commonly excluded from evolutionary processes and tested a method that included an additional restabilization step. The esterase from the thermophilic bacterium Alicyclobacillus acidocaldarius was used as a model protein, and its activity at its optimum temperature of 65 °C was improved by evolutionary experiments using random mutations by error-prone PCR. After restabilization using low-stability variants with low-temperature (37 °C) activity, several re-stabilizing variants were obtained from a large number of variant libraries. Some of the restabilized variants achieved by removing the destabilizing mutations showed higher activity than that of the wild-type protein. This implies that low-stability variants with low-temperature activity can be re-evolved for future use. This method will enable further diversification of evolutionary pathways.

Survey on Actual Management of Osteoporosis with the Japanese Medical Data Vision Database in Elderly Patients Undergoing Spinal Fusion.

Background: No actual data on spinal fusion and management of osteoporosis in Japan have been reported. The aim of the survey was to investigate pre- and post-operative management of osteoporosis, including testing and prescription, in elderly patients undergoing spinal fusion in Japan. Methods: Medical data on patients aged 65 years or older undergoing spinal fusion from April 2018 to March 2022 were extracted from the medical data vision (MDV) database containing health insurance claims data from Japanese acute care hospitals to investigate fusion area, pre- and post-operative osteoporosis tests (bone mineral density and osteoporosis markers), prescriptions of osteoporosis medications, and other information. Results: The analysis set consisted of 26,959 patients. Annual pre-operative BMD testing rates and osteoporosis markers testing rates were higher than the post-operative rates without significant annual changes. The post-operative prescription rate of osteoporosis medications throughout the target period was approximately two times higher than the preoperative rate. The drug with highest pre- and post-operative prescription rates was teriparatide (TPTD) followed by bisphosphonates, showing that the prescription rate of TPTD proportionally increased with the length of fusion area. Conclusions: It was suggested that patients aged 65 years or older undergoing spinal fusion might receive insufficient osteoporosis tests. Despite no trend in the testing rate with the length of fusion area, some tendency was observed in the selection of osteoporosis medications. In patients with osteoporosis undergoing spinal fusion, early examination, diagnosis, and therapeutic intervention may improve the prognoses, and solid testing and prescriptions are therefore expected.

Engineering a monobody specific to monomeric Cu/Zn-superoxide dismutase associated with amyotrophic lateral sclerosis.

Misfolding of mutant Cu/Zn-superoxide dismutase (SOD1) has been implicated in familial form of amyotrophic lateral sclerosis (ALS). A natively folded SOD1 forms a tight homodimer, and the dimer dissociation has been proposed to trigger the oligomerization/aggregation of SOD1. Besides increasing demand for probes allowing the detection of monomerized forms of SOD1 in various applications, the development of probes has been limited to conventional antibodies. Here, we have developed Mb(S4) monobody, a small synthetic binding protein based on the fibronectin type III scaffold, that recognizes a monomeric but not dimeric form of SOD1 by performing combinatorial library selections using phage and yeast-surface display methods. Although Mb(S4) was characterized by its excellent selectivity to the monomeric conformation of SOD1, the monomeric SOD1/Mb(S4) complex was not so stable (apparent Kd ~ µM) as to be detected in conventional pull-down experiments. Instead, the complex of Mb(S4) with monomeric but not dimeric SOD1 was successfully trapped by proximity-enabled chemical crosslinking even when reacted in the cell lysates. We thus anticipate that Mb(S4) binding followed by chemical crosslinking would be a useful strategy for in vitro and also ex vivo detection of the monomeric SOD1 proteins.

Induced Circular Dichroism Analysis of Thermally Induced Conformational Changes on Protein Binding Sites Under a Crowding Environment.

Protein-ligand interactions in crowded cellular environments play a crucial role in biological functions. The crowded environment can perturb the overall protein structure and local conformation, thereby influencing the binding pathway of protein-ligand reactions within the cellular milieu. Therefore, a detailed understanding of the local conformation is crucial for elucidating the intricacies of protein-ligand interactions in crowded cellular environments. In this study, we investigated the feasibility of induced circular dichroism (ICD) using 8-anilinonaphthalene-1-sulfonic acid (ANS) for local conformational analysis at the binding site in a crowding environment. Bovine serum albumin (BSA) concentration-dependent measurements were performed to assess the feasibility of ANS-ICD for analyzing protein interior binding sites. The results showed distinct changes in the ANS-ICD spectra of BSA solutions, indicating their potential for analyzing the internal conformation of proteins. Moreover, temperature-dependent measurements were performed in dilute and crowding environments, revealing distinct denaturation pathways of BSA binding sites. Principal component analysis of ANS-ICD spectral changes revealed lower temperature pre-denaturation in the crowded solution than that in the diluted solution, suggesting destabilization of binding sites owing to self-crowding repulsive interactions. The established ANS-ICD method can provide valuable conformational insights into protein-ligand interactions in crowded cellular environments.

Conformation state-specific monobodies regulate the functions of flexible proteins through conformation trapping.

Synthetic binding proteins have emerged as modulators of protein functions through protein-protein interactions (PPIs). Because PPIs are influenced by the structural dynamics of targeted proteins, investigating whether the synthetic-binders-based strategy is applicable for proteins with large conformational changes is important. This study demonstrates the applicability of monobodies (fibronectin type-III domain-based synthetic binding proteins) in regulating the functions of proteins that undergo tens-of-angstroms-scale conformational changes, using an example of the A55C/C77S/V169C triple mutant (Adktm ; a phosphoryl transfer-catalyzing enzyme with a conformational change between OPEN/CLOSED forms). Phage display successfully developed monobodies that recognize the OPEN form (substrate-unbound form), but not the CLOSED form of Adktm . Two OPEN form-specific clones (OP-2 and OP-4) inhibited Adktm kinase activity. Epitope mapping with a yeast-surface display/flow cytometry indicated that OP-2 binds to the substrate-entry side of Adktm , whereas OP-4 binding occurs at another site. Small angle X-ray scattering  coupled with size-exclusion chromatography (SEC-SAXS) indicated that OP-4 binds to the hinge side opposite to the substrate-binding site of Adktm , retaining the whole OPEN-form structure of Adktm . Titration of the OP-4-Adktm complex with Ap5 A, a transition-state analog of Adktm , showed that the conformational shift to the CLOSED form was suppressed although Adktm retained the OPEN-form (i.e., substrate-binding ready form). These results show that OP-4 captures and stabilizes the OPEN-form state, thereby affecting the hinge motion. These experimental results indicate that monobody-based modulators can regulate the functions of proteins that show tens-of-angstroms-scale conformational changes, by trapping specific conformational states generated during large conformational change process that is essential for function exertion.

Structures of a FtsZ single protofilament and a double-helical tube in complex with a monobody.

FtsZ polymerizes into protofilaments to form the Z-ring that acts as a scaffold for accessory proteins during cell division. Structures of FtsZ have been previously solved, but detailed mechanistic insights are lacking. Here, we determine the cryoEM structure of a single protofilament of FtsZ from Klebsiella pneumoniae (KpFtsZ) in a polymerization-preferred conformation. We also develop a monobody (Mb) that binds to KpFtsZ and FtsZ from Escherichia coli without affecting their GTPase activity. Crystal structures of the FtsZ-Mb complexes reveal the Mb binding mode, while addition of Mb in vivo inhibits cell division. A cryoEM structure of a double-helical tube of KpFtsZ-Mb at 2.7 Å resolution shows two parallel protofilaments. Our present study highlights the physiological roles of the conformational changes of FtsZ in treadmilling that regulate cell division.

An efficient selenium transport pathway of selenoprotein P utilizing a high-affinity ApoER2 receptor variant and being independent of selenocysteine lyase.

Selenoprotein P (SeP, encoded by the SELENOP gene) is a plasma protein that contains selenium in the form of selenocysteine residues (Sec, a cysteine analog containing selenium instead of sulfur). SeP functions for the transport of selenium to specific tissues in a receptor-dependent manner. Apolipoprotein E receptor 2 (ApoER2) has been identified as a SeP receptor. However, diverse variants of ApoER2 have been reported, and the details of its tissue specificity and the molecular mechanism of its efficiency remain unclear. In the present study, we found that human T lymphoma Jurkat cells have a high ability to utilize selenium via SeP, while this ability was low in human rhabdomyosarcoma cells. We identified an ApoER2 variant with a high affinity for SeP in Jurkat cells. This variant had a dissociation constant value of 0.67 nM and a highly glycosylated O-linked sugar domain. Moreover, the acidification of intracellular vesicles was necessary for selenium transport via SeP in both cell types. In rhabdomyosarcoma cells, SeP underwent proteolytic degradation in lysosomes and transported selenium in a Sec lyase-dependent manner. However, in Jurkat cells, SeP transported selenium in Sec lyase-independent manner. These findings indicate a preferential selenium transport pathway involving SeP and high-affinity ApoER2 in a Sec lyase-independent manner. Herein, we provide a novel dynamic transport pathway for selenium via SeP.

Sleep duration and food intake in people with type 2 diabetes mellitus and factors affecting confectionery intake.

AIMS/INTRODUCTION: We carried out a cross-sectional study of people with type 2 diabetes mellitus to elucidate the association between sleep duration and food intake. MATERIALS AND METHODS: Overall, 2,887 participants with type 2 diabetes mellitus (mean age 63.0 years; 61.1% men; mean glycated hemoglobin level 7.5%) were included in this study. The participants' self-reported dietary habits and sleep duration were evaluated using a brief self-administered dietary history questionnaire and Pittsburgh Sleep Quality Index, respectively. The participants were categorized into the following four groups based on sleep duration: <6, 6-6.9, 7-7.9 (reference) and ≥8 h. RESULTS: No significant differences were observed between the groups regarding energy intake (kcal/day), absolute intake (g/day) or relative intake (% energy) of carbohydrates, total fat, proteins and fibers. However, confectionery intake was higher in the <6 h group and lower in the ≥8 h group than in the reference group after adjustment for confounding factors. In multivariate analysis, sleep durations <6 h and ≥8 h significantly correlated with increased (95% confidence interval 0.55 to 3.6; P = 0.0078) and decreased (95% confidence interval -4.0 to -0.32; P = 0.021) confectionery intake, respectively. Confectionery intake was positively correlated with female sex, glycated hemoglobin level and dyslipidemia, whereas it was negatively correlated with alcohol consumption and current smoking status. CONCLUSIONS: Short sleep duration is associated with high confectionery intake in people with type 2 diabetes mellitus; this might disturb their glycemic control. Therefore, short sleepers with type 2 diabetes mellitus could improve their glycemic control by avoiding confectionery intake and maintaining adequate sleep duration.

Association of continuous positive airway pressure therapy on cardiac hypertrophy in patients with sleep apnea comorbid with type 2 diabetes mellitus.

Previous reports indicated the therapeutic effect of chronic continuous positive airway pressure (CPAP) therapy on cardiac hypertrophy due to sleep apnea syndrome. However, little is known for cases involving diabetic complications. This retrospective observational study examined the effects of CPAP therapy on left ventricular hypertrophy (LVH) in patients with obstructive sleep apnea syndrome (OSAS) and type 2 diabetes mellitus (T2DM). For all cases, the observation period was 3 years from the time when the patient was introduced to CPAP therapy. Overall, 123 patients were divided into a good CPAP group (CPAP ≥4 h/day, n = 63) and non-adherence group (CPAP <4 h/day, n = 60). The mean CPAP usage times were 5.58 ± 1.23 and 1.03 ± 1.17 h/day in the good CPAP and non-adherence groups, respectively. Regression tendencies of the thickness of the left ventricular posterior (-0.30 ± 1.19 mm) and interventricular septal walls (-0.48 ± 1.22 mm) were observed in the good CPAP group. Hypertrophic tendencies of the left ventricular posterior wall (+0.59 ± 1.44 mm) and interventricular septal wall thickness (+0.59 ± 1.43) were observed in the non-adherence group. Left ventricular posterior wall thickness (coefficient: -0.254, p = 0.0376) and interventricular septal wall thickness (coefficient: -0.426, p = 0.0006) were more likely to be greater in the non-adherence group than in the good CPAP group. Patients in the non-adherence group with an apnea hypopnea index ≥30 had increased left ventricular posterior wall thickness (coefficient: -0.263, p = 0.0673) and interventricular septal wall thickness (coefficient: -0.450, p = 0.0011). In conclusion, appropriate CPAP therapy is an effective treatment for LVH in patients with T2DM and OSAS, especially for severe cases.

Dispersion Effect of Molecular Crowding on Ligand-Protein Surface Binding Sites of Escherichia coli RNase HI.

The molecular crowding effect on ligand-protein interactions, which plays several crucial roles in life processes, has been investigated using various models by adding crowding agents to mimic the intracellular environment. Several studies evaluating this effect have focused on the ligand-protein binding reaction of well-structured binding sites with rigid conformations. However, the crowding effect on flexible binding sites is not well-understood, especially in terms of the conformations. In this work, to elucidate the detailed molecular mechanism underlying the ligand-protein interactions with flexible binding sites on a protein surface, we studied the interaction between the basic protrusion of Escherichia coli ribonuclease HI (RNase HI) and 8-anilinonaphthalene-1-sulfonic acid (ANS). The RNase HI concentration-dependent measurement of ANS fluorescence combined with the multivariate analysis and the fluorescence vibronic structure analysis revealed an increase in the heterogeneous species with an increase in the protein concentration, which is a different behavior from that of proteins with rigid binding sites. This result indicates that ANS molecules bind to the additional binding sites because of the destabilization of the main sites by the excluded volume effect in a crowded environment. The fluorescence vibronic structure analysis yields a detailed molecular picture, indicating that the main species of ANS can have a distorted structure. On the other hand, some ANS molecules move to the minor binding sites of a different microenvironment to secure a stabilized structure. These spectroscopic analyses may show a hypothesis, suggesting that the decrease in the ΔG difference between the main and minor sites due to destabilization of the main binding site could lower the potential barrier between them, inducing the dispersion of binding pathways.

A cross-sectional study of the relationship between quality of life and sleep quality in Japanese patients with type 1 diabetes mellitus.

This study aimed to reveal the relationship between quality of life (QOL) and sleep quality in patients with type 1 diabetes mellitus (T1DM). Overall, 202 patients with T1DM were registered in our study, and 192 were eligible for analysis. Baseline characteristics and laboratory values were determined. Patients completed the Japanese versions of the Pittsburgh Sleep Quality Index (PSQI) and Diabetes Therapy-Related QOL (DTR-QOL) questionnaires. We investigated the relationship between the global PSQI and DTR-QOL total scores by using linear regression analysis. In univariate regression analysis, DTR-QOL total scores were associated with body mass index, alcohol consumption, hypertension, hemoglobin A1c (HbA1c), and global PSQI score (all p-value <0.05) but not with sleep duration. When the association between PSQI subscales and DTR-QOL total scores was examined, DTR-QOL total scores were significantly related to subjective sleep quality and daytime dysfunction. In a multivariate regression analysis, the global PSQI score was negatively related to DTR-QOL total scores. Patients with an HbA1c concentration ≥8.0% had significantly lower DTR-QOL total scores. We revealed a relationship between QOL and sleep quality in T1DM patients and showed that the relationship between QOL and PSQI subscales in T1DM patients may be different from that in patients with type 2 diabetes mellitus. Assessing and managing sleep quality may be necessary for patients with diabetes to improve QOL.

Multicolor imaging of calcium-binding proteins in human kidney stones for elucidating the effects of proteins on crystal growth.

The pathogenesis of kidney stone formation includes multi-step processes involving complex interactions between mineral components and protein matrix. Calcium-binding proteins in kidney stones have great influences on the stone formation. The spatial distributions of these proteins in kidney stones are essential for evaluating the in vivo effects of proteins on the stone formation, although the actual distribution of these proteins is still unclear. We reveal micro-scale distributions of three different proteins, namely osteopontin (OPN), renal prothrombin fragment 1 (RPTF-1), and calgranulin A (Cal-A), in human kidney stones retaining original mineral phases and textures: calcium oxalate monohydrate (COM) and calcium oxalate dihydrate (COD). OPN and RPTF-1 were distributed inside of both COM and COD crystals, whereas Cal-A was distributed outside of crystals. OPN and RPTF-1 showed homogeneous distributions in COM crystals with mosaic texture, and periodically distributions parallel to specific crystal faces in COD crystals. The unique distributions of these proteins enable us to interpret the different in vivo effects of each protein on CaOx crystal growth based on their physico-chemical properties and the complex physical environment changes of each protein. This method will further allow us to elucidate in vivo effects of different proteins on kidney stone formation.

Dynamic Assembly/Disassembly of Staphylococcus aureus FtsZ Visualized by High-Speed Atomic Force Microscopy.

FtsZ is a key protein in bacterial cell division and is assembled into filamentous architectures. FtsZ filaments are thought to regulate bacterial cell division and have been investigated using many types of imaging techniques such as atomic force microscopy (AFM), but the time scale of the method was too long to trace the filament formation process. Development of high-speed AFM enables us to achieve sub-second time resolution and visualize the formation and dissociation process of FtsZ filaments. The analysis of the growth and dissociation rates of the C-terminal truncated FtsZ (FtsZt) filaments indicate the net growth and dissociation of FtsZt filaments in the growth and dissociation conditions, respectively. We also analyzed the curvatures of the full-length FtsZ (FtsZf) and FtsZt filaments, and the comparative analysis indicated the straight-shape preference of the FtsZt filaments than those of FtsZf. These findings provide insights into the fundamental dynamic behavior of FtsZ protofilaments and bacterial cell division.

Structure-function studies of ultrahigh molecular weight isoprenes provide key insights into their biosynthesis.

Some plant trans-1,4-prenyltransferases (TPTs) produce ultrahigh molecular weight trans-1,4-polyisoprene (TPI) with a molecular weight of over 1.0 million. Although plant-derived TPI has been utilized in various industries, its biosynthesis and physiological function(s) are unclear. Here, we identified three novel Eucommia ulmoides TPT isoforms-EuTPT1, 3, and 5, which synthesized TPI in vitro without other components. Crystal structure analysis of EuTPT3 revealed a dimeric architecture with a central hydrophobic tunnel. Mutation of Cys94 and Ala95 on the central hydrophobic tunnel no longer synthesizd TPI, indicating that Cys94 and Ala95 were essential for forming the dimeric architecture of ultralong-chain TPTs and TPI biosynthesis. A spatiotemporal analysis of the physiological function of TPI in E. ulmoides suggested that it is involved in seed development and maturation. Thus, our analysis provides functional and mechanistic insights into TPI biosynthesis and uncovers biological roles of TPI in plants.

Revisiting the Rate-Limiting Step of the ANS-Protein Binding at the Protein Surface and Inside the Hydrophobic Cavity.

8-Anilino-1-naphthalenesulfonic acid (ANS) is used as a hydrophobic fluorescence probe due to its high intensity in hydrophobic environments, and also as a microenvironment probe because of its unique ability to exhibit peak shift and intensity change depending on the surrounding solvent environment. The difference in fluorescence can not only be caused by the microenvironment but can also be affected by the binding affinity, which is represented by the binding constant (K). However, the overall binding process considering the binding constant is not fully understood, which requires the ANS fluorescence binding mechanism to be examined. In this study, to reveal the rate-limiting step of the ANS-protein binding process, protein concentration-dependent measurements of the ANS fluorescence of lysozyme and bovine serum albumin were performed, and the binding constants were analyzed. The results suggest that the main factor of the binding process is the microenvironment at the binding site, which restricts the attached ANS molecule, rather than the attractive diffusion-limited association. The molecular mechanism of ANS-protein binding will help us to interpret the molecular motions of ANS molecules at the binding site in detail, especially with respect to an equilibrium perspective.

Insertion loop-mediated folding propagation governs efficient maturation of hyperthermophilic Tk-subtilisin at high temperatures.

The serine protease Tk-subtilisin from the hyperthermophilic archaeon Thermococcus kodakarensis possesses three insertion loops (IS1-IS3) on its surface, as compared to its mesophilic counterparts. Although IS1 and IS2 are required for maturation of Tk-subtilisin at high temperatures, the role of IS3 remains unknown. Here, CD spectroscopy revealed that IS3 deletion arrested Tk-subtilisin folding at an intermediate state, in which the central nucleus was formed, but the subsequent folding propagation into terminal subdomains did not occur. Alanine substitution of the aspartate residue in IS3 disturbed the intraloop hydrogen-bonding network, as evidenced by crystallographic analysis, resulting in compromised folding at high temperatures. Taking into account the high conservation of IS3 across hyperthermophilic homologues, we propose that the presence of IS3 is important for folding of hyperthermophilic subtilisins in high-temperature environments.

Exploring mutable conserved sites and fatal non-conserved sites by random mutation of esterase from Sulfolobus tokodaii and subtilisin from Thermococcus kodakarensis.

Homologous proteins differ in their amino acid sequences at several positions. Generally, conserved sites are recognized as not suitable for amino acid substitution, and thus in evolutionary protein engineering, non-conserved sites are often selected as mutation sites. However, there have also been reports of possible mutations in conserved sites. In this study, we explored mutable conserved sites and immutable non-conserved sites by testing random mutations of two thermostable proteins, an esterase from Sulfolobus tokodaii (Sto-Est) and a subtilisin from Thermococcus kodakarensis (Tko-Sub). The subtilisin domain of Tko-Sub needs Ca2+ ions and the propeptide domain for stability, folding and maturation. The results from the two proteins showed that about one-third of the mutable sites were detected in conserved sites and some non-conserved sites lost enzymatic activity at high temperatures due to mutation. Of the conserved sites in Sto-Est, the sites on the loop, on the surface, and far from the active site are more resistant to mutation. In Tko-Sub, the sites flanking Ca2+-binding sites and propeptide were undesirable for mutation. The results presented here serve as an index for selecting mutation sites and contribute to the expansion of available sequence range by introducing mutations at conserved sites.

Spectroscopic Evidence of the Salt-Induced Conformational Change around the Localized Electric Charges on the Protein Surface of Fibronectin Type III.

The effect of salt on the electrostatic interaction of a protein is an important issue, because addition of salt affects protein stability and association/aggregation. Although adding salt is a generally recognized strategy to improve protein stability, this improvement does not necessarily occur. The lack of an effect upon the addition of salt was previously confirmed for the tenth fibronectin type III domain from human fibronectin (FN3) by thermal stability analysis. However, the detailed molecular mechanism is unknown. In the present study, by employing the negatively charged carboxyl triad on the surface of FN3 as a case study, the molecular mechanism of the inefficient NaCl effect on protein stability was experimentally addressed using spectroscopic methods. Complementary analysis using Raman spectroscopy and 8-anilino-1-naphthalenesulfonic acid fluorescence revealed the three-phase behavior of the salt-protein interaction between NaCl and FN3 over a wide salt concentration range from 100 mM to 4.0 M, suggesting that the Na+-specific binding to the negatively charged carboxyl triad causes a local conformational change around the binding site with an accompanying structural change in the overall protein, which contributes to the protein's structural destabilization. This spectroscopic evidence clarifies the molecular understanding of the inefficiency of salt to improve protein stability. The findings will inform the optimization of formulation conditions.

Microflow system promotes acetaminophen crystal nucleation.

It is known that interfaces have various impacts on crystallization from a solution. Here, we describe crystallization of acetaminophen using a microflow channel, in which two liquids meet and form a liquid-liquid interface due to laminar flow, resulting in uniform mixing of solvents on the molecular scale. In the anti-solvent method, the microflow mixing promoted the crystallization more than bulk mixing. Furthermore, increased flow rate encouraged crystal formation, and a metastable form appeared under a certain flow condition. This means that interface management by the microchannel could be a beneficial tool for crystallization and polymorph control.

Highly active enzymes produced by directed evolution with stability-based selection.

In a directed evolution aimed at improving enzymatic activity, a situation occurs where highly active variants can no longer be obtained from a template protein because the template is already located at a peak (local maximum) in the fitness landscape of activity for the sequence space. To overcome this situation, the template needs to descend the mountain (lose activity) once and climb another higher mountain. However, there is no solid guideline of how the template should go down. Here, we propose a stability index. Previous studies have shown that protein evolution is potentially governed by stability, and that proteins with low activity but high stability are more favorable templates for producing highly active variants. In our earlier works on conventional directed evolution by random mutagenesis of an esterase from Sulfolobus tokodaii, we identified variants with 3-fold higher activity than the wild-type as the highest activity variants. In this work, as a first step, stability-keeping variants were selected by five rounds of random mutagenesis and screening based on halo formation assay using the substrate tributyrin at 70 °C after heat treatment for 30 min at 90 °C. These variants are likely to be scattered at the feet of various mountains in the fitness landscape. Next, these variants were pooled and used as parental proteins for a conventional experiment with activity-based selection, where the activity of variants was assayed using their cell-free extracts on the substrate p-nitrophenyl butyrate at 75 °C. After two rounds of random mutagenesis, we successfully obtained a variant with 9-fold higher activity than the wild-type. These results indicate that the two-step selection by stability and activity enables us more easily to produce markedly activity-improving variants.