Modeling of solvent-dependent conformational transitions in Burkholderia cepacia lipase
© Trodler et al; licensee BioMed Central Ltd. 2009
Received: 17 October 2008
Accepted: 28 May 2009
Published: 28 May 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 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.
Simulation of open BCL in toluene
Simulation of closed BCL in toluene
Additional file 2: Interfacial activation of BCL. Lid opening of closed BCL during interfacial activation in a 30 ns simulation of BCL in toluene. (WMV 5 MB)
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 Å).
Reversed movement of the lid
Conformation of BCL in toluene after [ns]
Distance C α 138-C α 250 before simulation [Å]
Distance C α 138-C α 250 after simulation [Å]
Δ distance C α 138-C α 250 [Å]
Coordination of Asp130
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 .
The work was carried out in the framework of the Sonderforschungsbereich 716 funded by the Deutsche Forschungsgemeinschaft (DFG).
- Kazlauskas RJ, Weber HK: Improving hydrolases for organic synthesis. Curr Opin Chem Biol 1998, 2: 121–126.View ArticlePubMedGoogle Scholar
- Bornscheuer UT, Kazlauskas RJ: Hydrolases in organic synthesis: regio- and stereolelective biotransformations. Weinheim: Wiley-VCH; 1999.Google Scholar
- Schmid RD, Verger R: Lipases: Interfacial enzymes with attractive applications. Angew Chem, Int Edit 1998, 37: 1609–1633.View ArticleGoogle Scholar
- Jaeger KE, Eggert T: Lipases for biotechnology. Curr. Opin. Biotechnol 2002, 13: 390–397.View ArticlePubMedGoogle Scholar
- Schrag JD, Cygler M: Lipases and alpha/beta hydrolase fold. Lipases, Part A 1997, 284: 85–107.View ArticleGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- Warshel A, Narayszabo G, Sussman F, Hwang JK: How Do Serine Proteases Really Work. Biochemistry 1989, 28: 3629–3637.View ArticlePubMedGoogle Scholar
- Verger R: "Interfacial activation" of lipases: facts and artefacts. Trends Biotechnol 1997, 15: 32–38.View ArticleGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- Sarda L, Desnuelle P: Action De La Lipase Pancreatique Sur Les Esters En Emulsion. Biochim Biophys Acta 1958, 30: 513–521.View ArticlePubMedGoogle Scholar
- Verger R: Enzyme kinetics of lipolysis. Methods Enzymol 1980, 64: 340–392.View ArticlePubMedGoogle Scholar
- 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.View ArticleGoogle Scholar
- Derewenda ZS: Structure and Function of Lipases. Adv Protein Chem 1994, 45: 1–52.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- Thuren T: A Model for the Molecular Mechanism of Interfacial Activation of Phospholipase-A2 Supporting the Substrate Theory. FEBS Lett 1988, 229: 95–99.View ArticlePubMedGoogle Scholar
- 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.PubMedGoogle Scholar
- Grochulski P, Li Y, Schrag JD, Cygler M: 2 Conformational States of Candida-Rugosa Lipase. Protein Sci 1994, 3: 82–91.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- Peters GH, Bywater RP: Computational analysis of chain flexibility and fluctuations in Rhizomucor miehei lipase. Protein Eng 1999, 12: 747–754.View ArticlePubMedGoogle Scholar
- Lee J, Suh SW, Shin S: Computational studies of essential dynamics of Pseudomonas cepacia lipase. J Biomol Struct Dyn 2000, 18: 297–309.View ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- Cygler M, Schrag JD: Structure as basis for understanding interfacial properties of lipases. Lipases, Part A 1997, 284: 3–27.View ArticleGoogle Scholar
- Noble ME, Cleasby A, Johnson LN, Egmond MR, Frenken LG: Analysis of the structure of Pseudomonas glumae lipase. Protein Eng 1994, 7: 559–562.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- Cajal Y, Svendsen A, Girona V, Patkar SA, Alsina MA: Interfacial control of lid opening in Thermomyces lanuginosa lipase. Biochemistry 2000, 39: 413–423.View ArticlePubMedGoogle Scholar
- 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.PubMedGoogle Scholar
- Peters GH, Bywater RP: Influence of a lipid interface on protein dynamics in a fungal lipase. Biophys J 2001, 81: 3052–3065.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- Zaks A, Klibanov AM: Enzyme-catalyzed processes in organic solvents. Proc Natl Acad Sci USA 1985, 82: 3192–3196.PubMed CentralView ArticlePubMedGoogle Scholar
- Colombo G, Ottolina G, Carrea G: Modelling of enzyme properties in organic solvents. Monatsh Chem 2000, 131: 527–547.View ArticleGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- Micaelo NM, Soares CM: Modeling hydration mechanisms of enzymes in nonpolar and polar organic solvents. FEBS J 2007, 274: 2424–2436.View ArticlePubMedGoogle Scholar
- Trodler P, Pleiss J: Modeling structure and flexibility of Candida antarctica lipase B in organic solvents. BMC Struct Biol 2008, 8: 9.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.PubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- Martel PJ, Baptista A, Petersen SB: Protein electrostatics. Biotechnol Annu Rev 1996, 2: 315–372.View ArticleGoogle Scholar
- Havranek JJ, Harbury PB: Tanford-Kirkwood electrostatics for protein modeling. Proc Natl Acad Sci USA 1999, 96: 11145–11150.PubMed CentralView ArticlePubMedGoogle Scholar
- Petersen MT, Martel P, Petersen EI, Drablos F, Petersen SB: Surface and electrostatics of cutinases. Methods Enzymol 1997, 284: 130–154.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- Honig B, Nicholls A: Classical electrostatics in biology and chemistry. Science 1995, 268: 1144–1149.View ArticlePubMedGoogle Scholar
- Sitkoff D, Sharp KA, Honig B: Accurate Calculation of Hydration Free-Energies Using Macroscopic Solvent Models. J Phys Chem 1994, 98: 1978–1988.View ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- Vriend G: What If – a Molecular Modeling and Drug Design Program. J Mol Graphics 1990, 8: 52.View ArticleGoogle Scholar
- Kabsch W, Sander C: Dictionary of Protein Secondary Structure – Pattern-Recognition of Hydrogen-Bonded and Geometrical Features. Biopolymers 1983, 22: 2577–2637.View ArticlePubMedGoogle Scholar
- DeLano WL: The PyMOL Molecular Graphics System. DeLano Scientific, San Carlos, CA, USA; 2002.Google Scholar
- Humphrey W, Dalke A, Schulten K: VMD: Visual molecular dynamics. J Mol Graphics 1996, 14: 33.View ArticleGoogle Scholar
- Eisenberg D, Schwarz E, Komarony M, Wall R: Amino acid scale: Normalized consensus hydrophobicity scale. J Mol Biol 1984, 179: 125–142.View ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar