PubMed Health⌕ Search

Biomedical subjects

R M Knegtel

Publications and source records attributed to R M Knegtel.

At least 19 recordsLinked to original sources
PubMed2000

Similarity-driven flexible ligand docking.

A similarity-driven approach to flexible ligand docking is presented. Given a reference ligand or a pharmacophore positioned in the protein active site, the method allows inclusion of a similarity term during docking. Two different algorithms have been implemented, namely, a similarity-penalized docking (SP-DOCK) and a similarity-guided docking (SG-DOCK). The basic idea is to maximally exploit the structural information about the ligand binding mode present in cases where ligand-bound protein structures are available, information that is usually ignored in standard docking procedures. SP-DOCK and SG-DOCK have been derived as modified versions of the program DOCK 4.0, where the similarity program MIMIC acts as a module for the calculation of similarity indices that correct docking energy scores at certain steps of the calculation. SP-DOCK applies similarity corrections to the set of ligand orientations at the end of the ligand incremental construction process, penalizing the docking energy and, thus, having only an effect on the relative ordering of the final solutions. SG-DOCK applies similarity corrections throughout the entire ligand incremental construction process, thus affecting not only the relative ordering of solutions but also actively guiding the ligand docking. The performance of SP-DOCK and SG-DOCK for binding mode assessment and molecular database screening is discussed. When applied to a set of 32 thrombin ligands for which crystal structures are available, SG-DOCK improves the average RMSD by ca. 1 A when compared with DOCK. When those 32 thrombin ligands are included into a set of 1,000 diverse molecules from the ACD, DIV, and WDI databases, SP-DOCK significantly improves the retrieval of thrombin ligands within the first 10% of each of the three databases with respect to DOCK, with minimal additional computational cost. In all cases, comparison of SP-DOCK and SG-DOCK results with those obtained by DOCK and MIMIC is performed.

Crystallography, X-Ray↗
PubMed2000

Sequence annotation of nuclear receptor ligand-binding domains by automated homology modeling.

The quality of three-dimensional homology models derived from protein sequences provides an independent measure of the suitability of a protein sequence for a certain fold. We have used automated homology modeling and model assessment tools to identify putative nuclear hormone receptor ligand-binding domains in the genome of Caenorhabditis elegans. Our results indicate that the availability of multiple crystal structures is crucial to obtaining useful models in this receptor family. The majority of annotated mammalian nuclear hormone receptors could be assigned to a ligand-binding domain fold by using the best model derived from any of four template structures. This strategy also assigned the ligand-binding domain fold to a number of C.elegans. sequences without prior annotation. Interestingly, the retinoic acid receptor crystal structure contributed most to the number of sequences that could be assigned to a ligand-binding domain fold. Several causes for this can be suggested, including the high quality of this protein structure in terms of our assessment tools, similarity between the biological function or ligand of this receptor and the modeled genes and gene duplication in C.elegans.

Animals↗
PubMed1999

Efficacy and selectivity in flexible database docking.

Flexible database docking with DOCK 4.0 has been evaluated for its ability to retrieve biologically active molecules from a database of approximately 1,000 compounds with known activities against thrombin and the progesterone receptor. The retrieval of known actives and chemically similar but inactive molecules was monitored as a function of conformational and orientational sampling. The largest enrichment of actives among the 10% highest ranking molecules is obtained when only five conformations are used to seed the next round of ligand reconstruction and limited sampling is applied to place the base fragment in the binding site. The performance of energy and chemical scoring, as implemented in DOCK 4.0, was found to depend on the protein used for docking. For the progesterone receptor, energy scoring yields the largest enrichments (64%) in terms of actives retrieved among the 10% top scoring molecules, while chemical scoring performs best for thrombin (94%). With the exception of the application of energy scoring to the progesterone receptor, both energy-based scoring schemes applied in this study do not discriminate well between true actives and chemically similar but inactive compounds. In conclusion, flexible docking is able to effectively prioritize high-throughput screening databases, using less conformational sampling than normally required for appropriate reconstruction of protein-ligand complexes. The more subtle discrimination between chemically similar classes of active and inactive compounds remains, however, problematic.

Algorithms↗
PubMed1999

Comparison of two implementations of the incremental construction algorithm in flexible docking of thrombin inhibitors.

