Search PubMed⌕ Search

Biomedical subjects

Martin Stahl

Publications and source records attributed to Martin Stahl.

At least 19 recordsLinked to original sources

Dose-dependent antithrombotic activity of an orally active tissue factor/factor VIIa inhibitor without concomitant enhancement of bleeding propensity.

The discovery of a highly potent and selective tissue factor/factor VIIa inhibitor is described. Upon oral administration of its double prodrug in the guinea pig, a dose-dependent antithrombotic effect is observed in an established model of arterial thrombosis without prolonging bleeding time. The pharmacodynamic properties of this selective inhibitor are compared to the behaviour of a mixed factor VIIa/factor Xa inhibitor.

Administration, Oral↗

Integrating molecular design resources within modern drug discovery research: the Roche experience.

Various computational disciplines, such as cheminformatics, ADME modeling, virtual screening, chemogenomics search strategies and classic structure-based design, should be seen as one multifaceted discipline contributing to the early drug discovery process. Although significant resources enabling these activities have been established, their true integration into daily research should not be taken for granted. This article reviews value-adding activities from target assessment to lead optimization, and highlights the technical and process-related aspects that can be considered essential for performance and alignment within the research organization.

Animals↗

Comparison of the binding pockets of two chemically unrelated allosteric antagonists of the mGlu5 receptor and identification of crucial residues involved in the inverse agonism of MPEP.

