Long-read sequencing identifies an SVA_D retrotransposon insertion deep within the intron of ATP7A as a novel cause of occipital horn syndrome.

BACKGROUND: SINE-VNTR-Alu (SVA) retrotransposons move from one genomic location to another in a 'copy-and-paste' manner. They continue to move actively and cause monogenic diseases through various mechanisms. Currently, disease-causing SVA retrotransposons are classified into human-specific young SVA_E or SVA_F subfamilies. In this study, we identified an evolutionarily old SVA_D retrotransposon as a novel cause of occipital horn syndrome (OHS). OHS is an X-linked, copper metabolism disorder caused by dysfunction of the copper transporter, ATP7A. METHODS: We investigated a 16-year-old boy with OHS whose pathogenic variant could not be detected via routine molecular genetic analyses. RESULTS: A 2.8 kb insertion was detected deep within the intron of the patient's ATP7A gene. This insertion caused aberrant mRNA splicing activated by a new donor splice site located within it. Long-read circular consensus sequencing enabled us to accurately read the entire insertion sequence, which contained highly repetitive and GC-rich segments. Consequently, the insertion was identified as an SVA_D retrotransposon. Antisense oligonucleotides (AOs) targeting the new splice site restored the expression of normal transcripts and functional ATP7A proteins. AO treatment alleviated excessive accumulation of copper in patient fibroblasts in a dose-dependent manner. Pedigree analysis revealed that the retrotransposon had moved into the OHS-causing position two generations ago. CONCLUSION: This is the first report of a human monogenic disease caused by the SVA_D retrotransposon. The fact that the evolutionarily old SVA_D is still actively transposed, leading to increased copy numbers may make a notable impact on rare genetic disease research.

Achieving unprecedented stability in lyophilized recombinase polymerase amplification with thermostable pyruvate kinase from Thermotoga maritima.

Recombinase polymerase amplification (RPA) is an isothermal DNA amplification reaction at around 41 °C using recombinase (Rec), single-stranded DNA-binding protein (SSB), strand-displacing DNA polymerase (Pol), and an ATP-regenerating enzyme. Considering the onsite use of RPA reagents, lyophilized RPA reagents with long storage stability are highly desired. In this study, as one of the approaches to solve this problem, we attempted to use a thermostable pyruvate kinase (PK). PK gene was isolated from a thermophilic bacterium Thermotoga maritima (Tma-PK). Tma-PK was expressed in Escherichia coli and purified from the cells. Tma-PK exhibited higher thermostability than human PK. The purified Tma-PK preparation was applied to RPA as an ATP-regenerating enzyme. Liquid RPA reagent with Tma-PK exhibited the same performance as that with human PK. Lyophilized RPA reagent with Tma-PK exhibited higher performance than that with human PK. Combined with our previous results of RPA reagents of thermostable Pol from a thermophilic bacterium, Aeribacillus pallidus, the results in this study suggest that thermostable enzymes are preferable to mesophilic ones as a component in lyophilized RPA reagents.

Producing Conventional and Transient Amino(mercapto)methyl Radicals by Addition of Muonium to a Crystalline Thioformamide (Mes*NHCH=S).

Muonium (Mu=µ+e-) is composed of a muon of light isotope of proton (µ+) and electron (e-) and can be used as a light surrogate for a hydrogen atom. In this paper, we investigated addition of muonium to a newly synthesized Mes*-substituted thioformamide (Mes*NHCH=S, Mes*=2,4,6-tBu3C6H2). Transverse-field muon spin rotation (TF-µSR) of a solution sample of the thioformamide confirmed addition of muonium to the sulfur atom leading to the corresponding C-centered radical [Mes*NHC(H)⋅-SMu]. Density functional theory (DFT) calculations assigned a conventional amino(mercapto)methyl radical, in which both nitrogen and carbon were slightly pyramidalized, and the calculated muon hyperfine coupling constant (hfcc) including the muon isotope effect was compatible with the experimentally determined parameter. However, the muon level-crossing resonance (µLCR) spectrum of an anisotropic crystalline sample indicated two paramagnetic species, and the major product showed the considerably larger muon hfcc compared with the conventional structure of the amino(mercapto)methyl radical. The unusual transient muoniated thioformamide with the larger muon hfcc that showed rapid relaxation could be only explained by a transient structure including planarization of the nitrogen and carbon atoms in Mes*NHC(H)⋅-SMu.

