A Real-World Perspective on Molecular Design.

We present a series of small molecule drug discovery case studies where computational methods were prospectively employed to impact Roche research projects, with the aim of highlighting those methods that provide real added value. Our brief accounts encompass a broad range of methods and techniques applied to a variety of enzymes and receptors. Most of these are based on judicious application of knowledge about molecular conformations and interactions: filling of lipophilic pockets to gain affinity or selectivity, addition of polar substituents, scaffold hopping, transfer of SAR, conformation analysis, and molecular overlays.

Challenges and Rewards in Medicinal Chemistry Targeting Cardiovascular and Metabolic Diseases.

Medicinal chemistry has been transformed by major technological and conceptual innovations over the last three decades: structural biology and bioinformatics, structure and property based molecular design, the concepts of multidimensional optimization (MDO), in silico and experimental high-throughput molecular property analysis. The novel technologies advanced gradually and in synergy with biology and Roche has been at the forefront. Applications in drug discovery programs towards new medicines in cardiovascular and metabolic diseases are highlighted to show impact and advancement: the early discovery of endothelin antagonists for endothelial dysfunction (Bosentan), 11-beta hydroxysteroid dehydrogenase (11β-HSD1) inhibitors for dysregulated cellular glucocorticoid tonus (type 2 diabetes and metabolic syndrome) and non-covalent hormone sensitive lipase (HSL) inhibitors to study the scope of direct inhibition of lipolysis in the conceptual frame of lipotoxicity and type 2 diabetes.

Challenges and Rewards in Medicinal Chemistry Targeting Cardiovascular and Metabolic Diseases.

Medicinal chemistry has been transformed by major technological and conceptual innovations over the last three decades: structural biology and bioinformatics, structure and property based molecular design, the concepts of multidimensional optimization (MDO), in silico and experimental high-throughput molecular property analysis. The novel technologies advanced gradually and in synergy with biology and Roche has been at the forefront. Applications in drug discovery programs towards new medicines in cardiovascular and metabolic diseases are highlighted to show impact and advancement: the early discovery of endothelin antagonists for endothelial dysfunction (Bosentan), 11-beta hydroxysteroid dehydrogenase (11β-HSD1) inhibitors for dysregulated cellular glucocorticoid tonus (type 2 diabetes and metabolic syndrome) and non-covalent hormone sensitive lipase (HSL) inhibitors to study the scope of direct inhibition of lipolysis in the conceptual frame of lipotoxicity and type 2 diabetes.

Novel Triazolopyrimidine-Derived Cannabinoid Receptor 2 Agonists as Potential Treatment for Inflammatory Kidney Diseases.

The cannabinoid receptor 2 (CB2) system is described to modulate various pathological conditions, including inflammation and fibrosis. A series of new heterocyclic small-molecule CB2 receptor agonists were identified from a high-throughput screen. Lead optimization gave access to novel, highly potent, and selective (over CB1) triazolopyrimidine derivatives.

3-Amido-3-aryl-piperidines: A Novel Class of Potent, Selective, and Orally Active GlyT1 Inhibitors.

3-Amido-3-aryl-piperidines were discovered as a novel structural class of GlyT1 inhibitors. The structure-activity relationship, which was developed, led to the identification of highly potent compounds exhibiting excellent selectivity against the GlyT2 isoform, drug-like properties, and in vivo activity after oral administration.

CONFECT: conformations from an expert collection of torsion patterns.

The generation of sets of low-energy conformations for a given molecule is a central task in drug design. Herein we present a new conformation generator called CONFECT that builds on our previously published library of torsion patterns. Conformations are generated as well as ranked by means of normalized frequency distributions derived from the Cambridge Structural Database (CSD).

Isothermal titration calorimetry with micelles: Thermodynamics of inhibitor binding to carnitine palmitoyltransferase 2 membrane protein.

Carnitine palmitoyl transferase 2 (CPT-2) is a key enzyme in the mitochondrial fatty acid metabolism. The active site is comprised of a Y-shaped tunnel with distinct binding sites for the substrate acylcarnitine and the cofactor CoA. We investigated the thermodynamics of binding of four inhibitors directed against either the CoA or the acylcarnitine binding sites using isothermal titration calorimetry (ITC).

Highly potent and selective cannabinoid receptor 2 agonists: initial hit optimization of an adamantyl hit series identified from high-through-put screening.

