- Research article
- Open Access
Modeling of solvent-dependent conformational transitions in Burkholderia cepacia lipase
BMC Structural Biologyvolume 9, Article number: 38 (2009)
The characteristic of most lipases is the interfacial activation at a lipid interface or in non-polar solvents. Interfacial activation is linked to a large conformational change of a lid, from a closed to an open conformation which makes the active site accessible for substrates. While for many lipases crystal structures of the closed and open conformation have been determined, the pathway of the conformational transition and possible bottlenecks are unknown. Therefore, molecular dynamics simulations of a closed homology model and an open crystal structure of Burkholderia cepacia lipase in water and toluene were performed to investigate the influence of solvents on structure, dynamics, and the conformational transition of the lid.
The conformational transition of B. cepacia lipase was dependent on the solvent. In simulations of closed B. cepacia lipase in water no conformational transition was observed, while in three independent simulations of the closed lipase in toluene the lid gradually opened during the first 10–15 ns. The pathway of conformational transition was accessible and a barrier was identified, where a helix prevented the lid from opening to the completely open conformation. The open structure in toluene was stabilized by the formation of hydrogen bonds.
In simulations of open lipase in water, the lid closed slowly during 30 ns nearly reaching its position in the closed crystal structure, while a further lid opening compared to the crystal structure was observed in toluene. While the helical structure of the lid was intact during opening in toluene, it partially unfolded upon closing in water. The closing of the lid in water was also observed, when with eight intermediate structures between the closed and the open conformation as derived from the simulations in toluene were taken as starting structures. A hydrophobic β-hairpin was moving away from the lid in all simulations in water, which was not observed in simulations in toluene. The conformational transition of the lid was not correlated to the motions of the β-hairpin structure.
Conformational transitions between the experimentally observed closed and open conformation of the lid were observed by multiple molecular dynamics simulations of B. cepacia lipase. Transitions in both directions occurred without applying restraints or external forces. The opening and closing were driven by the solvent and independent of a bound substrate molecule.
Burkholderia cepacia lipase (BCL) formerly known as Pseudomonas cepacia lipase has been shown to be a useful enzyme which catalyzes a broad range of different reactions in water and non-polar solvents under mild conditions. BCL is a highly selective catalyst for a broad range of substrates , including the kinetic resolution of racemic mixtures of secondary alcohols by hydrolysis in water or esterification in organic solvents [2–4]. All lipases have a similar architecture, the α/β hydrolase fold [5, 6], which consists of parallel β-strands flanked by α-helices, and an active site with a catalytic triad consisting of Ser, His, and Asp/Glu [6, 7]. The characteristic of most lipases is their activation upon binding to a hydrophobic interface . Most lipases consist of a mobile element at the surface, a lid, which covers the active site . The lid is opening at a hydrophobic interface, making the active site accessible for substrates and enhancing the activity of the lipase [10, 11]. Because the exterior of the lid is hydrophilic and its interior is hydrophobic, the hydrophobic surface of lipases increases upon lid opening . Therefore, in non-polar solvents the open conformation was expected to be thermodynamically favored .
The crystal structure of BCL has been determined without a bound inhibitor in an open conformation of the lid. The lid consists of a helix-loop-helix motif with helix α4 (residues 118–127) and helix α5 (residues 134–150) which are linked to helix α6 (residues 160–166) . BCL contains an essential Ca2+-ion site  which stabilizes a β-hairpin (residues 214–228) . While no experimental structure of BCL in the closed conformation is available, a homologous lipase from Burkholderia glumae (BGL) [PDB: 1QGE]  with a sequence identity of 84% has been crystallized in the closed conformation. The cores of both homologous lipases have an identical structure with a root-mean squared deviation (RMSD) of backbone atoms of 0.7 Å, while the main structural difference between the two conformations is the position of the lid (residues 118–150) with a RMSD of 5 Å. The length of helix α5 in BCL (residues 134–150) was identified to be longer compared to BGL (residues 136–150), and the length of helix α6 in BCL (residues 160–166) was shorter in BGL (residues 156–166). Thus, the X-ray structures of open BCL and closed BGL suggests that the interfacial activation of BCL involves a motion of helices α4 and α5  and a partial unfolding of helix α6 in BGL. However, the mechanism of lid opening is still discussed. Several mechanisms of interfacial activation of lipases have been supposed. The substrate theory suggests that the lid opening is induced by binding of a substrate [18, 19], while the enzyme theory suggests that the lid opens by adsorption to a hydrophobic interface . Despite the fact that structures of the closed [9, 13, 17, 20] and the open conformation [20–25] of several lipases are known, a complete description of the mechanism of lid opening and closing is still lacking. The identification of mobile elements was possible by the availability of open and closed structures, but attempts to obtain any intermediates between the closed and open conformation  and conformational transitions have not been successful in experiment. Therefore, short molecular dynamics simulations of different lipases have been reported to get more insight into the conformational transition of the lid [26–31]. However, the time-scale of conformational transition of the lid from the closed to the fully open conformation was not reached in short molecular dynamics simulations. Therefore, methods like essential dynamics , normal mode analysis , and Brownian dynamics simulations  were used to study the lid opening in non-polar solvents.
To investigate the solvent-dependent motions of the mobile lid (residues 118–150) and the β-hairpin (residues 214–228) of BCL, multiple molecular dynamics simulations of the open crystal structure and the closed homology model of BCL were performed in water and toluene without using restraints or external forces. Toluene was used as a typical non-polar solvent which is widely used in lipase-catalyzed conversions. It has been shown that BCL is stable in toluene without significant loss of activity. The details of the conformational pathway were examined, to distinguish whether the transition occurs via a rigid-body movement or via a partial or complete unfolding, and bottlenecks in the lid opening and closing were identified.
Three simulations of 15–30 ns with different starting velocity distributions were performed for four enzyme-solvent systems, the closed homology model and the open crystal structure of BCL in water and toluene. The observed conformational transitions of the lid were dependent on the solvent, while all three simulations of the same enzyme-solvent system showed essentially the same transition. In eleven simulations, the Ca2+-binding site was stable, while in one simulation of closed BCL in water the Ca2+-ion left the binding pocket after 26 ns of simulation time.
Analysis of crystal structure and homology model
All available crystal structures of BCL have been obtained in their open conformation [14, 23, 24, 35, 36] showing only small root mean squared deviations (RMSD) of backbone atoms in the range of 0.2–0.5 Å. In all crystal structures there are crystal contacts between helix α4 (residues 118–127) and helix α9 (residues 243–256) of a second monomer, as well as between residues 217–218 and 223–225 of the β-hairpin and residues 130–138 of helix α5 of the lid of a second monomer (additional file 1). The crystal contacts are stabilized by electrostatic interactions between the negative potential of the lid and the positive potential of the β-hairpin (data not shown). In addition, hydrophobic contacts were observed between helix α5 and its equivalent of a second monomer at residues 142 to 149. In contrast, there are no crystal contacts between the lid and the β-hairpin in the crystal structure in the closed conformation of the homologous Burkholderia glumae lipase (BGL) [PDB: 1QGE] , which was used as template for the homology model of closed BCL. The sequences of BCL and BGL are highly similar with a sequence identity of 84%. In the lid, only residues 125, 128, 129, 148, and 150 are different.
The homology model of closed BCL showed no significant differences to the core of the open crystal structure, but was different in the conformation of the lid. The major movement was expected for helix α5 and the loop between helices α4 and α5, while helix α4 slightly changed its orientation. Helix α5 of BCL increased in length in the open conformation (residues 134–150) as compared to the closed conformation (residues 136–150), while helix α6 decreased in length in the open conformation (residues 160–166) as compared to the closed conformation (residue 156–166) of BCL (figure 1), as already been observed in the closed conformation of BGL. Furthermore, small differences can be observed in the position of the β-hairpin and in other loops (residues 17–27 and residues 50–59). The quality of the closed homology model of BCL was analyzed by ProSA showing a Z-score of -5.9 and a local model quality similar to the template BGL with a Z-score of -6.7. In contrast, the Z-scores of both closed structures were slightly better as compared to the open crystal structure of BCL showing a Z-score of -4.9 and a local model quality with high energy at helix α5 of the lid.
The conformation of the lid was characterized by two distances: the movement of the whole lid was described by the distance between atom Cα of residue 138 in helix α5 to atom Cα of residue 250 in helix α9, which was 6.1 Å in the closed homology model and 24.5 Å in the open crystal structure. The movement of the loop between helices α4 and α5 was described by the distance of atom Cα of residue 132 from the closed to the open conformation, which was 11.1 Å between the homology model and the crystal structure.
Simulation of closed BCL in water
In all three simulations of closed BCL in water no significant conformational change was observed. The backbone RMSD of the conformation after 30 ns to the initial structure was 1.3 Å of the core and 2.0 Å of the lid. The core was stable, deviations were mainly observed at surface loops, especially at the β-hairpin (residues 214–228) and the loop between helices α4 and α5 of the lid (residues 129–133). The lid was slightly closing as indicated by the distance Cα 138-Cα 250 which decreased from 6.4 Å to 6.1 Å and the movement of the loop between helices α4 and α5 of the lid by a distance of 5.7 Å between atom Cα of residue 132 from the simulated to the initial structure. The position of helix α5 (residues 134–150) was unchanged, helix α4 (residues 118–127) showed a movement by 2.7 Å at its C-cap, near the loop between helices α4 and α5.
The calculated B-factors per residue identified helix α5, the loop between helix α4 and α5, the β-hairpin, helix α6 (residues 156–159), helix α7 (residues 170–178) and five further surface loops (loop 23–28, loop 48–61, loop 184–193, loop 198–203, loop 232–240) as the most flexible regions (figure 2). While the lid showed only a small movement, the flexible hydrophobic β-hairpin moved away from the lid by up to 9.7 Å (average of three simulations: 7.6 Å), measured by the distance of atom Cα 221 to its initial structure. Thus, the movement of the β-hairpin was independent of the movement of the lid.
Simulation of open BCL in toluene
In all three simulations of open BCL in toluene, the core was stable, while the lid opened slightly more (figure 3), as compared to the crystal structure. The backbone RMSD of the conformation after 30 ns to the initial structure was 1.0 Å of the core and 2.4 Å of the lid. The distance Cα 138-Cα 250 increased by 3 Å from 24.3 Å to 27.3 Å, indicating a further lid opening. The distance of atom Cα of residue 132 in the loop between helix α4 and α5 from the initial structure to the simulated structure was 10.7 Å, indicating a major movement of the loop. In contrast to the simulations of the closed structure in water, the flexible β-hairpin did not shift. During the further lid opening, the hydrophobic surface increased by 80 Å2 as compared to the crystal structure. The regions with highest flexibility as indicated by the calculated B-factors per residue were the lid, the β-hairpin, and further surface loops (residues 23–28, 185–189, 198–203, 232–240) (figure 4).
Simulation of closed BCL in toluene
In all three simulations of closed BCL in toluene, the lid gradually opened [additional file 2], while the core of the protein was stable. The backbone RMSD of the conformation after 30 ns to the initial structure was 1.3 Å of the core and 4.8 Å of the lid. The RMSD to the initial structure continuously increased during the first 13 ns of simulation (figure 5) by a fast opening of the lid (figure 6), while the distance Cα 138-Cα 250 increased by 17 Å. Between 13 and 22 ns the lid slightly opened further by a distance Cα 138-Cα 250 of 2 Å, approaching its position in the nearly open conformation of the crystal structure until reaching a barrier (additional file 3). After gradually opening, the movement of the lid was blocked by helix α6 at a distance Cα 138-Cα 250 of 23.2 Å, near to its position in the crystal structure with a distance Cα 138-Cα 250 of 24.5 Å. The movement of helix α5 was blocked by the N-cap of helix α6 (residues 156 to 159) (figure 7) and it was not possible to overcome this barrier until the end of simulations. While in the simulation residues 156 to 159 were helical as in the homology model of the closed structure, helix α6 (residues 160–166) in the crystal structure of the open conformation was partially unfolded near its N-cap, leaving residues 156 to 159 in a loop structure (additional file 4), which would allow the further movement of helix α5.
In all simulations in toluene, a salt bridge was formed between Asp121 at the N-cap of helix α4 and Arg258, which was not observed in the crystal structures of BCL and the homology model derived from BGL. This salt bridge prevented helix α4 from a movement in toluene. The opening of the lid was a hinge-type motion of a single helix α5 and the concomitant stretching of the flexible loop between helix α4 and α5. During the conformational transition the secondary structure of helix α5 was conserved. Thus, the lid opening involved conformational rearrangements of several secondary structure elements: a rigid body movement of helix α5 and a movement of the loop between helix α4 and α5, while helix α4 showed no significant movement (figure 8). The solvent accessible surface area (SASA) of the open crystal structure (13309 Å2) was 300 Å2 larger than the SASA of the closed conformation of the homology model (12978 Å2). During lid opening, the SASA increased by 800 Å2 to 13803 Å2. While the hydrophilic SASA decreased by 300 Å2, the hydrophobic SASA increased by 1100 Å2 exposing a large hydrophobic patch of the lid to the non-polar solvent (residues Phe119, Val123, Val126, Val138, Ile139, Phe142, Val143, Val145, Phe146, and Leu149). The increase in hydrophobic SASA mainly occurred during the first 13 ns in correlation with lid opening (figure 9). While in the closed structure this hydrophobic patch was shielded from the solvent, making the protein surface more hydrophilic (additional file 5), in the open conformation the hydrophobic interior of the lid was exposed to the solvent.
The β-hairpin was very flexible, but showed no directed movement in contrast to the simulations in water.
Simulation of open BCL in water
In the three simulations of open BCL in water, a partial lid closing was observed, while the core of the protein was stable. The backbone RMSD of the conformation after 30 ns to the initial structure of open BCL was 0.7 Å of the core and 3.1 Å of the lid. The distance Cα 138-Cα 250 decreased by 5.8 Å, from 24.4 Å to 18.6 Å. The distance of atom Cα of residue 132 in the loop between helix α4 and α5 to its position in the open crystal structure increased to 12.3 Å. However, the fully closed conformation was not reached. Predominantly helix α5 was moving, while helices α4, α6, and α9 showed no movement. The movement of helix α5 was very slow and was not blocked by other structural elements. In contrast to the opening in toluene, helix α5 partially unfolded at its N-cap (residues 134–138) and changed from a straight to a more curved structure. As in the simulations of closed BCL in water, the β-hairpin showed a movement away from the lid by up to 7.2 Å (average of three simulations: 5.3 Å).
To investigate whether the solvent is the only driving force for the observed conformational transition, eight simulations in water were performed starting with intermediate structures taken from the simulation of the lid opening in toluene. Conformers were taken at 6, 9, 12, 15, 18, 21, 24, and 27 ns and simulated for 6 ns after changing the solvent from toluene to water (figure 10A). In all simulations, the lid was partially closing and the distance Cα 138-Cα 250 changed between 1.9 to 14.2 Å (table 1), indicating that the conformational transition can be reversed by changing the solvent. In one simulation, using the conformer after 15 ns as starting structure, the lid closed completely after 6 ns. There was no blocking of the lid by other structural elements in the simulations. It is quite remarkable, that for all simulations in water four preferred distances Cα 138-Cα 250 were identified at 6, 10, 14 and 16.5 Å (figure 10B).
Coordination of Asp130
In the crystal structures of BCL and BGL, the carboxylic group of the side chain of Asp130, located in the loop between helix α4 and α5, is coordinated by hydrogen bonds between the backbone and the side chain of Thr132, and the backbone of Ser135, stabilizing the conformation of the loop. During the simulations of open and closed BCL in water, the side chain of Asp130 was also coordinated by the backbone and the side chain of Thr132, and the backbone of Ser135, but these hydrogen bonds were partially lost during the simulation and replaced by hydrogen bonds to water molecules of the solvent. However, during simulations of closed BCL in toluene, Asp130 formed additional hydrogen bonds during lid opening and became fourfold coordinated by the backbones and side chains of Thr132 and Ser135. Furthermore, in the simulations of open BCL in toluene, Asp130 was sixfold coordinated forming two additional hydrogen bonds to backbone and side chain of Thr136 (figure 11). Thus, in toluene the loop between helix α4 and α5 was stabilized in the open conformation by the coordination of Asp130 in an extended hydrogen bonding to other amino acids of the loop.
The interfacial activation of lipases is associated with a conformational change of the lid [11, 37]. The detailed analysis of the pathway of the conformational transition is not accessible in experiment. To get insight into the interfacial activation of BCL, multiple molecular dynamics simulations of the open and closed conformation of BCL in explicit solvent models of water and toluene were carried out without using restraints or external forces. An implicit organic solvent model was not used, because essential water molecules are supposed to be important to retain the flexibility of the lipase. Reproducibility of the observed properties was confirmed by performing three simulations for each system. The primary purpose was to apply molecular dynamics simulations examine the motion and conformational transition of the lid (residues 118–150) and the β-hairpin (residues 214–228) of BCL, which is not yet accessible to an experimental investigation.
For each system, a similar behaviour was observed: a fast opening of the lid in toluene starting with the closed structure, a closing of the open structure in water, and no major conformational changes in the simulation of the closed and open structure in water and toluene, respectively. Our major observations, the solvent-dependent opening and closing of the lid and the decoupling of the motions of lid and β-hairpin, were confirmed in an independent work in the group of P. Monsan (INSA, Toulouse) on simulation of conformational transitions of BCL: while lid opening occured at a water-hexane interface, closing of the lid was observed in water (I. André, personal communication).
Simulation of closed BCL in water
In simulations of closed BCL in water, the lid stayed in its closed conformation. A slight further lid closing was observed which indicates that the homology model of closed BCL might not represent the closed conformation in solution. Furthermore, the lid showed a high conformational mobility in the absence of an oil-water interface which is supported by experimental observations. High B-factors of the lid were observed in all lipases in closed conformation [38–40], which led initially to the hypothesis that the lid of BGL can undergo a conformational transition. In crystallographic studies of closed Thermomyces lanuginosa lipase (TLL) in water, the lid was disordered , and in Rhizopus delemar lipase two different conformations of the closed lid were observed .
Simulation of open BCL in toluene
In simulations of open BCL in toluene, the lid further opened. This was also observed in Brownian dynamics simulations of Rhizomucor miehei lipase (RML) . At a lipid binding zone, which is comparable to a non-polar solvent, the lid further opened starting from the open crystal structure . The further lid opening might be due to a different stabilization of the open structure in non-polar solvents as compared to the crystal structure. While in non-polar solvents, the lid is stabilized by an increase of hydrophobic surface, the conformation of the lid in the crystal is the result of energetic compromises to form stable contact interactions . Upon molecular dynamics simulations of monomers in solution, starting with crystal structures, the structure is released from its crystal contacts and relaxed to a conformation close to that in solution . The overall flexibility in simulations of open BCL in toluene was about half the flexibility of closed BCL in water, which is in agreement with a general observation of a decreased flexibility of proteins in non-polar solvents in experiment [46, 47] and in molecular dynamics simulations [26, 48–51].
Simulation of closed BCL in toluene
In our simulations of closed BCL in toluene, the lid opening by 18 Å was observed by a rigid-body movement of helix α5 without unfolding of helices α4 and α5. A similar rigid-body hinge-type motion of a single helix was also concluded from the existence of a closed and open crystal structure of RML . The lid opening of RML was observed in a low-dielectric medium, whereas in a high-dielectric medium no opening was observed . In our simulations the hydrophobic surface of BCL increased upon lid opening, which is in agreement with a general observation of most lipases, where lid opening is associated with an increase in hydrophobic surface [12, 20, 52, 53].
Previously, a partial opening of the lid has been investigated in simulations for different lipases. In molecular dynamics simulations in a continuum with different dielectric constants, lid opening was enhanced in a medium with a low dielectric constant, while the opening time of the lid increased for media with increasing dielectric constants . In previous Brownian dynamics simulations, the lid of RML opened during 100 ns in a low-dielectric medium, while in a high-dielectric medium the lid stayed in its closed conformation for 900 ns . In a 1 ns molecular dynamics simulation of RML in explicit water, the lid was highly flexible and a lid opening was concluded from essential dynamics analysis . The partial lid opening was also observed for RML , Humicola lanuginosa lipase (HLL) , and Candida rugosa lipase (CRL)  in explicit solvents. In short molecular dynamics simulations of open BCL in water and vacuum, no conformational transition of the lid and a similar behaviour in different solvent environments was observed . It was suggested that lid opening of HLL is a multi-step process involving more than two conformational transitions before a fully activated conformation is assumed . This is in agreement with our simulations of BCL, where a barrier was identified, blocking the movement of helix α5 by the N-cap of helix α6. While the generated homology model of closed BCL and the template of closed BGL show a long helix α6 (residues 156 to 166), helix α6 is shorter in the open structure of BCL (residues 160 to 166), where Asp157 is replaced by Asn in the sequence of BCL. Thus, the unfolding of the N-cap of helix α6 (residues 156 to 159) would be necessary to make the complete movement of the lid towards the fully open structure feasible. However, this unfolding was not observed in the time scale in our simulations of 30 ns. The blocking of helix α5 by helix α6 seems to result from a different secondary structure of helix α6 in the open and closed structure.
A correlation of the lid movement with a surface loop was previously observed in restrained simulations of Pseudomonas aeruginosa lipase, which was suggested to trigger the lid opening . In all our simulations of closed and open BCL in water, the β-hairpin drifted away from the lid, while it did not move in toluene, independent of the conformation of the lid. The conformational transition of the lid of BCL was not coupled to the β-hairpin or other surface loops, which was also not observed for other lipases [28, 30, 43, 55].
Different mechanisms of interfacial activation have been already proposed. While in the 'substrate theory', changes in the conformation of the lid are dependent on a bound substrate molecule [18, 19], in the 'enzyme theory' conformational changes on the enzyme are based on the adsorption to a lipid interface . The lid opening in our simulations was only driven by the hydrophobicity of the solvent, where a hydrophobic solvent favored the opening of the lid, which supports the 'enzyme theory' and experimental observations. However, the solvent might not be the only reason leading to a complete opening of the lid, and we cannot exclude a possible influence of a substrate intermediate covalently bound to the catalytic serine. Several lipases have been crystallized in an open conformation in the absence of a bound inhibitor [57, 14, 36]. The lid opening of Human pancreatic lipase by water-miscible organic solvents without a bound substrate was also observed in experiment using antibodies , and a hydrophobic patch led to a conformational transition to the open conformation of CRL .
Simulation of open BCL in water
The open conformation of BCL in water showed a partial lid closing during the simulations. However, lid closing in water was much slower than lid opening in toluene, and a complete movement of the lid to the closed conformation of the homology model was not observed. For CRL it was observed in experiment that the conversion between the open and closed conformation was very slow in aqueous solution in the absence of interfaces, which made the separation of the two conformations possible, while the exposure of CRL to a hydrophobic patch accelerated the transition . However, a transition state between the closed and the open conformation could not be isolated. Short molecular dynamics simulations of the open conformation of BCL in water and vacuum demonstrated the high flexibility of the lid and the β-hairpin . It was concluded from its high flexibility that the open structure of BCL becomes unstable in water leading to lid closing, however, a conformational transition from the open to the closed conformation was not observed.
The conformational transition of BCL from the closed to the open conformation of the lid in toluene could be reversed after changing the solvent from toluene to water. The analysis of eight simulations starting with intermediate positions of lid opening in toluene indicated that there are at least three local energy minima in addition to the closed crystal structure. In one simulation the lid closed completely, while in all other simulations a partially closing was observed. In contrast to the lid opening in non-polar solvents, the lid closing in water was a stepwise rather than a continuous process, associated with a partial unfolding and deformation of the helical structures in the lid.
Coordination of Asp130
The solvent was the driving force of the opening and closing of the lid. Local unfolding by breaking hydrogen bonds was mediated by water molecules, while the formation of hydrogen bonds was driven by the non-polar solvent. The formation of a local hydrogen bond network was most prominent in the loop between helix α4 and helix α5, which seems to direct the conformational transition. Asp 130 plays a pivotal role in the hydrogen bond formation and stabilizes the open conformation of the lid in toluene. While Asp130 was sixfold coordinated by backbone and side chain atoms after simulation in toluene, which was closely related to the movement of the loop, these hydrogen bonds were lost in water by coordination of water molecules of the solvent, leading to a destabilization of the open structure. A pivotal role of single amino acids directing the lid opening might not be restricted to BCL, as observed in previous simulations of TLL, where a solvent-dependent switch of Arg84 in the lid was observed .
The observed conformational transitions of the lid were dependent on the solvent. In simulations of closed BCL in water, no significant conformational change was observed, while in simulations of closed BCL in toluene the lid gradually opened. In simulations of open BCL in toluene, the lid opened slightly more than in the crystal structure, while in simulations of open BCL in water a partial closing of the lid was observed. In simulations of closed and open BCL in water, the hydrophobic β-hairpin moved away from the lid, while it showed no directed movement in simulations in toluene. Thus, the movement of the β-hairpin was independent of the movement of the lid.
The crystal structure of open BCL was taken from the Protein Data Bank [PDB: 3LIP]  with a resolution of 2.0 Å by X-ray diffraction as initial structure for the simulations. There was no crystal structure of the closed conformation of BCL available. A homology model of closed BCL was built by the protein modeling server SWISS-MODEL  using the closed conformation of the homologous lipase from Burkholderia glumae (BGL) [PDB: 1QGE]  as template.
Calculation of protonation states
pKa values and protonation states of titratable sites Arg, Lys, Asp, Glu and His at pH 7 were calculated using TITRA , based on the Tanford-Kirkwood model , using standard parameters. The solvent accessible surface area (SASA) of each residue was calculated by the program acc_run . In toluene BCL was protonated assuming pH memory from the protonation states in water [49, 63]. The electrostatic potential at pH 7 was calculated using DELPHI V. 4  with a dielectric constant of the solvent of 80 and a molecule interior of 4. Atomic charges and radii were taken from the PARSE charge and radii files .
Simulations were set up in XLEAP of the Amber 7.0 program package. Hydrogens were added as calculated by TITRA, His was protonated at Hδ-position. Disulfide bridges were built from information in the crystal structure. The BCL structure, including 193 crystal water molecules, was solvated using the explicit TIP3 water model (dielectric constant 78.5)  and a non-polar solvent model of toluene (dielectric constant 2.4)  in a truncated octahedral box with a minimal distance of 14 Å between the box boundary and the protein. Three Na+-ions were added in XLEAP as counter ions to neutralize BCL for simulations.
Molecular dynamics simulations
Multiple molecular dynamics simulations of the protein-solvent systems were performed using the AMBER 7 program package . The all-atom AMBER force-field ff99  was used to present the protein system and Ca2+-ions. The simulations were done in a truncated octahedral box under periodic boundary conditions. The Sander tool of AMBER 7 was used for minimization and simulations. Non-bonded interactions were calculated using the Particle-Mesh Ewald method to a cutoff distance of 10 Å and the SHAKE algorithm  was applied to constrain all bonds involving hydrogen atoms. The initial structures were energy minimized by applying 500 steps steepest descent and 50 steps conjugate gradient to relax clashes in the system. After the minimization steps the temperature was set to 300 K and the pressure to 1 bar under restrained conditions using a harmonic potential for all backbone atoms C, O and CA. The force constants were gradually decreased every 50 ps steps from 10, 5, 1 to 0.1 kcal/mol force constant, followed by a non-restrained simulation without additional external forces under periodic boundary conditions. Simulations were performed at 300 K using a time step of 1 fs. Temperature and pressure of the system were controlled using a weak coupling to an external heat bath  with a temperature coupling constant of 1.0 and a pressure coupling constant of 1.2. Three simulations, one up to 30 ns and two up to 15 ns, for the open and closed system in water and toluene were carried out using different starting atomic velocities from a Maxwell distribution to sample the conformational space. Each distribution was generated by a random number generator. The resulting trajectories were analyzed using PTRAJ of AMBER 7 after fitting the backbone atoms of the core (except residues 118–150 and residues 214–228) of each conformer to the initial structure. In PTRAJ the root-mean squared deviation (RMSD) of the backbone atoms between each conformer and the initial structure and between all conformers (2D-RMSD), atomic positional fluctuations relating to B-factors, and atom-atom distances were calculated.
Crystals were built up using SWISS-PDB Viewer and crystal contacts were analyzed by the WHAT IF web interface , including a 5.0 Å shell of symmetry related residues around the molecule, pairs of atoms were analyzed only between different asymmetric units. Secondary structure elements and the solvent accessible surface area (SASA) using a probe radius of 1.4 Å were calculated by DSSP . Protein structures were visualized using PyMol 0.98  and VMD . Electrostatic potentials and hydrophobicity, using the hydrophobicity scale of Eisenberg , were mapped on the surface in PyMol 0.98. The quality of the homology model was analyzed using the ProSA-web interface .
Kazlauskas RJ, Weber HK: Improving hydrolases for organic synthesis. Curr Opin Chem Biol 1998, 2: 121–126.
Bornscheuer UT, Kazlauskas RJ: Hydrolases in organic synthesis: regio- and stereolelective biotransformations. Weinheim: Wiley-VCH; 1999.
Schmid RD, Verger R: Lipases: Interfacial enzymes with attractive applications. Angew Chem, Int Edit 1998, 37: 1609–1633.
Jaeger KE, Eggert T: Lipases for biotechnology. Curr. Opin. Biotechnol 2002, 13: 390–397.
Schrag JD, Cygler M: Lipases and alpha/beta hydrolase fold. Lipases, Part A 1997, 284: 85–107.
Ollis DL, Cheah E, Cygler M, Dijkstra B, Frolow F, Franken SM, Harel M, Remington SJ, Silman I, Schrag J, et al.: The Alpha/Beta-Hydrolase Fold. Protein Eng 1992, 5: 197–211.
Warshel A, Narayszabo G, Sussman F, Hwang JK: How Do Serine Proteases Really Work. Biochemistry 1989, 28: 3629–3637.
Verger R: "Interfacial activation" of lipases: facts and artefacts. Trends Biotechnol 1997, 15: 32–38.
Brady L, Brzozowski AM, Derewenda ZS, Dodson E, Dodson G, Tolley S, Turkenburg JP, Christiansen L, Hugejensen B, Norskov L, et al.: A Serine Protease Triad Forms the Catalytic Center of a Triacylglycerol Lipase. Nature 1990, 343: 767–770.
Sarda L, Desnuelle P: Action De La Lipase Pancreatique Sur Les Esters En Emulsion. Biochim Biophys Acta 1958, 30: 513–521.
Verger R: Enzyme kinetics of lipolysis. Methods Enzymol 1980, 64: 340–392.
Vantilbeurgh H, Egloff MP, Martinez C, Rugani N, Verger R, Cambillau C: Interfacial Activation of the Lipase Procolipase Complex by Mixed Micelles Revealed by X-Ray Crystallography. Nature 1993, 362: 814–820.
Derewenda ZS: Structure and Function of Lipases. Adv Protein Chem 1994, 45: 1–52.
Schrag JD, Li Y, Cygler M, Lang D, Burgdorf T, Hecht HJ, Schmid R, Schomburg D, Rydel TJ, Oliver JD, et al.: The open conformation of a Pseudomonas lipase. Structure 1997, 5: 187–202.
El Khattabi M, Van Gelder P, Bitter W, Tommassen J: Role of the calcium ion and the disulfide bond in the Burkholderia glumae lipase. J Mol Catal B:Enzym 2003, 22: 329–338.
Tanaka A, Sugimoto H, Muta Y, Mizuno T, Senoo K, Obata H, Inouye K: Differential scanning Calorimetry of the effects of Ca2+ on the thermal unfolding of Pseudomonas cepacia lipase. Biosci Biotechnol Biochem 2003, 67: 207–210.
Lang D, Hofmann B, Haalck L, Hecht HJ, Spener F, Schmid RD, Schomburg D: Crystal structure of a bacterial lipase from Chromobacterium viscosum ATCC 6918 refined at 1.6 angstrom resolution. J Mol Biol 1996, 259: 704–717.
Thuren T: A Model for the Molecular Mechanism of Interfacial Activation of Phospholipase-A2 Supporting the Substrate Theory. FEBS Lett 1988, 229: 95–99.
Brockman HL, Law JH, Kezdy FJ: Catalysis by Adsorbed Enzymes – Hydrolysis of Tripropionin by Pancreatic Lipase Adsorbed to Siliconized Glass Beads. J Biol Chem 1973, 248: 4965–4970.
Grochulski P, Li Y, Schrag JD, Cygler M: 2 Conformational States of Candida-Rugosa Lipase. Protein Sci 1994, 3: 82–91.
Brzozowski AM, Derewenda U, Derewenda ZS, Dodson GG, Lawson DM, Turkenburg JP, Bjorkling F, Hugejensen B, Patkar SA, Thim L: A Model for Interfacial Activation in Lipases from the Structure of a Fungal Lipase-Inhibitor Complex. Nature 1991, 351: 491–494.
Nardini M, Lang DA, Liebeton K, Jaeger KE, Dijkstra BM: Crystal structure of Pseudomonas aeruginosa lipase in the open conformation – The prototype for family I.1 of bacterial lipases. J Biol Chem 2000, 275: 31219–31225.
Luic M, Tomic S, Lescic I, Ljubovic E, Sepac D, Sunjic V, Vitale L, Saenger W, Kojic-Prodic B: Complex of Burkholderia cepacia lipase with transition state analogue of 1-phenoxy-2-acetoxybutane – Biocatalytic, structural and modelling study. Eur J Biochem 2001, 268: 3964–3973.
Mezzetti A, Schrag JD, Cheong CS, Kazlauskas RJ: Mirror-image packing in enantiomer discrimination: Molecular basis for the enantioselectivity of B-cepacia lipase toward 2-methyl-3-phenyl-1-propanol. Chem Biol 2005, 12: 427–437.
Brzozowski AM, Savage H, Verma CS, Turkenburg JP, Lawson DM, Svendsen A, Patkar S: Structural origins of the interfacial activation in Thermomyces (Humicola) lanuginosa lipase. Biochemistry 2000, 39: 15071–15082.
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: 548–559.
Peters GH, Bywater RP: Computational analysis of chain flexibility and fluctuations in Rhizomucor miehei lipase. Protein Eng 1999, 12: 747–754.
Lee J, Suh SW, Shin S: Computational studies of essential dynamics of Pseudomonas cepacia lipase. J Biomol Struct Dyn 2000, 18: 297–309.
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 Comput Biol 2005, 1: e28.
Jensen MO, Jensen TR, Kjaer K, Bjornholm T, Mouritsen OG, Peters GH: Orientation and conformation of a lipase at an interface studied by molecular dynamics simulations. Biophys J 2002, 83: 98–111.
Secundo F, Carrea G, Tarabiono C, Gatti-Lafranconi P, Brocca S, Lotti M, Jaeger KE, Puls M, Eggert T: The lid is a structural and functional determinant of lipase activity and selectivity. J Mol Catal B:Enzym 2006, 39: 166–170.
Peters GH, Toxvaerd S, Olsen OH, Svendsen A: Computational studies of the activation of lipases and the effect of a hydrophobic environment. Protein Eng 1997, 10: 137–147.
Jaaskelainen S, Verma CS, Hubbard RE, Linko P, Caves LSD: Conformational change in the activation of lipase: An analysis in terms of low-frequency normal modes. Protein Sci 1998, 7: 1359–1367.
Peters GH, Olsen OH, Svendsen A, Wade RC: Theoretical investigation of the dynamics of the active site lid in Rhizomucor miehei lipase. Biophys J 1996, 71: 119–129.
Lang DA, Mannesse MLM, De Haas GH, Verheij HM, Dijkstra BW: Structural basis of the chiral selectivity of Pseudomonas cepacia lipase. Eur J Biochem 1998, 254: 333–340.
Kim KK, Song HK, Shin DH, Hwang KY, Suh SW: The crystal structure of a triacylglycerol lipase from Pseudomonas cepacia reveals a highly open conformation in the absence of a bound inhibitor. Structure 1997, 5: 173–185.
Cygler M, Schrag JD: Structure as basis for understanding interfacial properties of lipases. Lipases, Part A 1997, 284: 3–27.
Noble ME, Cleasby A, Johnson LN, Egmond MR, Frenken LG: Analysis of the structure of Pseudomonas glumae lipase. Protein Eng 1994, 7: 559–562.
Mancheno JM, Pernas MA, Martinez MJ, Ochoa B, Rua ML, Hermoso JA: Structural insights into the lipase/esterase behavior in the Candida rugosa lipases family: Crystal structure of the lipase 2 isoenzyme at 1.97 angstrom resolution. J Mol Biol 2003, 332: 1059–1069.
Angkawidjaja C, You DJ, Matsumura H, Kuwahara K, Koga Y, Takano K, Kanaya S: Crystal structure of a family I.3 lipase from Pseudomonas sp MIS38 in a closed conformation. FEBS Lett 2007, 581: 5060–5064.
Cajal Y, Svendsen A, Girona V, Patkar SA, Alsina MA: Interfacial control of lid opening in Thermomyces lanuginosa lipase. Biochemistry 2000, 39: 413–423.
Derewenda U, Swenson L, Wei YY, Green R, Kobos PM, Joerger R, Haas MJ, Derewenda ZS: Conformational Lability of Lipases Observed in the Absence of an Oil-Water Interface – Crystallographic Studies of Enzymes from the Fungi Humicola-Lanuginosa and Rhizopus-Delemar. J Lipid Res 1994, 35: 524–534.
Peters GH, Bywater RP: Influence of a lipid interface on protein dynamics in a fungal lipase. Biophys J 2001, 81: 3052–3065.
Kossiakoff AA, Randal M, Guenot J, Eigenbrot C: Variability of Conformations at Crystal Contacts in Bpti Represent True Low-Energy Structures – Correspondence among Lattice Packing and Molecular-Dynamics Structures. Proteins 1992, 14: 65–74.
Koumanov A, Karshikoff A, Friis EP, Borchert TV: Conformational averaging in pK calculations: Improvement and limitations in prediction of ionization properties of proteins. J Phys Chem B 2001, 105: 9339–9344.
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: 12657–12663.
Zaks A, Klibanov AM: Enzyme-catalyzed processes in organic solvents. Proc Natl Acad Sci USA 1985, 82: 3192–3196.
Colombo G, Ottolina G, Carrea G: Modelling of enzyme properties in organic solvents. Monatsh Chem 2000, 131: 527–547.
Soares CM, Teixeira VH, Baptista AM: Protein structure and dynamics in nonaqueous solvents: insights from molecular dynamics simulation studies. Biophys J 2003, 84: 1628–1641.
Micaelo NM, Soares CM: Modeling hydration mechanisms of enzymes in nonpolar and polar organic solvents. FEBS J 2007, 274: 2424–2436.
Trodler P, Pleiss J: Modeling structure and flexibility of Candida antarctica lipase B in organic solvents. BMC Struct Biol 2008, 8: 9.
Derewenda U, Brzozowski AM, Lawson DM, Derewenda ZS: Catalysis at the Interface – the Anatomy of a Conformational Change in a Triglyceride Lipase. Biochemistry 1992, 31: 1532–1541.
Miled N, De Caro A, De Caro J, Verger R: A conformational transition between an open and closed form of human pancreatic lipase revealed by a monoclonal antibody. Bba-Protein Struct M 2000, 1476: 165–172.
Peters GH, vanAalten DMF, Edholm O, Toxvaerd S, Bywater R: Dynamics of proteins in different solvent systems: Analysis of essential motion in lipases. Biophys J 1996, 71: 2245–2255.
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: 4377–4383.
Gunasekaran K, Ma BY, Nussinov R: Triggering loops and enzyme function: Identification of loops that trigger and modulate movements. J Mol Biol 2003, 332: 143–159.
Grochulski P, Li YG, Schrag JD, Bouthillier F, Smith P, Harrison D, Rubin B, Cygler M: Insights into Interfacial Activation from an Open Structure of Candida-Rugosa Lipase. J Biol Chem 1993, 268: 12843–12847.
Turner NA, Needs EC, Khan JA, Vulfson EN: Analysis of conformational states of Candida rugosa lipase in solution: implications for mechanism of interfacial activation and separation of open and closed forms. Biotechnol Bioeng 2001, 72: 108–118.
Arnold K, Bordoli L, Kopp J, Schwede T: The SWISS-MODEL workspace: a web-based environment for protein structure homology modelling. Bioinformatics 2006, 22: 195–201.
Martel PJ, Baptista A, Petersen SB: Protein electrostatics. Biotechnol Annu Rev 1996, 2: 315–372.
Havranek JJ, Harbury PB: Tanford-Kirkwood electrostatics for protein modeling. Proc Natl Acad Sci USA 1999, 96: 11145–11150.
Petersen MT, Martel P, Petersen EI, Drablos F, Petersen SB: Surface and electrostatics of cutinases. Methods Enzymol 1997, 284: 130–154.
Costantino HR, Griebenow K, Langer R, Klibanov AM: On the pH memory of lyophilized compounds containing protein functional groups. Biotechnol Bioeng 1997, 53: 345–348.
Honig B, Nicholls A: Classical electrostatics in biology and chemistry. Science 1995, 268: 1144–1149.
Sitkoff D, Sharp KA, Honig B: Accurate Calculation of Hydration Free-Energies Using Macroscopic Solvent Models. J Phys Chem 1994, 98: 1978–1988.
Jorgensen WL, Chandrasekhar J, Madura JD, Impey RW, Klein ML: Comparison of Simple Potential Functions for Simulating Liquid Water. J Chem Phys 1983, 79: 926–935.
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. J Comput Chem 2005, 26: 1668–1688.
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: 5179–5197.
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: 327–341.
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.
Vriend G: What If – a Molecular Modeling and Drug Design Program. J Mol Graphics 1990, 8: 52.
Kabsch W, Sander C: Dictionary of Protein Secondary Structure – Pattern-Recognition of Hydrogen-Bonded and Geometrical Features. Biopolymers 1983, 22: 2577–2637.
DeLano WL: The PyMOL Molecular Graphics System. DeLano Scientific, San Carlos, CA, USA; 2002.
Humphrey W, Dalke A, Schulten K: VMD: Visual molecular dynamics. J Mol Graphics 1996, 14: 33.
Eisenberg D, Schwarz E, Komarony M, Wall R: Amino acid scale: Normalized consensus hydrophobicity scale. J Mol Biol 1984, 179: 125–142.
Wiederstein M, Sippl MJ: ProSA-web: interactive web service for the recognition of errors in three-dimensional structures of proteins. Nucleic Acids Res 2007, 35: W407-W410.
The work was carried out in the framework of the Sonderforschungsbereich 716 funded by the Deutsche Forschungsgemeinschaft (DFG).
PT carried out the simulations, RDS contributed to the discussion of the results, JP was the principal investigator and directed the research. All authors read and approved the final manuscript.
Electronic supplementary material
Additional file 1: Crystal contacts of BCL. In the crystal structure of BCL [PDB: 3LIP]  crystal contacts were observed. The β-hairpin (residues 214–228) (red) is in contact to the lid (residues 118–150) (yellow) of the next monomer. (DOC 470 KB)
Additional file 3: 2D-RMSD in simulations of BCL. The root mean squared deviation of every conformation to all other conformations as a function of time during simulation of closed BCL in 30 ns simulation in toluene is shown in the 2D-RMSD. A stable conformation after lid opening is indicated by the green area. (DOC 978 KB)
Additional file 4: Barrier during lid opening. The lid opening of helix α5 was blocked by helix α6, indicated by the conformation at the end of the simulation of closed BCL in toluene (red), between the open conformation of the crystal structure (green) and the closed conformation of the homology model (blue). (DOC 501 KB)
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.