PubMed Health⌕ Search

Biomedical subjects

Jürgen Pleiss

Publications and source records attributed to Jürgen Pleiss.

16 recordsLinked to original sources

DWARF--a data warehouse system for analyzing protein families.

BACKGROUND: The emerging field of integrative bioinformatics provides the tools to organize and systematically analyze vast amounts of highly diverse biological data and thus allows to gain a novel understanding of complex biological systems. The data warehouse DWARF applies integrative bioinformatics approaches to the analysis of large protein families. DESCRIPTION: The data warehouse system DWARF integrates data on sequence, structure, and functional annotation for protein fold families. The underlying relational data model consists of three major sections representing entities related to the protein (biochemical function, source organism, classification to homologous families and superfamilies), the protein sequence (position-specific annotation, mutant information), and the protein structure (secondary structure information, superimposed tertiary structure). Tools for extracting, transforming and loading data from public available resources (ExPDB, GenBank, DSSP) are provided to populate the database. The data can be accessed by an interface for searching and browsing, and by analysis tools that operate on annotation, sequence, or structure. We applied DWARF to the family of alpha/beta-hydrolases to host the Lipase Engineering database. Release 2.3 contains 6138 sequences and 167 experimentally determined protein structures, which are assigned to 37 superfamilies 103 homologous families. CONCLUSION: DWARF has been designed for constructing databases of large structurally related protein families and for evaluating their sequence-structure-function relationships by a systematic analysis of sequence, structure and functional annotation. It has been applied to predict biochemical properties from sequence, and serves as a valuable tool for protein engineering.

Amino Acid Sequence↗

Multiple molecular dynamics simulations of human p450 monooxygenase CYP2C9: the molecular basis of substrate binding and regioselectivity toward warfarin.

To examine the molecular basis of activity and regioselectivity of the clinically important human microsomal cytochrome P450 (CYP) monooxygenase 2C9 toward its substrate warfarin, 22 molecular dynamics simulations (3-5 ns each) were performed in the presence and absence of warfarin. The resulting trajectories revealed a stable protein core and mobile surface elements. This mobility leads to the formation of two surface channels in the region between F-G loop, B' helix/B-B' loop, beta(1)-sheet, and between helices F and I and the turn in the C-terminal antiparallel beta-sheet in the presence of warfarin. Besides the nonproductive state of the CYP2C9 warfarin complex captured in the crystal structure, three additional states were observed. These states differ in the shape of the substrate binding cavity and the position of the warfarin molecule relative to heme. In one of these states, the 7- and 6-positions of warfarin contact the heme with a marked geometrical preference for position 7 over position 6. This modeling result is consistent with experimentally determined regioselectivity (71 and 22% hydroxylation in positions 7 and 6, respectively). Access to the heme group is limited by the core amino acids Ala297, Leu362, Leu366, and Thr301, which therefore are expected to have a major impact on regioselectivity. In addition, modeling predicts that autoactivation of warfarin is sterically hindered. Our study demonstrates how the combination of mobile surface and rigid core leads to interesting properties: a broad substrate profile and simultaneously a high regioselectivity.

Aryl Hydrocarbon Hydroxylases↗

Factors mediating activity, selectivity, and substrate specificity for the thiamin diphosphate-dependent enzymes benzaldehyde lyase and benzoylformate decarboxylase.

Benzaldehyde lyase from Pseudomonas fluorescens and benzoylformate decarboxylase from Pseudomonas putida are homologous thiamin diphosphate-dependent enzymes that catalyze carboligase and carbolyase reactions. Both enzymes catalyze the formation of chiral 2-hydroxy ketones from aldehydes. However, the reverse reaction has only been observed with benzaldehyde lyase. Whereas benzaldehyde lyase is strictly R specific, the stereoselectivity of benzoylformate decarboxylase from P. putida is dependent on the structure and orientation of the substrate aldehydes. In this study, the binding sites of both enzymes were investigated by using molecular modelling studies to explain the experimentally observed differences in the activity, stereo- and enantioselectivity and substrate specificity of both enzymes. We designed a detailed illustration that describes the shape of the binding site of both enzymes and sufficiently explains the experimental effects observed with the wild-type enzymes and different variants. These findings demonstrate that steric reasons are predominantly responsible for the differences observed in the (R)-benzoin cleavage and in the formation of chiral 2-hydroxy ketones.

Aldehyde-Lyases↗

How to find soluble proteins: a comprehensive analysis of alpha/beta hydrolases for recombinant expression in E. coli.

