Biomedical subjects

P D Grootenhuis

Publications and source records attributed to P D Grootenhuis.

At least 19 recordsLinked to original sources

PubMed2000

A rapid computational method for lead evolution: description and application to alpha(1)-adrenergic antagonists.

The high failure rate of drugs in the development phase requires a strategy to reduce risks by generating lead candidates from different chemical classes. We describe a new three-dimensional computational approach for lead evolution, based on multiple pharmacophore hypotheses. Using full conformational models for both active and inactive compounds, a large number of pharmacophore hypotheses are analyzed to select the set or "ensemble" of hypotheses that, when combined, is most able to discriminate between active and inactive molecules. The ensemble hypothesis is then used to search virtual chemical libraries to identify compounds for synthesis. This method is very rapid, allowing very large virtual libraries on the order of a million compounds to be filtered efficiently. In applying this method to alpha(1)-adrenergic receptor ligands, we have demonstrated lead evolution from heterocyclic alpha(1)-adrenergic receptor ligands to highly dissimilar active N-substituted glycine compounds. Our results also show that the active N-substituted glycines are part of our smaller filtered library and thus could have been identified by synthesizing only a portion of the N-substituted glycine library.

Adrenergic alpha-Antagonists↗

PubMed1999

Computer simulations of the dynamics of human choriogonadotropin and its alpha subunit.

Human choriogonadotropin (hCG) belongs to a family of heterodimeric glycoprotein hormones involved in reproduction. Over 75 ns of molecular dynamics simulations of this heterodimer and the free alpha subunit were performed and validated by experimental information to arrive at a qualitative dynamical description of these molecules. A number of 5-ns simulations at 400 degrees K describe a sufficiently stable heterodimer structure, whereas the free alpha subunit shows the experimentally observed partial unfolding. From the main collective fluctuations of the free alpha subunit, it can be derived that residues alpha35-55 form a domain that is highly flexible with respect to the other domain, which contains all five disulfide bonds. The apparent loss of secondary structure in the region alpha33-58 may very well be induced by this. Dynamic domains can also be determined from the hCG heterodimer simulations. The most important collective mode of motion shows that the flexibility of the alpha subunit is reduced by concerted rotation with both the long loop and the determinant loop of the beta subunit. The motion of the free alpha subunit does not differ significantly from the motion it has in the hCG heterodimer, but the amplitudes along the most important eigenvectors are larger.

Carbohydrates↗

PubMed1999

Comparative molecular field analysis and energy interaction studies of thrombin-inhibitor complexes.

A Comparative Molecular Field Analysis (CoMFA) and an interaction energy-based method were applied on a database holding the 3D structures of 29 thrombin-inhibitor complexes. Several parameters were optimized in both methods in order to obtain the best correlation between theoretical and experimentally determined binding (Ki) data. CoMFA, which only uses the information of the inhibitors, performed best (r = 0.99, q2 = 0.46, N = 29) when HF 6-31G charges were used in combination with a pharmacophore-based alignment. Inclusion of hydrophobic fields did not lead to improvements. The interaction energy-based approach uses the information of the whole thrombin-inhibitor complex. A statistically significant correlation (r = 0.74, N = 14) could only be obtained for a subset of the database containing the high resolution structures. Geometry optimization of the ligand only in combination with downscaled electrostatics performed best.

Antithrombins↗

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↗

PubMed1999

Polar molecular surface as a dominating determinant for oral absorption and brain penetration of drugs.

PURPOSE: To study oral absorption and brain penetration as a function of polar molecular surface area. METHODS: Measured brain penetration data of 45 drug molecules were investigated. The dynamic polar surface areas were calculated and correlated with the brain penetration data. Also the static polar surface areas of 776 orally administered CNS drugs that have reached at least Phase II efficacy studies were calculated. The same was done for a series of 1590 orally administered non-CNS drugs that have reached at least Phase II efficacy studies. RESULTS: A linear relationship between brain penetration and dynamic polar surface area (A2) was found (n = 45, R = 0.917, F1,43 = 229). Brain penetration decreases with increasing polar surface area. A clear difference between the distribution of the polar surface area of the 776 CNS and 1590 non-CNS drugs was found. It was deduced that orally active drugs that are transported passively by the transcellular route should not exceed a polar surface area of about 120 A2. They can be tailored to brain penetration by decreasing the polar surface to <60-70 A2. This conclusion is supported by the inverse linear relationship between experimental brain penetration data and the dynamic polar surface area of 45 drug molecules. CONCLUSIONS: The polar molecular surface area is a dominating determinant for oral absorption and brain penetration of drugs that are transported by the transcellular route. This property should be considered in the early phase of drug screening.

