Email updates

Keep up to date with the latest news and content from BMC Structural Biology and BioMed Central.

Open Access Highly Accessed Research article

Modeling structure and flexibility of Candida antarctica lipase B in organic solvents

Peter Trodler and Jürgen Pleiss*

Author Affiliations

Institute of Technical Biochemistry, University of Stuttgart, Allmandring 31, D-70569 Stuttgart, Germany

For all author emails, please log on.

BMC Structural Biology 2008, 8:9  doi:10.1186/1472-6807-8-9

The electronic version of this article is the complete one and can be found online at: http://www.biomedcentral.com/1472-6807/8/9


Received:28 September 2007
Accepted:6 February 2008
Published:6 February 2008

© 2008 Trodler and Pleiss; licensee BioMed Central Ltd.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

Background

The structure and flexibility of Candida antarctica lipase B in water and five different organic solvent models was investigated using multiple molecular dynamics simulations to describe the effect of solvents on structure and dynamics. Interactions of the solvents with the protein and the distribution of water molecules at the protein surface were examined.

Results

The simulated structure was independent of the solvent, and had a low deviation from the crystal structure. However, the hydrophilic surface of CALB in non-polar solvents decreased by 10% in comparison to water, while the hydrophobic surface is slightly increased by 1%. There is a large influence on the flexibility depending on the dielectric constant of the solvent, with a high flexibility in water and a low flexibility in organic solvents. With decreasing dielectric constant, the number of surface bound water molecules significantly increased and a spanning water network with an increasing size was formed.

Conclusion

The reduced flexibility of Candida antarctica lipase B in organic solvents is caused by a spanning water network resulting from less mobile and slowly exchanging water molecules at the protein-surface. The reduced flexibility of Candida antarctica lipase B in organic solvent is not only caused by the interactions between solvent-protein, but mainly by the formation of a spanning water network.

Background

Candida antarctica lipase B (CALB) is an efficient catalyst for hydrolysis in water and esterification in organic solvents [1]. It is used in many industrial applications because of its high enantioselectivity, wide range of substrates, thermal stability, and stability in organic solvents [2]. CALB belongs to the α/β hydrolase fold family with a conserved catalytic triad consisting of Ser, His, and Asp/Glu [3]. The binding pocket for the substrates consists of an acyl binding pocket, a large and a medium binding pocket for the small and large moiety of secondary alcohols, respectively. In contrast to most lipases, CALB has no lid covering the entrance to the active site and shows no interfacial activation [4].

It has been shown that many enzymes retain activity in organic solvents [5] and have interesting catalytic properties such as higher thermostability and altered stereoselectivity [6,5,8], which was also observed for CALB [9,10]. As compared to water, organic solvents offer new possibilities in biocatalysis. In organic solvents, enzyme-catalyzed esterifications become feasible, and can be efficiently used due to the high solubility of hydrophobic substrates. Despite many advantages of enzymatic reactions in organic solvents, in most cases the catalytic activity in organic solvents is orders of magnitude lower than in aqueous systems [11-13], because of diffusional limitations, changes in protein flexibility, or destabilization of the enzyme [14,15]. It has been shown that the enzymatic activity in organic solvents can be correlated with the solvent polarity as indicated by the octanol-water partition coefficient logP [5]. However, it has also been reported that activity might be affected by the solvent without correlation to the logP [16,17]. Enzymes that are active in organic solvents retain their native structure upon transfer from water to organic solvents [18]. However, it is essential to add small amounts of water to maintain stability and flexibility of enzymes in organic solvent. Thus enzyme-bound water is essential for catalysis and serves as a lubricant for the enzyme [18]. In contrast, fully dry enzymes are inactive and enzymes in organic solvents with high amounts of water show denaturation [19].

Computer simulations of enzymes in different solvents are a valuable tool to investigate the effect on the structure and dynamics of proteins [20-23]. In most simulations of lipases, organic solvents are implicitly modeled [24]. In simulations with explicit organic solvent models, structural conservation and reduced flexibility in organic solvents as compared to water simulations have already been observed [25,24,31].

The aim of this work was to get insight in the structure and flexibility of CALB in water and five different organic solvents using multiple molecular dynamics simulations. CALB was selected in this study, because of its industrial relevance and its stability in different organic solvents. Using explicit solvent models and a low amount of protein-bound water molecules allows to analyze the position of single water molecules, their exchange during the simulation, and the effect of molecular properties of the organic solvents.

Results

Structure of CALB

In the simulations of CALB in water (WAT) (Figure 1) and in the polar solvent methanol (MET), the root mean square deviation (RMSD) of the backbone atoms to the initial structure increased gradually during 1 ns. After this equilibration phase the RMSD in WAT and MET were in the range from 0.7 to 0.9 Å and CALB was in equilibrium without showing any significant conformational changes. The simulations of CALB in the non-polar solvents chloroform (CL3), isopentane (ISO), toluene (TOL), and cyclohexane (CHE) reached the equilibrium already after 500 ps (Figure 1). After equilibration, the RMSD to the initial structure were in the range from 0.7 to 0.9 Å, which is slightly smaller than simulations in WAT and MET. No significant structural differences of CALB backbone in all solvents were observed, and the RMSD between the conformations of CALB in all solvents were in a range from 0.4 to 0.8 Å. The 2D-RMSD plot, where the root mean square deviation of every conformation to all other conformations of a simulation is shown, demonstrated that the conformational space sampled by CALB in the simulations in WAT and MET was larger than in non-polar solvents (Figure 2). In CALB, hydrophobic and hydrophilic patches are equally distributed over the surface (Additional file 1) and the solvent accessible surface area of hydrophilic and hydrophobic surface residues is similar (Table 1). In the simulations in WAT the total surface increased by 600 Å2 compared to the crystal structure, while in non-polar solvents the total surface is similar to the crystal structure. After simulations in non-polar organic solvents the hydrophilic surface is decreased by about 600 Å2 as compared to simulations in WAT.