A set of 32 known thrombin inhibitors representing different chemical classes has been used to evaluate the performance of two implementations of incremental construction algorithms for flexible molecular docking: DOCK 4.0 and FlexX 1.5. Both docking tools are able to dock 10-35% of our test set within 2 A of their known, bound conformations using default sampling and scoring parameters. Although flexible docking with DOCK or FlexX is not able to reconstruct all native complexes, it does offer a significant improvement over rigid body docking of single, rule-based conformations, which is still often used for docking of large databases. Docking of sets of multiple conformers of each inhibitor, obtained with a novel protocol for diverse conformer generation and selection, yielded results comparable to those obtained by flexible docking. Chemical scoring, which is an empirically modified force field scoring method implemented in DOCK 4.0, outperforms both interaction energy scoring by DOCK and the Böhm scoring function used by FlexX in rigid and flexible docking of thrombin inhibitors. Our results indicate that for reliable docking of flexible ligands the selection of anchor fragments, conformational sampling and currently available scoring methods still require improvement.

Algorithms↗
PubMed1998

Practical evaluation of comparative modelling and threading methods.

Six protein pairs, all with known 3D-structures, were used to evaluate different protein structure prediction tools. Firstly, alignments between a target sequence and a template sequence or structure were obtained by sequence alignment with QUANTA or by threading with THREADER, 123D and PHD Topits. Secondly, protein structure models were generated using MODELLER. The two protein structure assessment tools used were the root mean square deviation (RMSD) compared with the experimental target structure and the total 3D profile score. Also the accuracy of the active sites of models built in the absence and presence of ligands was investigated. Our study confirms that threading methods are able to yield more accurate models than comparative modelling in cases of low sequence identity (< 30%). However, a gap of 2 A (RMSD) exists between the theoretically best model and the models obtained by threading methods. For high sequence identities (> 30%) comparative modelling using MODELLER resulted in accurate models. Furthermore, the total 3D profile score was not always able to distinguish correct from incorrect folds when different alignment methods were used. Finally, we found it to be important to include possible ligands in the model-building process in order to prevent unrealistic filling of active site areas.

Algorithms↗
PubMed1998

Rational design of novel antimicrobials: blocking purine salvage in a parasitic protozoan.

All parasitic protozoa obtain purine nucleotides solely by salvaging purine bases and/or nucleosides from their host. This observation suggests that inhibiting purine salvage may be a good way of killing these organisms. To explore this idea, we attempted to block the purine salvage pathway of the parasitic protozoan Tritrichomonas foetus. T. foetus is a good organism to study because its purine salvage depends primarily on a single enzyme, hypoxanthine-guanine-xanthine phosphoribosyltransferase (HGXPRTase), and could provide a good model for rational drug design through specific enzyme inhibition. Guided by the crystal structure of T. foetus HGXPRTase, we used structure-based drug design to identify several non-purine compounds that inhibited this enzyme without any detectable effect on human HGPRTase. One of these compounds, 4-[N-(3, 4-dichlorophenyl)carbamoyl]phthalic anhydride (referred to as TF1), was selected for further characterization. TF1 was shown to be a competitive inhibitor of T. foetus HGXPRTase with respect to both guanine (in the forward reaction; Ki = 13 microM) and GMP (in the reverse reaction; Ki = 10 microM), but showed no effect on the homologous human enzyme at concentrations of up to 1 mM. TF1 inhibited the in vitro growth of T. foetus with an EC50 of approximately 40 microM. This inhibitory effect was associated with a decrease in the incorporation of exogenous guanine into nucleic acids, and could be reversed by supplementing the growth medium with excess exogenous hypoxanthine or guanine. Thus, rationally targeting an essential enzyme in a parasitic organism has yielded specific enzyme inhibitors capable of suppressing that parasite's growth.

Animals↗
PubMed1997

Molecular docking to ensembles of protein structures.

Until recently, applications of molecular docking assumed that the macromolecular receptor exists in a single, rigid conformation. However, structural studies involving different ligands bound to the same target biomolecule frequently reveal modest but significant conformational changes in the target. In this paper, two related methods for molecular docking are described that utilize information on conformational variability from ensembles of experimental receptor structures. One method combines the information into an "energy-weighted average" of the interaction energy between a ligand and each receptor structure. The other method performs the averaging on a structural level, producing a "geometry-weighted average" of the inter-molecular force field score used in DOCK 3.5. Both methods have been applied in docking small molecules to ensembles of crystal and solution structures, and we show that experimentally determined binding orientations and computed energies of known ligands can be reproduced accurately. The use of composite grids, when conformationally different protein structures are available, yields an improvement in computational speed for database searches in proportion to the number of structures.

