Towards novel biocatalysts via protein design: the case of lipases

Abstract: During the past 3 years, the tertiary structures of several lipases have been solved by X-ray analysis. The structures revealed unique features such as hydrophobic 'patches' on the surface, presumably involved in lipid supersubstrate binding, and a lid structure which covers the active site in the absence of substrate. Only very recently the first X-ray structure of a bacterial lipase has been solved, and further structural features different from lipases of eukaryotic origin became apparent. Many lipase genes have been cloned and sequenced recently, and expression systems for the preparation of recombinant enzymes in good yields are available. As an example, the lipase from Rhizopus oryzae has been successfully expressed by us in Escherichia coli , and the resulting inclusion bodies were renatured in high yields. Consequently, the mechanism of action of lipases is now being studied via site-directed mutagenesis, and the rational design of lipases for the selective transformation of substrates is presently addressed in several laboratories.     

X-ray structure determination and comparison of two crystal forms of a variant (Asn115Arg) of the alkaline protease from Bacillus alcalophilus refined at 1.85 A resolution.

The X-ray structure determination, refinement and comparison of two crystal forms of a variant (Asn115Arg) of the alkaline protease from Bacillus alcalophilus is described. Under identical conditions crystals were obtained in the orthorhombic space group P2(1)2(1)2(1) (form I) and the rhombohedral space group R32 (form II). For both space groups the structures of the protease were solved by molecular replacement and refined at 1.85 A resolution. The final R-factors are 17.9% and 17.1% for form I and form II, respectively. The root-mean-square deviation between the two forms is 0.48 A and 0.86 A for main-chain and side-chain atoms, respectively. Due to differences in crystal lattice contacts and packing, the structures of the two crystal forms differ in intermolecular interaction affecting the local conformation of three flexible polypeptide sequences (Ser50-Glu55, Ser99-Gly102, Gly258-Ser259) at the surface of the protein. While the two overall structures are very similar, the differences are significantly larger than the errors inherent in the structure determination. As expected, the differences in the temperature factors in form I and II are correlated with the solvent accessibility of the corresponding amino acid residues. In form II, two symmetry-related substrate binding sites face each other, forming a tight intermolecular interaction. Some residues contributing to this intermolecular interaction are also found to be involved in the formation of the complex between subtilisin Carlsberg and the proteinaceous inhibitor eglin C. This demonstrates that the two symmetry-related molecules interact with each other at the same molecular surface area that is used for binding of substrates and inhibitors.     

An Irreversible and Kinetically Controlled Process: Thermal Induced Denaturation of <Emphasis Type="SmallCaps">L</Emphasis>-2-Hydroxyisocaproate Dehydrogenase from <Emphasis Type="Italic">Lactobacillus confusus</Emphasis>

The thermal denaturation of Lactobacillus confusus l-2-Hydroxyisocaproate Dehydrogenase (l-HicDH) has been studied by Differential Scanning Calorimetry (DSC). The stability of this enzyme has been investigated at different pH conditions. The results of this study indicate that the thermal denaturation of this enzyme is irreversible and the T _m is dependent on the scan-rate, which suggests that the denaturation process of l-HicDH is kinetically determined. The heat capacity function of l-HicDH shows a single peak with the T _m values between 52.14°C and 55.89°C at pH 7.0 at different scan rates. These results indicate that the whole l-HicDH could unfold as a single cooperative unit, and intersubunit interactions of this homotetrameric enzyme must play a significant role in the stabilization of the whole enzyme. The rate constant of the unfolding is analyzed as a first order kinetic constant with the Arrhenius equation, and the activation energy has been calculated. The variation of the activation energy values obtained with different methods does not support the validity of the one-step irreversible model. The denaturation pathway was described by a three-state model, N → U → F, in which the dissociation of the tetramer takes place as an irreversible step before the irreversible unfolding of the monomers. The calorimetric enthalpy associated with the irreversible dissociation and the calorimetric enthalpy associated with the unfolding of the monomer were obtained from the best fitting procedure. Thermal unfolding of l-HicDH was also studied using Circular Dichroism (CD) spectroscopy. Both methods yielded comparable values.     

GTP plus water mimic ATP in the active site of protein kinase CK2

The structures of the catalytic subunit of protein kinase CK2 from Zea mays complexed with Mg2+ and with analogs of ATP or GTP were determined to 2.2 A resolution. Unlike most other protein kinases, CK2 from various sources shows 'dual-cosubstrate specificity', that is, the ability to efficiently use either ATP or GTP as a cosubstrate. The structures of these complexes demonstrate that water molecules are critical to switch the active site of CK2 from an ATP- to a GTP-compatible state. An understanding of the structural basis of dual-cosubstrate specificity may help in the design of drugs that target CK2 or other kinases with this property.     

Efficient methods for filtering and ranking fragments for the prediction of structurally variable regions in proteins