Additional file 1. Hydrophobicity of CALB. Surface of CALB colored by hydrophobicity (hydrophobic parts red; hydrophilic parts white)

Format: PDF Size: 63KB Download file

This file can be viewed with: Adobe Acrobat ReaderOpen Data

Table 1. Total and hydrophilic surface of CALB in the crystal structure and averaged over the last 1 ns of each simulation in six solvents. Hydrophilic residues are by negative Eisenberg

thumbnailFigure 1. RMS-Deviation in simulations of CALB. Root mean squared deviations of CALB backbone atoms from the X-ray structure as a function of time during simulations of 2 ns in (A) water and (B) in toluene showing three different simulations of the same system using different initial velocity distributions, indicated by different colors

thumbnailFigure 2. 2D-RMSD in simulations of CALB. The root mean square deviation of every conformation to all other conformations of a simulation, as a function of time of CALB backbone atoms during a 2 ns simulation in (A) water and (B) cyclohexane, is shown in the 2D-RMSD

Flexibility of CALB

To evaluate the flexibility of CALB in different solvents, B-factors per residue were calculated from the simulations. In all solvents CALB showed lower flexibility in the core and the active site, and higher flexibility of residues located on the surface. In the WAT, the overall flexibility, as estimated from the sum of all B-factors, is significantly higher than in the organic solvents (Table 2). The overall flexibility of CALB decreases in organic solvents in the order methanol, isopentane, chloroform, toluene, and cyclohexane. The regions with higher flexibility are the same in every solvent, but the degree of decreased flexibility is different. While the flexibility of the rigid core residues is similar in all solvents, the solvent-dependent changes of flexibility are limited to only 5 surface elements (Figure 3): a short α-helix (residues 139–150), a long α-helix (residues 266–289) which form the entrance to the active site and three surface loops (residues 26–30, and 92–97, 215–222), showing higher difference in flexibility between water and non-polar solvents. The medium binding pocket for secondary alcohols (residues 39, 42, 47, 104, 225) and the acyl binding pocket (residues 134, 138, 157, 189, 190) showed low flexibility, while the large binding pocket for secondary alcohols (residues 141,144,154,285,289,290) showed higher flexibility (Figure 4).

Table 2. Total B-factors (sum of B-factors per residue) of CALB averaged over the last 1 ns of each simulation in six solvents

thumbnailFigure 3. Flexibility of CALB in different solvents. Flexibility (indicated by B-factors per residue) of CALB calculated from the last 1 ns of simulations (A) in water (black), in the polar solvent methanol (green), and in the non-polar organic solvents toluene (red) and cyclohexane (blue) and (B) the non-polar organic solvents chloroform (red) and isopentane (blue)

thumbnailFigure 4. Flexible regions in CALB. Flexibility (indicated by B-factors per residue) of CALB from simulations in water mapped on the surface of the crystal structure (blue: low, red: high flexibility); highly flexible including regions are labeled

The flexibility in different solvents was correlated to logP and the dielectric constant of the solvents. The simulations showed a decrease of CALB flexibility by increasing logP and a decreasing dielectric constant (Figure 5). The simulations in ISO are an exception, which will be discussed later. For all twenty amino acids, the highest flexibility was observed in the WAT. However, the extent of the decrease of flexibility in non-polar organic solvents depends on the type of amino acid. While in non-polar solvents the flexibility of glutamine and the hydrophobic amino acids V, Y, T, F, L, I, and W decreased only slightly to 80, 73, 74, 76, 80, 81, 82 and 98%, respectively, of their flexibility in water, the flexibility of alanine and the hydrophilic amino acids D, R, E and N was significantly decreased in non-polar organic solvents to 54, 44, 55, 56 and 62%.

thumbnailFigure 5. Dependence of CALB flexibility on logP and dielectric constant. Flexibility of CALB (sum of B-factors) calculated from the last 1 ns of simulations as a function of (A) logP and (B) dielectric constant from simulations in different solvents. logP values and dielectric constants were taken from the CRC Handbook of Chemistry and Physics [75]

Water at the surface

During the simulations in WAT and in organic solvents, the nine buried water molecules which have been identified in the crystal structure [3] were not exchanged. In contrast, all water molecules that were initially bound at the surface were exchanged by other water molecules during the simulations in the WAT, either by water molecules from the bulk or by other surface-bound water molecules. In non-polar solvents the water exchange at the surface is less frequent and the bound water molecules were less mobile depending on the type of solvent (Figure 6), using a B-factor of 40 Å2 as cutoff. While in WAT there were only 13 ± 2 bound water molecules (averaged over three simulations) with a low B-factor less than 40 Å2 identified, this number increased with increasing logP of the solvent to 26 ± 2 in MET, 31 ± 7 in ISO, 50 ± 17 in CL3, 63 ± 9 in TOL, and 103 ± 6 in CHE out of a total of 286 water molecules included in the organic solvent simulations. However, while the number of less mobile water molecules at the protein surface was similar in the three simulations of the same system, the binding sites of these water-molecules were different. In CHE simulations the 103 less mobile surface water molecules formed a spanning water network covering the whole surface of CALB (Figure 7). This spanning water network consists of water molecules bound to polar side chains or to the backbone of the enzyme, as well as of water molecules bridging these bound waters (Figure 8). The spanning water network is variable and differs slightly in each simulation of the same protein-solvent system. Due to the lower number of less mobile water molecules at the CALB surface in TOL, ISO, and CL3, the spanning water network decreased in size with an decreasing logP. A spanning water network was not observed in MET and WAT simulations. Water clusters include water molecules bound to the protein and water bridges between two water molecules. In CHE eight water clusters of less mobile water molecules were identified by visual inspection (Figure 8). The largest cluster consisted of eight water molecules, where five water molecules were directly bound to polar side chains or to backbone atoms, while three water molecules are bridging these protein-bound waters, but did not interact with the protein itself (Additional file 2, 3, 4, 5). Visual inspection indicated that the active site is accessible to organic solvents. In simulations in MET five methanol and two water molecules were located in the substrate binding site after 2 ns of simulation, in CL3 one chloroform and five water molecules, in ISO one isopentane and six water molecules, in TOL one toluene and five water molecules, and in CHE one cyclohexane and six water molecules. In MET four water molecules that were initially present in the substrate binding site were replaced by methanol, which was not observed by non-polar solvent molecules.

Additional file 2. Cluster I – Structure. Structure of cluster I in the simulation of CALB in cyclohexane, water molecules are displayed as red and white spheres

Format: PDF Size: 66KB Download file

This file can be viewed with: Adobe Acrobat ReaderOpen Data

Additional file 3. Cluster I – Ligands. Coordination of water molecules in cluster I in the simulation of CALB in cyclohexane

Format: PDF Size: 15KB Download file

This file can be viewed with: Adobe Acrobat ReaderOpen Data

Additional file 4. Cluster II – Structure. Structure of cluster II in the simulation of CALB in cyclohexane, water molecules are displayed as red and white spheres

Format: PDF Size: 57KB Download file

This file can be viewed with: Adobe Acrobat ReaderOpen Data

Additional file 5. Cluster II – Ligands. Coordination of water molecules in cluster II in the simulation of CALB in cyclohexane

Format: PDF Size: 12KB Download file

This file can be viewed with: Adobe Acrobat ReaderOpen Data

thumbnailFigure 6. Dependence of slow exchanged water molecules on dielectric constant. Number of less mobile water molecules at the surface of CALB, with a B-factor lower than 40 Å2 calculated from the last 1 ns of simulations, as a function of the dielectric constant from simulations in different solvents

thumbnailFigure 7. Less mobile water molecules at the surface of CALB. Location of less mobile water molecules, with a B-factor lower than 40 Å2, at the surface of the last snapshot of CALB during the simulation (A) in water and (B) in cyclohexane; water molecules are displayed as red and white spheres

thumbnailFigure 8. Water network and water clusters on CALB. Snapshot of the spanning water network during the simulation of CALB in cyclohexane with representative clusters I and II, which consist of surface-bound and bridging water molecules; water molecules are displayed as red and white spheres

Discussion

Three simulations with different initial velocity distributions instead of a single trajectory were used to sample the conformational space in the protein-solvent system, because a single trajectory samples only a small fraction of the conformational space than multiple short trajectories [32-34].

The structure of CALB shows a high stability in all solvents and is therefore a useful system to examine the effect of different solvents on structure and flexibility. Multiple MD simulations of each protein-solvent system confirmed that the structures of CALB in different solvent deviated from each other by less than 0.8 Å. This structural difference is within the deviation of 0.6–0.8 Å of the three simulations of the same protein-solvent system. Our observation that the structure of CALB is independent of its environment is supported by the fact that CALB does not undergo conformational transitions [4] and by a comparison of different crystal structures obtained under different crystallization conditions [3,35] which show backbone RMSD values below 0.3 Å. Also in most molecular dynamics simulations structures showed no significant changes in different solvents [20,22,31] with the exception of Rhizomucor miehei lipase for which a solvent-induced conformational change was observed [30]. In addition, circular dichroism measurements of CALB in different solvents showed that its secondary structure did not change [36]. This is also confirmed by X-ray structures of various serine proteases, crystallized in the presence of small amounts of different organic solvents, which are nearly indistinguishable from their structure in water [37-40]. Flexible and rigid regions of CALB identified in the simulations were similar to regions of high and low B-factors reported in the crystal structure [3]. In contrast to structure, the flexibility is solvent-dependent. In organic solvents, the flexibility of proteins is decreased, which has been confirmed by different experimental techniques such as time-resolved fluorescence anisotropy [41], ESR [42], and dielectric relaxation spectroscopy [43]. This was also observed in simulations of lipases [30] and subtilisin [20], while no significant differences between water and organic solvents have been observed in simulations of subtilisin [44].

It is observed that solvents with a lower dielectric constant lead to a decreased protein flexibility as shown by EPR [45,12], which is in general consistent with the results of our simulations. However, there is one outlier. In ISO the flexibility is higher than expected from its dielectric constant. Interestingly, isopentane is the only solvent in our simulations with one freely rotatable, bond which is not considered in the rigid solvent model. We suppose that a flexible solvent model would be necessary to properly treat the effects of this solvent. It has also been suggested that the size of organic solvent molecules correlates with the protein flexibility [46], but no correlation was found in our simulations. A further effect of organic solvents in our simulations is a decrease in solvent-accessible surface, especially of the hydrophilic surface. In our simulations, the hydrophilic surface decreased by 600 Å2 from simulations in WAT to CHE, while the hydrophobic surface increased slightly by 50 Å2. It has been suggested that in water polar side chains orient toward the surface, thus increasing the hydrophilic surface and decreasing the polar intra-molecular interactions that mediate the rigidity of the protein [47], while in organic solvents the surface area is reduced which leads to improved packing and increased stability [42].

