Search PubMed⌕ Search

Biomedical subjects

Gert Vriend

Publications and source records attributed to Gert Vriend.

At least 19 recordsLinked to original sources

Correcting ligands, metabolites, and pathways.

BACKGROUND: A wide range of research areas in bioinformatics, molecular biology and medicinal chemistry require precise chemical structure information about molecules and reactions, e.g. drug design, ligand docking, metabolic network reconstruction, and systems biology. Most available databases, however, treat chemical structures more as illustrations than as a datafield in its own right. Lack of chemical accuracy impedes progress in the areas mentioned above. We present a database of metabolites called BioMeta that augments the existing pathway databases by explicitly assessing the validity, correctness, and completeness of chemical structure and reaction information. DESCRIPTION: The main bulk of the data in BioMeta were obtained from the KEGG Ligand database. We developed a tool for chemical structure validation which assesses the chemical validity and stereochemical completeness of a molecule description. The validation tool was used to examine the compounds in BioMeta, showing that a relatively small number of compounds had an incorrect constitution (connectivity only, not considering stereochemistry) and that a considerable number (about one third) had incomplete or even incorrect stereochemistry. We made a large effort to correct the errors and to complete the structural descriptions. A total of 1468 structures were corrected and/or completed. We also established the reaction balance of the reactions in BioMeta and corrected 55% of the unbalanced (stoichiometrically incorrect) reactions in an automatic procedure. The BioMeta database was implemented in PostgreSQL and provided with a web-based interface. CONCLUSION: We demonstrate that the validation of metabolite structures and reactions is a feasible and worthwhile undertaking, and that the validation results can be used to trigger corrections and improvements to BioMeta, our metabolite database. BioMeta provides some tools for rational drug design, reaction searches, and visualization. It is freely available at http://www.cmbi.ru.nl/biometa/ provided that the copyright notice of all original data is cited. The database will be useful for querying and browsing biochemical pathways, and to obtain reference information for identifying compounds. However, these applications require that the underlying data be correct, and that is the focus of BioMeta.

Computational Biology↗

Costs and benefits of processivity in enzymatic degradation of recalcitrant polysaccharides.

Many enzymes that hydrolyze insoluble crystalline polysaccharides such as cellulose and chitin guide detached single-polymer chains through long and deep active-site clefts, leading to processive (stepwise) degradation of the polysaccharide. We have studied the links between enzyme efficiency and processivity by analyzing the effects of mutating aromatic residues in the substrate-binding groove of a processive chitobiohydrolase, chitinase B from Serratia marcescens. Mutation of two tryptophan residues (Trp-97 and Trp-220) close to the catalytic center (subsites +1 and +2) led to reduced processivity and a reduced ability to degrade crystalline chitin, suggesting that these two properties are linked. Most remarkably, the loss of processivity in the W97A mutant was accompanied by a 29-fold increase in the degradation rate for single-polymer chains as present in the soluble chitin-derivative chitosan. The properties of the W220A mutant showed a similar trend, although mutational effects were less dramatic. Processivity is thought to contribute to the degradation of crystalline polysaccharides because detached single-polymer chains are kept from reassociating with the solid material. The present results show that this processivity comes at a large cost in terms of enzyme speed. Thus, in some cases, it might be better to focus strategies for enzymatic depolymerization of polysaccharide biomass on improving substrate accessibility for nonprocessive enzymes rather than on improving the properties of processive enzymes.

Chitinases↗

Negative constraints underlie the ErbB specificity of epidermal growth factor-like ligands.