Helenus and Ajax, Two Groups of Non-Autonomous LTR Retrotransposons, Represent a New Type of Small RNA Gene-Derived Mobile Elements.

Terminal repeat retrotransposons in miniature (TRIMs) are short non-autonomous long terminal repeat (LTR) retrotransposons found from various eukaryotes. Cassandra is a unique TRIM lineage which contains a 5S rRNA-derived sequence in its LTRs. Here, two new groups of TRIMs, designated Helenus and Ajax, are reported based on bioinformatics analysis and the usage of Repbase. Helenus is found from fungi, animals, and plants, and its LTRs contain a tRNA-like sequence. It includes two LTRs and between them, a primer-binding site (PBS) and polypurine tract (PPT) exist. Fungal and plant Helenus generate 5 bp target site duplications (TSDs) upon integration, while animal Helenus generates 4 bp TSDs. Ajax includes a 5S rRNA-derived sequence in its LTR and is found from two nemertean genomes. Ajax generates 5 bp TSDs upon integration. These results suggest that despite their unique promoters, Helenus and Ajax are TRIMs whose transposition is dependent on autonomous LTR retrotransposon. These TRIMs can originate through an insertion of SINE in an LTR of TRIM. The discovery of Helenus and Ajax suggests the presence of TRIMs with a promoter for RNA polymerase III derived from a small RNA gene, which is here collectively termed TRIMp3.

Increase in the solubility of uvsY using a site saturation mutagenesis library for application in a lyophilized reagent for recombinase polymerase amplification.

BACKGROUND: Recombinase uvsY from bacteriophage T4, along with uvsX, is a key enzyme for recombinase polymerase amplification (RPA), which is used to amplify a target DNA sequence at a constant temperature. uvsY, though essential, poses solubility challenges, complicating the lyophilization of RPA reagents. This study aimed to enhance uvsY solubility. METHODS: Our hypothesis centered on the C-terminal region of uvsY influencing solubility. To test this, we generated a site-saturation mutagenesis library for amino acid residues Lys91-Glu134 of the N-terminal (His)6-tagged uvsY. RESULTS: Screening 480 clones identified A116H as the variant with superior solubility. Lyophilized RPA reagents featuring the uvsY variant A116H demonstrated enhanced performance compared to those with wild-type uvsY. CONCLUSIONS: The uvsY variant A116H emerges as an appealing choice for RPA applications, offering improved solubility and heightened lyophilization feasibility.

Application of recombinant human pyruvate kinase in recombinase polymerase amplification.

Recombinase polymerase amplification (RPA) is an isothermal DNA amplification reaction at around 41°C using recombinase (Rec), single-stranded DNA-binding protein (SSB), strand-displacing DNA polymerase (Pol), and an ATP-regenerating enzyme. In this study, we attempted to use pyruvate kinase instead of creatine kinase (CK) that has been consistently used as an ATP-regenerating enzyme in RPA. Human pyruvate kinase M1 (PKM) was expressed in Escherichia coli and purified from the cells. RPA with PKM was performed at 41°C with the in vitro synthesized urease subunit ß (ureB) DNA from Ureaplasma parvum serovar 3 as a standard DNA. The optimal concentrations of PKM and phosphoenolpyruvate were 20 ng/µL and 10 mM, respectively. The RPA reaction with PKM was more sensitive than that with CK. PKM exhibited higher thermostability than CK, suggesting that the RPA reagents with PKM are preferable to those with CK for onsite use.

Daidara: A gigantic Gypsy LTR retrotransposon lineage in the springtail Allacma fusca genome.

Long terminal repeat (LTR) retrotransposons are the major contributor to genome size expansion, as in the cases of the maize genome or the axolotl genome. Despite their impact on the genome size, the length of each retrotransposon is limited, compared to DNA transposons, which sometimes exceed over 100 kb. The longest LTR retrotransposon known to date is Burro-1 from the planarian Schmidtea medierranea, which is around 35.7 kb long. Here through bioinformatics analysis, a new lineage of gigantic LTR retrotransposons, designated Daidara, is reported from the springtail Allacma fusca genome. Their entire length (25-33 kb) rivals Burro families, while their LTRs are shorter than 1.5 kb, in contrast to other gigantic LTR retrotransposon lineages Burro and Ogre, whose LTRs are around 5 kb long. Daidara encodes three core proteins corresponding to gag, pol, and an additional protein of unknown function. The phylogenetic analysis supports the independent gigantification of Daidara from Burro or Ogre.

