Biochemical analysis and the preliminary crystallographic characterization of d-tagatose 3-epimerase from Rhodobacter sphaeroides
© The Author(s) 2017
Received: 12 August 2017
Accepted: 1 November 2017
Published: 9 November 2017
d-Tagatose 3-epimerase epimerizes d-fructose to yield d-psicose, which is a rare sugar that exists in small quantities in nature and is difficult to synthesize chemically. We aim to explore potential industrial biocatalysts for commercial-scale manufacture of this rare sugar. A d-tagatose 3-epimerase from Rhodobacter sphaeroides (RsDTE) has recently been identified as a d-tagatose 3-epimerase that can epimerize d-fructose to yield d-psicose with a high conversion rate.
The purified RsDTE by Ni-affinity chromatography, ionic exchange chromatography and gel filtration forms a tetramer in solution. The maximal activity was in Tris–HCl buffer pH 8.5, and the optimal temperature was at 35 °C. The product, d-psicose, was confirmed using HPLC and NMR. Crystals of RsDTE were obtained using crystal kits and further refined under crystallization conditions such as 10% PEG 8000,0.1 M HEPES pH 7.5, and 8% ethylene glycol at 20 °C using the sitting-drop vapor diffusion method. The RsDTE homology model showed that it possessed the characteristic TIM-barrel fold. Four residues, Glu156, Asp189, Gln215 and Glu250, forms a hydrogen bond network with the active Mn(II) for the hydride transfer reaction. These residues may constitute the catalytic tetrad of RsDTE. The residues around O1, O2 and O3 of the substrates were conserved. However, the binding-site residues are different at O4, O5 and O6. Arg118 formed the unique hydrogen bond with O4 of d-fructose which indicates RsDTE’s preference of d-fructose more than any other family enzymes.
RsDTE possesses a different metal-binding site. Arg118, forming unique hydrogen bond with O4 of d-fructose, regulates the substrate recognition. The research on d-tagatose 3-epimerase or d-psicose 3-epimerase enzymes attracts enormous commercial interest and would be widely used for rare sugar production in the future.
d-Psicose, a rare sugar that exists in small quantities in nature, is a valuable low-calorie sweetener food additive [1, 2] due to its poor absorbance in the digestive tract and almost no energy. The supplemental d-psicose might be helpful for diabetic patients in preventing post prandial hyperglycemia [3, 4]. d-Psicose is commercially produced by alkali isomerization. However, biosynthesis of rare sugars is environmentally friendly and sustainable with moderate reaction conditions and high efficiency compared with chemical synthesis [5, 6]. Therefore, the most promising approach to produce d-psicose is to use an enzymatic reaction catalyzed by d-psicose-tagatose 3-epimerase (DTE) or d-psicose 3-epimerase (DPE) [6–8]. DTE and DPE were named based on the substrate specificities toward d-tagatose or d-psicose, respectively. Both of these enzymes catalyze the epimerization of various ketohexoses, such as d-tagatose and d-fructose to d-sorbose and d-psicose, respectively, by catalyzing the epimerization at carbon-3 (C3) position [9–11].
The characterization study of DTE/DPE from Pseudomonas cichorii (PcDTE, 290 amino acid residues, 32,615 Da) and Agrobacterium tumefaciens (AtDPE, 283 amino acid residues, 30,650 Da) was reported to efficiently catalyze epimerization of not only d-tagatose to d-sorbose, but also d-fructose to d-psicose [12, 13]. AtDPE has a sequence similarity of 39% with PcDTE. However, there are significant differences in enzymatic properties between them. For instance, PcDTE shows the highest epimerization activity toward d-tagatose, while AtDPE has higher bioconversion rate toward d-psicose than PcDTE . Furthermore, the crystal structures of AtDPE and PcDTE were reported [14, 15], suggesting that both enzymes showed a (β/α)8 TIM barrel fold with a Mn2+ metal ion in the active site. Two glutamate residues and a metal ion conduct the epimerization reaction analogous to the catalytic mechanism of d-ribose-5-phosphate 3-epimerase. Protein engineering was performed on PcDTE to improve its catalytic activity using its structural information . Based on the structural information [8, 14, 15, 17, 18], DTE and DPE shared the same catalytic mechanism of deprotonation/protonation at C3 by two Glu residues. One of the Glu residues removes a proton from C3 to generate a cis-enediolate intermediate, and then the other one protonates C3 on the opposite side.
A novel identified DTE from Rhodobacter sphaeroides (RsDTE) has shown the highest activity for d-fructose with a bioconversion rate to d-psicose of 23% at 40 °C and pH 9.0 . However, the reason why RsDTE showed a high specificity for d-fructose remains unclear. The structure analysis will provide more information for the catalytic reaction. In this study, the RsDTE was purified, crystallized and characterized. The substrate-binding site was analyzed based on RsDTE homology model. The research on DTE/DPE enzymes attracts enormous commercial interest and would be widely used for rare sugar production in the future.
Cloning and expression
Genomic DNA of DTE was prepared from R. sphaeroides strain as previously described  and was used as a template for genomic polymerase chain reaction (PCR) with the PrimeSTAR HS DNA polymerase (TaKaRa, Dalian, China) and a pair of specific primers, F: (5′-GGAATTCCATATGAAAAATCCTGTCGGCATCATCTCG-3′) and R: (5′-CCGGAATTCTCAGTGGGTCACCTCCGCC-3′), for the 5′- and 3′-untranslated regions, respectively. PCR was conducted using temperature settings of 95 °C for 5 min followed by 30 cycles of 95 °C for 15 s, 55 °C for 15 s and 72 °C for 45 s. The final step for extension was 72 °C for 2 min. The RsDTE gene (GenBank Accession No. NC_007494.2) was cloned into the vector pET28a(+) (Novagen, Madison, WI, USA) between the NdeI and EcoRI sites (Additional file 1: Figure S1) with 20 residues (MGSSHHHHHHSSGLVPRGSH) at the N-terminus, including His6 tag and thrombin digestion site. E. coli BL21(DE3) cells harboring the pET28a(+) plasmid with the RsDTE gene (pET28a-RsDTE) were transformed and grown in Lysogeny broth (LB) at 37 °C. Isopropyl β-d-1-thiogalactopyranoside (IPTG) was added at a final concentration of 0.5 mM when the OD600 value reached 0.6, and the cultures were further incubated at 25 °C overnight.
Purification of RsDTE enzyme
After harvesting by centrifugation at 5000g and 4 °C for 15 min, the cells were resuspended in 100 mL lysis buffer (20 mM Tris–HCl pH 8.0, 10 mM imidazole, 0.5 M NaCl, and 1 mM dithiothreitol), disrupted by sonication using an ultrasonicator, set at 1 s pulse, 1 s output and 50% duty cycle for 30 min, and the cell debris was removed by centrifugation at 40,000g at 4 °C for 30 min. Cleared lysate was trapped on 3 mL of Ni–NTA Superflow resin (Qiagen, Hilden, Germany). After washing with 10 mL lysis buffer, the protein was eluted with 15 mL elution buffer (20 mM Tris–HCl pH 8.0, 300 mM imidazole, 100 mM NaCl, and 1 mM dithiothreitol). The eluted solution was dialyzed against 20 mM Tris–HCl pH 8.0, 1 mM dithiothreitol and was further purified by ion exchange using Source Q 4.6/100 PE (volume: 1 mL, flow rate: 3 mL/min, GE Healthcare) and gel filtration chromatography using a Superdex200 column 10/300 GL (volume: 24 mL, flow rate: 0.5 mL/min, GE Healthcare) in running buffer of 20 mM Tris–HCl pH 8.0, 100 mM NaCl, and 1 mM dithiothreitol [19, 20]. The concentration of total protein after each purification step was determined by the BCA assay following to the manufacturer’s protocol. The amount and purity of target protein were analyzed by SDS–PAGE and densitometry of CBB-stained gels using Image Lab Software (Bio-Rad, Hercules, California, USA). The resulting solution containing RsDTE was used for activity assays.
Circular dichroism (CD) spectra (190–250 nm) were recorded using a MOS-450 CD spectropolarimeter (Biologic, Claix, Charente, France) with a 1 mm path-length cell at room temperature. The spectra were obtained as the average of four scans with a bandwidth of 0.1 nm, a step resolution of 0.1 nm and a scan rate of 1 nm/s. The CD spectra of RsDTE (0.08 mg/mL) were recorded in 20 mM Tris–HCl (pH 8.0), 0.1 M NaCl, and 1 mM DTT. Analysis of the protein secondary structure was performed using the program SELCON3 (http://www.dichroweb.cryst.bbk.ac) .
Confirmation of the product d-psicose using HPLC and NMR
The product d-psicose was analyzed on a high-performance liquid chromatography (HPLC) system equipped with a Prevail Carbohydrate ES column-W (5 μm, 4.6 × 250 mm, Agela Technologies, China), an Agilent (USA) multichannel interface, and a XWK-III pump. Acetonitrile (85%) was used as the eluent at a flow rate of 1 mL/min. The column temperature was kept at 40 °C. The product was detected using an evaporative light-scattering detector (Agilent 1260 Infinity, USA).
The 1H and 13CNMR spectra were recorded in D2O on a Bruker AV-500 spectrometer at working frequencies 400 and 100 MHz, respectively. Chemical shifts are expressed in ppm (δ) values and coupling constants (J) in Hz.
The activity was determined by measuring d-psicose formation using d-fructose as a substrate [22, 23]. The reactions were performed carried out with 5 μM RsDTE, 10 g/L d-fructose and 10 μM Mn2+ in 20 mM Tris–HCl buffer (pH 8.5) at 35 °C for 10 min. After the reaction, d-psicose was measured using HPLC.
The optimal pH of the purified enzyme was determined in the following buffers: 20 mM MES buffer (pH 5.5–6.5), 20 mM PBS buffer (pH 7.0–8.0), 20 mM Tris–HCl buffer (pH 8.5–9.0) and 20 mM CAPS buffer (pH 9.5–11.0) at 35 °C for 10 min. The optimal temperature was determined by incubating RsDTE at different temperatures ranging from 20 to 70 °C. To determine the effect of metal ions on the activity of RsDTE, 0.1 mM final concentrations of Ca2+, Mg2+, Mn2+, Ni2+, Co2+, Zn2+, Cu2+, Fe2+ Fe3+ and EDTA were added to the reaction system. The activity was then measured under standard reaction conditions. All assays were repeated three times, and the data are shown as mean ± SD.
Crystallization of RsDTE
Crystallization of RsDTE was performed with the sparse-matrix screening kits Crystal Screen HT (Hampton Research), Index HT (Hampton Research), Wizard I and II (Emerald BioSystems) and JCSG+ (Qiagen) in 96-well plates using the sitting-drop vapour-diffusion method. For refinement of the crystallization conditions, 0.5 μL of protein was mixed with an equal volume of reservoir solution and equilibrated against 0.5 μL of reservoir solution at 20 °C in 96-well plates. The X-ray diffraction data for RsDTE crystals were collected using an in-house X-ray diffractometer (Rigaku FR-E rotating-anode X-ray generator with R-AXIS VII imaging-plate detector).
Structure modeling of RsDTE
The three-dimensional (3D) homology model of RsDTE was generated using Modeller 9.9 . The crystal structure of PcDTE (PDB ID: 2QUL, 1.79 Å), which has 31% sequence identity to the target protein RsDTE, was chosen as the template. The align2d command was used to automatically generate a sequence alignment between the template and RsDTE. Subsequently, homology modeling was performed by the automodel command. Thereafter, each model was first optimized by the variable target function method with conjugate gradients. Simulated annealing MD simulations were used to refine the structure. Finally, the best model was chosen from the values of the Modeller objective function and the DOPE assessment scores. The PyMol molecular Graphics System (http://www.pymol.org)  was used to visualize and analyze the generated model structure.
The RsDTE mutants were generated using a one-step PCR method with plasmid RsDTE_pET28a wild-type as the template. The primers used for mutagenesis are summarized in Additional file 1: Table S1. After amplification, the PCR reaction mixtures were treated with DpnI to completely digest the template and then transformed into E. coli JM109 cells. All the mutations were confirmed by DNA sequencing. Genes of RsDTE mutants were expressed and proteins were purified using the method described above for wild-type enzyme.
Results and discussion
Gene cloning, overexpression of RsDTE and purification of enzyme
The DTE gene was amplified by PCR using the genomic DNA of R. sphaeroides as a template. The results showed that this gene contained 885 bp of a complete open reading frame encoding a protein (295 amino acids, 31.75 kDa). This gene was deposited in GenBank with the Accession Number NC_007494.2. The recombinant plasmid (pET28a-RsDTE) was confirmed by double-enzyme digestion with NdeI and EcoRI (Additional file 1: Figure S1) and gene sequencing. The results show that the RsDTE gene was subcloned into pET28a plasmid, and the expression vector was constructed successfully.
Secondary-structure contents of RsDTE determined by CD spectroscopy in the wavelength region from 190 to 260 nm
Secondary structure (%)
Characterization of RsDTE and product confirmation
Kinetic parameters of RsDTE wild-type on four substrates
K m (mM)
k cat (s−1)
k cat/K m (s−1 mM−1)
Relative activity (%)
78 ± 2.3
42.60 ± 0.89
0.55 ± 0.008
215 ± 4.8
20.47 ± 0.63
0.10 ± 0.005
23.2 ± 1.8
138 ± 0.8
48.31 ± 0.52
0.35 ± 0.006
69.1 ± 2.3
359 ± 3.5
8.30 ± 0.34
0.02 ± 0.001
4.5 ± 0.4
162 ± 4.2
41.05 ± 0.43
0.25 ± 0.008
58.46 ± 3.9
193 ± 6.2
26.04 ± 0.58
0.14 ± 0.008
22.78 ± 2.7
98 ± 3.5
51.25 ± 0.54
0.52 ± 0.013
148 ± 2.8
15.36 ± 0.44
0.10 ± 0.001
16.28 ± 1.6
Crystallization of RsDTE
We obtained the apo-RsDTE crystals (Additional file 1: Figure S5) in crystallization conditions as shown in Additional file 1: Table S2 after 3 days. The crystallization conditions were further refined in different pHs and precipitants as followings: (A) 10% PEG 8000, 0.1 M HEPES pH 7.5, 8% ethylene glycol; (B) 12% PEG 20,000, 0.1 M MES pH 6.5; and (C) 12% PEG 8000, 0.1 M HEPES pH 7.5, 8% ethylene glycol at 20 °C in 96-well plates (Additional file 1: Figure S6). X-ray diffraction data were collected to about 5–6 Å resolution. Additional file 1: Figure S7 shows an X-ray diffraction image of the RsDTE crystal. The crystal needs to be further refined to obtain higher resolution for structure determination.
Phylogenetic tree analysis of RsDTE
A search for enzymes with different amino acid sequence identities above 40% was performed using BLAST. These sequential and homologous proteins were compared with RsDTE (Additional file 1: Figure S8). RsDTE was very similar in sequence to its homologs in Pannonibacter phragmitetus (WP_050473007.1), Paracoccus alcaliphilus (SEN75178.1) and Poseidonocella pacifica (SFB17294.1), which shares the sequence identities of 84, 82 and 77% with these family members, respectively. The crystal structures of DTE from P. cichorii, A. fabrum str. C58, C. cellulolyticum H10 and T. maritime were determined, however, they only shared 28, 30, 36 and 20% sequence identity with RsDTE, respectively.
Homology modeling of RsDTE
Characteristics of the active site
The modeling structure of RsDTE shows that there were four residues, Glu156, Asp189, Gln215 and Glu250, around the metal-binding site, which form a hydrogen bond network for the hydride transfer reaction (Fig. 5b). This result strongly suggests that these residues constitute the catalytic tetrad of RsDTE, which is partly conserved in the DTE/DPE superfamily (Fig. 5c). Three residues of Glu156, Asp189 and Glu250 were completely conserved. However, Gln215 was not conserved, in which histidine is positioned, instead of glutamate, in P. cichorii, A. tumefaciens, C. cellulolyticum and T. maritime. The substrate d-fructose coordinated Mn(II) in a bidentate manner using its O2 and O3 groups, which form a distorted octahedral coordination geometry complex (Fig. 5b). This modeling structure supports the previously proposed mechanism of deprotonation/protonation at C3 of substrate by two Glu residues (Glu156 and Glu250) [14, 15, 18].
Proposal interactions between enzyme and substrates
The purified d-tagatose 3-epimerase from Rhodobacter sphaeroides catalyzes the epimerization of d-fructose to d-psicose at the C3 position. RsDTE shows optimal conditions in Tris–HCl pH 8.5 at 35 °C. RsDTE was crystallized under refined crystallization conditions. Arg118 around the substrate-binding site was investigated, and R118 W improves the substrate recognition and activity. Consistent with biochemical data, R118 forms a hydrogen bond with O4 of d-fructose. The research of RsDTE in its active form provides a valuable tool for solving this enzyme’s crystal structure in the future.
FL and H-MQ designed the research; ZQ, ZZ, SL and J-WW performed the experiments; ZQ, QG and PX analyzed data; H-MQ, FL and ZQ wrote the paper. All authors read and approved the final manuscript.
We thank Stewart Dods, Ph.D., from Liwen Bianji, Edanz Editing China (http://www.liwenbianji.cn/ac), for editing the English text of a draft of this manuscript.
The authors declare that they have no competing interests.
Availability of data and materials
The data collected upon which this article is based upon are all included in this manuscript and the additional files associated with it.
Consent for publication
Our manuscript does not contain any individual data in any form.
Ethics approval and consent to participate
No animals or human subjects were used in the above research.
This work was supported by National Key Research and Development Project (2016YFD0400803), National Natural Science Foundation of China (31771911), Natural Science Foundation of Tianjin (16JCQNJC09200) and the Overseas High-level Talents Program of Tianjin University of Science and Technology to H.-M. Qin, China, National High Technology Research and Development Program (SS2015AA020601), Science Foundation of Tianjin University of Science and Technology (2015LG22).
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Matsuo T, Izumori K. Effects of dietary d-psicose on diurnal variation in plasma glucose and insulin concentrations of rats. Biosci Biotechnol Biochem. 2006;70:2081–5.View ArticleGoogle Scholar
- Matsuo T, Suzuki H, Hashiguchi M, Izumori K. d-Psicose is a rare sugar that provides no energy to growing rats. J Nutr Sci Vitaminol. 2002;48:77–80.View ArticleGoogle Scholar
- Lu Y, Levin GV, Donner TW. Tagatose, a new antidiabetic and obesity control drug. Diabetes Obes Metab. 2008;10:109–34.View ArticleGoogle Scholar
- Matsuo T, Baba Y, Hashiguchi M, Takeshita K, Izumori K, Suzuki H. Dietary d-psicose, a C-3 epimer of d-fructose, suppresses the activity of hepatic lipogenic enzymes in rats. Asia Pac J Clin Nutr. 2001;10:233–7.View ArticleGoogle Scholar
- Emmadi M, Kulkarni SS. Recent advances in synthesis of bacterial rare sugar building blocks and their applications. Nat Prod Rep. 2014;31:870–9.View ArticleGoogle Scholar
- Zhang WL, Zhang T, Jiang B, Mu WM. Enzymatic approaches to rare sugar production. Biotechnol Adv. 2017;35:267–74.View ArticleGoogle Scholar
- Granstrom TB, Takata G, Tokuda M, Izumori K. Izumoring: a novel and complete strategy for bioproduction of rare sugars. J Biosci Bioeng. 2004;97:89–94.View ArticleGoogle Scholar
- van Overtveldt S, Verhaeghe T, Joosten HJ, van den Bergh T, Beerens K, Desmet T. A structural classification of carbohydrate epimerases: from mechanistic insights to practical applications. Biotechnol Adv. 2015;33:1814–28.View ArticleGoogle Scholar
- Ishida Y, Kamiya T, Izumori K. Production of d-tagatose 3-epimerase of Pseudomonas cichorii ST-24 using recombinant Escherichia coli. J Ferment Bioeng. 1997;84:348–50.View ArticleGoogle Scholar
- Itoh H, Sato T, Izumori K. Preparation of d-psicose from d-fructose by immobilized d-tagatose 3-epimerase. J Ferment Bioeng. 1995;80:101–3.View ArticleGoogle Scholar
- Zhang LT, Mu WM, Jiang B, Zhang T. Characterization of d-tagatose-3-epimerase from Rhodobacter sphaeroides that converts d-fructose into d-psicose. Biotechnol Lett. 2009;31:857–62.View ArticleGoogle Scholar
- Itoh H, Okaya H, Khan AR, Tajima S, Hayakawa S, Izumori K. Purification and characterization of d-tagatose 3-epimerase from Pseudomonas sp. ST-24. Biosci Biotechnol Biochem. 1994;58:2168–71.View ArticleGoogle Scholar
- Kim HJ, Hyun EK, Kim YS, Lee YJ, Oh DK. Characterization of d-psicose 3-epimerase from Agrobacterium tumefaciens converting d-fructose into d-psicose. Appl Environ Microbiol. 2006;72:981–5.View ArticleGoogle Scholar
- Kim K, Kim HJ, Oh DK, Cha SS, Rhee S. Crystal structure of d-psicose 3-epimerase from Agrobacterium tumefaciens and its complex with true substrate d-fructose: a pivotal role of metal in catalysis, an active site for the non-phosphorylated substrate, and its conformational changes. J Mol Biol. 2006;361:920–31.View ArticleGoogle Scholar
- Yoshida H, Yamada M, Nishitani T, Takada G, Izumori K, Karnitori S. Crystal structures of d-tagatose 3-epimerase from Pseudomonas cichorii and its complexes with d-tagatose and d-fructose. J Mol Biol. 2007;374:443–53.View ArticleGoogle Scholar
- Bosshart A, Hee CS, Bechtold M, Schirmer T, Panke S. Directed divergent evolution of a thermostable d-tagatose epimerase towards improved activity for two hexose substrates. ChemBioChem. 2015;16:592–601.View ArticleGoogle Scholar
- Sakuraba H, Yoneda K, Satomura T, Kawakami R, Ohshima T. Structure of a d-tagatose 3-epimerase-related protein from the hyperthermophilic bacterium Thermotoga maritima. Acta Crystallogr Sect F Struct Biol Cryst Commun. 2009;65:199–203.View ArticleGoogle Scholar
- Chan HC, Zhu Y, Hu Y, Ko TP, Huang CH, Ren F, Chen CC, Ma Y, Guo RT, Sun Y. Crystal structures of d-psicose 3-epimerase from Clostridium cellulolyticum H10 and its complex with ketohexose sugars. Protein Cell. 2012;3:123–31.View ArticleGoogle Scholar
- Qin HM, Miyakawa T, Jia MZ, Nakamura A, Ohtsuka J, Xue YL, Kawashima T, Kasahara T, Hibi M, Ogawa J, Tanokura M. Crystal structure of a novel N-substituted l-amino acid dioxygenase from Burkholderia ambifaria AMMD. PLoS ONE. 2013;8:e63996.View ArticleGoogle Scholar
- Qin HM, Yamamura A, Miyakawa T, Kataoka M, Nagai T, Kitamura N, Urano N, Maruoka S, Ohtsuka J, Nagata K, Shimizu S, Tanokura M. Structure of conjugated polyketone reductase from Candida parapsilosis IFO 0708 reveals conformational changes for substrate recognition upon NADPH binding. Appl Microbiol Biotechnol. 2014;98:243–9.View ArticleGoogle Scholar
- Whitmore L, Wallace BA. Protein secondary structure analyses from circular dichroism spectroscopy: methods and reference databases. Biopolymers. 2008;89:392–400.View ArticleGoogle Scholar
- Mu WM, Chu FF, Xing QC, Yu SH, Zhou L, Jiang B. Cloning, expression, and characterization of a d-psicose 3-epimerase from Clostridium cellulolyticum H10. J Agric Food Chem. 2011;59:7785–92.View ArticleGoogle Scholar
- Ishida Y, Kamiya T, Itoh H, Kimura Y, Izumori K. Cloning and characterization of the d-tagatose 3-epimerase gene from Pseudomonas cichorii ST-24. J Ferment Bioeng. 1997;83:529–34.View ArticleGoogle Scholar
- Sali A, Blundell TL. Comparative protein modelling by satisfaction of spatial restraints. J Mol Biol. 1993;234:779–815.View ArticleGoogle Scholar
- Humphrey W, Dalke A, Schulten K. VMD: visual molecular dynamics. J Mol Graph. 1996;14:33–8.View ArticleGoogle Scholar