Epidermal growth factor (EGF)-like growth factors bind their ErbB receptors in a highly selective manner, but the molecular basis for this specificity is poorly understood. We have previously shown that certain residues in human EGF (Ser(2)-Asp(3)) and TGFalpha (Glu(26)) are not essential for their binding to ErbB1 but prevent binding to ErbB3 and ErbB4. In the present study, we have used a phage display approach to affinity-optimize the C-terminal linear region of EGF-like growth factors for binding to each ErbB receptor and thereby shown that Arg(45) in EGF impairs binding to both ErbB3 and ErbB4. By omitting all these so-called negative constraints from EGF, we designed a ligand designated panerbin that binds ErbB1, ErbB3, and ErbB4 with similarly high affinity as their wild-type ligands. Homology models, based on the known crystal structure of TGFalpha-bound ErbB1, showed that panerbin is able to bind ErbB1, ErbB3, and ErbB4 in a highly similar manner with respect to position and number of interaction sites. Upon in silico introduction of the experimentally known negative constraints into panerbin, we found that Arg(45) induced local charge repulsion and Glu(26) induced steric hindrance in a receptor-specific manner, whereas Ser(2)-Asp(3) impaired binding due to a disordered conformation. Furthermore, radiolabeled panerbin was used to quantify the level of all three receptors on human breast cancer cells in a single radioreceptor assay. It is concluded that the ErbB specificity of EGF-like growth factors primarily results from the presence of a limited number of residues that impair the unintended interaction with other ErbB receptors.

Amino Acid Sequence↗

Distinct clinical phenotypes associated with a mutation in the mitochondrial translation elongation factor EFTs.

The 13 polypeptides encoded in mitochondrial DNA (mtDNA) are synthesized in the mitochondrial matrix on a dedicated protein-translation apparatus that resembles that found in prokaryotes. Here, we have investigated the genetic basis for a mitochondrial protein-synthesis defect associated with a combined oxidative phosphorylation enzyme deficiency in two patients, one of whom presented with encephalomyopathy and the other with hypertrophic cardiomyopathy. Sequencing of candidate genes revealed the same homozygous mutation (C997T) in both patients in TSFM, a gene coding for the mitochondrial translation elongation factor EFTs. EFTs functions as a guanine nucleotide exchange factor for EFTu, another translation elongation factor that brings aminoacylated transfer RNAs to the ribosomal A site as a ternary complex with guanosine triphosphate. The mutation predicts an Arg333Trp substitution at an evolutionarily conserved site in a subdomain of EFTs that interacts with EFTu. Molecular modeling showed that the substitution disrupts local subdomain structure and the dimerization interface. The steady-state levels of EFTs and EFTu in patient fibroblasts were reduced by 75% and 60%, respectively, and the amounts of assembled complexes I, IV, and V were reduced by 35%-91% compared with the amounts in controls. These phenotypes and the translation defect were rescued by retroviral expression of either EFTs or EFTu. These data clearly establish mutant EFTs as the cause of disease in these patients. The fact that the same mutation is associated with distinct clinical phenotypes suggests the presence of genetic modifiers of the mitochondrial translation apparatus.

Amino Acid Sequence↗

Fast empirical pKa prediction by Ewald summation.

pK(a) calculations for macromolecules are normally performed by solving the Poisson-Boltzmann equation, accounting for the different dielectric constants of solvent and solute, as well as the ionic strength. Despite the large number of successful applications, there are some situations where the current algorithms are not suitable: (1) large scale, high-throughput analysis which requires calculations to be completed within a fraction of a second, e.g. when permanently monitoring pK(a) shifts during a molecular dynamics simulation; (2) prediction of pK(a)s in periodic boundaries, e.g. when reconstructing entire protein crystal unit cells from PDB files, including the correct protonation patterns at experimental pH. Such in silico crystals are needed by 'self-parameterizing' molecular dynamics force fields like YASARA YAMBER, that optimize their parameters while energy-minimizing high-resolution protein crystals. To address both problems, we define an empirical equation that expresses the pK(a) as a function of electrostatic potential, hydrogen bonds and accessible surface area. The electrostatic potential is evaluated by Ewald summation, which captures periodic crystal environments and the uncertainty in atom positions using Gaussian charge densities. The empirical proportionality constants are derived from 217 experimentally determined pK(a)s, and despite its simplicity, this pK(a) calculation method reaches a high overall jack-knifed accuracy, and is fast enough to be used during a molecular dynamics simulation. A reliable null-model to judge pK(a) prediction accuracies is also presented.

Algorithms↗

GRIS: glycoprotein-hormone receptor information system.

