Prediction of Small-Molecule Developability Using Large-Scale In Silico ADMET Models

Early in silico assessment of the potential of a series of compounds to deliver a drug is one of the major challenges in computer-assisted drug design. The goal is to identify the right chemical series of compounds out of a large chemical space to then subsequently prioritize the molecules with the highest potential to become a drug. Although multiple approaches to assess compounds have been developed over decades, the quality of these predictors is often not good enough and compounds that agree with the respective estimates are not necessarily druglike. Here, we report a novel deep learning approach that leverages large-scale predictions of ∼100 ADMET assays to assess the potential of a compound to become a relevant drug candidate. The resulting score, which we termed bPK score, substantially outperforms previous approaches and showed strong discriminative performance on data sets where previous approaches did not.

PREFER: A New Predictive Modeling Framework for Molecular Discovery

Machine-learning and deep-learning models have been extensively used in cheminformatics to predict molecular properties, to reduce the need for direct measurements, and to accelerate compound prioritization. However, different setups and frameworks and the large number of molecular representations make it difficult to properly evaluate, reproduce, and compare them. Here we present a new PREdictive modeling FramEwoRk for molecular discovery (PREFER), written in Python (version 3.7.7) and based on AutoSklearn (version 0.14.7), that allows comparison between different molecular representations and common machine-learning models. We provide an overview of the design of our framework and show exemplary use cases and results of several representation-model combinations on diverse data sets, both public and in-house. Finally, we discuss the use of PREFER on small data sets. The code of the framework is freely available on GitHub.

Chemical reactivity prediction: current methods and different application areas

The ability to predict chemical reactivity of a molecule is highly desirable in drug discovery, both ex vivo (synthetic route planning, formulation, stability) and in vivo: metabolic reactions determine pharmacodynamics, pharmacokinetics and potential toxic effects, and early assessment of liabilities is vital to reduce attrition rates in later stages of development. Quantum mechanics offer a precise description of the interactions between electrons and orbitals in the breaking and forming of new bonds. Modern algorithms and faster computers have allowed the study of more complex systems in a punctual and accurate fashion, and answers for chemical questions around stability and reactivity can now be provided. Through machine learning, predictive models can be built out of descriptors derived from quantum mechanics and cheminformatics, even in the absence of experimental data to train on. In this article, current progress on computational reactivity prediction is reviewed: applications to problems in drug design, such as modelling of metabolism and covalent inhibition, are highlighted and unmet challenges are posed.

Discovery of 4-((2S,4S)-4-Ethoxy-1-((5-methoxy-7-methyl-1H-indol-4-yl)methyl)piperidin-2-yl)benzoic Acid (LNP023), a Factor B Inhibitor Specifically Designed To Be Applicable to Treating a Diverse Array of Complement Mediated Diseases

The alternative pathway (AP) of the complement system is a key contributor to the pathogenesis of several human diseases including age-related macular degeneration, paroxysmal nocturnal hemoglobinuria (PNH), atypical hemolytic uremic syndrome (aHUS), and various glomerular diseases. The serine protease factor B (FB) is a key node in the AP and is integral to the formation of C3 and C5 convertase. Despite the prominent role of FB in the AP, selective orally bioavailable inhibitors, beyond our own efforts, have not been reported previously. Herein we describe in more detail our efforts to identify FB inhibitors by high-throughput screening (HTS) and leveraging insights from several X-ray cocrystal structures during optimization efforts. This work culminated in the discovery of LNP023 (41), which is currently being evaluated clinically in several diverse AP mediated indications.

Medicinal Chemistry Aware Database GDBMedChem

The generated database GDB17 enumerates 166.4 billion possible molecules up to 17 atoms of C, N, O, S and halogens following simple chemical stability and synthetic feasibility rules, however medicinal chemistry criteria are not taken into account. Here we applied rules inspired by medicinal chemistry to exclude problematic functional groups and complex molecules from GDB17, and sampled the resulting subset uniformly across molecular size, stereochemistry and polarity to form GDBMedChem as a compact collection of 10 million small molecules. This collection has reduced complexity and better synthetic accessibility than the entire GDB17 but retains higher sp³ -carbon fraction and natural product likeness scores compared to known drugs. GDBMedChem molecules are more diverse and very different from known molecules in terms of substructures and represent an unprecedented source of diversity for drug design. GDBMedChem is available for 3D-visualization, similarity searching and for download at http://gdb.unibe.ch.

Structure-based design and synthesis of macrocyclic human rhinovirus 3C protease inhibitors

The design and synthesis of macrocyclic inhibitors of human rhinovirus 3C protease is described. A macrocyclic linkage of the P1 and P3 residues, and the subsequent structure-based optimization of the macrocycle conformation and size led to the identification of a potent biochemical inhibitor 10 with sub-micromolar antiviral activity.

Building Compound Archives for the Future

Will the targets of the future be covered by the compound libraries of today? This communication will cover a critical review of past strategies before turning to a new measure of diversity, protein pockets. A fingerprint descriptor for pockets will be described.

Discovery and structure-based optimization of adenain inhibitors

The cysteine protease adenain is the essential protease of adenovirus and, as such, represents a promising target for the treatment of ocular and other adenoviral infections. Through a concise two-pronged hit discovery approach we identified tetrapeptide nitrile 1 and pyrimidine nitrile 2 as complementary starting points for adenain inhibition. These hits enabled the first high-resolution X-ray cocrystal structures of adenain with inhibitors bound and revealed the binding mode of 1 and 2. The screening hits were optimized by a structure-guided medicinal chemistry strategy into low nanomolar drug-like inhibitors of adenain.