A series of highly potent & selective adamantane derived CB2 agonists was identified in a high-throughput screen. A SAR was established and physicochemical properties were significantly improved. This was accompanied by potency of the compounds on the Q63R variant and varying β-arrestin data which will support the insight into their relevance for the in vivo situation.

Torsion angle preferences in druglike chemical space: a comprehensive guide.

Crystal structure databases offer ample opportunities to derive small molecule conformation preferences, but the derived knowledge is not systematically applied in drug discovery research. We address this gap by a comprehensive and extendable expert system enabling quick assessment of the probability of a given conformation to occur. It is based on a hierarchical system of torsion patterns that cover a large part of druglike chemical space.

TFD: Torsion Fingerprints as a new measure to compare small molecule conformations.

Advantages like intuitive interpretation, objectivity, general applicability, and its easy, automated calculation make the rmsd (root-mean-squared deviation) the measure of choice for the investigation of the accuracy of conformational model generators. For comparing conformations of a single molecule this is a clearly superior method. Single molecule analysis is, however, a rare scenario.

NAOMI: on the almost trivial task of reading molecules from different file formats.

In most cheminformatics workflows, chemical information is stored in files which provide the necessary data for subsequent calculations. The correct interpretation of the file formats is an important prerequisite to obtain meaningful results. Consistent reading of molecules from files, however, is not an easy task.

De novo design by pharmacophore-based searches in fragment spaces.

De novo ligand design supports the search for novel molecular scaffolds in medicinal chemistry projects. This search can either be based on structural information of the targeted active site (structure-based approach) or on similarity to known binders (ligand-based approach). In the absence of structural information on the target, pharmacophores provide a way to find topologically novel scaffolds.

Structure-assisted discovery of an aminothiazole derivative as a lead molecule for inhibition of bacterial fatty-acid synthesis.

Fatty-acid synthesis in bacteria is of great interest as a target for the discovery of antibacterial compounds. The addition of a new acetyl moiety to the growing fatty-acid chain, an essential step in this process, is catalyzed by beta-ketoacyl-ACP synthase (KAS). It is inhibited by natural antibiotics such as cerulenin and thiolactomycin; however, these lack the requirements for optimal drug development.

Recore: a fast and versatile method for scaffold hopping based on small molecule crystal structure conformations.

Replacing central elements of known active structures is a common procedure to enter new compound classes. Different computational methods have already been developed to help with this task, varying in the description of possible replacements, the query input, and the similarity measure used. In this paper, a novel approach for scaffold replacement and a corresponding software tool, called Recore, is introduced.

Access of the substrate to the active site of yeast oxidosqualene cyclase: an inhibition and site-directed mutagenesis approach.

A structural model of Saccharomyces cerevisiae oxidosqualene cyclase (SceOSC) suggests that some residues of the conserved sequence Pro-Ala-Glu-Val-Phe-Gly (residues 524-529) belong to a channel constriction that gives access to the active-site cavity. Starting from the SceOSC C457D mutant, which lacks the cysteine residue next to the catalytic Asp456 residue Cys457 has been replaced but Asp456 is still there, we prepared two further mutants where the wild-type residues Ala525 and Glu526 were individually replaced by cysteine. These mutants, especially E526C, were very sensitive to the thiol-reacting agent dodecyl-maleimide.

Enzyme redesign: two mutations cooperate to convert cycloartenol synthase into an accurate lanosterol synthase.

Efforts to modify the catalytic specificity of enzymes consistently show that it is easier to broaden the substrate or product specificity of an accurate enzyme than to restrict the selectivity of one that is promiscuous. Described herein are experiments in which cycloartenol synthase was redesigned to become a highly accurate lanosterol synthase. Several single mutants have been described that modify the catalytic specificity of cycloartenol to form some lanosterol.

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.

Scoring functions for protein-ligand interactions: a critical perspective.

Scoring functions play an essential role in structure-based virtual screening. They are required to guide the docking of candidate compounds to structures of receptor binding sites, to select probable binding modes, and to discriminate binders from non-binders. Although many scoring functions have successfully been used to identify novel ligands for a wide variety of targets, much work remains to be done to avoid incorrect prediction of binding modes and high numbers of false positives.

Ensemble methods for classification in cheminformatics.

We describe the application of ensemble methods to binary classification problems on two pharmaceutical compound data sets. Several variants of single and ensembles models of k-nearest neighbors classifiers, support vector machines (SVMs), and single ridge regression models are compared. All methods exhibit robust classification even when more features are given than observations.

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.

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.