Trends in small molecule drug properties: A developability molecule assessment perspective

Developability molecule assessment is a key interfacial capability across the biopharmaceutical industry, screening and staging molecules discovered by medicinal chemists for successful chemistry manufacturing controls (CMC) development and launch. The breadth of responsibility and expertise such teams possess puts them in a unique position to understand the impact of the physicochemical properties of a drug during its initial discovery and subsequent development. However, most of the publications describing trends in physicochemical properties are written from a medicinal chemistry perspective with the aim to identify molecules with better ADMET profiles that are either lead-like or drug-like, failing to describe the impact these properties have on CMC development. To systematically uncover knowledge obtained from recent trends in physicochemical properties and the corresponding impact on CMC development, a comprehensive analysis was conducted on molecules in the drug repurposing hub dataset. The only physicochemical property that seems to have been preserved in FDA-approved oral molecules over the decades (1900-2020) is a constant H-bond donor count, highlighting the importance this property has on cell permeability and lattice energy. Pharmaceutical attrition analysis suggests that partition-distribution coefficient, H-bond acceptors, polar surface area and the fraction of sp³ carbons are properties that are associated with compound attrition. Looking at pharmaceutical attrition asynchronously with the temporal analysis of FDA-approved oral molecules highlights the opposing trends, risks and diminishing effects some of these physiochemical properties (cLogP, cLogD and Fsp3) have on describing compound attrition during the past decade. Trellising the dataset by target class suggests that certain formulation and drug delivery strategies can be anticipated or put into place based on target class of a molecule. For example, molecules binding to nuclear hormone receptors are amenable to lipid-based drug delivery systems with proven commercial success. Although the poor solubility of kinase inhibitors is a combination of hydrophobicity (due to aromaticity) required to bind to its target and high lattice energy (melting point), they are a challenging target class to formulate. The influence of drug targets on physicochemical properties and the temporal nature of these properties is highlighted when comparing molecules in the drug repurposing dataset to those developed at Amgen. An improved understanding of the impact of molecular properties on performance attributes can accelerate decisions and facilitate risk assessments during candidate selection and development.

In Silico Resources of Drug-Likeness as a Mirror: What Are We Lacking in Pesticide-Likeness?

Unfavorable bioavailability is an important aspect underlying the failure of drug candidates. Computational approaches for evaluating drug-likeness can minimize these risks. Over the past decades, computational approaches for evaluating drug-likeness have sped up the process of drug development and were also quickly derived to pesticide-likeness. As a result of many critical differences between drugs and pesticides, many kinds of methods for drug-likeness cannot be used for pesticide-likeness. Therefore, it is crucial to comprehensively compare and analyze the differences between drug-likeness and pesticide-likeness, which may provide a basis for solving the problems encountered during the evaluation of pesticide-likeness. Here, we systematically collected the recent advances of drug-likeness and pesticide-likeness and compared their characteristics. We also evaluated the current lack of studies on pesticide-likeness, the molecular descriptors and parameters adopted, the pesticide-likeness model on pesticide target organisms, and comprehensive analysis tools. This work may guide researchers to use appropriate methods for developing pesticide-likeness models. It may also aid non-specialists to understand some important concepts in drug-likeness and pesticide-likeness.

Physicochemical Profiling and Comparison of Research Antiplasmodials and Advanced Stage Antimalarials with Oral Drugs