The glycoprotein-hormone receptor information system (GRIS) presents a comprehensive view on all available molecular data for the lutropin/choriogonadotropin receptor, follitropin receptor, and thyrotropin receptor G protein-coupled receptors. It features a mutation database presently containing 696 point mutations, combined with all sequences and the associated homology models. The mutation information was automatically extracted from the literature and manually augmented with respect to constitutivity, surface expression, sensitivity to hormones, and binding affinity. All information in this integrated system is presented in a G protein-coupled receptor specialist-friendly way. A series of interactive tools such as rotamer analysis, mutation prediction, or cavity visualization aids with the design and interpretation of experiments. A universal residue numbering system has been introduced to ease database searches as well as the use of the information in conjunction with literature data from diverse origins. Users can upload new mutations. GRIS is freely accessible at http://gris.ulb.ac.be/.

Amino Acid Sequence↗

Traditional biomolecular structure determination by NMR spectroscopy allows for major errors.

One of the major goals of structural genomics projects is to determine the three-dimensional structure of representative members of as many different fold families as possible. Comparative modeling is expected to fill the remaining gaps by providing structural models of homologs of the experimentally determined proteins. However, for such an approach to be successful it is essential that the quality of the experimentally determined structures is adequate. In an attempt to build a homology model for the protein dynein light chain 2A (DLC2A) we found two potential templates, both experimentally determined nuclear magnetic resonance (NMR) structures originating from structural genomics efforts. Despite their high sequence identity (96%), the folds of the two structures are markedly different. This urged us to perform in-depth analyses of both structure ensembles and the deposited experimental data, the results of which clearly identify one of the two models as largely incorrect. Next, we analyzed the quality of a large set of recent NMR-derived structure ensembles originating from both structural genomics projects and individual structure determination groups. Unfortunately, a visual inspection of structures exhibiting lower quality scores than DLC2A reveals that the seriously flawed DLC2A structure is not an isolated incident. Overall, our results illustrate that the quality of NMR structures cannot be reliably evaluated using only traditional experimental input data and overall quality indicators as a reference and clearly demonstrate the urgent need for a tight integration of more sophisticated structure validation tools in NMR structure determination projects. In contrast to common methodologies where structures are typically evaluated as a whole, such tools should preferentially operate on a per-residue basis.

Amino Acid Sequence↗

Integrins of Anopheles gambiae and a putative role of a new beta integrin, BINT2, in phagocytosis of E. coli.

Mosquitoes use effective immune responses, including phagocytosis, to fight microbial infection. Here we show that in an Anopheles gambiae immune responsive cell line, RGD recognizing receptors play an important role in the phagocytic response, suggesting overlap between molecular components implicated in adhesion and phagocytosis. Integrins are a major class of adhesive receptors that recognize ligands containing an RGD motif. We have cloned a gene encoding a new beta integrin, BINT2, and demonstrated its involvement in Escherichia coli engulfment. Based on molecular modeling, we propose a structural reason for the role of BINT2, but not BINT1, on phagocytosis of Gram-negative bacteria. Using bioinformatic tools, we have identified and compared the complete A. gambiae integrin repertoire as a prelude to a future systematic functional study.

Animals↗

A text-mining analysis of the human phenome.

A number of large-scale efforts are underway to define the relationships between genes and proteins in various species. But, few attempts have been made to systematically classify all such relationships at the phenotype level. Also, it is unknown whether such a phenotype map would carry biologically meaningful information. We have used text mining to classify over 5000 human phenotypes contained in the Online Mendelian Inheritance in Man database. We find that similarity between phenotypes reflects biological modules of interacting functionally related genes. These similarities are positively correlated with a number of measures of gene function, including relatedness at the level of protein sequence, protein motifs, functional annotation, and direct protein-protein interaction. Phenotype grouping reflects the modular nature of human disease genetics. Thus, phenotype mapping may be used to predict candidate genes for diseases as well as functional relations between genes and proteins. Such predictions will further improve if a unified system of phenotype descriptors is developed. The phenotype similarity data are accessible through a web interface at http://www.cmbi.ru.nl/MimMiner/.