Animals↗

PubMed1999

Towards minimized gonadotropins with full bioactivity.

Gonadotropins are highly complex glycoprotein hormones consisting of two noncovalently associated subunits, which are heavily glycosylated. Using the X-ray structure of human choriogonadotropin and structure/activity relationships we aimed to design 'minimized' gonadotropins of reduced complexity. Our results show that it is possible to reduce the size of natural human choriogonadotropin by one-third of its molecular weight while retaining its wild-type biopotency. To our knowledge, such 'mini'-human choriogonadotropins represent the smallest gonadotropins described so far with an lutropin/choriogonadotropin receptor affinity and in vitro biological activity comparable with that of natural human choriogonadotropin. It provides an important step towards the structure/function-based design of small molecule drugs to the human gonadotropin receptors.

Animals↗

PubMed1999

Random mutagenesis and screening of complex glycoproteins: expression of human gonadotropins in Dictyostelium discoideum.

The soil amoeba Dictyostelium discoideum is a host cell that provides simple genetics in combination with complex protein synthesis. We show that the complex human heterodimeric gonadotropins can be produced and secreted by this organism. Furthermore, both follicle stimulation hormone and choriogonadotropin produced by D. dictyostelium bind to their human receptors and elicit a biological response comparable to the wild-type hormones. We also show that structure-function analysis using random mutagenesis and screening of recombinant glycoprotein hormones is feasible. Thus, expression of gonadotropins in D. dictyostelium opens the way to the engineering of potential new therapeutic analogues.

Animals↗

PubMed1998

Novel acylguanidine containing thrombin inhibitors with reduced basicity at the P1 moiety.

Replacement of the noragmatine group in thrombin inhibitors with a beta-alanyl-guanidine group resulted in a nearly equipotent and more selective compound 8 despite the fact that the pKa of this P1 moiety is five orders of magnitude lower. Further modification resulted in a nonpeptide inhibitor with this beta-alanyl-guanidine group, compound 28. This is an active and selective thrombin inhibitor and in view of its nonpeptide/low basicity structure selected for further pharmacological studies.

Antithrombins↗

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

Expression of a bioactive, single-chain choriogonadotropin in Dictyostelium discoideum.

Human choriogonadotropin (hCG) is a highly complex glycoprotein consisting of two non-covalently associated subunits. We aimed for the expression of a single-chain hCG in the soil amoebae Dictyostelium discoideum, a host which, in principle, provides simple genetics in combination with complex protein synthesis. To limit anticipated problems in mRNA translation, the first 30 bases of the coding sequence were altered to conform to the Dictyostelium preferred codon usage. We show that, immunologically, active single-chain hCG is indeed produced by Dictyostelium. Furthermore, this single-chain hCG is able to bind to the human luteinizing hormone/CG receptor and elicit a biological response. Its receptor-binding affinity is comparable to single-chain hCG produced by mammalian cells. We conclude that Dictyostelium is able to express bioactive highly complex heterologous glycoproteins.

Amino Acid Sequence↗

PubMed1998

Partially deglycosylated human choriogonadotropin, stabilized by intersubunit disulfide bonds, shows full bioactivity.

Several studies indicate that in human choriogonadotropin the N-linked oligosaccharide at position 52 of the alpha-subunit is important for bioactivity. We have generated choriogonadotropin mutants in which the alpha52 glycosylation site is removed and the alpha and beta subunits are covalently linked by intersubunit disulfide bonds. These mutants display wild-type receptor binding and bioactivity. Furthermore, we show that removal of the alpha52 sugar leads to instability of heterodimeric choriogonadotropin. Therefore, we conclude that the alpha52 oligosaccharide of choriogonadotropin is not involved in signal transduction, but in the stability of the heterodimer.

Animals↗

PubMed1997

X-ray structure of antistasin at 1.9 A resolution and its modelled complex with blood coagulation factor Xa.

