IADE: a system for intelligent automatic design of bioisosteric analogs

Molecular mechanisms of caramel-like odorant-olfactory receptor interactions based on a computational chemistry approach

Molecular mechanisms of caramel-like odorant-olfactory receptor interactions were investigated based on molecular docking and molecular dynamics simulations. The transmembrane regions TM-3, TM-5 and TM-6 of receptors were main contributors of amino acid residues in the docking. Molecular docking results showed that hydrogen bonding and pi-pi stacking were the key forces for the stabilization of caramel-like odorants. The binding energies were positively correlated with the molecular weight of caramel-like odorants. Residues Asn155 (84%, OR2W1), Asn206 (86%, OR8D1), Ser155 (77%, OR8D1), Asp179 (87%, OR5M3), Val182 (84%, OR2J2) and Tyr260 (94%, OR2J2) with high frequencies played an important role in the complexes formation. Odorants 4-hydroxy-5-methylfuran-3(2H)-one (16#) and methylglyoxal (128#) were screened by molecular field-based similarity analysis, which tended to bind to the receptors OR1G1 and OR52H1 respectively, resulting a caramel-like aroma perception. The obtained results are useful for better understanding the perception of caramel-like odorants and their high-throughput screening.

Individual and collective human intelligence in drug design: evaluating the search strategy

In recent years, individual and collective human intelligence, defined as the knowledge, skills, reasoning and intuition of individuals and groups, have been used in combination with computer algorithms to solve complex scientific problems. Such approach was successfully used in different research fields such as: structural biology, comparative genomics, macromolecular crystallography and RNA design. Herein we describe an attempt to use a similar approach in small-molecule drug discovery, specifically to drive search strategies of de novo drug design. This is assessed with a case study that consists of a series of public experiments in which participants had to explore the huge chemical space in silico to find predefined compounds by designing molecules and analyzing the score associate with them. Such a process may be seen as an instantaneous surrogate of the classical design-make-test cycles carried out by medicinal chemists during the drug discovery hit to lead phase but not hindered by long synthesis and testing times. We present first findings on (1) assessing human intelligence in chemical space exploration, (2) comparing individual and collective human intelligence performance in this task and (3) contrasting some human and artificial intelligence achievements in de novo drug design.

MolOpt: A Web Server for Drug Design using Bioisosteric Transformation

BACKGROUND

Bioisosteric replacement is widely used in drug design for lead optimization. However, the identification of a suitable bioisosteric group is not an easy task.

METHODS

In this work, we present MolOpt, a web server for in silico drug design using bioisosteric transformation. Potential bioisosteric transformation rules were derived from data mining, deep generative machine learning and similarity comparison. MolOpt tries to assist the medicinal chemist in his/her search for what to make next.

RESULTS AND DISCUSSION

By replacing molecular substructures with similar chemical groups, MolOpt automatically generates lists of analogues. MolOpt also evaluates forty important pharmacokinetic and toxic properties for each newly designed molecule. The transformed analogues can be assessed for possible future study.

CONCLUSION

MolOpt is useful for the identification of suitable lead optimization ideas. The MolOpt Server is freely available for use on the web at http://xundrug.cn/molopt.

Molecular modeling approaches to address drug-metabolizing enzymes (DMEs) mediated chemoresistance: a review

Resistance against clinically approved anticancer drugs is the main roadblock in cancer treatment. Drug metabolizing enzymes (DMEs) that are capable of metabolizing a variety of xenobiotic get overexpressed in malignant cells, therefore, catalyzing drug inactivation. As evident from the literature reports, the levels of DMEs increase in cancer cells that ultimately lead to drug inactivation followed by drug resistance. To puzzle out this issue, several strategies inclusive of analog designing, prodrug designing, and inhibitor designing have been forged. On that front, the implementation of computational tools can be considered a fascinating approach to address the problem of chemoresistance. Various research groups have adopted different molecular modeling tools for the investigation of DMEs mediated toxicity problems. However, the utilization of these in-silico tools in maneuvering the DME mediated chemoresistance is least considered and yet to be explored. These tools can be employed in the designing of such chemotherapeutic agents that are devoid of the resistance problem. The current review canvasses various molecular modeling approaches that can be implemented to address this issue. Special focus was laid on the development of specific inhibitors of DMEs. Additionally, the strategies to bypass the DMEs mediated drug metabolism were also contemplated in this report that includes analogs and pro-drugs designing. Different strategies discussed in the review will be beneficial in designing novel chemotherapeutic agents that depreciate the resistance problem.

The Derivation of a Matched Molecular Pairs Based ADME/Tox Knowledge Base for Compound Optimization

Matched Molecular Pairs (MMP) analysis is a well-established technique for Structure Activity and Property Analysis (SAR and SPR). Summarizing multiple MMPs that describe the same structural change into a single chemical transform can be a powerful tool for prediction (termed Transform from here on). This is particularly useful in the area of Absorption, Distribution, Metabolism, and Elimination (ADME) analysis that is less influenced by 3D structural binding effects. The creation of a knowledge database containing many of these Transforms across typical ADME assays promises to be a powerful approach to aid multidimensional optimization. We present a detailed workflow for the derivation of such a database. We include details of an MMP fragmentation algorithm with associated statistical summarization methods for the derivation of Transforms. This is made freely available as part of the LillyMol software package. We describe the application of this method to several ADME/Tox (Toxicity) assay data sets and highlight multiple cases where the impact of traditional medicinal chemistry Transforms is contradicted by MMP data. We also describe the internal software interface used by medicinal chemists to aid the design of new compounds via automated suggestion. This approach utilizes the matched pairs database to "suggest" improved compounds in an automated design scenario. A nonvisual script-based version of the automated suggestions code with an associated set of described chemical Transforms is also made freely available along with this paper and as part of the LillyMol software package. Finally, we contrast this knowledge database against a larger database of all MMPs derived from a 2 million compound diversity set and a subset of MMPs seen in historical discovery projects. The comparison against all transforms in the diversity collection highlights the very low coverage of the transform database as compared to all possible transforms involving 15 atom fragments. The comparison against a smaller subset of Transforms seen on internal Medicinal Chemistry projects shows better coverage of the transform database for a small set of common medicinal chemistry strategies. Within the context of all possible transforms available to a medicinal chemistry project team, the challenge remains to move beyond mere idea generation from past projects toward high quality prediction for novel ADME/Tox modulating Transforms.

The In Silico Drug Discovery Toolbox: Applications in Lead Discovery and Optimization

BACKGROUND

Discovery and development of a new drug is a long lasting and expensive journey that takes around 20 years from starting idea to approval and marketing of new medication. Despite R&D expenditures have been constantly increasing in the last few years, the number of new drugs introduced into market has been steadily declining. This is mainly due to preclinical and clinical safety issues, which still represent about 40% of drug discontinuation. To cope with this issue, a number of in silico techniques are currently being used for an early stage evaluation/prediction of potential safety issues, allowing to increase the drug-discovery success rate and reduce costs associated with the development of a new drug.

METHODS

In the present review, we will analyse the early steps of the drug-discovery pipeline, describing the sequence of steps from disease selection to lead optimization and focusing on the most common in silico tools used to assess attrition risks and build a mitigation plan.

RESULTS

A comprehensive list of widely used in silico tools, databases, and public initiatives that can be effectively implemented and used in the drug discovery pipeline has been provided. A few examples of how these tools can be problem-solving and how they may increase the success rate of a drug discovery and development program have been also provided. Finally, selected examples where the application of in silico tools had effectively contributed to the development of marketed drugs or clinical candidates will be given.

CONCLUSION

The in silico toolbox finds great application in every step of early drug discovery: (i) target identification and validation; (ii) hit identification; (iii) hit-to-lead; and (iv) lead optimization. Each of these steps has been described in details, providing a useful overview on the role played by in silico tools in the decision-making process to speed-up the discovery of new drugs.

GuacaMol: Benchmarking Models for de Novo Molecular Design

De novo design seeks to generate molecules with required property profiles by virtual design-make-test cycles. With the emergence of deep learning and neural generative models in many application areas, models for molecular design based on neural networks appeared recently and show promising results. However, the new models have not been profiled on consistent tasks, and comparative studies to well-established algorithms have only seldom been performed. To standardize the assessment of both classical and neural models for de novo molecular design, we propose an evaluation framework, GuacaMol, based on a suite of standardized benchmarks. The benchmark tasks encompass measuring the fidelity of the models to reproduce the property distribution of the training sets, the ability to generate novel molecules, the exploration and exploitation of chemical space, and a variety of single and multiobjective optimization tasks. The benchmarking open-source Python code and a leaderboard can be found on https://benevolent.ai/guacamol .

Evaluation of a semi-automated workflow for fragment growing

Intelligent Automatic Design (IADE) is an expert system developed at Novartis to identify nonclassical bioisosteres. In addition to bioisostere searching, one could also use IADE to grow a fragment bound to a protein. Here we report an evaluation of IADE as a tool for fragment growing. Three examples from the literature served as test cases. In all three cases, IADE generated close analogues of the published compounds and reproduced their crystallographic binding modes. This exercise validated the use of the IADE system for fragment growing. We have also gained experience in optimizing the performance of IADE for this type of application.

In silico design of anti-atherogenic biomaterials

Atherogenesis, the uncontrolled deposition of modified lipoproteins in inflamed arteries, serves as a focal trigger of cardiovascular disease (CVD). Polymeric biomaterials have been envisioned to counteract atherogenesis based on their ability to repress scavenger mediated uptake of oxidized lipoprotein (oxLDL) in macrophages. Following the conceptualization in our laboratories of a new library of amphiphilic macromolecules (AMs), assembled from sugar backbones, aliphatic chains and poly(ethylene glycol) tails, a more rational approach is necessary to parse the diverse features such as charge, hydrophobicity, sugar composition and stereochemistry. In this study, we advance a computational biomaterials design approach to screen and elucidate anti-atherogenic biomaterials with high efficacy. AMs were quantified in terms of not only 1D (molecular formula) and 2D (molecular connectivity) descriptors, but also new 3D (molecular geometry) descriptors of AMs modeled by coarse-grained molecular dynamics (MD) followed by all-atom MD simulations. Quantitative structure-activity relationship (QSAR) models for anti-atherogenic activity were then constructed by screening a total of 1164 descriptors against the corresponding, experimentally measured potency of AM inhibition of oxLDL uptake in human monocyte-derived macrophages. Five key descriptors were identified to provide a strong linear correlation between the predicted and observed anti-atherogenic activity values, and were then used to correctly forecast the efficacy of three newly designed AMs. Thus, a new ligand-based drug design framework was successfully adapted to computationally screen and design biomaterials with cardiovascular therapeutic properties.