BACKGROUND: In screening of libraries derived by expression cloning, expression of active proteins in E. coli can be limited by formation of inclusion bodies. In these cases it would be desirable to enrich gene libraries for coding sequences with soluble gene products in E. coli and thus to improve the efficiency of screening. Previously Wilkinson and Harrison showed that solubility can be predicted from amino acid composition (Biotechnology 1991, 9(5):443-448). We have applied this analysis to members of the alpha/beta hydrolase fold family to predict their solubility in E. coli. alpha/beta hydrolases are a highly diverse family with more than 1800 proteins which have been grouped into homologous families and superfamilies. RESULTS: The predicted solubility in E. coli depends on hydrolase size, phylogenetic origin of the host organism, the homologous family and the superfamily, to which the hydrolase belongs. In general small hydrolases are predicted to be more soluble than large hydrolases, and eukaryotic hydrolases are predicted to be less soluble in E. coli than prokaryotic ones. However, combining phylogenetic origin and size leads to more complex conclusions. Hydrolases from prokaryotic, fungal and metazoan origin are predicted to be most soluble if they are of small, medium and large size, respectively. We observed large variations of predicted solubility between hydrolases from different homologous families and from different taxa. CONCLUSION: A comprehensive analysis of all alpha/beta hydrolase sequences allows more efficient screenings for new soluble alpha/beta hydrolases by the use of libraries which contain more soluble gene products. Screening of hydrolases from families whose members are hard to express as soluble proteins in E. coli should first be done in coding sequences of organisms from phylogenetic groups with the highest average of predicted solubility for proteins of this family. The tools developed here can be used to identify attractive target genes for expression using protein sequences published in databases. This analysis also directs the design of degenerate, family-specific primers to amplify new members from homologous families or superfamilies with a high probability of soluble alpha/beta hydrolases.

Animals↗

Congenital adrenal hyperplasia: the molecular basis of 21-hydroxylase deficiency in H-2(aw18) mice.

The mouse strain H-2(aw18) shows typical characteristics of 21-hydroxylase deficiency (21-OHD). A deletion of the active Cyp21a1 gene has been postulated; however, the changes on the nucleotide level are still unknown. To investigate whether this animal model, the only one available, is suitable for studying congenital adrenal hyperplasia in man, a detailed analysis of the Cyp21 locus has been performed to ascertain the genetic cause of 21-OHD in H-2(aw18) mice. We demonstrate that 21-OHD is caused by unequal crossing over between the active Cyp21a1 gene and the pseudogene resulting in a hybrid Cyp21a1-Cyp21a2-p gene including a partial deletion of Cyp21a1. Next to several pseudogene-specific point mutations, various novel missense mutations and a nonsense mutation are present. Enzyme activity for each point mutation has been determined in vitro and the structure-function relationship has been studied by sequence conservation analysis and a three-dimensional murine 21-hydroxylase protein (Cyp21) structure model. The mutations are classified in three classes: I, no or minor decrease in enzyme activity: R238Q, P465L, R361K, A362V, P458L; II, loss of enzyme activity caused by inefficient electron flux: R346H, R400C; III, loss of activity due to deficient substrate binding: I462F, L464F. The combination of in vitro protein expression and three-dimensional structure modeling provides a valuable tool to understand the role of the different mutations and polymorphisms on the resulting enzyme activity. The underlying genetic mechanisms are also known to be responsible for 21-OHD in humans, so rodent 21-OHD turns out to be an excellent genetic model for studying the human disease.

Adrenal Hyperplasia, Congenital↗

Sequence and structure of epoxide hydrolases: a systematic analysis.

Epoxide hydrolases (EC 3.3.2.3) are ubiquitous enzymes that catalyze the hydrolysis of epoxides to the corresponding vicinal diols. More than 100 epoxide hydrolases (EH) have been identified or predicted, and 3 structures are available. Although they catalyze the same chemical reaction, sequence similarity is low. To identify conserved regions, all EHs were aligned. Phylogenetic analysis identified 12 homologous families, which were grouped into 2 major superfamilies: the microsomal EH superfamily, which includes the homologous families of Mammalian, Insect, Fungal, and Bacterial EHs, and the cytosolic EH superfamily, which includes Mammalian, Plant, and Bacterial EHs. Bacterial EHs show a high sequence diversity. Based on structure comparison of three known structures from Agrobacterium radiobacter AD1 (cytosolic EH), Aspergillus niger (microsomal EH), Mus musculus (cytosolic EH), and multisequence alignment and phylogenetic analysis of 95 EHs, the modular architecture of this enzyme family was analyzed. Although core and cap domain are highly conserved, the structural differences between the EHs are restricted to only two loops: the NC-loop connecting the core and the cap and the cap-loop, which is inserted into the cap domain. EHs were assigned to either of three clusters based on loop length. By using this classification, core and cap region of all EHs, NC-loops and cap-loops of 78% and 89% of all EHs, respectively, could be modeled. Representative models are available from the Lipase Engineering Database, http://www.led.uni-stuttgart.de.