There are in general two types of water molecules observed in organic solvents. The 'inside class' water molecules that are bound in the interior of the protein and can play an important role for active conformation in organic solvents by a hydrogen bond network and the 'contact class' water molecules that are weakly bound to the surface of the enzyme and can be rapidly exchanged by other water molecules [48]. The dynamical properties of surface-bound water molecules differ considerably from the bulk water as shown by X-ray crystallography [49] and NMR experiments [50]. The residence times of most surface-bound water molecules are between 10 and 100 ps, bound and free water molecules at the surface of a protein are in a dynamic equilibrium [51]. In agreement with these experiments all water molecules at the surface were rapidly exchanged during the simulation time of 2 ns. In correlation to our simulations, where an increased number of less mobile water molecules was observed at higher logP at the surface of CALB, the number of water molecules with high B-factor was decreased by an increasing hydrophobicity of alcohols at lysozyme studied by X-ray [52]. In organic solvents the amount of water molecules in the organic solvent phase is low and therefore the probability for an exchange of bound and structured water at the surface is low. A high rate of exchange of water molecules at the surface, observed in our simulations, might be the reason for an increased flexibility, in agreement to previous observations where the protein mobility increased with an increasing amount of water [22].

In agreement with experimental results [14] in the simulations organic solvents strip just a few water molecules from the enzymes surface, while polar solvents strip most water from the surface. This stripping of water from the surface by polar solvents has been already shown in experiments [12] and in simulations [31]. The solvent dielectric constant correlated with stripping of water molecules in simulations with increasing polarity of the solvent [44]. It was shown in experiments, that desorption is independent of the kind of the protein, increasing with the dielectric constant of the solvent [53]. Like in our simulations of CALB, a spanning water network consisting of several small water clusters was observed in previous simulations [54,55,22,31] and experiments [49], in which non-polar solvents enhanced the formation of clusters [31,54]. In agreement to our simulations, the water clusters in a cutinase consisted of 2 to 8 water molecules at a hydration level of 15% (w/w) [31], depending on the solvent, in the crystal structure of CALB the largest cluster of water molecules with B-factors lower 40 Å2 consists of 14 water molecules at most. The spanning water network resulted by a slow exchange of water molecules at the surface in organic solvents. Polar groups favor direct interactions with water molecules and form hydrogen bonds, while non-polar groups enhance interactions among water molecules and enhance the local structure of neighboring water molecules. The concept of hydrophobic hydration and the freezing of water to clusters around hydrophobic surface was previously suggested [56] and is supported by the solvent dependent flexibility observed in simulations.

From the results it can be concluded, that the reduced flexibility of CALB in non-polar solvents is not only a consequence of the interaction between organic solvent molecules and the protein, but also due to the interaction with the enzyme-bound water and its exchange on the surface. Despite the higher fluidity of organic solvents, the flexibility of CALB is decreased, because the water exchange at the surface is restricted.

Methods

Parametrization of solvent models

The organic solvent molecules of cyclohexane, isopentane, and toluene were parameterized. The geometric parameters were derived by ab initio geometry optimization on the HF/6-31G* level using Gaussian98 [57] in the gas phase. The partial charges were derived by fitting partial charges using the RESP program [58] of AMBER 7.0 [59] to the electrostatic potential (Additional file 6). For each molecule a periodic solvent box was built in XLEAP (Additional file 7). The boxes were equilibrated by molecular dynamics simulations using AMBER 7.0 and the all-atom AMBER force-field ff99 [60] by pressure coupling, all solvent molecules were treated as rigid. The calculated densities of the boxes in equilibrium were in good agreement with experimental data (Additional file 7). After minimization (500 steps steepest descent followed by 200 steps conjugate gradient) the systems were heated during 30 ps to 300 K using a temperature coupling constant of 0.8 ps at a pressure of 1 bar using a pressure coupling constant of 1.0 ps [61]. Molecular dynamics simulations of the systems were performed for 2 ns, applying the SHAKE algorithm [62] to constrain the bond lengths. Electrostatic interactions were calculated using Ewald summation [59], Van der Waals interactions were calculated using a 16 Å cut-off.

Additional file 6. Parameters of solvent models. The parameters of the solvent models cyclohexane, isopentane and toluene parametrized by ab initio geometry optimization

Format: PDF Size: 64KB Download file

This file can be viewed with: Adobe Acrobat ReaderOpen Data

Additional file 7. Data of solvent Boxes. Equilibration of solvent boxes by molecular dynamics simulations

Format: PDF Size: 362KB Download file

This file can be viewed with: Adobe Acrobat ReaderOpen Data

System setup

The crystal structure of CALB [PDB: 1TCA] [3] with a resolution of 1.55 Å was taken from the Protein Data Bank as initial structure for the simulations.