The three-dimensional structure of antistasin, a potent inhibitor of blood coagulation factor Xa, from the Mexican leech Haementeria officinalis was determined at 1.9 A resolution by X-ray crystallography. The structure reveals a novel protein fold composed of two homologous domains, each resembling the structure of hirustasin, a related 55-residue protease inhibitor. However, hirustasin has a different overall shape than the individual antistasin domains, it contains four rather than two beta-strands, and does not inhibit factor Xa. The two antistasin domains can be subdivided into two similarly sized subdomains with different relative orientations. Consequently, the domain shapes are different, the N-terminal domain being wedge-shaped and the C-terminal domain flat. Docking studies suggest that differences in domain shape enable the N-terminal, but not C-terminal, domain of antistasin to bind and inhibit factor Xa, even though both have a very similar reactive site. Furthermore, a putative exosite binding region could be defined in the N-terminal domain of antistasin, comprising residues 15-17, which is likely to interact with a cluster of positively charged residues on the factor Xa surface (Arg222/Lys223/Lys224). This exosite binding region explains the specificity and inhibitory potency of antistasin towards factor Xa. In the C-terminal domain of antistasin, these exosite interactions are prevented due to the different overall shape of this domain.

Amino Acid Sequence↗

PubMed1997

The second Kunitz domain of human tissue factor pathway inhibitor: cloning, structure determination and interaction with factor Xa.

Tissue Factor Pathway Inhibitor (TFPI) is a 36 kDa glycoprotein that helps maintain haemostasis by inhibiting Factor Xa and the Factor VIIa/Tissue Factor (TF) complex. TFPI contains three tandemly linked Kunitz inhibitor domains, of which the second inhibits factor Xa. We have undertaken a multidisciplinary approach to study the structure and function of the second Kunitz domain of TFPI, with a view towards the rational design of factor Xa inhibitors. Amino acid residues 93 to 154 of the mature TFPI protein, corresponding to the second Kunitz domain (TFPI-kII), were expressed in Escherichia coli. The protein was purified to near homogeneity by ion exchange, hydrophobic interaction, and size exclusion chromatography, respectively. TFPI-kII is a potent factor Xa inhibitor with a Ki of 1.5 x 10(-10) M, a value that does not differ significantly from that of intact TFPI. The three-dimensional structure of TFPI-kII in aqueous solution was determined by 1H nuclear magnetic resonance spectroscopy (NMR). A set of 30 conformers was calculated with the program DIANA using 906 distance constraints derived from nuclear Overhauser effects and 23 dihedral angle constraints. This set, representing the solution structure of TFPI-kII, has an average root-mean-square deviation of 0.78 A for the backbone atoms and 1.38 A for all heavy atoms of residues 1 to 58. The structure of TFPI-kII has also been determined in complex with porcine trypsin using X-ray crystallographic techniques. The complex has been solved to a resolution of 2.6 A, with a final R-factor of 16.2%. Comparison of the NMR derived structure with that of TFPI-kII in complex with trypsin reveals little divergence of the two structures, with the exception of residue Tyr17. Superposition of the trypsin:TFPI-kII complex on factor Xa provides insights into macromolecular determinants for the inhibition of factor Xa. Complexation would require a degree of reorganisation of factor Xa residues, in particular of TyrF99, but also perhaps of the F148-loop. The interaction was further investigated using restrained molecular dynamics. Electrostatic interactions would appear to play a major role. The reorganisation of factor Xa is in contrast to the proposed factor Xa:TAP interaction, where TAP would bind to the "ground state" structure of factor Xa.

Amino Acid Sequence↗

PubMed1997

Evaluation of subunit truncation and the nature of the spacer for single chain human gonadotropins.

Three single chain gonadotropins were designed based on the three-dimensional-structure of human choriogonadotropin and structure/activity relationships of the glycoprotein hormones. In each single chain, the C-terminal end of the human choriogonadotropin beta subunit is connected via Ser-Gly repeats to the N-terminal end of the alpha subunit. In addition, two of the single chains have truncated subunits. The three mutants were expressed in CHO cells. In vitro binding of two of the three mutants to the human lutropin/choriogonadotropin receptor was found to be comparable to wild-type lutropin. In contrast, both the receptor binding and the in vitro bioactivity of the mutant with truncated alpha and beta subunits in which the beta:26-110 disulphide bond cannot be formed, are lowered relative to wild-type lutropin. The fact that this mutant still displays biological activity shows that the seat-belt arrangement proposed by Isaacs and coworkers [Lapthorn, A. J., Harris, D. C., Littlejohn, A., Lustbader, J. W., Canfield, R. E., Machin, K. J.. Morgan, F. J. & Isaacs, N. W. (1994) Nature 369, 455-461] is important but not essential for receptor binding and biological activity in the context of single chain gonadotropins. Single chains in which Ser-Gly spacers are combined with truncated subunits, provide an attractive approach towards the design and generation of novel, biologically active gonadotropins.