Amino Acid Sequence↗

The database of epoxide hydrolases and haloalkane dehalogenases: one structure, many functions.

UNLABELLED: The epoxide hydrolases and haloalkane dehalogenases database (EH/HD) integrates sequence and structure of a highly diverse protein family, including mainly the Asp-hydrolases of EHs and HDs but also proteins, such as Ser-hydrolases non-heme peroxidases, prolyl iminopetidases and 2-hydroxymuconic semialdehyde hydrolases. These proteins have a highly conserved structure, but display a remarkable diversity in sequence and function. A total of 305 protein entries were assigned to 14 homologous families, forming two superfamilies. Annotated multisequence alignments and phylogenetic trees are provided for each homologous family and superfamily. Experimentally derived structures of 19 proteins are superposed and consistently annotated. Sequence and structure of all 305 proteins were systematically analysed. Thus, deeper insight is gained into the role of a highly conserved sequence motifs and structural elements. AVAILABILITY: The EH/HD database is available at http://www.led.uni-stuttgart.de

Amino Acid Motifs↗

Molecular modeling of family GH16 glycoside hydrolases: potential roles for xyloglucan transglucosylases/hydrolases in cell wall modification in the poaceae.

Family GH16 glycoside hydrolases can be assigned to five subgroups according to their substrate specificities, including xyloglucan transglucosylases/hydrolases (XTHs), (1,3)-beta-galactanases, (1,4)-beta-galactanases/kappa-carrageenases, "nonspecific" (1,3/1,3;1,4)-beta-D-glucan endohydrolases, and (1,3;1,4)-beta-D-glucan endohydrolases. A structured family GH16 glycoside hydrolase database has been constructed (http://www.ghdb.uni-stuttgart.de) and provides multiple sequence alignments with functionally annotated amino acid residues and phylogenetic trees. The database has been used for homology modeling of seven glycoside hydrolases from the GH16 family with various substrate specificities, based on structural coordinates for (1,3;1,4)-beta-D-glucan endohydrolases and a kappa-carrageenase. In combination with multiple sequence alignments, the models predict the three-dimensional (3D) dispositions of amino acid residues in the substrate-binding and catalytic sites of XTHs and (1,3/1,3;1,4)-beta-d-glucan endohydrolases; there is no structural information available in the databases for the latter group of enzymes. Models of the XTHs, compared with the recently determined structure of a Populus tremulos x tremuloides XTH, reveal similarities with the active sites of family GH11 (1,4)-beta-D-xylan endohydrolases. From a biological viewpoint, the classification, molecular modeling and a new 3D structure of the P. tremulos x tremuloides XTH establish structural and evolutionary connections between XTHs, (1,3;1,4)-beta-D-glucan endohydrolases and xylan endohydrolases. These findings raise the possibility that XTHs from higher plants could be active not only on cell wall xyloglucans, but also on (1,3;1,4)-beta-D-glucans and arabinoxylans, which are major components of walls in grasses. A role for XTHs in (1,3;1,4)-beta-D-glucan and arabinoxylan modification would be consistent with the apparent overrepresentation of XTH sequences in cereal expressed sequence tags databases.

Amino Acid Sequence↗

Potent mechanism-based inhibition of human CYP2B6 by clopidogrel and ticlopidine.

The thienopyridine derivatives ticlopidine and clopidogrel are inhibitors of ADP-induced platelet aggregation. Pharmacological activity of these prodrugs depends on cytochrome P450 (P450)-dependent oxidation to the active antithrombotic agent. In this study, we investigated the interaction potential of clopidogrel and ticlopidine by using human liver microsomes and recombinantly expressed P450 isoforms. Both clopidogrel and ticlopidine inhibited CYP2B6 with highest potency and CYP2C19 with lower potency. Clopidogrel also inhibited CYP2C9, and ticlopidine also inhibited CYP1A2, with lower potency. Inhibition of CYP2B6 was time- and concentration-dependent, and as shown by dialysis experiments, it was irreversible and dependent on NADPH, suggesting a mechanism-based mode of action. Inactivation was of nonpseudo-firstorder type with maximal rates of inactivation (K(inact)) for clopidogrel and ticlopidine in microsomes (recombinant CYP2B6) of 0.35 (1.5 min(-1)) and 0.5 min(-1) (0.8 min(-1)), respectively, and half-maximal inactivator concentrations (KI) were 0.5 microM (1.1 microM) for clopidogrel and 0.2 microM (0.8 microM) for ticlopidine. Inhibition was attenuated by the presence of alternative active site ligands but not by nucleophilic trapping agents or reactive oxygen scavengers, further supporting mechanism-based action. A chemical mechanism is discussed based on the known metabolic activation of clopidogrel and on the finding that hemoprotein integrity of recombinant CYP2B6 was not affected by irreversible inhibition. These results suggest the possibility of drug interactions between thienopyridine derivates and drug substrates of CYP2B6 and CYP2C19.

Aryl Hydrocarbon Hydroxylases↗

A molecular mechanism of enantiorecognition of tertiary alcohols by carboxylesterases.

Carboxylesterases containing the sequence motif GGGX catalyze the hydrolysis of esters of chiral tertiary alcohols, albeit with only low to moderate enantioselectivity, for three model substrates (linalyl acetate, methyl-1-pentin-1-yl acetate, 2-phenyl-3-butin-2-yl acetate). In order to understand the molecular mechanism of enantiorecognition and to improve enantioselectivity for this interesting substrate class, the interaction of both enantiomers with the substrate binding sites of acetylcholinesterases and p-nitrobenzyl esterase from Bacillus subtilis was modeled and correlated to experimental enantioselectivity. For all substrate-enzyme pairs, enantiopreference and ranking by enantioselectivity could be predicted by the model. In p-nitrobenzyl esterase, one of the key residues in determining enantioselectivity was G105: exchange of this amino acid for an alanine residue led to a sixfold increase of enantioselectivity (E = 19) towards 2-phenyl-3-butin-2-yl acetate. However, the effect of this mutation is specific: the same mutant had the opposite enantiopreference towards the substrate linalyl acetate. Thus, depending on the substrate structure, the same mutant has either increased enantioselectivity or opposite enantiopreference compared to the wild-type enzyme.

Alcohols↗

The Lipase Engineering Database: a navigation and analysis tool for protein families.

The Lipase Engineering Database (LED) (http://www.led.uni-stuttgart.de) integrates information on sequence, structure, and function of lipases, esterases, and related proteins. Sequence data on 806 protein entries are assigned to 38 homologous families, which are grouped into 16 superfamilies with no global sequence similarity between each other. For each family, multisequence alignments are provided with functionally relevant residues annotated. Pre-calculated phylogenetic trees allow navigation inside superfamilies. Experimental structures of 45 proteins are superposed and consistently annotated. The LED has been applied to systematically analyze sequence-structure-function relationships of this vast and diverse enzyme class. It is a useful tool to identify functionally relevant residues apart from the active site residues, and to design mutants with desired substrate specificity.

Amino Acid Sequence↗

Insight into the mechanism of the IMP-1 metallo-beta-lactamase by molecular dynamics simulations.

Two models, a purely nonbonded model and a cationic dummy atom approach, were examined for the modeling of the binuclear zinc-containing IMP-1 metallo-beta-lactamase in complex with a mercaptocarboxylate inhibitor. The cationic dummy atom approach had substantial advantages as it maintained the initial, experimentally determined geometry of the metal-containing active site during molecular dynamics simulations in water. The method was extended to the modeling of the free enzyme and the enzyme in complex with a cephalosporin substrate docked in an intermediate structure. For all three systems, the modeled complexes and the tetrahedral coordination of the zinc ions were stable. The average zinc-zinc distance increased by approximately 1 A in the substrate complex compared with the inhibitor complex and the free enzyme in which a hydroxide ion acts as a bridging ligand. Thus, the zinc ions are predicted to undergo a back and forth movement upon the cycle of hydrolysis. In contrast to previous assumptions, no interaction of the Asn167 side chain with the bound cephalosporin substrate was observed. Our observations are in agreement with quantum-mechanical calculations and experimental data and indicate that the cationic dummy atom approach is useful to model zinc-containing metallo-beta-lactamases as free proteins, in complex with inhibitors and in complex with substrates.

Binding Sites↗

Regioselectivity of CYP2B6: homology modeling, molecular dynamics simulation, docking.

Human cytochrome P450 (CYP) 2B6 activates the anticancer prodrug cyclophosphamide (CPA) by 4-hydroxylation. In contrast, the same enzyme catalyzes N-deethylation of a structural isomer, the prodrug ifosfamide (IFA), thus causing severe adverse drug effects. To model the molecular interactions leading to a switch in regioselectivity, the structure of CYP2B6 was modeled based on the structure of rabbit CYP2C5. We modeled the missing 22-residue loop in CYP2C5 between helices F and G (the F-G loop), which is not resolved in the X-ray structure, by molecular dynamics (MD) simulations using a simulated annealing protocol. The modeled conformation of the loop was validated by unconstrained MD simulations of the complete enzymes (CYP2C5 and CYP2B6) in water for 70 and 120 ps, respectively. The simulations were stable and led to a backbone r.m.s. deviation of 1.7 A between the two CYPs. The shape of the substrate binding site of CYP2B6 was further analyzed. It consists of three well-defined hydrophobic binding pockets adjacent to the catalytic heme. Size, shape and hydrophobicity of these pockets were compared to the shapes of the two structurally isomeric substrates. In their preferred orientation in the binding site, both substrates fill all three binding pockets without repulsive interactions. The distance to the heme iron is short enough for 4-hydroxylation and N-deethylation to occur for CPA and IFA, respectively. However, if the substrates are docked in the non-preferred orientation (such that 4-hydroxylation and N-deethylation would occur for IFA and CPA, respectively), one pocket is left empty, and clashes were observed between the substrates and the enzyme.

Amino Acid Sequence↗

Mapping of a hapten-binding site: molecular modeling and site-directed mutagenesis study of an anti-atrazine antibody.

A three-dimensional model of the variable domain of the atrazine-specific Fab fragment K411B was constructed by molecular modeling using known structures of highly homologous immunoglobulins as templates. Molecular dynamic simulations and cross-reactivity data were used to predict residues responsible for the binding of the hapten 4-chloro-6-(isopropylamino)-1,3,5-triazine-2-(6-aminohexanecarboxylic acid) (iPr/Cl/C6) instead of atrazine. Specific binding pockets could be defined for the chlorine, the isopropylamino group and the C6-spacer of the hapten. The influence of various amino acids on hapten binding was investigated by site-directed mutagenesis, and the effect of these mutations was analyzed by capture ELISA using the hapten iPr/Cl/C6 and 4-amino-6-chloro-1,3,5-triazine-2-(6-aminohexanecarboxylic acid) (H/Cl/C6). GlyH100a seems to be important in determining the conformation of the heavy-chain complementarity determining region H3; replacing it with any other residue prevented the binding of the hapten. Altering residues responsible for the binding of the chlorine atom (TrpH33, GluH50 and TyrL96) decreased the affinity significantly. Hapten-spacer recognition can be attributed to the interaction with PheL32; replacing PheL32 by leucine reduced the affinity towards iPr/Cl/C6. A triple mutant Fab fragment (GlnL89Glu, ValH37Ile and GluL3Val) showed an affinity 5-fold greater towards iPr/Cl/C6 compared to the wild-type K411B, as a result of better recognition of the isopropylamino group of iPr/Cl/C6.

Amino Acid Sequence↗

Blocking the tunnel: engineering of Candida rugosa lipase mutants with short chain length specificity.

The molecular basis of chain length specificity of Candida rugosa lipase 1 was investigated by molecular modeling and site-directed mutagenesis. The synthetic lip1 gene and the lipase mutants were expressed in Pichia pastoris and assayed for their chain length specificity in single substrate assays using triglycerides as well as in a competitive substrate assay using a randomized oil. Mutation of amino acids at different locations inside the tunnel (P246F, L413F, L410W, L410F/S300E, L410F/S365L) resulted in mutants with a different chain length specificity. Mutants P246F and L413F have a strong preference for short chain lengths whereas substrates longer than C10 are hardly hydrolyzed. Increasing the bulkiness of the amino acid at position 410 led to mutants that show a strong discrimination of chain lengths longer than C14. The results obtained can be explained by a simple mechanical model: the activity for a fatty acid sharply decreases as it becomes long enough to reach the mutated site. In contrast, a mutation at the entrance of the tunnel (L304F) has a strong impact on C4 and C6 substrates. This mutant is nevertheless capable of hydrolyzing chain lengths longer than C8.

Binding, Competitive↗