Difluoromethylborates and Muonium for the Study of Isonitrile Insertion Affording Phenanthridines via Imidoyl Radicals.

The 6-(difluoromethyl)phenanthridine unit is a highly attractive fluoroalkyl-substituted planar nitrogen heterocycle in pharmaceutical and agrochemical research. In this paper, we report that difluoromethylborates can be used as a source of difluoromethyl radicals for isonitrile insertion, leading to 6-(difluoromethyl)phenanthridines. Tuning the aryl substituents in the difluoromethylborates and oxidizing reagents enabled the synthesis of 6-(difluoromethyl)phenanthridines through the generation of difluoromethyl radical and spontaneous intramolecular cyclization of the CF2H-imidoyl radical intermediates. The presence of difluoromethyl radicals was experimentally confirmed, and the reaction mechanisms including imidoyl radical and prompt cyclization reactions could be supported theoretically. Furthermore, we obtained valuable information about the imidoyl radical intermediate by performing transverse-field muon spin rotation (TF-µSR) measurements of 2-isocyano-4'-methoxy-1,1'-biphenyl and using density functional theory (DFT) calculations to interpret the spectra. Muonium, a simple free radical, preferentially adds to the carbon atom of the isonitrile unit, yielding the corresponding imidoyl radical. The temperature dependence of the muon hyperfine coupling constant and the spin relaxation of the muoniated radical signal are compatible with the intramolecular cyclization of biaryl-substituted imidoyl radicals on the µs time scale.

Two superconducting states with broken time-reversal symmetry in FeSe1-xSx.

Iron-chalcogenide superconductors FeSe1-xSx possess unique electronic properties such as nonmagnetic nematic order and its quantum critical point. The nature of superconductivity with such nematicity is important for understanding the mechanism of unconventional superconductivity. A recent theory suggested the possible emergence of a fundamentally new class of superconductivity with the so-called Bogoliubov Fermi surfaces (BFSs) in this system. However, such an ultranodal pair state requires broken time-reversal symmetry (TRS) in the superconducting state, which has not been observed experimentally. Here, we report muon spin relaxation (µSR) measurements in FeSe1-xSx superconductors for 0 ≤ x ≤ 0.22 covering both orthorhombic (nematic) and tetragonal phases. We find that the zero-field muon relaxation rate is enhanced below the superconducting transition temperature Tc for all compositions, indicating that the superconducting state breaks TRS both in the nematic and tetragonal phases. Moreover, the transverse-field µSR measurements reveal that the superfluid density shows an unexpected and substantial reduction in the tetragonal phase (x > 0.17). This implies that a significant fraction of electrons remain unpaired in the zero-temperature limit, which cannot be explained by the known unconventional superconducting states with point or line nodes. The TRS breaking and the suppressed superfluid density in the tetragonal phase, together with the reported enhanced zero-energy excitations, are consistent with the ultranodal pair state with BFSs. The present results reveal two different superconducting states with broken TRS separated by the nematic critical point in FeSe1-xSx, which calls for the theory of microscopic origins that account for the relation between nematicity and superconductivity.

Base-excision restriction enzymes: expanding the world of epigenetic immune systems.

The restriction enzymes examined so far are phosphodiesterases, which cleave DNA strands by hydrolysing phosphodiester bonds. Based on the mobility of restriction-modification systems, recent studies have identified a family of restriction enzymes that excise a base in their recognition sequence to generate an abasic (AP) site unless the base is properly methylated. These restriction glycosylases also show intrinsic but uncoupled AP lyase activity at the AP site, generating an atypical strand break. Action of an AP endonuclease at the AP site may generate another atypical break, rejoining/repairing of which is difficult. This PabI family of restriction enzymes contain a novel fold (HALFPIPE) and show unusual properties, such as non-requirement of divalent cations for cleavage. These enzymes are present in Helicobacteraceae/Campylobacteraceae and in few hyperthermophilic archaeal species. In Helicobacter genomes, their recognition sites are strongly avoided, and the encoding genes are often inactivated by mutations or replacement, indicating that their expression is toxic for the cells. The discovery of restriction glycosylases generalizes the concept of restriction-modification systems to epigenetic immune systems, which may use any mode of damage to DNA that are considered 'non-self' based on epigenetic modifications. This concept will add to our understanding of immunity and epigenetics.