Animals↗

PubMed1997

Structure-based design and protein engineering of intersubunit disulfide bonds in gonadotropins.

Pairs of cystine residues were introduced in the alpha- and beta-subunits of human choriogonadotropin at positions with optimal geometries for the formation of disulfide bonds. Using the homology with luteinizing hormone and follicle stimulating hormone, similar mutations were carried out in these glycoprotein hormones. In nearly all mutants the corresponding disulfide bonds were formed leading to a non-natural, covalent linkage between the alpha- and beta-subunits. The mutants typically display wild-type receptor binding and bioactivity. The mutants with non-natural intersubunit disulfide bonds display enhanced thermostabilities relative to the corresponding heterodimeric glycoprotein hormones, rendering them candidates for long acting gonadotropins with enhanced shelf lives.

Animals↗

PubMed1996

Expression of biologically active fusion genes encoding the common alpha subunit and either the CG beta or FSH beta subunits: role of a linker sequence.

The gonadotropin/thyrotropin hormone family is characterized by a heterodimeric structure composed of a common alpha subunit non-covalently linked to a hormone-specific beta subunit. The conformation of the heterodimer is essential for controlling secretion, hormone-specific post-translational modifications and signal transduction. Structure-function studies of FSH and the other glycoprotein hormones are often hampered by mutagenesis induced defects in subunit combination. Thus, the ability to overcome the limitation of subunit assembly would expand the range of structure activity relationships that can be performed on these hormones. Here we converted the FSH heterodimer to a single chain by genetically fusing the carboxyl end of the FSH beta subunit to the amino end of the alpha subunit in the presence or absence of a natural linker sequence. In the absence of the CTP linker, the secretion rate was decreased over three fold. (The CTP sequence is the last 28 amino acids of the CG beta sequence and contains four serine-linked oligosaccharides). Unexpectedly however receptor binding/signal transduction was unaffected by absence of the linker. Molecular modelling of the tethers lacking the linker sequence show that the alignment of the alpha/beta domains in the single chain differ substantially from that seen in the heterodimer. These data show that the single chain FSH was secreted efficiently and is biologically active and that the conformation determinants required for secretion and biologic activity are not the same.

Animals↗

PubMed1996

Functionality map analysis of the active site cleft of human thrombin.

The Multiple Copy Simultaneous Search methodology has been used to construct functionality maps for an extended region of human thrombin, including the active site. This method allows the determination of energetically favorable positions and orientations for functional groups defined by the user on the three-dimensional surface of a protein. The positions of 10 functional group sites are compared with those of corresponding groups of four thrombin-inhibitor complexes. Many, but not all features, of known thrombin inhibitors are reproduced by the method. The results indicate that certain aspects of the binding modes of these inhibitors are not optimal. In addition, suggestions are made for improving binding by interaction with functional group sites on the thrombin surface that are not used by the thrombin inhibitors.

Amino Acid Chloromethyl Ketones↗

PubMed1995

Correlation of binding affinities with non-bonded interaction energies of thrombin-inhibitor complexes.

Several empirical modeling protocols are evaluated allowing a quantification of the interaction between an enzyme and a series of inhibitors. The evaluation and optimization of the modeling protocols used a database of 35 non-covalently bound inhibitors of human thrombin. Intermolecular interaction energies were calculated with CHARMm after energy minimization of the modeled complexes using various dielectric functions and constraining strategies. These calculated binding energies were correlated with the experimentally obtained binding constants of the inhibitors. The best protocols resulted in linear correlations with correlation coefficients > 0.80. In the best protocols the enzyme was fully constrained, the ligand was allowed to move freely and electrostatics were scaled down drastically or fully neglected during the energy minimization. For the interaction energy evaluation step, a distance-dependent dielectric function epsilon = R proved to be optimal. This simple empirical protocol, that neglects solvation or entropy effects, can be implemented readily in other force field packages and may be applied on various enzyme-inhibitor complexes, providing a tool for the evaluation and rank-ordering of newly designed inhibitors.

Journal Article↗