Computer Simulation↗
PubMed1997

Alkylation of a catalytic aspartate group of the SIV protease by an epoxide inhibitor.

Specific irreversible inhibition of the SIV protease by FMOC-protected piperidine epoxide 1 involves alkylation of the protein. Tryptic digestion of the alkylated protein and mass spectrometric analysis of the peptides identify an active site aspartic acid (Asp-25) as the single residue that is alkylated. Computer modeling of 1 bound in the crystal structure of the SIV protease using DOCK 3.5 indicates that 1 has appropriate access to the active site. It is able to align in an orientation that allows a proton to be transferred to the epoxide from one of the catalytic aspartic acid groups in conjunction with nucleophilic attack on the epoxide of the carboxylate moiety of the second catalytic aspartic acid residue. Hydrophobic interactions are not optimal for this process due, in part, to the rigidity of the inhibitor ring system and the planar conformation of the amide. The combination of modeling with protein alkylation can provide insights into structural modifications of the inhibitor that may lead to improved inhibitory activity.

Alkylation↗
PubMed1996

The raw starch binding domain of cyclodextrin glycosyltransferase from Bacillus circulans strain 251.

The E-domain of cyclodextrin glycosyltransferase (CGTase) (EC 2.4.1.19) from Bacillus circulans strain 251 is a putative raw starch binding domain. Analysis of the maltose-dependent CGTase crystal structure revealed that each enzyme molecule contained three maltose molecules, situated at contact points between protein molecules. Two of these maltoses were bound to specific sites in the E-domain, the third maltose was bound at the C-domain. To delineate the roles in raw starch binding and cyclization reaction kinetics of the two maltose binding sites in the E-domain, we replaced Trp-616 and Trp-662 of maltose binding site 1 and Tyr-633 of maltose binding site 2 by alanines using site-directed mutagenesis. Purified mutant CGTases were characterized with respect to raw starch binding and cyclization reaction kinetics on both soluble and raw starch. The results show that maltose binding site 1 is most important for raw starch binding, whereas maltose binding site 2 is involved in guiding linear starch chains into the active site. beta-Cyclodextrin causes product inhibition by interfering with catalysis in the active site and the function of maltose binding site 2 in the E-domain. CGTase mutants in the E-domain maltose binding site 1 could no longer be crystallized as maltose-dependent monomers. Instead, the W616A mutant CGTase protein was successfully crystallized as a carbohydrate-independent dimer; its structure has been refined to 2.2 A resolution. The three-dimensional structure shows that, within the error limits, neither the absence of carbohydrates nor the W616A mutation caused significant further conformational changes. The modified starch binding and cyclization kinetic properties observed with the mutant CGTase proteins thus can be directly related to the amino acid replacements.

Amino Acid Sequence↗
PubMed1996

Structure of cyclodextrin glycosyltransferase complexed with a maltononaose inhibitor at 2.6 angstrom resolution. Implications for product specificity.

Crystals of the Y195F mutant of cyclodextrin glycosyltransferase from Bacillus circulans strain 251 were subjected to a double soaking procedure, in which they were first soaked in a solution containing the inhibitor acarbose and subsequently in a solution containing maltohexaose. The refined structure of the resulting protein-carbohydrate complex has final crystallographic and free R-factors for data in the 8-2.6 angstrom resolution range of 15.0% and 21.5%, respectively, and reveals that a new inhibitor, composed of nine saccharide residues, is bound in the active site. The first four residues correspond to acarbose and occupy the same subsites near the catalytic residues as observed in the previously reported acarbose-enzyme complex [Strokopytov et al. (1995) Biochemistry 34, 2234-2240]. An oliogosaccharide consisting of five glucose residues has been coupled to the nonreducing end of acarbose. At the fifth residue the polysaccharide chain makes a sharp turn, allowing it to interact with residues Tyr89, Phe195, and Asn193 and a flexible loop formed by residues 145-148. On the basis of the refined model of the complex an explanation is given for the product specificity of CGTases.