IS481EU Shows a New Connection between Eukaryotic and Prokaryotic DNA Transposons.

DDD/E transposase gene is the most abundant gene in nature and many DNA transposons in all three domains of life use it for their transposition. A substantial number of eukaryotic DNA transposons show similarity to prokaryotic insertion sequences (ISs). The presence of IS481-like DNA transposons was indicated in the genome of Trichomonas vaginalis. Here, we surveyed IS481-like eukaryotic sequences using a bioinformatics approach and report a group of eukaryotic IS481-like DNA transposons, designated IS481EU, from parabasalids including T. vaginalis. The lengths of target site duplications (TSDs) of IS481EU are around 4 bps, around 15 bps, or around 25 bps, and strikingly, these discrete lengths of TSDs can be observed even in a single IS481EU family. Phylogenetic analysis indicated the close relationships of IS481EU with some of the prokaryotic IS481 family members. IS481EU was not well separated from IS3EU/GingerRoot in the phylogenetic analysis, but was distinct from other eukaryotic DNA transposons including Ginger1 and Ginger2. The unique characteristics of IS481EU in protein sequences and the distribution of TSD lengths support its placement as a new superfamily of eukaryotic DNA transposons.

Recombinase polymerase amplification using novel thermostable strand-displacing DNA polymerases from Aeribacillus pallidus and Geobacillus zalihae.

Recombinase polymerase amplification (RPA) is an isothermal DNA amplification reaction at around 41 °C using recombinase (Rec), single-stranded DNA-binding protein (SSB), and strand-displacing DNA polymerase (Pol). Component instability and the need to store commercial kits in a deep freezer until use are some limitations of RPA. In a previous study, Bacillus stearothermophilus Pol (Bst-Pol) was used as a thermostable strand-displacing DNA polymerase in RPA. Here, we attempted to optimize the lyophilization conditions for RPA with newly isolated thermostable DNA polymerases for storage at room temperature. We isolated novel two thermostable strand-displacing DNA polymerases, one from a thermophilic bacterium Aeribacillus pallidus (H1) and the other from Geobacillus zalihae (C1), and evaluated their performances in RPA reaction. Urease subunit ß (UreB) DNA from Ureaplasma parvum serovar 3 was used as a model target for evaluation. The RPA reaction with H1-Pol or C1-Pol was performed at 41 °C with the in vitro synthesized standard UreB DNA. The minimal initial copy numbers of standard DNA from which the amplified products were observed were 600, 600, and 6000 copies for RPA with H1-Pol, C1-Pol, and Bst-Pol, respectively. Optimization was carried out using RPA components, showing that the lyophilized RPA reagents containing H1-Pol exhibited the same performance as the corresponding liquid RPA reagents. In addition, lyophilized RPA reagents with H1-Pol showed almost the same activity after two weeks of storage at room temperature as the freshly prepared liquid RPA reagents. These results suggest that lyophilized RPA reagents with H1-Pol are preferable to liquid RPA reagents for onsite use.

A unique eukaryotic lineage of composite-like DNA transposons encoding a DDD/E transposase and a His-Me finger homing endonuclease.

BACKGROUND: DNA transposons are ubiquitous components of eukaryotic genomes. A major group of them encode a DDD/E transposase and contain terminal inverted repeats (TIRs) of varying lengths. The Kolobok superfamily of DNA transposons has been found in a wide spectrum of organisms. RESULTS: Here we report a new Kolobok lineage, designated KolobokP. They were identified in 7 animal phyla (Mollusca, Phoronida, Annelida, Nemertea, Bryozoa, Chordata, and Echinodermata), and are especially rich in bivalves. Unlike other Kolobok families, KolobokP adopts a composite-like architecture: an internal region (INT) flanked by two long terminal direct repeats (LTDRs), which exhibit their own short terminal inverted repeats ranging up to 18 bps. The excision of LTDRs was strongly suggested. The LTDR lengths seem to be constrained to be either around 450-bp or around 660-bp. The internal region encodes a DDD/E transposase and a small His-Me finger nuclease, which likely originated from the homing endonuclease encoded by a group I intron from a eukaryotic species. The architecture of KolobokP resembles composite DNA transposons, usually observed in bacterial genomes, and long terminal repeat (LTR) retrotransposons. In addition to this monomeric LTDR-INT-LTDR structure, plenty of solo LTDRs and multimers represented as (LTDR-INT)n-LTDR are also observed. Our structural and phylogenetic analysis supported the birth of KolobokP in the late stage of the Kolobok evolution. We propose KolobokP families propagate themselves in two ways: the canonical transposition catalyzed by their transposase and the sequence-specific cleavage by their endonuclease followed by the multimerization through the unequal crossover. CONCLUSIONS: The presence of homing endonuclease and long terminal direct repeats of KolobokP families suggest their unique dual replication mechanisms: transposition and induced unequal crossover.