Discovery of C-(1-aryl-cyclohexyl)-methylamines as selective, orally available inhibitors of dipeptidyl peptidase IV

The successful launches of dipeptidyl peptidase IV (DPP IV) inhibitors as oral anti-diabetics warrant and spur the further quest for additional chemical entities in this promising class of therapeutics. Numerous pharmaceutical companies have pursued their proprietary candidates towards the clinic, resulting in a large body of published chemical structures associated with DPP IV. Herein, we report the discovery of a novel chemotype for DPP IV inhibition based on the C-(1-aryl-cyclohexyl)-methylamine scaffold and its optimization to compounds which selectively inhibit DPP IV at low-nM potency and exhibit an excellent oral pharmacokinetic profile in the rat.

Design, synthesis, and functional evaluation of leukocyte function associated antigen-1 antagonists in early and late stages of cancer development

The integrin leukocyte function associated antigen 1 (LFA-1) binds the intercellular adhesion molecule 1 (ICAM-1) by its α(L)-chain inserted domain (I-domain). This interaction plays a key role in cancer and other diseases. We report the structure-based design, small-scale synthesis, and biological activity evaluation of a novel family of LFA-1 antagonists. The design led to the synthesis of a family of highly substituted homochiral pyrrolidines with antiproliferative and antimetastatic activity in a murine model of colon carcinoma, as well as potent antiadhesive properties in several cancer cell lines in the low micromolar range. NMR analysis of their binding to the isolated I-domain shows that they bind to the I-domain allosteric site (IDAS), the binding site of other allosteric LFA-1 inhibitors. These results provide evidence of the potential therapeutic value of a new set of LFA-1 inhibitors, whose further development is facilitated by a synthetic strategy that is versatile and fully stereocontrolled.

Magnetic Effects of Disulfide Bridges: A Density Functional and Semiempirical Study

Density functional chemical shielding calculations are reported for methane and hydrogen disulfide dimers. The calculations show that the contributions of disulfide bridges to the chemical shielding of neighboring protons is sizable at distances that are frequently sampled in protein structures. A semiempirical model of the quantum chemical data is developed. It is shown that magnetic anisotropy effects of disulfide are poorly described by the McConnell equation, both qualitatively and quantitatively. In particular, the ratio of magnetic anisotropy contributions to shielding along and perpendicular to the magnetic anisotropy principal axis do not conform to the predictions of the McConnell equation, and magnetic anisotropy effects are not null along the magic angle axis. A sulfur-based model of the magnetic anisotropy of the disulfide is developed and shown to give much better agreement with the quantum chemical data.

Functionality maps of the ATP binding site of DNA gyrase B: generation of a consensus model of ligand binding

The Multiple Copy Simultaneous Search method (MCSS) was used to construct consensus functionality maps for functional group binding in the ATP binding site of DNA gyrase B. To account for the conformational flexibility of the protein active site, which involves small side chain fluctuations as well as large-scale loop motions, the calculations were done for three different conformations of the 24 kDa subdomain of DNA gyrase B. A postprocessing procedure that employs a continuum dielectric model to include solvent effects was used to calculate the binding free energy for every functional group. These results were ranked according to their affinity for DNA gyrase B and clustered using a new procedure based on van der Waals contacts that is better adapted for cases where multiple conformations are being considered. A total of 23 different functional groups were tested. The results gave consensus maps that indicate those functional group binding sites that are insensitive to the specific protein conformation. The maps also demonstrate that functional groups other than those found in the known ligands may bind competitively in the binding sites of known ligands. This suggests numerous scaffolds that can be used in the development of new ligands for the ATP and coumarinic binding sites in DNA gyrase B. Finally, the calculations show the existence of alternative binding sites near the known binding sites that could be targeted in the rational design for new inhibitors.

Structure activity relationship by NMR and by computer: a comparative study

There has recently been considerable interest in using NMR spectroscopy to identify ligand binding sites of macromolecules. In particular, a modular approach has been put forward by Fesik et al. (Shuker, S. B.; Hajduk, P. J.; Meadows, R. P.; Fesik, S. W. Science 1996, 274, 1531-1534) in which small ligands that bind to a particular target are identified in a first round of screening and subsequently linked together to form ligands of higher affinity. Similar strategies have also been proposed for in silico drug design, where the binding sites of small chemical groups are identified, and complete ligands are subsequently assembled from different groups that have favorable interactions with the macromolecular target. In this paper, we compare experimental and computational results on a selected target (FKBP12). The binding sites of three small ligands ((2S)1-acetylprolinemethylester, 1-formylpiperidine, 1-piperidinecarboxamide) in FKBP12 were identified independently by NMR and by computational methods. The subsequent comparison of the experimental and computational data showed that the computational method identified and ranked favorably ligand positions that satisfy the experimental NOE constraints.

Structural and dynamic characterization of omega-conotoxin MVIIA: the binding loop exhibits slow conformational exchange

omega-Conotoxin MVIIA is a 25-residue, disulfide-bridged polypeptide from the venom of the sea snail Conus magus that binds to neuronal N-type calcium channels. It forms a compact folded structure, presenting a loop between Cys8 and Cys15 that contains a set of residues critical for its binding. The loop does not have a unique defined structure, nor is it intrinsically flexible. Broadening of a subset of resonances in the NMR spectrum at low temperature, anomalous temperature dependence of the chemical shifts of some resonances, and exchange contributions to J(0) from (13)C relaxation measurements reveal that conformational exchange affects the residues in this loop. The effects of this exchange on the calculated structure of omega-conotoxin MVIIA are discussed. The exchange appears to be associated with a change in the conformation of the disulfide bridge Cys8-Cys20. The implications for the use of the omega-conotoxins as a scaffold for carrying other functions is discussed.