The prediction of protein 3D structures close to insertions and deletions or, more generally, loop prediction, is still one of the major challenges in homology modeling projects. In this article, we developed ranking criteria and selection filters to improve knowledge-based loop predictions. These criteria were developed and optimized for a test data set containing 678 insertions and deletions. The examples are, in principle, predictable from the used loop database with an RMSD < 1 A and represent realistic modeling situations. Four noncorrelated criteria for the selection of fragments are evaluated. A fast prefilter compares the distance between the anchor groups in the template protein with the stems of the fragments. The RMSD of the anchor groups is used for fitting and ranking of the selected loop candidates. After fitting, repulsive close contacts of loop candidates with the template protein are used for filtering, and fragments with backbone torsion angles, which are unfavorable according to a knowledge-based potential, are eliminated. By the combined application of these filter criteria to the test set, it was possible to increase the percentage of predictions with a global RMSD < 1 A to over 50% among the first five ranks, with average global RMSD values for the first rank candidate that are between 1.3 and 2.2 A for different loop lengths. Compared to other examples described in the literature, our large numbers of test cases are not self-predictions, where loops are placed in a protein after a peptide loop has been cut out, but are attempts to predict structural changes that occur in evolution when a protein is affected by insertions and deletions.     

The Helmholtz Network for Bioinformatics: an integrative web portal for bioinformatics resources

SUMMARY: The Helmholtz Network for Bioinformatics (HNB) is a joint venture of eleven German bioinformatics research groups that offers convenient access to numerous bioinformatics resources through a single web portal. The 'Guided Solution Finder' which is available through the HNB portal helps users to locate the appropriate resources to answer their queries by employing a detailed, tree-like questionnaire. Furthermore, automated complex tool cascades ('tasks'), involving resources located on different servers, have been implemented, allowing users to perform comprehensive data analyses without the requirement of further manual intervention for data transfer and re-formatting. Currently, automated cascades for the analysis of regulatory DNA segments as well as for the prediction of protein functional properties are provided. AVAILABILITY: The HNB portal is available at http://www.hnbioinfo.de     

Expression, purification, and aggregation studies of His-tagged thermoalkalophilic lipase from Bacillus thermocatenulatus

The His-tagged lipase BTL2 from Bacillus thermocatenulatus was expressed in Escherichia coli and purified to homogeneity by a simple, one-step purification protocol using immobilized metal affinity chromatography. The success of protein separation and purification was pH-dependent and increased with decreasing pH. The purified BTL2 lipase showed a strong tendency to aggregate upon concentration, which prevented a reproducible crystallization. Aggregation studies using dynamic light-scattering (DLS) analysis were performed to improve the purification and concentration of BTL2 lipase. Different chemical classes of additives were tested to manipulate the aggregation behaviour of BTL2 lipase with the aim of obtaining a monodisperse sample to use for crystallization. For the process of concentration of BTL2 lipase in monomeric form, the alcohol 2-propanol and the ionic detergent dodecyl dimethylamine-N-oxide (LDAO) were found to be necessary. For the concentrated lipase, the availability of 5% 2-propanol was sufficient to hold the lipase in monomeric form and no additional detergent was needed.     

The crystal structure of R-specific alcohol dehydrogenase from Lactobacillus brevis suggests the structural basis of its metal dependency

The crystal structure of the apo-form of an R-specific alcohol dehydrogenase from Lactobacillus brevis (LB-RADH) was solved and refined to 1.8A resolution. LB-RADH is a member of the short-chain dehydrogenase/reductase (SDR) enyzme superfamily. It is a homotetramer with 251 amino acid residues per subunit and uses NADP(H) as co-enzyme. NADPH and the substrate acetophenone were modelled into the active site. The enantiospecificity of the enzyme can be explained on the basis of the resulting hypothetical ternary complex. In contrast to most other SDR enzymes, the catalytic activity of LB-RADH depends strongly on the binding of Mg(2+). Mg(2+) removal by EDTA inactivates the enzyme completely. In the crystal structure, the Mg(2+)-binding site is well defined. The ion has a perfect octahedral coordination sphere and occupies a special position concerning crystallographic and molecular point symmetry, meaning that each RADH tetramer contains two magnesium ions. The magnesium ion is no direct catalytic cofactor. However, it is structurally coupled to the putative C-terminal hinge of the substrate-binding loop and, via an extended hydrogen bonding network, to some side-chains forming the substrate binding region. Therefore, the presented structure of apo-RADH provides plausible explanations for the metal dependence of the enzyme.     

Positioning of anchor groups in protein loop prediction: the importance of solvent accessibility and secondary structure elements

The prediction of loop regions in the process of protein structure prediction by homology is still an unsolved problem. In an earlier publication, we could show that the correct placement of the amino acids serving as an anchor group to be connected by a loop fragment with a predicted geometry is a highly important step and an essential requirement within the process (Lessel and Schomburg, Proteins 1999; 37:56-64). In this article, we present an analysis of the quality of possible loop predictions with respect to gap length, fragment length, amino acid type, secondary structure, and solvent accessibility. For 550 insertions and 544 deletions, we test all possible positions for anchor groups with an inserted loop of a length between 3 and 12 amino acids. We could show that approximately 80% of the indel regions could be predicted within 1.5 A RMSD from a knowledge-based loop data base if criteria for the correct localization of anchor groups could be found and the loops can be sorted correctly. From our analysis, several conclusions regarding the optimal placement of anchor groups become obvious: (1) The correct placement of anchor groups is even more important for longer gap lengths, (2) medium length fragments (length 5-8) perform better than short or long ones, (3) the placement of anchor groups at hydrophobic amino acids gives a higher chance to include the best possible loop, (4) anchor groups within secondary structure elements, in particular beta-sheets are suitable, (5) amino acids with lower solvent accessibility are better anchor group. A preliminary test using a combination of the anchor group positioning criteria deduced from our analysis shows very promising results.