Diversity and Evolution of DNA Transposons Targeting Multicopy Small RNA Genes from Actinopterygian Fish.

Dada is a unique superfamily of DNA transposons, inserted specifically in multicopy RNA genes. The zebrafish genome harbors five families of Dada transposons, whose targets are U6 and U1 snRNA genes, and tRNA-Ala and tRNA-Leu genes. Dada-U6, which is inserted specifically in U6 snRNA genes, is found in four animal phyla, but other target-specific lineages have been reported only from one or two species. Here, vertebrate genomes and transcriptomes were surveyed to characterize Dada families with new target specificities, and over 120 Dada families were characterized from the genomes of actinopterygian fish. They were classified into 12 groups with confirmed target specificities. Newly characterized Dada families target tRNA genes for Asp, Asn, Arg, Gly, Lys, Ser, Tyr, and Val, and 5S rRNA genes. Targeted positions inside of tRNA genes are concentrated in two regions: around the anticodon and the A box of RNA polymerase III promoter. Phylogenetic analysis revealed the relationships among actinopterygian Dada families, and one domestication event in the common ancestor of carps and minnows belonging to Cyprinoidei, Cypriniformes. Sequences targeted by phylogenetically related Dada families show sequence similarities, indicating that the target specificity of Dada is accomplished through the recognition of primary nucleotide sequences.

Modified uvsY by N-terminal hexahistidine tag addition enhances efficiency of recombinase polymerase amplification to detect SARS-CoV-2 DNA.

BACKGROUND: Recombinase (uvsY and uvsX) from bacteriophage T4 is a key enzyme for recombinase polymerase amplification (RPA) that amplifies a target DNA sequence at a constant temperature with a single-stranded DNA-binding protein and a strand-displacing polymerase. The present study was conducted to examine the effects of the N- and C-terminal tags of uvsY on its function in RPA to detect SARS-CoV-2 DNA. METHODS: Untagged uvsY (uvsY-Δhis), N-terminal tagged uvsY (uvsY-Nhis), C-terminal tagged uvsY (uvsY-Chis), and N- and C-terminal tagged uvsY (uvsY-NChis) were expressed in Escherichia coli and purified. RPA reaction was carried out with the in vitro synthesized standard DNA at 41 °C. The amplified products were separated on agarose gels. RESULTS: The minimal initial copy numbers of standard DNA from which the amplified products were observed were 6 × 105, 60, 600, and 600 copies for the RPA with uvsY-Δhis, uvsY-Nhis, uvsY-Chis, and uvsY-NChis, respectively. The minimal reaction time at which the amplified products were observed were 20, 20, 30, and 20 min for the RPA with uvsY-Δhis, uvsY-Nhis, uvsY-Chis, and uvsY-NChis, respectively. The RPA with uvsY-Nhis exhibited clearer bands than that with either of other three uvsYs. CONCLUSIONS: The reaction efficiency of RPA with uvsY-Nhis was the highest, suggesting that uvsY-Nhis is suitable for use in RPA.

Development of robust isothermal RNA amplification assay for lab-free testing of RNA viruses.