To understand the property space of antimalarials, we collated a large dataset of research antiplasmodial (RAP) molecules with known in vitro potencies and advanced stage antimalarials (ASAMs) with established oral bioavailability. While RAP molecules are "non-druglike", ASAM molecules display properties closer to Lipinski's and Veber's thresholds. Comparison within the different potency groups of RAP molecules indicates that the in vitro potency is positively correlated to the molecular weight, the calculated octanol-water partition coefficient (clog P), aromatic ring counts (#Ar), and hydrogen bond acceptors. Despite both categories being bioavailable, the ASAM molecules are relatively larger and more lipophilic, have a lower polar surface area, and possess a higher count of heteroaromatic rings than oral drugs. Also, antimalarials are found to have a higher proportion of aromatic (#ArN) and basic nitrogen (#BaN) counts, features implicitly used in the design of antimalarial molecules but not well studied hitherto. We also propose using descriptors scaled by the sum of #ArN and #BaN (SBAN) to define an antimalarial property space. Together, these results may have important applications in the identification and optimization of future antimalarials.

meta-Non-flat substituents: a novel molecular design to improve aqueous solubility in small molecule drug discovery

Aqueous solubility is a key requirement for small-molecule drug candidates. Here, we investigated the regioisomer-physicochemical property relationships of disubstituted benzenes. We found that meta-isomers bearing non-flat substituents tend to possess the lowest melting point and the highest thermodynamic aqueous solubility among the regioisomers. The examination of pharmaceutical compounds containing a disubstituted benzene moiety supported the idea that the introduction of a non-flat substituent at the meta position of a benzene substructure would be a promising approach for medicinal chemists aiming to improve the thermodynamic aqueous solubility of drug candidates, even though it might not be universally effective.

Facts, Patterns, and Principles in Drug Discovery: Appraising the Rule of 5 with Measured Physicochemical Data

The rule of 5 was designed to estimate the likelihood of poor absorption or permeation, noting the impact of poor solubility. This Perspective explores the impact of various physicochemical descriptors and contemporary lipophilicity measurements on permeability and solubility, showing that the distribution coefficient log D_7.4 (rather than log P) is the most impactful parameter. Molecular weight, almost invariably the defining characteristic of "beyond the rule of 5" compounds, has little impact on solubility when log D_7.4 measurements and aromaticity are considered. Predicting permeation is more complex, given passive and carrier transport mechanisms; however, notable patterns of behavior are apparent, giving insight even "beyond the rule of 5". Recommended best practices should involve using the facts (measurements) and the patterns they reveal to establish informative principles rather than fastidious rules.

Squaramides as novel class I and IIB histone deacetylase inhibitors for topical treatment of cutaneous t-cell lymphoma

A series of squaramide-based hydroxamic acids were designed, synthesized and evaluated against human HDAC enzyme. Squaramides were found to be potent in the Hut78 cell line, but initially suffered from low solubility. Leads with improved solubility and metabolic profiles were shown to be class I, IIB and IV selective.

Matched Molecular Pair Analysis on Large Melting Point Datasets: A Big Data Perspective

A matched molecular pair (MMP) analysis was used to examine the change in melting point (MP) between pairs of similar molecules in a set of ∼275k compounds. We found many cases in which the change in MP (ΔMP) of compounds correlates with changes in functional groups. In line with the results of a previous study, correlations between ΔMP and simple molecular descriptors, such as the number of hydrogen bond donors, were identified. In using a larger dataset, covering a wider chemical space and range of melting points, we observed that this method remains stable and scales well with larger datasets. This MMP-based method could find use as a simple privacy-preserving technique to analyze large proprietary databases and share findings between participating research groups.

Novel bisnaphthalimidopropyl (BNIPs) derivatives as anticancer compounds targeting DNA in human breast cancer cells

Bisnaphthalimidopropyl (BNIP) derivatives are a family of compounds that exert anti-cancer activities in vitro and, according to previous studies, variations in the linker sequence have increased their DNA binding and cytotoxic activities. By modifying the linker sequence of bisnaphthalimidopropyl diaminodicyclohexylmethane (BNIPDaCHM), a previously synthesised BNIP derivative with anti-cancer properties, three novel BNIP derivatives were designed. Bisnaphthalimidopropyl-piperidylpropane (BNIPPiProp), a structural isomer of BNIPDaCHM, bisnaphthalimidopropyl ethylenedipiperidine dihydrobromide (BNIPPiEth), an isoform of BNIPDaCHM with a shorter linker chain, and (trans(trans))-bisnaphthalimidopropyl diaminodicyclohexylmethane (trans,trans-BNIPDaCHM), a stereoisomer of BNIPDaCHM, were successfully synthesised (72.3-29.5% yield) and characterised by nuclear magnetic resonance spectroscopy (NMR) and mass spectrometry (MS). Competitive displacement of ethidium bromide (EtBr) and UV binding studies were used to study the interactions of BNIP derivatives with Calf Thymus DNA. The cytotoxicity of these derivatives was assessed against human breast cancer MDA-MB-231 and SKBR-3 cells by MTT assay. Propidium iodide (PI) flow cytometry was conducted in order to evaluate the cellular DNA content in both breast cancer cell lines before and after treatment with BNIPs. The results showed that all novel BNIPs exhibit strong DNA binding properties in vitro, and strong cytotoxicity, with IC₅₀ values in the range of 0.2-3.3 μM after 24 hours drug treatment. Two of the novel BNIP derivatives, BNIPPiEth and trans,trans-BNIPDaCHM, exhibited greater cytotoxicity against the two breast cancer cell lines studied, compared to BNIPDaCHM. By synthesising enantiopures and reducing the length of the linker sequence, the cytotoxicity of the BNIP derivatives was significantly improved compared to BNIPDaCHM, while maintaining DNA binding and bis-intercalating properties. In addition, cell cycle studies indicated that trans,trans-BNIPDaCHM, the most cytotoxic BNIP derivative, induced sub-G1 cell cycle arrest, indicative of apoptotic cell death. Based on these findings, further investigation is under way to assess the potential efficacy of trans,trans-BNIPDaCHM and BNIPPiEth in treating human breast cancer.

Discovery and process development of a novel TACE inhibitor for the topical treatment of psoriasis

Targeting the TNFα pathway is a validated approach to the treatment of psoriasis. In this pathway, TACE stands out as a druggable target and has been the focus of in-house research programs. In this article, we present the discovery of clinical candidate 26a. Starting from hits plagued with poor solubility or genotoxicity, 26a was identified through thorough multiparameter optimisation. Showing robust in vivo activity in an oxazolone-mediated inflammation model, the compound was selected for development. Following a polymorph screen, the hydrochloride salt was selected and the synthesis was efficiently developed to yield the API in 47% overall yield.

Insights into the crystallisation process from anhydrous, hydrated and solvated crystal forms of diatrizoic acid

Diatrizoic acid (DTA), a clinically used X-ray contrast agent, crystallises in two hydrated, three anhydrous and nine solvated solid forms, all of which have been characterised by X-ray crystallography. Single-crystal neutron structures of DTA dihydrate and monosodium DTA tetrahydrate have been determined. All of the solid-state structures have been analysed using partial atomic charges and hardness algorithm (PACHA) calculations. Even though in general all DTA crystal forms reveal similar intermolecular interactions, the overall crystal packing differs considerably from form to form. The water of the dihydrate is encapsulated between a pair of host molecules, which calculations reveal to be an extraordinarily stable motif. DTA presents functionalities that enable hydrogen and halogen bonding, and whilst an extended hydrogen-bonding network is realised in all crystal forms, halogen bonding is not present in the hydrated crystal forms. This is due to the formation of a hydrogen-bonding network based on individual enclosed water squares, which is not amenable to the concomitant formation of halogen bonds. The main interaction in the solvates involves the carboxylic acid, which corroborates the hypothesis that this strong interaction is the last one to be broken during the crystal desolvation and nucleation process.