pKa values and protonation states of titratable groups Arg, Lys, Asp, Glu, and His were calculated at pH 7 using MEAD [63] and TITRA [64]. The online tool PCE [65] of MEAD was used to solve numerically the Poisson-Boltzmann equation using its MULTIFLEX program with a protein internal dielectric constant ε = 20, solvent dielectric constant ε = 80, ionic strength 0.145 M, grid spacing of 1 Å in a cubic box and the PARSE parameters for radii and charges [66]. The REDTI program was used to compute the protonation states [67]. TITRA is based on the Tanford-Kirkwood model, with default parameters and ε = 20 for the protein dielectric constant. The solvent accessible surface area of each residue used in TITRA was calculated by the program acc_run [68]. Both methods resulted in the same protonation states at pH 7. CALB was calculated to be neutral, the same protonation states were used for the simulations in water and in organic solvents, assuming pH memory from the protonation in aqueous solution after lyophilization [14]. In XLEAP of the Amber 7.0 program package hydrogens were added as calculated by TITRA. The CALB crystal structure including 286 crystal water molecules, corresponding to a hydration level of 15.4% (w/w), was solvated in six different solvent boxes using a minimal distance of 14 Å between the box boundary and the protein. The equilibrated boxes of cyclohexane (CHE), isopentane (ISO), and toluene (TOL) were used as parametrized, boxes of TIP3 water (WAT) [69], methanol (MET) [70] and chloroform (CL3) [71] were used as given in the AMBER package.

Molecular dynamics simulations

Multiple molecular dynamics simulations of the protein-solvent systems were performed using the Amber 7 program package [59] and the all-atom AMBER force field ff99 [60]. The simulations were done in a truncated octahedral box under periodic boundary conditions. Non-bonded interactions were calculated at a cutoff distance of 10 Å. The SHAKE algorithm [62] was applied to all bonds. The simulations were performed at 300 K and 1 bar using a time step of 1 fs. Temperature and pressure of the system were controlled using a weak coupling to an external heat bath [61] with a temperature coupling constant of 1.0 ps and a pressure coupling constant of 1.2 ps. The initial structures were energy minimized (500 steepest descent and 50 conjugate gradient) and followed by a simulation at 300 K and 1 bar by restraining the position of all Cα atoms using a harmonic potential. The force constant was gradually decreased every 50 ps from 10 to 5, 1 and 0.1 kcal/mol followed by an unrestrained simulation of 2 ns. Three simulations of each system were performed using different initial random velocity distributions. All snapshots from the resulting trajectories were fitted the backbone atoms to the initial structure. The root mean square deviation of the backbone atoms between each conformer and the initial structure (RMSD) and between all conformers (2D-RMSD) and B-factors were analyzed using the PTRAJ of AMBER 7.0 [59]. The solvent accessible surface area was calculated by DSSP [72] using a probe radius of 1.4 Å. The protein structures were visualized and the hydrophobicity [73] was mapped on the surface using PyMol 0.98 [74]. Hydrophilic and hydrophobic residues were identified by their negative or positive hydrophobicity index [73], respectively.

Authors' contributions

PT performed the simulations, JP was the principal investigator, conceived the project and guided its development. All authors read and approved the final manuscript.