Amino Acid Sequence↗
PubMed1996

Crystal structure at 2.3 A resolution and revised nucleotide sequence of the thermostable cyclodextrin glycosyltransferase from Thermonanaerobacterium thermosulfurigenes EM1.

The crystal structure of the cyclodextrin glycosyltransferase (CGTase) from the thermophilic microorganism Thermoanaerobacterium thermosulfurigenes EM1 has been elucidated at 2.3 A resolution. The final model consists of all 683 amino acid residues, two calcium ions and 343 water molecules, and has a crystallographic R-factor of 17.9% (Rfree 24.9%) with excellent stereochemistry. The overall fold of the enzyme is highly similar to that reported for mesophilic CGTases and differences are observed only at surface loop regions. Closer inspection of these loop regions and comparison with other CGTase structures reveals that especially loops 88-95, 335-339 and 534-539 possibly contribute with novel hydrogen bonds and apolar contacts to the stabilization of the enzyme. Other structural features that might confer thermostability to the T. thermosulfurigenes EM1 CGTase are the introduction of five new salt-bridges and three Gly to Ala/Pro substitutions. The abundance of Ser, Thr and Tyr residues near the active site and oligosaccharide binding sites might explain the increased thermostability of CGTase in the presence of starch, by allowing amylose chains to bind non-specifically to the protein. Additional stabilization of the A/E domain interface through apolar contacts involves residues Phe273 and Tyr187. No additional or improved calcium binding is observed in the structure, suggesting that the observed stabilization in the presence of calcium ions is caused by the reduced exchange of calcium from the protein to the solvent, rendering it less susceptible to unfolding. The 50% decrease in cyclization activity of the T. thermosulfurigenes EM1 CGTase compared with that of B. circulans strain 251 appears to be caused by the changes in the conformation and amino acid composition of the 88-95 loop. In the T. thermosulfurigenes EM1 CGTase there is no residue homologous to Tyr89, which was observed to take part in stacking interactions with bound substrate in the case of the B. circulans strain 251 CGTase. The lack of this interaction in the enzyme-substrate complex is expected to destabilize bound substrates prior to cyclization. Apparently, some catalytic functionality of CGTase has been sacrificed for the sake of structural stability by modifying loop regions near the active site.

Amino Acid Sequence↗
PubMed1995

Crystallographic studies of the interaction of cyclodextrin glycosyltransferase from Bacillus circulans strain 251 with natural substrates and products.

Asp-229, Glu-257, and Asp-328 constitute the catalytic residues in cyclodextrin glycosyl transferase from Bacillus circulans strain 251. Via site-directed mutagenesis constructed D229N, E257Q, and D328N mutant proteins showed a 4,000-60,000-fold reduction of cyclization activity. A D229N/E257Q double mutant showed a 700,000-fold reduction and was crystallized for use in soaking experiments with alpha-cyclodextrin. Crystal structures were determined of wild type CGTase soaked at elevated pH with alpha-cyclodextrin (resolution, 2.1 A) and maltoheptaose (2.4 A). In addition, structures at cryogenic temperature were solved of the unliganded enzyme (2.2 A) and of the D229N/E257Q mutant after soaking with alpha-cyclodextrin (2.6 A). In the crystals soaked in alpha-cyclodextrin and maltoheptaose, a maltotetraose molecule is observed to bind in the active site. Residue 229 is at hydrogen bonding distance from the C-6 hydroxyl group of the sugar, which after cleavage will contain the new reducing end. In the D229N/E257Q double mutant structure, two alpha-cyclodextrins are observed to replace two maltoses at the E-domain, thus providing structural information on product inhibition via binding to the enzyme's raw starch binding domain.

Bacillus↗
PubMed1995

A model for the LexA repressor DNA complex.

A structural model for the interaction of the LexA repressor DNA binding domain (DBD) with operator DNA is derived by means of Monte Carlo docking. Protein-DNA complexes were generated by docking the LexA repressor DBD NMR solution structure onto both rigid and bent B-DNA structures while giving energy bonuses for contacts in agreement with experimental data. In the resulting complexes, helix III of the LexA repressor DBD is located in the major groove of the DNA and residues Asn-41, Glu-44, and Glu-45 form specific hydrogen bonds with bases of the CTGT DNA sequence. Ser-39, Ala-42, and Asn-41 are involved in a hydrophobic interaction with the methyl group of the first thymine base. Residues in the loop region connecting the two beta-sheet strands are involved in nonspecific contacts near the dyad axis of the operator. The contacts observed in the docked complexes cover the entire consensus CTGT half-site DNA operator, thus explaining the specificity of the LexA repressor for such sequences. In addition, a large number of nonspecific interactions between protein and DNA is observed. The agreement between the derived model for the LexA repressor DBD/DNA complex and experimental biochemical results is discussed.