Simple tests of infectiousness that return results in minutes and directly from samples even with low viral loads could be a potential game-changer in the fight against COVID-19. Here, we describe an improved isothermal nucleic acid amplification assay, termed the RICCA (RNA Isothermal Co-assisted and Coupled Amplification) reaction, that consists of a simple one-pot format of 'sample-in and result-out' with a primary focus on the detection of low copy numbers of RNA virus directly from saliva without the need for laboratory processing. We demonstrate our assay by detecting 16S rRNA directly from E. coli cells with a sensitivity as low as 8 CFU/µL and RNA fragments from a synthetic template of SARS-CoV-2 with a sensitivity as low as 1740 copies/µL. We further demonstrate the applicability of our assay for real-time testing at the point of care by designing a closed format for paper-based lateral flow assay and detecting heat-inactivated SARS-COV-2 virus in human saliva at concentrations ranging from 28,000 to 2.8 copies/µL with a total assay time of 15-30 min.

Muonium Addition to a peri-Trifluoromethylated 9-Phosphaanthracene Producing a High-Energy Paramagnetic π-Conjugated Fused Heterocycle.

In this communication, we report muon spin rotation/resonance (µSR) studies for understanding radical reactivity of 10-mesityl-1,8-bis(trifluoromethyl)-9-phosphaanthracene. Transverse-field muon spin rotation (TF-µSR) and muon avoided level-crossing resonance (µLCR) measurements successfully visualized a paramagnetic species produced by regioselective addition of muonium (Mu) to the skeletal phosphorus atom. Density functional theory (DFT) calculations for the P-muoniation product suggested two possible isomers. Whereas the most stable isomer including the envelope-type phosphorus heterocycle shows considerably different hyperfine coupling constants (hfcs) from those of the TF-µSR and µLCR, the metastable structure accompanying the almost planar tricyclic π-conjugated skeleton could simulate the experimentally determined hfcs. The metastable planar π-conjugated paramagnetic tricyclic-fused skeleton is promoted by the larger zero-point energy due to the light muon (µ+ ), one ninth of the proton mass.

Dimensional reduction by geometrical frustration in a cubic antiferromagnet composed of tetrahedral clusters.

Dimensionality is a critical factor in determining the properties of solids and is an apparent built-in character of the crystal structure. However, it can be an emergent and tunable property in geometrically frustrated spin systems. Here, we study the spin dynamics of the tetrahedral cluster antiferromagnet, pharmacosiderite, via muon spin resonance and neutron scattering. We find that the spin correlation exhibits a two-dimensional characteristic despite the isotropic connectivity of tetrahedral clusters made of spin 5/2 Fe3+ ions in the three-dimensional cubic crystal, which we ascribe to two-dimensionalisation by geometrical frustration based on spin wave calculations. Moreover, we suggest that even one-dimensionalisation occurs in the decoupled layers, generating low-energy and one-dimensional excitation modes, causing large spin fluctuation in the classical spin system. Pharmacosiderite facilitates studying the emergence of low-dimensionality and manipulating anisotropic responses arising from the dimensionality using an external magnetic field.

Optimization of reaction condition of recombinase polymerase amplification to detect SARS-CoV-2 DNA and RNA using a statistical method.

Recombinase polymerase amplification (RPA) is an isothermal reaction that amplifies a target DNA sequence with a recombinase, a single-stranded DNA-binding protein (SSB), and a strand-displacing DNA polymerase. In this study, we optimized the reaction conditions of RPA to detect SARS-CoV-2 DNA and RNA using a statistical method to enhance the sensitivity. In vitro synthesized SARS-CoV-2 DNA and RNA were used as targets. After evaluating the concentration of each component, the uvsY, gp32, and ATP concentrations appeared to be rate-determining factors. In particular, the balance between the binding and dissociation of uvsX and DNA primer was precisely adjusted. Under the optimized condition, 60 copies of the target DNA were specifically detected. Detection of 60 copies of RNA was also achieved. Our results prove the fabrication flexibility of RPA reagents, leading to an expansion of the use of RPA in various fields.

Improvement of Moloney murine leukemia virus reverse transcriptase thermostability by introducing a disulfide bridge in the ribonuclease H region.

Moloney murine leukemia virus (MMLV) reverse transcriptase (RT) is widely used in research and clinical diagnosis. Improvement of MMLV RT thermostability has been an important topic of research for increasing the efficiency of cDNA synthesis. In this study, we attempted to increase MMLV RT thermostability by introducing a disulfide bridge in its RNase H region using site-directed mutagenesis. Five variants were designed, focusing on the distance between the two residues to be mutated into cysteine. The variants were expressed in Escherichia coli and purified. A551C/T662C was determined to be the most thermostable variant.