Fenobam [N-(3-chlorophenyl)-N'-(4,5-dihydro-1-methyl-4-oxo-1H-imidazole-2-yl)urea], a clinically validated non-benzodiazepine anxiolytic, has been shown to be a potent and non-competitive metabotropic glutamate (mGlu)-5 receptor antagonist. In the present study, we have used the site-directed mutagenesis coupled with three-dimensional receptor-based pharmacophore modelling to elucidate the interacting mode of fenobam within the seven-transmembrane domain (7TMD) of mGlu5 receptor and its comparison with that of 2-methyl-6-(phenylethynyl)pyridine (MPEP), the prototype antagonist. The common residues involved in the recognition of MPEP and fenobam include Pro654(3.36), Tyr658(3.40), Thr780(6.44), Trp784(6.48), Phe787(6.51), Tyr791(6.55) and Ala809(7.47). The differentiating residues between both modulators' interacting modes are Arg647(3.29), Ser657(3.39) and Leu743(5.47). Our data suggest that these chemically unrelated mGlu5 antagonists act similarly, probing a functionally unique region of the 7TMD. Using [3H]inositol phosphates accumulation assay, we have also identified the critical residues involved in the inverse agonist effect of MPEP. The mutation W784(6.48)A completely blocked the inverse agonist activity of MPEP; two mutations F787(6.51)A and Y791(6.55)A, caused a drastic decrease in the MPEP inverse agonism. Furthermore, these three mutations led to an increased efficacy of quisqualate without having any effect on its potency. The fact that the residues Trp784(6.48) and Phe787(6.51) are essential equally in antagonism and inverse agonism effects emphasizes again the key role of these residues and the involvement of a common transmembrane network in receptor inactivation by MPEP.

Binding Sites↗

The environment of amide groups in protein-ligand complexes: H-bonds and beyond.

A comprehensive structural analysis of interactions involving amide NH and C=O groups in protein-ligand complexes has been performed based on 3,275 published crystal structures (resolution < or =2.5 A). Most of the amide C=O and NH groups at the protein-ligand interface are highly buried within the binding site and involved in H-bonds with corresponding counter-groups. Small percentages of C=O and NH groups are solvated or embedded in hydrophobic environments. In particular, C=O groups show a higher propensity to be solvated or embedded in a hydrophobic environment than NH groups do. A small percentage of carbonyl groups is involved in weak hydrogen bonds with CH. Cases of dipolar interactions, involving carbonyl oxygen and electrophilic carbon atoms, such as amide, amidinium, guanidium groups, are also identified. A higher percentage of NH are in contact with aromatic carbons, interacting either through hydrogen bonds (preferably with the NH group pointing towards a ring carbon atom) or through stacking between amide plane and ring plane. Comprehensive studies such as the present one are thought to be important for future improvements in the molecular design area, in particular for the development of new scoring functions. [Figure: see text].

Amides↗

Focused library design in GPCR projects on the example of 5-HT(2c) agonists: comparison of structure-based virtual screening with ligand-based search methods.

The aim of this study was to investigate the usefulness of structure-based virtual screening (VS) for focused library design in G protein-coupled receptors (GPCR) projects on the example of 5-HT(2c) agonists. We compared the performance of structure-based VS against two different homology models using FRED for docking and ScreenScore, FlexX, and PMF for rescoring with the results of 12 ligand-based similarity searches using four different query compounds and three different similarity metrics (Daylight, FTree, Phacir). The result of the similarity search showed much variation, from an enrichment factor up to 3.2 to worse than random, whereas the structure-based VS gave a more stable result with a constant enrichment factor around 2. Additionally, actives retrieved by the structure-based approach were more diverse than the actives among the top scorers of the similarity searches. Based on these results, we suggest basing a focused library design for a GPCR project on a combination of a ligand-based similarity search and structure-based docking.

Cluster Analysis↗

A robust clustering method for chemical structures.

A clustering method based on finding the largest set of disconnected fragments that two chemical compounds have in common is shown to be able to group structures in a way that is ideally suited to medicinal chemistry programs. We describe how markedly improved results can be obtained by using a similarity metric that accounts not just for the size of the shared fragments but also on their relative arrangement in the two parent compounds. The use of a physiochemical atom typing scheme is also shown to provide significant contributions. Results from calculations using a test set consisting of actives from nine different important biological target proteins demonstrate the strengths of our clustering method and the advantages over other approaches that are widely used throughout the pharmaceutical industry.

Chemistry, Organic↗

Validation and use of the MM-PBSA approach for drug discovery.

The MM-PBSA approach has become a popular method for calculating binding affinities of biomolecular complexes. Published application examples focus on small test sets and few proteins and, hence, are of limited relevance in assessing the general validity of this method. To further characterize MM-PBSA, we report on a more extensive study involving a large number of ligands and eight different proteins. Our results show that applying the MM-PBSA energy function to a single, relaxed complex structure is an adequate and sometimes more accurate approach than the standard free energy averaging over molecular dynamics snapshots. The use of MM-PBSA on a single structure is shown to be valuable (a) as a postdocking filter in further enriching virtual screening results, (b) as a helpful tool to prioritize de novo design solutions, and (c) for distinguishing between good and weak binders (DeltapIC(50) > or = 2-3), but rarely to reproduce smaller free energy differences.

Binding Sites↗

Prediction of pKa shifts in proteins using a combination of molecular mechanical and continuum solvent calculations.

The prediction of pKa shifts of ionizable groups in proteins is of great relevance for a number of important biological phenomena. We present an implementation of the MM-GBSA approach, which combines molecular mechanical (MM) and generalized Born (GB) continuum solvent energy terms, to the calculation of pKa values of a panel of nine proteins, including 69 individual comparisons with experiment. While applied so far mainly to the calculation of biomolecular binding free energies, we show that this method can also be used for the estimation of protein pKa shifts, with an accuracy around 1 pKa unit, even for strongly shifted residues. Our analysis reveals that the nonelectrostatic terms that are part of the MM-GBSA free energy expression are important contributors to improved prediction accuracy. This suggests that most of the previous approaches that focus only on electrostatic interactions could be improved by adding other nonpolar energy terms to their free energy expression. Interestingly, our method yields best accuracy at protein dielectric constants of epsilonint = 2-4, which is in contrast to previous approaches that peak at higher epsilonint > or = 8. An important component of our procedure is an intermediate minimization step of each protonation state involving different rotamers and tautomers as a way to explicitly model protein relaxation upon (de)protonation.

Algorithms↗

Fluorine in medicinal chemistry.

Fluorinated compounds are synthesized in pharmaceutical research on a routine basis and many marketed compounds contain fluorine. The present review summarizes some of the most frequently employed strategies for using fluorine substituents in medicinal chemistry. Quite often, fluorine is introduced to improve the metabolic stability by blocking metabolically labile sites. However, fluorine can also be used to modulate the physicochemical properties, such as lipophilicity or basicity. It may exert a substantial effect on the conformation of a molecule. Increasingly, fluorine is used to enhance the binding affinity to the target protein. Recent 3D-structure determinations of protein complexes with bound fluorinated ligands have led to an improved understanding of the nonbonding protein-ligand interactions that involve fluorine.

Fluorine↗

Novel dihydrofolate reductase inhibitors. Structure-based versus diversity-based library design and high-throughput synthesis and screening.

Novel 2,4-diaminopyrimidines bearing N,N-disubstituted aminomethyl residues at the 5-position were designed as dihydrofolate reductase (DHFR) inhibitors. These compounds were obtained by treatment of 1-[(2,4-diamino-5-pyrimidinyl)methyl]pyridinium bromide with secondary amines in a polar solvent and in the presence of triethylamine at room temperature. The procedure was found to be very efficient and suitable for application in high-throughput synthesis. In addition, we found that high-throughput screening for enzymatic and in vitro antibacterial activity could be performed on crude reaction mixtures, thus avoiding any purification step. Over 1200 proprietary secondary amines were selected for high-throughput synthesis, based on structural and diversity-related criteria, and the resulting products were submitted to high-throughput screening. A greater number of hits, and significantly more active compounds, were obtained through structure-based library design than through diversity-based library design. Different classes of inhibitors of DHFR were identified in this way, including compounds derived from di-, tri-, and tetracyclic amines. In general, these products showed high activity against the enzymes derived from both TMP-sensitive and TMP-resistant Streptococcus pneumoniae. Some compounds possessed appreciable selectivity for the bacterial over the human enzyme, whereas other compounds were not at all selective. In most cases, active enzyme inhibitors also displayed antibacterial activity.

Amino Acid Sequence↗

Mechanistic insights into oxidosqualene cyclizations through homology modeling.

2,3-Oxidosqualene cyclases (OSC) are key enzymes in sterol biosynthesis. They catalyze the stereoselective cyclization and skeletal rearrangement of (3S)-2,3-oxidosqualene to lanosterol in mammals and fungi and to cycloartenol in algae and higher plants. Sequence information and proposed mechanism of 2,3-oxidosqualene cyclases are closely related to those of squalene-hopene cyclases (SHC), which represent functional analogs of OSCs in bacteria. SHCs catalyze the cationic cyclization cascade converting the linear triterpene squalene to fused ring compounds called hopanoids. High stereoselectivity and precision of the skeletal rearrangements has aroused the interest of researchers for nearly half a century, and valuable data on studying mechanistic details in the complex enzyme-catalyzed cyclization cascade has been collected. Today, interest in cyclases is still unbroken, because OSCs became targets for the development of antifungal and hypocholesterolemic drugs. However, due to the large size and membrane-bound nature of OSCs, three-dimensional structural information is still not available, thus preventing a complete understanding of the atomic details of the catalytic mechanism. In this work, we discuss results gained from homology modeling of human OSC based on structural information of SHC from Alicyclobacillus acidocaldarius and propose a structural model of human OSC. The model is in accordance with previously performed experimental studies with mechanism-based suicide inhibitors and mutagenesis experiments with altered activity and product specificity. Structural insight should strongly stimulate structure-based design of antifungal or cholesterol-lowering drugs.

Amino Acid Sequence↗

Binding site characteristics in structure-based virtual screening: evaluation of current docking tools.

Two new docking programs FRED (OpenEye Scientific Software) and Glide (Schrödinger, Inc.) in combination with various scoring functions implemented in these programs have been evaluated against a variety of seven protein targets (cyclooxygenase-2, estrogen receptor, p38 MAP kinase, gyrase B, thrombin, gelatinase A, neuraminidase) in order to assess their accuracy in virtual screening. Sets of known inhibitors were added to and ranked relative to a random library of drug-like compounds. Performance was compared in terms of enrichment factors and CPU time consumption. Results and specific features of the two new tools are discussed and compared to previously published results using FlexX (Tripos, Inc.) as a docking engine. In addition, general criteria for the selection of docking algorithms and scoring functions based on binding-site characteristics of specific protein targets are proposed. Figure Enrichment factors obtained with FlexX, Glide and FRED docking engines in combination with different scoring functions for seven selected targets with highly variable binding sites

Algorithms↗

A validation study on the practical use of automated de novo design.

The de novo design program Skelgen has been used to design inhibitor structures for four targets of pharmaceutical interest. The designed structures are compared to modeled binding modes of known inhibitors (i) visually and (ii) by means of a novel similarity measure considering the size and spatial proximity of the maximum common substructure of two small molecules. It is shown that the Skelgen algorithm generates representatives of many inhibitor classes within a very short time and that the new similarity measure is useful for comparing and clustering designed structures. The results demonstrate the necessity of properly defining search constraints in practical applications of de novo design.

Algorithms↗

Multiscale contrast enhancement for radiographies: Laplacian Pyramid versus fast wavelet transform.

Contrast enhancement of radiographies based on a multiscale decomposition of the images recently has proven to be a far more versatile and efficient method than regular unsharp-masking techniques, while containing these as a subset. In this paper, we compare the performance of two multiscale-methods, namely the Laplacian Pyramid and the fast wavelet transform (FWT). We find that enhancement based on the FWT suffers from one serious drawback-the introduction of visible artifacts when large structures are enhanced strongly. By contrast, the Laplacian Pyramid allows a smooth enhancement of large structures, such that visible artifacts can be avoided. Only for the enhancement of very small details, for denoising applications or compression of images, the FWT may have some advantages over the Laplacian Pyramid.

Algorithms↗

Assignment of Relative Configuration to Acyclic Compounds Based on (13)C NMR Shifts. A Density Functional and Molecular Mechanics Study.

1,3-Dimethylated hydrocarbon segments occur frequently as structural elements in polyketide natural products. The (13)C NMR chemical shifts of a series of model compounds containing such segments can be well reproduced by a combination of molecular mechanics and SOS-DFPT/IGLO calculations. (13)C NMR chemical shifts are calculated on MM3 geometries and are Boltzmann weighted according to the MM3 energies. On the basis of the resulting thermally averaged chemical shifts, all diastereomers of the model compounds can be unequivocally distinguished. Significant differences in chemical shifts occur at methyl groups and methylene groups that are adjacent to a single stereogenic center. The method is applied to predict the relative configuration of two stereocenters in the side chains of two natural products, sambutoxin and the bradykinin inhibitor L-755,897.

Journal Article↗

Prediction of UV and ESI-MS signal intensities.

All major pharmaceutical companies maintain large collections of compounds that are used either for screening against biological targets or as synthetic precursors. The quality assessment of these compounds is typically done by liquid chromatography combined with mass spectroscopy (LC/MS) and UV purity control. To facilitate the analysis of the analytical data, we have built computational models to predict UV and MS signal intensities under experimental LC/MS conditions. The discriminant partial-least-squares technique was used for classifying compounds into those most likely to yield a MS signal and others where the signal is below the detection limit (94% and 88% correct predictions, respectively). In the case of UV prediction, we compared this statistical linear-regression technique to a knowledge-based approach. A combination of both techniques proved to be the most reliable (96/98% correct predictions of UV-active/ UV-inactive compounds). Both models have been incorporated into the automated compound integrity profiling at F. Hoffmann-La Roche.

Chromogenic Compounds↗