Amino Acid Sequence↗
PubMed1995

NMR structural studies on the zinc finger domains of nuclear hormone receptors.

This chapter presents an overview of the application of modern NMR methods in structural studies of the DNA binding domains (DBDs) of nuclear hormone receptors. The DBDs studied so far comprise those of the glucocorticoid, estrogen, retinoic acid and retinoid X receptors. NMR spectroscopy has allowed the elucidation of the first structures of this family of C4-type zinc fingers, which led to a better understanding of their role in gene regulation. Crystallographic studies provided insight in protein-protein and protein-DNA interactions. Subsequent studies, applying NMR, have provided deeper insight into a diversity of issues concerning these proteins, ranging from backbone dynamics and metal coordination to the interaction of these domains with their DNA target sites. From this work a picture emerges of a class of closely related zinc-binding proteins which, despite their strong sequence homology, exhibit interesting structural and functional differences between members of different subfamilies.

Amino Acid Sequence↗
PubMed1994

MONTY: a Monte Carlo approach to protein-DNA recognition.

A Monte Carlo method is described for automated docking of proteins on DNA. The simulation program MONTY keeps the entire DNA and the protein backbone and core fixed while protein surface side-chains are allowed to rotate freely. The entire protein is rotated and translated by small random steps in order to find the best fit with the DNA. New configurations are accepted on basis of their Boltzmann probability. Protein-DNA interaction is represented by square well potentials for hydrogen bond and van der Waals interactions. The structure with the largest interaction energy encountered during the simulation is saved. The method is tested on complexes of the 434 Cro protein and its operator DNA where the protein is shifted up or down one or two base-pairs and is subsequently allowed to find back its native binding site. This protocol is performed for shifted complexes derived from the crystal structure, shifted complexes where the crystal structure DNA is replaced by standard B-DNA and shifted complexes where in addition the protein is replaced by protein from the uncomplexed crystal structure. In all three cases the six lowest energy structures correspond to complexes close to the native complex. The quality of sequence specific recognition diminishes, however, when the molecular surface complementarity between protein and DNA decreases.

Base Composition↗
PubMed1994

Monte Carlo docking of protein-DNA complexes: incorporation of DNA flexibility and experimental data.

A Monte Carlo simulation program (MONTY) has been developed to dock proteins onto DNA. Protein and DNA interact via square-well potentials for hydrogen bond and van der Waals interactions. The effect of the inclusion of DNA flexibility and experimentally derived restraints has been tested on members of the helix-turn-helix family of DNA binding proteins. Unwinding and bending the DNA double helix improves the number of correctly retrieved hydrogen bonds in simulations starting from the 434 cro protein monomer complexed with a standard B-DNA OR1 half-site. Agreement with phosphate ethylation interference and mutagenesis data is rewarded with energy bonuses. This protocol was tested on protein-DNA complexes of 434 cro, lac headpiece and a mutant lac headpiece resembling the gal repressor headpiece with the recognition helices in correct and reversed orientations in the DNA major groove. The inclusion of experimental data gives an improved convergence of the correctly oriented structures and allows for an easier discrimination between correctly and incorrectly docked complexes.

Algorithms↗
PubMed1993

Heteronuclear 113Cd-1H NMR study of metal coordination in the human retinoic acid receptor-beta DNA binding domain.

The two zinc fingers of the DNA binding domain of the human retinoic acid receptor-beta were labelled with 113Cd. Two- and three-dimensional heteronuclear nuclear magnetic resonance (NMR) experiments show that the first eight conserved cysteine residues coordinate the two zinc ions tetrahedrally. The ninth conserved cysteine is not involved in metal coordination. In each finger one cysteine exhibits a heteronuclear 113Cd-1H coupling constant substantially smaller than those of the other metal binding cysteines.

Amino Acid Sequence↗