References

  1. Anderson EM, Larsson KM, Kirk O: One biocatalyst - Many applications: The use of Candida antarctica B-lipase in organic synthesis.

    Biocatal Biotransfor 1998, 16(3):181-204. OpenURL

  2. Rotticci D, Rotticci-Mulder JC, Denman S, Norin T, Hult K: Improved enantioselectivity of a lipase by rational protein engineering.

    Chembiochem 2001, 2(10):766-770. PubMed Abstract | Publisher Full Text OpenURL

  3. Uppenberg J, Hansen MT, Patkar S, Jones TA: The sequence, crystal structure determination and refinement of two crystal forms of lipase B from Candida antarctica.

    Structure 1994, 2(4):293-308. PubMed Abstract | Publisher Full Text OpenURL

  4. Martinelle M, Holmquist M, Hult K: On the interfacial activation of Candida antarctica lipase A and B as compared with Humicola lanuginosa lipase.

    Biochim Biophys Acta 1995, 1258(3):272-276. PubMed Abstract OpenURL

  5. Zaks A, Klibanov AM: Enzyme-catalyzed processes in organic solvents.

    Proc Natl Acad Sci U S A 1985, 82(10):3192-3196. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  6. Wescott CR, Klibanov AM: The Solvent Dependence of Enzyme Specificity.

    Bba-Protein Struct M 1994, 1206(1):1-9. OpenURL

  7. Russell AJ, Klibanov AM: Enzymes in Organic-Solvents.

    Biochem Soc T 1989, 17(6):1145-1145. OpenURL

  8. Stahl M, Jeppssonwistrand U, Mansson MO, Mosbach K: Induced Stereoselectivity and Substrate Selectivity of Bio-Imprinted Alpha-Chymotrypsin in Anhydrous Organic Media.

    J Am Chem Soc 1991, 113(24):9366-9368. Publisher Full Text OpenURL

  9. Kirk O, Würtz Christensen M: Lipases from Candida antarctica: Unique Biocatalysts from a Unique Origin.

    Organic Process Research & Development 2002, 6:446-451. Publisher Full Text OpenURL

  10. Magnusson AO, Rotticci-Mulder JC, Santagostino A, Hult K: Creating space for large secondary alcohols by rational redesign of Candida antarctica lipase B.

    Chembiochem 2005, 6(6):1051-1056. PubMed Abstract | Publisher Full Text OpenURL

  11. Kitaguchi H, Fitzpatrick PA, Huber JE, Klibanov AM: Enzymatic Resolution of Racemic Amines - Crucial Role of the Solvent.

    J Am Chem Soc 1989, 111(8):3094-3095. Publisher Full Text OpenURL

  12. Clark DS: Characteristics of nearly dry enzymes in organic solvents: implications for biocatalysis in the absence of water.

    Philos T Roy Soc B 2004, 359(1448):1299-1307. Publisher Full Text OpenURL

  13. Overbeeke PL, Govardhan C, Khalaf N, Jongejan JA, Heijnen JJ: Influence of lid conformation on lipase enantioselectivity.

    journal of Molecular Catalysis B: Enzymatic 2000, 10(4):385-393. Publisher Full Text OpenURL

  14. Klibanov AM: Why are enzymes less active in organic solvents than in water?

    Trends Biotechnol 1997, 15(3):97-101. PubMed Abstract | Publisher Full Text OpenURL

  15. Dodson G, Verma CS: Protein flexibility: its role in structure and mechanism revealed by molecular simulations.

    Cellular and Molecular Life Sciences 2006, 63(2):207-219. Publisher Full Text OpenURL

  16. Lee YS, Hong JH, Jeon NY, Won K, Kim BT: Highly enantioselective acylation of rac-alkyl lactates using Candida antarctica lipase B.

    Organic Process Research & Development 2004, 8(6):948-951. Publisher Full Text OpenURL

  17. Degn P, Zimmermann W: Optimization of carbohydrate fatty acid ester synthesis in organic media by a lipase from Candida antarctica.

    Biotechnology and Bioengineering 2001, 74(6):483-491. Publisher Full Text OpenURL

  18. Zaks A, Klibanov AM: Enzymatic catalysis in nonaqueous solvents.

    J Biol Chem 1988, 263(7):3194-3201. PubMed Abstract | Publisher Full Text OpenURL

  19. Griebenow K, Klibanov AM: On protein denaturation in aqueous-organic mixtures but not in pure organic solvents.

    J Am Chem Soc 1996, 118(47):11695-11700. Publisher Full Text OpenURL

  20. Colombo G, Ottolina G, Carrea G: Modelling of enzyme properties in organic solvents.

    Monatsh Chem 2000, 131(6):527-547. OpenURL

  21. Colombo G, Carrea G: Modeling enzyme reactivity in organic solvents and water through computer simulations.

    J Biotechnol 2002, 96(1):23-33. PubMed Abstract | Publisher Full Text OpenURL

  22. Soares CM, Teixeira VH, Baptista AM: Protein structure and dynamics in nonaqueous solvents: insights from molecular dynamics simulation studies.

    Biophys J 2003, 84(3):1628-1641. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  23. Affleck R, Xu ZF, Suzawa V, Focht K, Clark DS, Dordick JS: Enzymatic catalysis and dynamics in low-water environments.

    Proc Natl Acad Sci U S A 1992, 89(3):1100-1104. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  24. Peters GH, Olsen OH, Svendsen A, Wade RC: Theoretical investigation of the dynamics of the active site lid in Rhizomucor miehei lipase.

    Biophysical Journal 1996, 71(1):119-129. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  25. Tejo BA, Salleh AB, Pleiss J: Structure and dynamics of Candida rugosa lipase: the role of organic solvent.

    Journal of Molecular Modeling 2004, 10(5-6):358-366. PubMed Abstract | Publisher Full Text OpenURL

  26. Kahlow UH, Schmid RD, Pleiss J: A model of the pressure dependence of the enantioselectivity of Candida rugosalipase towards (+/-)-menthol.

    Protein Sci 2001, 10(10):1942-1952. PubMed Abstract | Publisher Full Text OpenURL

  27. Lee J, Suh SW, Shin S: Computational studies of essential dynamics of Pseudomonas cepacia lipase.

    J Biomol Struct Dyn 2000, 18(2):297-309. PubMed Abstract | Publisher Full Text OpenURL

  28. Cherukuvada SL, Seshasayee ASN, Raghunathan K, Anishetty S, Pennathur G: Evidence of a Double-Lid Movement in Pseudomonas aeruginosa Lipase: Insights from Molecular Dynamics Simulations.

    PLOS computational biology 2005., 1(3) PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  29. James JJ, Lakshmi BS, Seshasayee AS, Gautam P: Activation of Candida rugosa lipase at alkane-aqueous interfaces: A molecular dynamics study.

    FEBS Lett 2007, 581(23):4377-4383. PubMed Abstract | Publisher Full Text OpenURL

  30. Norin M, Haeffner F, Hult K, Edholm O: Molecular dynamics simulations of an enzyme surrounded by vacuum, water, or a hydrophobic solvent.

    Biophys J 1994, 67(2):548-559. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  31. Micaelo NM, Soares CM: Modeling hydration mechanisms of enzymes in nonpolar and polar organic solvents.

    Febs J 2007, 274(9):2424-2436. PubMed Abstract | Publisher Full Text OpenURL

  32. Auffinger P, Louisemay S, Westhof E: Multiple Molecular-Dynamics Simulations of the Anticodon Loop of Trna(Asp) in Aqueous-Solution with Counterions.

    J Am Chem Soc 1995, 117(25):6720-6726. Publisher Full Text OpenURL

  33. Caves LSD, Evanseck JD, Karplus M: Locally accessible conformations of proteins: Multiple molecular dynamics simulations of crambin.

    Protein Science 1998, 7(3):649-666. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  34. Worth GA, Nardi F, Wade RC: Use of multiple molecular dynamics trajectories to study biomolecules in solution: The YTGP peptide.

    J Phys Chem B J Phys Chem B 1998, 102(32):6260-6272. OpenURL

  35. Uppenberg J, Ohrner N, Norin M, Hult K, Kleywegt GJ, Patkar S, Waagen V, Anthonsen T, Jones TA: Crystallographic and molecular-modeling studies of lipase B from Candida antarctica reveal a stereospecificity pocket for secondary alcohols.

    Biochemistry 1995, 34(51):16838-16851. PubMed Abstract | Publisher Full Text OpenURL

  36. McCabe RW, Rodger A, Taylor A: A study of the secondary structure of Candida antarctica lipase B using synchrotron radiation circular dichroism measurements.

    Enzyme and Microbial Technology 2005, 36(1):70-74. Publisher Full Text OpenURL

  37. Schmitke JL, Stern LJ, Klibanov AM: The crystal structure of subtilisin Carlsberg in anhydrous dioxane and its comparison with those in water and acetonitrile.

    P Natl Acad Sci USA 1997, 94(9):4250-4255. Publisher Full Text OpenURL

  38. Yennawar NH, Yennawar HP, Farber GK: X-ray crystal structure of gamma-chymotrypsin in hexane.

    Biochemistry 1994, 33(23):7326-7336. PubMed Abstract | Publisher Full Text OpenURL

  39. Fitzpatrick PA, Ringe D, Klibanov AM: X-ray crystal structure of cross-linked subtilisin Carlsberg in water vs. acetonitrile.

    Biochem Biophys Res Commun 1994, 198(2):675-681. PubMed Abstract | Publisher Full Text OpenURL

  40. Wangikar PP, Michels PC, Clark DS, Dordick JS: Structure and function of subtilisin BPN' solubilized in organic solvents.

    J Am Chem Soc 1997, 119(1):70-76. Publisher Full Text OpenURL

  41. Broos J, Visser AJWG, Engbersen JFJ, Verboom W, vanHoek A, Reinhoudt DN: Flexibility of enzymes suspended in organic solvents probed by time-resolved fluorescence anisotropy. Evidence that enzyme activity and enantioselectivity are directly related to enzyme flexibility.

    J Am Chem Soc 1995, 117(51):12657-12663. Publisher Full Text OpenURL

  42. Toba S, Merz KM: The Concept of Solvent Compatibility and Its Impact on Protein Stability and Activity Enhancement in Nonaqueous Solvents.

    J Am Chem Soc 1997, 119(42):9939-9948. Publisher Full Text OpenURL

  43. Mijovic J, Bian Y, Gross RA, Chen B: Dynamics of proteins in hydrated state and in solution as studied by dielectric relaxation spectroscopy.

    Macromolecules 2005, 38(26):10812-10819. Publisher Full Text OpenURL

  44. Yang L, Dordick JS, Garde S: Hydration of enzyme in nonaqueous media is consistent with solvent dependence of its activity.

    Biophysical Journal 2004, 87(2):812-821. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  45. Affleck R, Haynes CA, Clark DS: Solvent Dielectric Effects on Protein Dynamics.

    P Natl Acad Sci USA 1992, 89(11):5167-5170. Publisher Full Text OpenURL

  46. Ottosson J, Fransson L, King JW, Hult K: Size as a parameter for solvent effects on Candida antarctica lipase B enantioselectivity.

    Biochim Biophys Acta 2002, 1594(2):325-334. PubMed Abstract | Publisher Full Text OpenURL

  47. Hartsough DS, Merz KM: Protein Dynamics and Solvation in Aqueous and Nonaqueous Environments.

    J Am Chem Soc 1993, 115(15):6529-6537. Publisher Full Text OpenURL

  48. Nakasako M: Water-protein interactions from high-resolution protein crystallography.

    Philos T Roy Soc B 2004, 359(1448):1191-1204. Publisher Full Text OpenURL

  49. Teeter MM: Water-Protein Interactions - Theory and Experiment.

    Annu Rev Biophys Bio 1991, 20:577-600. Publisher Full Text OpenURL

  50. Brunne RM, Liepinsh E, Otting G, Wuthrich K, van Gunsteren WF: Hydration of Proteins - a Comparison of Experimental Residence Times of Water-Molecules Solvating the Bovine Pancreatic Trypsin-Inhibitor with Theoretical-Model Calculations.

    J Mol Biol J Mol Biol 1993, 231(4):1040-1048. Publisher Full Text OpenURL

  51. Bhattacharyya SM, Wang ZG, Zewail AH: Dynamics of water near a protein surface.

    Journal of Physical Chemistry B 2003, 107(47):13218-13228. Publisher Full Text OpenURL

  52. Deshpande A, Nimsadkar S, Mande SC: Effect of alcohols on protein hydration: crystallographic analysis of hen egg-white lysozyme in the presence of alcohols.

    Acta Crystallogr D 2005, 61:1005-1008. PubMed Abstract | Publisher Full Text OpenURL

  53. Gorman LAS, Dordick JS: Organic-Solvents Strip Water Off Enzymes.

    Biotechnology and Bioengineering 1992, 39(4):392-397. Publisher Full Text OpenURL

  54. Smolin N, Oleinikova A, Brovchenko I, Geiger A, Winter R: Properties of spanning water networks at protein surfaces.

    Journal of Physical Chemistry B 2005, 109(21):10995-11005. Publisher Full Text OpenURL

  55. Oleinikova A, Smolin N, Brovchenko I, Geiger A, Winter R: Formation of spanning water networks on protein surfaces via 2D percolation transition.

    Journal of Physical Chemistry B 2005, 109(5):1988-1998. Publisher Full Text OpenURL

  56. Franks F, Eagland D: The role of solvent interactions in protein conformation.

    Critical Reviews in Biochemistry 1975, 3(2):165-219. PubMed Abstract | Publisher Full Text OpenURL

  57. Frisch MJ, Trucks GW, Schlegel HB, Scuseria GE, Robb MA, Cheeseman JR, Jr JAM, Vreven T, Kudin KN, Burant JC, Millam JM, Iyengar SS, Tomasi J, Barone V, Mennucci B, Cossi M, Scalmani G, Rega N, Petersson GA, Nakatsuji H, Hada M, Ehara M, Toyota K, Fukuda R, Hasegawa J, Ishida M, Nakajima T, Honda Y, Kitao O, Nakai H, Klene M, Li X, Knox JE, Hratchian HP, Cross JB, Adamo C, Jaramillo J, Gomperts R, Stratmann RE, Yazyev O, Austin AJ, Cammi R, Pomelli C, Ochterski JW, Ayala PY, Morokuma K, Voth GA, Salvador P, Dannenberg JJ, Zakrzewski VG, Dapprich S, Daniels AD, Strain MC, Farkas O, Malick DK, Rabuck AD, Raghavachari K, Foresman JB, Ortiz JV, Cui Q, Baboul AG, Clifford S, Cioslowski J, Stefanov BB, Liu G, Liashenko A, Piskorz P, Komaromi I, Martin RL, Fox DJ, Keith T, Al-Laham MA, Peng CY, Nanayakkara A, Challacombe M, Gill PMW, Johnson B, Chen W, Wong MW, Gonzalez C, Pople JA: GAUSSIAN 98 (Revision A.7). Gaussian, Inc., Pittsburgh PA; 1998.

  58. Bayly CI, Cieplak P, Cornell WD, Kollman PA: A Well-Behaved Electrostatic Potential Based Method Using Charge Restraints for Deriving Atomic Charges - the Resp Model.

    J Phys Chem-Us 1993, 97(40):10269-10280. Publisher Full Text OpenURL

  59. Case DA, Cheatham TE, Darden T, Gohlke H, Luo R, Merz KM, Onufriev A, Simmerling C, Wang B, Woods RJ: The Amber biomolecular simulation programs.

    Journal of Computational Chemistry 2005, 26(16):1668-1688. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  60. Cornell WD, Cieplak P, Bayly CI, Gould IR, Merz KM, Ferguson DM, Spellmeyer DC, Fox T, Caldwell JW, Kollman PA: A 2nd Generation Force-Field for the Simulation of Proteins, Nucleic-Acids, and Organic-Molecules.

    J Am Chem Soc 1995, 117(19):5179-5197. Publisher Full Text OpenURL

  61. Berendsen HJC, Postma JP, van Gunsteren WF, DiNola A, Haak JR: Molecular dynamics with coupling to an external bath.

    J Chem Phys 1984, 81:3684-3690. Publisher Full Text OpenURL

  62. Ryckaert JP, Ciccotti G, Berendsen HJC: Numerical-Integration of Cartesian Equations of Motion of a System with Constraints - Molecular-Dynamics of N-Alkanes.

    J Comput Phys 1977, 23(3):327-341. Publisher Full Text OpenURL

  63. You TJ, Bashford D: Conformation and hydrogen ion titration of proteins: A continuum electrostatic model with conformational flexibility.

    Biophysical Journal 1995, 69(5):1721-1733. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  64. Martel PJ, Baptista A, Petersen SB: Protein electrostatics.

    Biotechnology Annual Review 1996, 2:315-372. OpenURL

  65. Miteva MA, Tuffery P, Villoutreix BO: PCE: web tools to compute protein continuum electrostatics.

    Nucleic Acids Res 2005, 33(Web Server issue):W372-5. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  66. Sitkoff D, Sharp KA, Honig B: Accurate Calculation of Hydration Free-Energies Using Macroscopic Solvent Models.

    J Phys Chem-Us 1994, 98(7):1978-1988. Publisher Full Text OpenURL

  67. Bashford D, Karplus M: Multiple-Site Titration Curves of Proteins - an Analysis of Exact and Approximate Methods for Their Calculation.

    J Phys Chem-Us 1991, 95(23):9556-9561. Publisher Full Text OpenURL

  68. Petersen SB, Jonson PH, Fojan P, Petersen EI, Petersen MTN, Hansen S, Ishak RJ, Hough E: Protein engineering the surface of enzymes.

    Journal of Biotechnology 1998, 66(1):11-26. PubMed Abstract | Publisher Full Text OpenURL

  69. Jorgensen WL, Chandrasekhar J, Madura JD, Impey RW, Klein ML: Comparison of Simple Potential Functions for Simulating Liquid Water.

    Journal of Chemical Physics 1983, 79(2):926-935. Publisher Full Text OpenURL

  70. Caldwell JW, Kollman PA: Structure and Properties of Neat Liquids Using Nonadditive Molecular-Dynamics - Water, Methanol, and N-Methylacetamide.

    J Phys Chem-Us 1995, 99(16):6208-6219. Publisher Full Text OpenURL

  71. Jorgensen WL, Briggs JM, Contreras ML: Relative Partition-Coefficients for Organic Solutes from Fluid Simulations.

    J Phys Chem-Us 1990, 94(4):1683-1686. Publisher Full Text OpenURL

  72. Kabsch W, Sander C: Dictionary of Protein Secondary Structure - Pattern-Recognition of Hydrogen-Bonded and Geometrical Features.

    Biopolymers 1983, 22(12):2577-2637. PubMed Abstract | Publisher Full Text OpenURL

  73. Eisenberg D., Schwarz E., Komarony M., R. W: Amino acid scale: Normalized consensus hydrophobicity scale.

    J Mol Biol 1984, 179:125-142. PubMed Abstract | Publisher Full Text OpenURL

  74. DeLano WL: The PyMOL Molecular Graphics System. Palo Alto, CA, USA , DeLano Scientific; 2002.

  75. Lide DR: CRC Handbook of Chemistry and Physics. 86th edition. London , CRC Press; 2005.