1Biotransformations Group and 3Structural Chemistry, Schering-Plough Research Institute, U-13-3000, 1011 Morris Avenue, Union, NJ 07083, USA
2 To whom correspondence should be addressed. e-mail: wen-chen.suen{at}spcorp.com
![]() |
Abstract |
---|
![]() ![]() ![]() ![]() ![]() ![]() |
---|
Keywords: Candida antarctica lipaseB/enantioselectivity/family shuffling/lipase/thermostability
![]() |
Introduction |
---|
![]() ![]() ![]() ![]() ![]() ![]() |
---|
|
We now report the application of DNA family shuffling coupled with a high-throughput screening assay to create chimeric thermostable lipase B with enhanced activity toward hydrolysis of DDG.
![]() |
Materials and methods |
---|
![]() ![]() ![]() ![]() ![]() ![]() |
---|
An EcoRI/SalI fragment of pWS52 (Zhang et al., 2003) containing the lipase B gene from C.antarctica ATCC 32657 was cloned into the EcoRI/XhoI site of the YEpFLAG-1 vector resulting in a recombinant plasmid designated pWS52a. The homologous lipase B gene from Hyphozyma sp. CBS 648.91 (Hoegh et al., 1995
; Hashida et al., 1999
) was amplified by PCR with forward (5'-AGTACGAATTCACACCCTTC CCCACGGG-3') and reverse (5'-AGTACCTCGAGTCATCC AGTGATGACGCCC-3') primers. EcoRI and XhoI restriction enzyme recognition sites in these primers are shown in italics. A recombinant plasmid designated pWS60 was constructed by cloning this PCR fragment, containing the Hyphozyma lipase B gene, into the EcoRI/XhoI site of YEpFLAG-1. The homologous lipase B gene from Crytococcus tsukubaensis ATCC 24555 was cloned by PCR amplification of an
700 bp PCR fragment using forward (5'-CCGCCGTTCATGCTCAAC GAC-3') and reverse (5'-GCGTAGGGCATGAGGTCAGG-3') primers. The full-length DNA sequence was completed using the Universal Genome Walker kit (Clontech, Palo Alto, CA). The sequences of the forward and reverse PCR primers were derived from internal conserved DNA sequence regions of the lipase B genes from C.antarctica strains ATCC 32657, ATCC 28323, ATCC 32189 and Hyphozyma sp. strain CBS648.91. The full-length gene of ATCC 24555 was subsequently amplified by PCR using forward (5'-AGTACGAATTCTTTACACCCTTCCCCACGGG-3') and reverse (5'-AGTACCTCGAGTCATCCAGTGATGACG CCC-3') primers. The resultant PCR fragment was then cloned into the EcoRI/XhoI site of YEpFLAG-1 resulting in the recombinant plasmid designated pWS57. Total genomic DNA isolation, PCR amplification and fragment purification, plasmid DNA isolation and nucleotide sequencing analysis were conducted as described previously (Zhang et al., 2003
). All three plasmids (pWS52a, pWS57 and pWS60) were transformed into Saccharomyces cerevisiae BJ3505 using the YEASTMAKER Yeast Transformation System according to the manufacturers instructions (Clontech).
DNA family shuffling library construction
The chimeric library was constructed using a modified DNA family shuffling method (Abecassis et al., 2000; Joern et al., 2002
). Three
1.2 kb DNA fragments containing the lipase B genes from pWS52a, pWS57 and pWS60 were amplified by using Y
N-21 (5'-AGCACAAATAACGGGTTATTG-3') and YcC-21 (5'-TACAGACGCGTGTACGCATGT) primers. The reaction mixture (100 µl) contained 1 µl of Pfu polymerase (2.5 units/µl) (Stratagene, La Jolla, CA), 250 ng of each primer, 50 ng of each plasmid DNA, 0.2 mM dNTP mix and 10 µl of 10x Pfu buffer. The resulting parent DNA fragments were then digested with bovine pancreas DNase I (Roche, Indianapolis, IN) as follows. A 25 µl volume of solution containing 1 µg of each parent DNA was mixed with 1.25 µl of 1 M TrisHCl (pH 7.5) and 1.25 µl of 0.2 M MnCl2. Digestion was initiated with the addition of 1 µl of DNase I (0.1 unit/µl) into the mixture at 15°C. Following incubation for 2030 s, the digestion was stopped by the addition of 12.5 µl of a solution containing 50 mM EDTA and 30% glycerol. The digested fragments were separated by gel electrophoresis. The desired 200450 bp DNA fragments were isolated and purified using the QIAquick gel extraction kit (Qiagen, Valencia, CA).
Reassembly of the DNase I digested fragments was conducted in a 50 µl reaction mixture containing 42 µl of fragment DNA, 5 µl of 10x Pfu buffer, 2 µl of dNTP (10 mM each) and 1 µl (2.5 units) of Pfu polymerase. A progressive hybridization PCR cycling program (Joern et al., 2003) was used for reassembly. The reassembled reaction mixture (1 µl) along with primers Y
N-21 and YcC-21 was used to amplify the full-length genes using a PCR cycling program described previously (Zhang et al., 2003
).
Purified YEpFLAG-1 (50 ng) cut with EcoRI/NotI was combined with purified PCR fragments (300 ng) containing the full-length reassembled genes. The resulting DNA mixture was transformed into freshly prepared S.cerevisiae BJ3505 (200 µl) using the YEASTMAKER Yeast Transformation System according to the manufacturers instructions. The transformed cells were plated on to SD-trp (QBIOgene, Carlsbad, CA) agar plates for isolation of the recombinant yeast clones.
Enzyme expression, library screening and characterization of positive clones
Saccharomyces cerevisiae transformants were inoculated into 96-deep well plates (Corning, Acton, MA) containing 0.6 ml of expression medium (1% yeast extract, 8% peptone, 1% glucose, 3% glycerol and 20 mM CaCl2). Cells were grown with agitation (200 r.p.m.) for 2 days at 30°C, followed by 7 days at 20°C. Supernatants containing the secreted enzymes were then clarified by centrifugation at 2500 g for 15 min.
Enzyme libraries were screened using the following protocol. Enzyme supernatant (10 µl) was transferred from each well of the 96-well growth plate into the corresponding well of the assay plate containing 180 µl of 4.45 mM MOPS (3-morpholinopropanesulfonic acid) buffer (pH 7.6), 2 µl of 0.5% bromothymol blue and 10 µl of a 120 mM solution of DDG dissolved in N,N-dimethylformamide (DMF). Addition of the enzyme solutions resulted in lowering of the pH from 7.6 to 7.0. The plates were then sealed with a ChemiSeal sealing film (USA Scientific, Ocala, FL) and incubated at 30°C with agitation at 200 r.p.m. for 3 h followed by centrifugation at 2500 g for 15 min. Supernatant (100 µl) was transferred to the corresponding well in the measurement plate for determination of A620 using a plate reader (Perkin-Elmer HTS7000). Enantioselectivity and DDG conversion were determined using a 96-deep well plate format as follows: culture supernatants (10 µl) were mixed with 180 µl of 4.45 mM MOPS buffer (pH 7.6) and 10 µl of 400 mM DDG in DMF. This reaction mixture was then agitated at 200 r.p.m. for 3 h at 30°C. The reaction was stopped by adjusting the pH to
3 with 10 µl of 0.2 M sulfuric acid. Samples were then extracted with ethyl acetate (420 µl) and analyzed by reversed-phase or chiral HPLC.
SDSPAGE was conducted using a 12% pre-cast gel (Ready Gel Cell system; Bio-Rad, Hercules, CA). The gel was stained with GelCode Blue Stain reagent (Pierce, Rockford, IL). The relative amount of recombinant lipase B protein in the culture supernatants was determined by comparing the intensity of protein bands among different samples using an Electrophoresis Documentation and Analysis System 120 (Eastman Kodak, Rochester, NY).
Analytical methods
Chiral HPLC was conducted using a Phenomenex Chirex NGLY&DNB column (150x4.6 mm i.d.) equilibrated with hexane2-propanoltrifluoroacetic acid (97:3:0.1). Reversed-phase HPLC was performed using a Phenomenex Luna 5µ C18 column (150x4.6 mm) equibrated with acetonitrilewatertrifluoroacetic acid (60:40:0.1).
The protein melting point (Tm) was determined using circular dichroism by fitting the temperature dependence of ellipticity into a two-state foldingunfolding model as described previously (Zhang et al., 2003).
Enzyme activity assays
The enzyme specific activity for the hydrolysis of DDG was determined by measuring the reaction rates at protein concentrations of 0.42 µg/ml and substrate concentrations of 350 mM in 50 mM MOPS buffer, pH 7.0, at 38°C. The reaction rates were independent of substrate concentrations under this condition. The product formation was quantified by HPLC. Data were analyzed using Grafit 5 software, version 5 (Erithacus Software, Horley, UK).
The enzyme specific activity for the hydrolysis of tributyrin was determined by measuring the initial rate of hydrolysis using volumetric titration at 25°C. The reaction mixture, containing 29 ml of 5 mM potassium phosphate buffer (pH 7) and 0.9 ml of tributyrin, was emulsified in a pH-stat cuvette using high-speed stirring. The reaction was initiated by the addition of 100 µl of purified enzyme (25 µg ). The pH was maintained at 7.0 by automatic addition of 25 mM NaOH using a pH-stat system (718 STAT Titrino; Metrohm, Herisau, Switzerland).
The enzyme specific activity for the hydrolysis of p-nitrophenyl butyrate (p-NB) was determined as described previously (Zhang et al., 2003).
The enzyme activity for the hydrolysis of Tween 80 was determined using Tween 80 agar plates as follows. Cells or culture supernatants were spotted onto Tween 80 agar plates (1% yeast extract, 4% peptone, 1% glucose, 3% glycerol, 2% Tween 80, 20 mM CaCl2 and 1.8% agar) and incubated at 30°C overnight. Active clones were identified by the formation of white fatty acid calcium salt precipitates in agar when Tween 80 was enzymatically hydrolyzed.
Protein purification and analysis
Recombinant lipase B from three parents and seven selected clones was purified by using an Anti-FLAG M1 monoclonal antibody affinity gel (Sigma, St Louis, MO) column and a BioCAD SPRINT perfusion chromatography system as described previously (Zhang et al., 2003). Protein purity was determined by SDSPAGE. Enzyme concentration was determined using Bio-Rad protein assay reagent with bovine serum albumin (BSA) as standard.
Native gel electrophoresis was conducted using the SDSPAGE procedure, except that SDS was excluded from the running buffer and native proteins were used for running samples. Following electrophoresis, the gels were stained with Gelcode Blue Stain reagent for 1 h and rinsed with water. The lipase activity was determined by incubating each protein band on Tween 80 agar plates overnight at 30°C.
Western blots were conducted as follows: purified enzymes (1µg) were separated on 12% SDSPAGE gels and electrophoretically transferred on to PVDF membranes (Bio-Rad). The membranes were soaked for 1 h in TBS buffer (20 mM Tris, 500 mM NaCl, pH 7.4) containing 1% (w/v) casein and then allowed to react with Anti-FLAG M1 primary antibodies (7 µg/ml) (Sigma) in TBS buffer containing 0.05% Tween 20 and 1 mM CaCl2 for 2 h. The membranes were washed with TBS buffer and allowed to react for 1 h with 1:2500 diluted secondary antibodies, conjugated with horseradish peroxidase (HRP) (Promega, Madison, WI). The membranes were then washed with TBS buffer and HRP activity was detected by the addition of the substrate, 3,3',5,5'-tetramethylbenzidine (Promega).
Protein deglycosylation with N-glycanase was conducted by following the protocol provided by the enzyme supplier (Prozyme, San Leandro, CA). Purified enzyme (5 µg) was denatured by boiling in a solution of 0.1% SDS, 50 mM ß-mercaptoethanol for 5 min. NP-40 (final concentration 0.75%) and N-glycanase (5 mU) were then added to the reaction mixture followed by incubation at 37°C for 16 h.
![]() |
Results and discussion |
---|
![]() ![]() ![]() ![]() ![]() ![]() |
---|
To take full advantage of natural diversity, we selected several homologous lipase B genes as candidates for DNA family shuffling. The genes from C.antarctica ATCC 32657 (Zhang et al., 2003) and Hyphozyma sp. CBS 648.91 (Hoegh et al., 1995
; Hashida et al., 1999
) were cloned by PCR as outlined in Materials and methods and determined to have a DNA and a protein sequence identity of 67.1 and 73.9%, respectively. Two additional homologous lipase B genes were cloned from C.antarctica (ATCC 28323 and ATCC 32189) by PCR using a primer pair derived from the beginning and the end of the CALB from strain ATCC 32657. Sequence analysis of these two cloned genes revealed that they have sequence identities with those of CALB from strain ATCC 32657 >91% at the DNA level and >97% at the protein level. To expand the pool of homologous lipase B genes, we decided to clone and sequence a relatively unknown lipase from C.tsukubaensis ATCC 24555. This strain was shown to produce a thermally stable enzyme with an immunological response similar to that of CALB from C.antarctica ATCC 32657 (Ishii, 1993
). Since attempts to clone this lipase using primers derived from the beginning and the end of CALB from strain ATCC 32657 were unsuccessful, we decided to utilize the primers derived from the internal conserved regions of four homologous lipase B genes described above. This strategy yielded an
700 bp fragment that had >70% amino acid sequence identity with the CALB from strain ATCC 32657. This fragment was then expanded to a full-length 951 bp gene by using the Universal Genome Walker kit. Comparison of this novel gene with that of CALB from ATCC 32657 revealed a sequence identity of 68.9 and 77.7% at the DNA and protein levels, respectively.
Of the five homologous lipase B genes cloned, three with the highest divergence, C.antarctica ATCC 32657 (P52), C.tsukubaensis ATCC 24555 (P57) and Hyphozyma sp. CBS 648.91 (P60), were selected for DNA family shuffling. Sequence analysis of these three genes revealed identities varying from 67.1 to 69.6% at the DNA level and from 73 to 81.1% at the protein level (Table I). All three genes were actively expressed in an S.cerevisiae Yeast FLAG expression system as secreted fusion proteins with an eight amino acid FLAG peptide at their N-terminus as described previously (Zhang et al., 2003).
|
Both in vitro and in vivo recombination of homologous genes in yeast were used to create a shuffled library (Abecassis et al., 2000; Joern et al., 2002
). Sequence analysis of 10 randomly picked clones revealed that nine out of 10 clones were hybrids. The number of crossovers per gene varied from one to 10 with an average of 3.4. Only two point mutations were detected in the 10 clones sequenced. This accounts for a mutational frequency of about 0.2 base substitutions per gene. The recombination efficiency and the point mutation rate were similar to those reported for family shuffling of three homologous dioxygenases (Joern et al., 2002
). Moreover, when 420 randomly picked clones were analyzed for Tween 80 and p-NB hydrolytic activity, distinct activity profiles characterized by a wide range of activity variation for each substrate were observed among different clones, indicating high functional diversity of the library (data not shown).
Assay development and the library screening
The hydrolysis of DDG results in the formation of an acid, prompting the use of a pH indicator-based assay for high-throughput screening. Although a number of high-throughput pH-based hydrolytic assays have been reported (Janes et al., 1998; Moris-Varas et al., 1999
), they had to be modified substantially to fit the specific assay needs. In the case of lipase-catalyzed hydrolysis of DDG, the MOPS buffer and the bromothymol blue indicator (both have pKa
7.1) were chosen to ensure proportionality between proton concentration and the absorbance. The concentrations of buffer, substrate, enzyme and indicator, in addition to the reaction time, were all optimized to achieve maximum assay sensitivity, while maintaining the ability to distinguish between clones with similar activity. The reaction progress was monitored by following the color transition from blue (alkaline) to yellow (acidic) visually and by measuring the change in absorbance at 620 nm.
Of the 2500 clones screened using the above assay, 69 were identified with a >2-fold improvement in the rate of DDG hydrolysis compared with the most active parent, P57. Even though the enantioselectivity of the 69 clones varied widely, 16 clones and all three parents catalyzed the hydrolysis of DDG to the desired (S)-monoglutarate product with >99% ee. Interestingly, two clones had reverse enantiomeric preference and produced (R)-monoglutarate with moderate enantiomeric purity (5560% ee).
In order to determine whether the enhanced activity of the 16 clones with desired enatioselectivity was due to their improved catalytic efficiency and/or their elevated protein expression level, we analyzed their crude enzyme preparations by SDSPAGE. The results showed a significant variation in lipase expression level among different clones (data not shown). We therefore selected the seven most active clones for further characterization. Four clones, 2A10, 2B3, 8C7 and 12E10, had an elevated lipase protein expression level compared with that of the most active parent, P57. The expression level of the remaining three clones, 3A4, 4D6 and 5D10, was similar to that of P57.
Characterization of lipase B parents and selected clones
Enzyme purification.
Three parent proteins, P52, P57 and P60, and the seven selected clones were purified using an Anti-FLAG M1 column as described in Materials and methods. Although the purification by affinity chromatography was expected to yield nearly homogeneous protein, SDSPAGE analysis revealed the presence of at least two protein populations in every sample except for the two parent proteins, P52 and P60. We hypothesized that the mixed protein population was probably due to a heterogeneous glycosylation pattern commonly observed during heterologous gene expression in yeast (Romanos et al., 1992; Morawski et al., 2000
). This hypothesis was supported by western blot analysis of the recombinant protein preparations with Anti-FLAG M1 monoclonal antibodies, which confirmed the presence of an N-terminal FLAG peptide in all protein populations detected by SDSPAGE. Hence the purified protein fractions contained completely processed mature recombinant lipase and were free of any other proteins (data not shown). To verify further that the heterogeneity of the purified protein fractions detected by SDSPAGE was caused by their distinct glycosylation pattern, we treated all protein fractions with N-glycanase, an enzyme that specifically removes the sugars from N-linked glycosylation sites. As expected, SDSPAGE of the deglycosylated protein fractions revealed a single band of decreased molecular weight (Figure 2), further supporting the variable glycosylation hypothesis. SDSPAGE of all seven selected clones exhibited similar patterns to those of P57 (data not shown). Consistent with earlier observations (Romanos et al., 1992
; Morawski et al., 2000
), this variable degree of glycosylation did not appear to have any noticeable effect on lipase activity. The discrete protein fractions of all native recombinant enzymes separated by native gel electrophoresis were found to be active in the Tween-80 assay (see Materials and methods). Consequently, no further attempts were made to separate the fractions with different glycosylation patterns.
|
|
|
|
|
|
|
In summary, DNA family shuffling was used to improve the activity of lipase B from C.antarctica towards the hydrolysis of DDG. An efficient high-throughput method was used to screen a library of 2500 clones leading to the identification of seven highly selective chimeras with 520-fold enhanced specific activity towards the substrate of interest. The thermal stability of the selected chimeras, characterized by an increase in Tm and a decrease in the rate of irreversible thermoinactivation at 45°C, was also significantly improved compared with that of CALB. To our knowledge, this is the first report describing the use of DNA family shuffling of homologous lipase B genes to generate chimeric enzymes with improved properties. Further study of these hybrid enzymes should expand our understanding of the structurefunction relationship for these valuable biocatalysts.
![]() |
Acknowledgements |
---|
![]() |
References |
---|
![]() ![]() ![]() ![]() ![]() ![]() |
---|
Anderson,E.M., Larsson,K.M. and Kirk,O. (1998) Biocatal. Biotransform., 16, 181204.[ISI]
Anthonsen,T. and Hoff,B.H. (1998) Chem. Phys. Lipids, 93, 199207.[CrossRef][ISI][Medline]
Cherry,J.R., Fidantsef L. (2003) Curr. Opin. Biotechnol., 14, 438443.[CrossRef][ISI][Medline]
Cherry,J.R., Lamsa,M.H., Schneider,P., Vind,J., Svendsen,A., Jones,A. and Pedersen,A.H. (1999) Nat. Biotechnol., 17, 379384.[CrossRef][ISI][Medline]
Crameri,A., Raillard,S.-A., Bermudez,E. and Stemmer,W.P.C. (1998) Nature, 391, 288291.[CrossRef][ISI][Medline]
Hashida,M., Abo,M., Takamura,Y., Kirk,O., Halkier,T., Pedersen,S., Patkar,S.A. and Hansen,M.T. (1999) US Patent, 5 856 163.
Hoegh,I., Shamkant, Halkier,T. and Hansen,M.T. (1995) Can. J. Bot., 73, S869S875.[ISI]
Homann,M.J., Vail,R., Morgan,B., Sabesan,V., Levy,C., Dodds,D.R. and Zaks,A. (2001) Adv. Synth. Catal., 343, 744749.[CrossRef][ISI]
Hult,K. and Berglund P. (2003) Curr. Opin. Biotechnol., 14, 395400.[CrossRef][ISI][Medline]
Ishii,M. (1993) US Patent, 5 273 898.
Jaeger,K.-E. and Eggert,T. (2002) Curr. Opin. Biotechnol., 13, 390397.[CrossRef][ISI][Medline]
Janes,L.E., Lowendahl,A.C. and Kazlauskas,R.J. (1998) Chem. Eur. J., 4, 23242331.[CrossRef][ISI]
Joern,J.M., Meihold,P. and Arnold,F.H. (2002) J. Mol. Biol., 316, 643656.[CrossRef][ISI][Medline]
Koops,B.C., Papadimou,E., Verheij,H.M., Slotboom,A.J. and Egmond,M.R. (1999) Appl. Microbiol. Biotechnol., 52, 791796.[CrossRef][ISI][Medline]
Lundhaug,K., Overbeeke,P.L.A., Jongejan,J.A. and Anthonsen,T. (1998) Tetrahedron: Asymmetry, 9, 28512856.[CrossRef][ISI]
Moore,J.C. and Arnold,F.H. (1996) Nat. Biotechnol., 14, 458467.[ISI][Medline]
Morawski,B., Lin,Z., Cirino,P., Joo,H., Bandara,G. and Arnold,F.H. (2000) Protein Eng., 13, 377384.[CrossRef][ISI][Medline]
Morgan,B., Dodds,D.R., Zaks,A., Andrew,D.R. and Klesse,R. (1997) J. Org. Chem., 62, 77367743.[CrossRef][ISI]
Moris-Varas,F., Amit Shah,J.A., Nadkarni,N.P., Rozzel,J.D. and Demirjian,D.C. (1999) Bioorg. Med. Chem., 7, 21832188.[CrossRef][ISI][Medline]
Ness,J.E., Welch,M., Giver,L., Bueno,M., Cherry,J.R., Borchert,T.V., Stemmer,W.P.C. and Minshull,J. (1999) Nat. Biotechnol., 17, 893896.[CrossRef][ISI][Medline]
Ness,J.E., Kim,S., Gottman,A., Pak,R., Krebber,A., Borchert,T.V., Govindarajan,S., Mundorff,E.C. and Minshull,J. (2002) Nat. Biotechnol., 20, 12511255.[CrossRef][ISI][Medline]
Park,S., Forro,E., Grewal,H., Fulop,F. and Kazlauskas,R.J. (2003) Adv. Synth. Catal., 345, 986995.[CrossRef][ISI]
Parker,M.-C., Brown,S.A., Robertson,L. and Turner,N.J. (1998) J Chem. Soc. Chem. Commun., 22472248.
Patkar,S.A., Svendsen,A., Kirk,O., Clausen,I.G. and Borch,K. (1997) J. Mol. Catal. B, 3, 5154.[ISI]
Patkar,S., Vind,J., Kelstrup,E., Christensen,M.W., Svendsen,A., Borch,K. and Kirk,O. (1998) Chem. Phys. Lipids, 93, 95101.[CrossRef][ISI][Medline]
Petrounia,L.P. and Arnold,F.H. (2000) Curr. Opin. Biotechnol., 11, 325330.[CrossRef][ISI][Medline]
Reetz,M.T. and Jaeger,K.-E. (2000) Chem. Eur. J., 6, 407412.[CrossRef][ISI]
Reetz,M.T., Tielmann,P., Wiesenhofer,W., Konen,W. and Zonta,A. (2003) Adv. Synth. Catal., 345, 717728.[CrossRef][ISI]
Romanos,M.A., Scorer,C.A. and Clare,J.J. (1992) Yeast, 8, 423488.[ISI][Medline]
Rotticci,D., Ottosson,J., Norin,T. and Hult,K. (2001a) Methods Biotechnol., 15, 261277.
Rotticci,D., Rotticci-Mulder,J.C., Denman,S., Norin,T. and Hult,K. (2001b) ChemBioChem, 2, 766770.[CrossRef][ISI][Medline]
Schmidt-Dannert,C. and Arnold,F.H. (1999) Trends Biotechnol., 17, 135136.[CrossRef][Medline]
Uppenberg,J., Hansen,M.T., Patkar,S. and Jones,T.A. (1994) Structure, 2, 293308.[ISI][Medline]
Zha,D., Eipper,A. and Reetz,M.T. (2003) ChemBioChem, 4, 3439.[CrossRef][ISI][Medline]
Zhang,N., Suen,W.-C., Windsor,W., Xiao,L., Madison,V. and Zaks,A. (2003) Protein Eng., 16, 599605.[CrossRef][ISI][Medline]
Received September 2, 2003; revised November 17, 2003; accepted December 10, 2003 Edited by Andrew Griffiths