Chromosome Mapping↗

GeneSeeker: extraction and integration of human disease-related information from web-based genetic databases.

The identification of genes underlying human genetic disorders requires the combination of data related to cytogenetic localization, phenotypes and expression patterns, to generate a list of candidate genes. In the field of human genetics, it is normal to perform this combination analysis by hand. We report on GeneSeeker (http://www.cmbi.ru.nl/GeneSeeker/), a web server that gathers and combines data from a series of databases. All database searches are performed via the web interfaces provided with the original databases, guaranteeing that the most recent data are queried, and obviating data warehousing. GeneSeeker makes the same selection of candidate genes as the human geneticists would have performed, and thus reducing the time-consuming process to a few minutes. GeneSeeker is particularly well suited for syndromes in which the disease gene displays altered expression patterns in the affected tissue(s).

Chromosome Mapping↗

Induced refolding of a temperature denatured llama heavy-chain antibody fragment by its antigen.

In a previous study we have shown that llama VHH antibody fragments are able to bind their antigen after a heat shock of 90 degrees C, in contrast to the murine monoclonal antibodies. However, the molecular mechanism by which antibody:antigen interaction occurs under these extreme conditions remains unclear. To examine in more detail the structural and thermodynamic aspects of the binding mechanism, an extensive CD, ITC, and NMR study was initiated. In this study the interaction between the llama VHH -R2 fragment and its antigen, the dye Reactive Red-6 (RR6) has been explored. The data show clearly that most of the VHH-R2 population at 80 degrees C is in an unfolded conformation. In contrast, CD spectra representing the complex between VHH-R2 and the dye remained the same up to 80 degrees C. Interestingly, addition of the dye to the denatured VHH-R2 at 80 degrees C yielded the spectrum of the native complex. These results suggest an induced refolding of denatured VHH-R2 by its antigen under these extreme conditions. This induced refolding showed some similarities with the well established "induced fit" mechanism of antibody-antigen interactions at ambient temperature. However, the main difference with the "induced fit" mechanism is that at the start of the addition of the antigen most of the VHH molecules are in an unfolded conformation. The refolding capability under these extreme conditions and the stable complex formation make VHHs useful in a wide variety of applications.

Animals↗

Definition of a new information-based per-residue quality parameter.

For biomolecular NMR structures typically only a poor correspondence is observed between statistics derived from the experimental input data and structural quality indicators obtained from the structure ensembles. Here, we investigate the relationship between the amount of available NMR data and structure quality. By generating datasets with a predetermined information content and evaluating the quality of the resulting structure ensembles we show that there is, in contrast to previous findings, a linear relation between the information contained in experimental data and structural quality. From this relation, a new quality parameter is derived that provides direct insight, on a per-residue basis, into the extent to which structural quality is governed by the experimental input data.

Animals↗

Isolation of llama antibody fragments for prevention of dandruff by phage display in shampoo.

As part of research exploring the feasibility of using antibody fragments to inhibit the growth of organisms implicated in dandruff, we isolated antibody fragments that bind to a cell surface protein of Malassezia furfur in the presence of shampoo. We found that phage display of llama single-domain antibody fragments (VHHs) can be extended to very harsh conditions, such as the presence of shampoo containing nonionic and anionic surfactants. We selected several VHHs that bind to the cell wall protein Malf1 of M. furfur, a fungus implicated in causing dandruff. In addition to high stability in the presence of shampoo, these VHHs are also stable under other denaturing conditions, such as high urea concentrations. Many of the stable VHHs were found to contain arginine at position 44. Replacement of the native amino acid at position 44 with arginine in the most stable VHH that lacked this arginine resulted in a dramatic further increase in the stability. The combination of the unique properties of VHHs together with applied phage display and protein engineering is a powerful method for obtaining highly stable VHHs that can be used in a wide range of applications.

Amino Acid Sequence↗

Making optimal use of empirical energy functions: force-field parameterization in crystal space.

Today's energy functions are not able yet to distinguish reliably between correct and almost correct protein models. Improving these near-native models is currently a major bottle-neck in homology modeling or experimental structure determination at low resolution. Increasingly accurate energy functions are required to complete the "last mile of the protein folding problem," for example during a molecular dynamics simulation. We present a new approach to reach this goal. For 50 high resolution X-ray structures, the complete unit cell was reconstructed, including disordered water molecules, counter ions, and hydrogen atoms. Simulations were then run at the pH at which the crystal was solved, while force-field parameters were iteratively adjusted so that the damage done to the structures was minimal. Starting with initial parameters from the AMBER force field, the optimization procedure converged at a new force field called YAMBER (Yet Another Model Building and Energy Refinement force field), which is shown to do significantly less damage to X-ray structures, often move homology models in the right direction, and occasionally make them look like experimental structures. Application of YAMBER during the CASP5 structure prediction experiment yielded a model for target 176 that was ranked first among 150 submissions. Due to its compatibility with the well-established AMBER format, YAMBER can be used by almost any molecular dynamics program. The parameters are freely available from www.yasara.org/yamber.

Computer Simulation↗

Rational engineering of enzyme stability.

During the past 15 years there has been a continuous flow of reports describing proteins stabilized by the introduction of mutations. These reports span a period from pioneering rational design work on small enzymes such as T4 lysozyme and barnase to protein design, and directed evolution. Concomitantly, the purification and characterization of naturally occurring hyperstable proteins has added to our understanding of protein stability. Along the way, many strategies for rational protein stabilization have been proposed, some of which (e.g. entropic stabilization by introduction of prolines or disulfide bridges) have reasonable success rates. On the other hand, comparative studies and efforts in directed evolution have revealed that there are many mutational strategies that lead to high stability, some of which are not easy to define and rationalize. Recent developments in the field include increasing awareness of the importance of the protein surface for stability, as well as the notion that normally a very limited number of mutations can yield a large increase in stability. Another development concerns the notion that there is a fundamental difference between the "laboratory stability" of small pure proteins that unfold reversibly and completely at high temperatures and "industrial stability", which is usually governed by partial unfolding processes followed by some kind of irreversible inactivation process (e.g. aggregation). Provided that one has sufficient knowledge of the mechanism of thermal inactivation, successful and efficient rational stabilization of enzymes can be achieved.

Biotechnology↗

DRESS: a database of REfined solution NMR structures.

Several studies have shown that biomolecular NMR structures are often of lower quality when compared to crystal structures, and consequently they are often excluded from structural analyses. We present a publicly available database of re-refined NMR structures, exhibiting significantly improved quality. This database (available at http://www.cmbi.kun.nl/dress/) presents a uniformly refined and validated set of structural models that improves the value of these NMR structures as input for experimental and theoretical studies in many fields of research.

Databases, Protein↗

Molecular determinants of the interaction of Mad with the PAH2 domain of mSin3.

The Sin3 co-repressor acts as a protein scaffold to recruit transcription factors via its four highly homologous paired amphipathic helix (PAH) domains. PAH2 has been shown to interact strongly with the Sin3 interacting domain (SID) of the tumor suppressor Mad. This PAH2/Mad complex has been studied extensively by NMR, but the molecular determinants that dictate the specificity of interaction remain to be elucidated. To uncover residues that convey the specificity of interaction between PAH2 and Mad, PAH2 residues contacted by the Mad-SID were introduced into the PAH1 domain of mSin3b and tested for gain-of-interaction in vivo in a yeast two-hybrid setting and further confirmed in a cell-free system. This approach led to the identification of PAH2-Phe-7 as a critical residue. Stabilization of the interaction between PAH1-Phe-7 and the Mad-SID was achieved by introducing Val-14 and Gln-39 into PAH1. Substitution of PAH2 residues contacted by the Mad-SID with their respective residues in PAH1 corroborated and extended the critical role of Phe-7 and the stabilizing role of Val-14 and Gln-39. We conclude that Phe-7 is the critical determinant and provides the molecular specificity for the association between Sin3 and Mad in regulating cell growth and differentiation.

Amino Acid Sequence↗