Past, Present, and Future Perspectives on Computer-Aided Drug Design Methodologies

Recent advances from computer-aided drug design to artificial intelligence drug design

Computer-aided drug design (CADD), a cornerstone of modern drug discovery, can predict how a molecular structure relates to its activity and interacts with its target using structure-based and ligand-based methods. Fueled by ever-increasing data availability and continuous model optimization, artificial intelligence drug design (AIDD), as an enhanced iteration of CADD, has thrived in the past decade. AIDD demonstrates unprecedented opportunities in protein folding, property prediction, and molecular generation. It can also facilitate target identification, high-throughput screening (HTS), and synthetic route prediction. With AIDD involved, the process of drug discovery is greatly accelerated. Notably, AIDD offers the potential to explore uncharted territories of chemical space beyond current knowledge. In this perspective, we began by briefly outlining the main workflows and components of CADD. Then through showcasing exemplary cases driven by AIDD in recent years, we describe the evolving role of artificial intelligence (AI) in drug discovery from three distinct stages, that is, chemical library screening, linker generation, and de novo molecular generation. In this process, we attempted to draw comparisons between the features of CADD and AIDD.

A computational workflow to determine drug candidates alternative to aminoglycosides targeting the decoding center of E. coli ribosome

The global antibiotic resistance problem necessitates fast and effective approaches to finding novel inhibitors to treat bacterial infections. In this study, we propose a computational workflow to identify plausible high-affinity compounds from FDA-approved, investigational, and experimental libraries for the decoding center on the small subunit 30S of the E. coli ribosome. The workflow basically consists of two molecular docking calculations on the intact 30S, followed by molecular dynamics (MD) simulations coupled with MM-GBSA calculations on a truncated ribosome structure. The parameters used in the molecular docking suits, Glide and AutoDock Vina, as well as in the MD simulations with Desmond were carefully adjusted to obtain expected interactions for the ligand-rRNA complexes. A filtering procedure was followed, considering a fingerprint based on aminoglycoside's binding site on the 30S to obtain seven hit compounds either with different clinical usages or aminoglycoside derivatives under investigation, suggested for in vitro studies. The detailed workflow developed in this study promises an effective and fast approach for the estimation of binding free energies of large protein-RNA and ligand complexes.

Machine Learning Guided AQFEP: A Fast and Efficient Absolute Free Energy Perturbation Solution for Virtual Screening

Structure-based methods in drug discovery have become an integral part of the modern drug discovery process. The power of virtual screening lies in its ability to rapidly and cost-effectively explore enormous chemical spaces to select promising ligands for further experimental investigation. Relative free energy perturbation (RFEP) and similar methods are the gold standard for binding affinity prediction in drug discovery hit-to-lead and lead optimization phases, but have high computational cost and the requirement of a structural analog with a known activity. Without a reference molecule requirement, absolute FEP (AFEP) has, in theory, better accuracy for hit ID, but in practice, the slow throughput is not compatible with VS, where fast docking and unreliable scoring functions are still the standard. Here, we present an integrated workflow to virtually screen large and diverse chemical libraries efficiently, combining active learning with a physics-based scoring function based on a fast absolute free energy perturbation method. We validated the performance of the approach in the ranking of structurally related ligands, virtual screening hit rate enrichment, and active learning chemical space exploration; disclosing the largest reported collection of free energy simulations to date.

Site Identification and Next Choice Protocol for Hit-to-Lead Optimization

Time efficiency and cost savings are major challenges in drug discovery and development. In this process, the hit-to-lead stage is expected to improve efficiency because it primarily exploits the trial-and-error approach of medicinal chemists. This study proposes a site identification and next choice (SINCHO) protocol to improve the hit-to-lead efficiency. This protocol selects an anchor atom and growth site pair, which is desirable for a hit-to-lead strategy starting from a 3D complex structure. We developed and fine-tuned the protocol using a training data set and assessed it using a test data set of the preceding hit-to-lead strategy. The protocol was tested for experimentally determined structures and molecular dynamics (MD) ensembles. The protocol had a high prediction accuracy for applying MD ensembles, owing to the consideration of protein flexibility. The SINCHO protocol enables medicinal chemists to visualize and modify functional groups in a hit-to-lead manner.

Another string to your bow: machine learning prediction of the pharmacokinetic properties of small molecules

INTRODUCTION

Prediction of pharmacokinetic (PK) properties is crucial for drug discovery and development. Machine-learning (ML) models, which use statistical pattern recognition to learn correlations between input features (such as chemical structures) and target variables (such as PK parameters), are being increasingly used for this purpose. To embed ML models for PK prediction into workflows and to guide future development, a solid understanding of their applicability, advantages, limitations, and synergies with other approaches is necessary.

AREAS COVERED

This narrative review discusses the design and application of ML models to predict PK parameters of small molecules, especially in light of established approaches including in vitro-in vivo extrapolation (IVIVE) and physiologically based pharmacokinetic (PBPK) models. The authors illustrate scenarios in which the three approaches are used and emphasize how they enhance and complement each other. In particular, they highlight achievements, the state of the art and potentials of applying machine learning for PK prediction through a comphrehensive literature review.

EXPERT OPINION

ML models, when carefully crafted, regularly updated, and appropriately used, empower users to prioritize molecules with favorable PK properties. Informed practitioners can leverage these models to improve the efficiency of drug discovery and development process.

Practical Three-Component Regioselective Synthesis of Drug-Like 3-Aryl(or heteroaryl)-5,6-dihydrobenzo[h]cinnolines as Potential Non-Covalent Multi-Targeting Inhibitors To Combat Neurodegenerative Diseases

Neurodegenerative diseases (NDs) are one of the prominent health challenges facing contemporary society, and many efforts have been made to overcome and (or) control it. In this research paper, we described a practical one-pot two-step three-component reaction between 3,4-dihydronaphthalen-1(2H)-one (1), aryl(or heteroaryl)glyoxal monohydrates (2a-h), and hydrazine monohydrate (NH2NH2•H2O) for the regioselective preparation of some 3-aryl(or heteroaryl)-5,6-dihydrobenzo[h]cinnoline derivatives (3a-h). After synthesis and characterization of the mentioned cinnolines (3a-h), the in silico multi-targeting inhibitory properties of these heterocyclic scaffolds have been investigated upon various Homo sapiens-type enzymes, including hMAO-A, hMAO-B, hAChE, hBChE, hBACE-1, hBACE-2, hNQO-1, hNQO-2, hnNOS, hiNOS, hPARP-1, hPARP-2, hLRRK-2(G2019S), hGSK-3β, hp38α MAPK, hJNK-3, hOGA, hNMDA receptor, hnSMase-2, hIDO-1, hCOMT, hLIMK-1, hLIMK-2, hRIPK-1, hUCH-L1, hPARK-7, and hDHODH, which have confirmed their functions and roles in the neurodegenerative diseases (NDs), based on molecular docking studies, and the obtained results were compared with a wide range of approved drugs and well-known (with IC50, EC50, etc.) compounds. In addition, in silico ADMET prediction analysis was performed to examine the prospective drug properties of the synthesized heterocyclic compounds (3a-h). The obtained results from the molecular docking studies and ADMET-related data demonstrated that these series of 3-aryl(or heteroaryl)-5,6-dihydrobenzo[h]cinnolines (3a-h), especially hit ones, can really be turned into the potent core of new drugs for the treatment of neurodegenerative diseases (NDs), and/or due to the having some reactionable locations, they are able to have further organic reactions (such as cross-coupling reactions), and expansion of these compounds (for example, with using other types of aryl(or heteroaryl)glyoxal monohydrates) makes a new avenue for designing novel and efficient drugs for this purpose.

In silico drug design strategies for discovering novel tuberculosis therapeutics

INTRODUCTION

Tuberculosis remains a significant concern in global public health due to its intricate biology and propensity for developing antibiotic resistance. Discovering new drugs is a protracted and expensive endeavor, often spanning over a decade and incurring costs in the billions. However, computer-aided drug design (CADD) has surfaced as a nimbler and more cost-effective alternative. CADD tools enable us to decipher the interactions between therapeutic targets and novel drugs, making them invaluable in the quest for new tuberculosis treatments.

AREAS COVERED

In this review, the authors explore recent advancements in tuberculosis drug discovery enabled by in silico tools. The main objectives of this review article are to highlight emerging drug candidates identified through in silico methods and to provide an update on the therapeutic targets associated with Mycobacterium tuberculosis.

EXPERT OPINION

These in silico methods have not only streamlined the drug discovery process but also opened up new horizons for finding novel drug candidates and repositioning existing ones. The continued advancements in these fields hold great promise for more efficient, ethical, and successful drug development in the future.

Ligand based-design of potential schistosomiasis inhibitors through QSAR, homology modeling, molecular dynamics, pharmacokinetics, and DFT studies

Objectives

Schistosomiasis, a neglected tropical disease, is a leading cause of mortality in affected geographic areas. Currently, because no vaccine for schistosomiasis is available, control measures rely on widespread administration of the drug praziquantel (PZQ). The mass administration of PZQ has prompted concerns regarding the emergence of drug resistance. Therefore, new therapeutic targets and potential compounds are necessary to combat schistosomiasis.

Methods

Twenty-four potent derivatives of PZQ were optimized via density functional theory (DFT) at the B3LYP/6-31G∗ level. Quantitative structureactivity relationship (QSAR) models were generated and statistically validated, and a lead candidate was selected to develop therapeutic options with improved efficacy against schistosomiasis. The biological and binding energies of the designed compounds were evaluated. In addition, molecular dynamics; drug-likeness; absorption, distribution, metabolism, excretion, and toxicity (ADMET); and DFT studies were performed on the newly designed compounds.

Results

Five QSAR models were generated, among which model 1 had favorable validation parameters (R2train: 0.957, R2adj: 0.941, LOF: 0.101, Q2cv: 0.906, and R2test: 0.783) and was chosen to identify a lead candidate. Other statistical parameters for the chosen model included variance inflation factor values ranging from 1.242 to 1.678, and a Y-scrambling coefficient (cRp2) of 0.747. Five new compounds were designed with improved predicted activity (ranging from 5.081 to 7.022) surpassing those of both the lead compound and PZQ (predicted pEC50 of 5.545). Molecular dynamics simulation revealed high binding affinity of the proposed compounds toward the target receptor. ADMET and drug-likeness assessments indicated adherence to Lipinski's rule of five criteria, thereby suggesting pharmacological and oral safety. In addition, DFT analysis indicated resistance to electronic alteration during chemical reactions.

Conclusion

The proposed compounds exhibited potential drug characteristics, thus indicating their suitability for further investigation to enhance schistosomiasis treatment options.

Rapid discovery of a new antifoulant: From in silico studies targeting barnacle chitin synthase to efficacy against barnacle settlement

Due to the adverse environmental impacts of toxic heavy metal-based antifoulants, the screening of environmentally friendly antifoulants has become important for the development of marine antifouling technology. Compared with the traditional lengthy and costly screening method, computer-aided drug design (CADD) offers a promising and efficient solution that can accelerate the screening process of green antifoulants. In this study, we selected barnacle chitin synthase (CHS, an important enzyme for barnacle settlement and development) as the target protein for docking screening. Three CHS genes were identified in the barnacle Amphibalanus amphitrite, and their encoded proteins were found to share a conserved glycosyltransferase domain. Molecular docking of 31,561 marine natural products with AaCHSs revealed that zoanthamine alkaloids had the best binding affinity (-11.8 to -12.6 kcal/mol) to AaCHSs. Considering that the low abundance of zoanthamine alkaloids in marine organisms would limit their application as antifoulants, a marine fungal-derived natural product, mycoepoxydiene (MED), which has a similar chemical structure to zoanthamine alkaloids and the potential for large-scale production by fermentation, was selected and validated for stable binding to AaCHS2L2 using molecular docking and molecular dynamics simulations. Finally, the efficacy of MED in inhibiting cyprid settlement of A. amphitrite was confirmed by a bioassay that demonstrated an EC50 of 1.97 μg/mL, suggesting its potential as an antifoulant candidate. Our research confirmed the reliability of using AaCHSs as antifouling targets and has provided insights for the efficient discovery of green antifoulants by CADD.

A hierarchical approach towards identification of novel inhibitors against L, D-transpeptidase YcbB as an anti-bacterial therapeutic target

The bacterial cell wall, being a vital component for cell viability, is regarded as a promising drug target. The L, D-Transpeptidase YcbB enzyme has been implicated for a significant role in cell wall polymers cross linking during typhoid toxin release, β-lactam resistance and outer membrane defect rescue. These observations have been recorded in different bacterial pathogens such as Salmonella Typhimurium, Citrobacter rodentium, and Salmonella typhi. In this work, we have shown structure based virtual screening of diverse natural and synthetic drug libraries against the enzyme and revealed three compounds as LAS_32135590, LAS_34036730 and LAS-51380924. These compounds showed highly stable energies and the findings are very competitive with the control molecule ((1RG or (4 R,5S)-3-({(3S,5S)-5-[(3-carboxyphenyl)carbamoyl]pyrrolidin-3-yl}sulfanyl)-5-[(1S,2R)-1-formyl-2-hydroxypropyl]-4-methyl-4,5-dihydro-1H-pyrrole-2-carboxylic acid or ertapenem)) used. Compared to control (which has binding energy score of -11.63 kcal/mol), the compounds showed better binding energy. The binding energy score of LAS_32135590, LAS_34036730 and LAS-51380924 is -12.63 kcal/mol, -12.22 kcal/mol and -12.10 kcal/mol, respectively. Further, the docked snapshot of the lead compounds and control were investigated for stability under time dependent dynamics environment. All the three leads complex and control system showed significant equilibrium (mean RMSD < 3 Å) both in term of intermolecular docked conformation and binding interactions network. Further validation on the complex's stability was acquired from the end-state MMPB/GBSA analysis that observed greater contribution from van der Waals forces and electrostatic energy while less contribution was noticed from solvation part. The compounds were also showed good drug-likeness and are non-toxic and non-mutagenic. In short, the compounds can be used in experimental testing's and might be subjected to structure modification to get better results.Communicated by Ramaswamy H. Sarma.

Shaping the Future of Obesity Treatment: In Silico Multi-Modeling of IP6K1 Inhibitors for Obesity and Metabolic Dysfunction

Recent research has uncovered a promising approach to addressing the growing global health concern of obesity and related disorders. The inhibition of inositol hexakisphosphate kinase 1 (IP6K1) has emerged as a potential therapeutic strategy. This study employs multiple ligand-based in silico modeling techniques to investigate the structural requirements for benzisoxazole derivatives as IP6K1 inhibitors. Firstly, we developed linear 2D Quantitative Structure-Activity Relationship (2D-QSAR) models to ensure both their mechanistic interpretability and predictive accuracy. Then, ligand-based pharmacophore modeling was performed to identify the essential features responsible for the compounds' high activity. To gain insights into the 3D requirements for enhanced potency against the IP6K1 enzyme, we employed multiple alignment techniques to set up 3D-QSAR models. Given the absence of an available X-ray crystal structure for IP6K1, a reliable homology model for the enzyme was developed and structurally validated in order to perform structure-based analyses on the selected dataset compounds. Finally, molecular dynamic simulations, using the docked poses of these compounds, provided further insights. Our findings consistently supported the mechanistic interpretations derived from both ligand-based and structure-based analyses. This study offers valuable guidance on the design of novel IP6K1 inhibitors. Importantly, our work exclusively relies on non-commercial software packages, ensuring accessibility for reproducing the reported models.

SSF-DDI: a deep learning method utilizing drug sequence and substructure features for drug-drug interaction prediction

BACKGROUND

Drug-drug interactions (DDI) are prevalent in combination therapy, necessitating the importance of identifying and predicting potential DDI. While various artificial intelligence methods can predict and identify potential DDI, they often overlook the sequence information of drug molecules and fail to comprehensively consider the contribution of molecular substructures to DDI.

RESULTS

In this paper, we proposed a novel model for DDI prediction based on sequence and substructure features (SSF-DDI) to address these issues. Our model integrates drug sequence features and structural features from the drug molecule graph, providing enhanced information for DDI prediction and enabling a more comprehensive and accurate representation of drug molecules.

CONCLUSION

The results of experiments and case studies have demonstrated that SSF-DDI significantly outperforms state-of-the-art DDI prediction models across multiple real datasets and settings. SSF-DDI performs better in predicting DDI involving unknown drugs, resulting in a 5.67% improvement in accuracy compared to state-of-the-art methods.

Exploring Citrus sinensis Phytochemicals as Potential Inhibitors for Breast Cancer Genes BRCA1 and BRCA2 Using Pharmacophore Modeling, Molecular Docking, MD Simulations, and DFT Analysis

BACKGROUND

Structure-activity relationship (SAR) is considered to be an effective in silico approach when discovering potential antagonists for breast cancer due to gene mutation. Major challenges are faced by conventional SAR in predicting novel antagonists due to the discovery of diverse antagonistic compounds. Methodologyand Results: In predicting breast cancer antagonists, a multistep screening of phytochemicals isolated from the seeds of the Citrus sinensis plant was applied using feasible complementary methodologies. A three-dimensional quantitative structure-activity relationship (3D-QSAR) model was developed through the Flare project, in which conformational analysis, pharmacophore generation, and compound alignment were done. Ten hit compounds were obtained through the development of the 3D-QSAR model. For exploring the mechanism of action of active compounds against cocrystal inhibitors, molecular docking analysis was done through Molegro software (MVD) to identify lead compounds. Three new proteins, namely, 1T15, 3EU7, and 1T29, displayed the best Moldock scores. The quality of the docking study was assessed by a molecular dynamics simulation. Based on binding affinities to the receptor in the docking studies, three lead compounds (stigmasterol P8, epoxybergamottin P28, and nobiletin P29) were obtained, and they passed through absorption, distribution, metabolism, and excretion (ADME) studies via the SwissADME online service, which proved that P28 and P29 were the most active allosteric inhibitors with the lowest toxicity level against breast cancer. Then, density functional theory (DFT) studies were performed to measure the active compound's reactivity, hardness, and softness with the help of Gaussian 09 software.

CONCLUSIONS

This multistep screening of phytochemicals revealed high-reliability antagonists of breast cancer by 3D-QSAR using flare, docking analysis, and DFT studies. The present study helps in providing a proper guideline for the development of novel inhibitors of BRCA1 and BRCA2.

Advances in Drug Discovery and Design using Computer-aided Molecular Modeling

Computer-aided molecular modeling is a rapidly emerging technology that is being used to accelerate the discovery and design of new drug therapies. It involves the use of computer algorithms and 3D structures of molecules to predict interactions between molecules and their behavior in the body. This has drastically improved the speed and accuracy of drug discovery and design. Additionally, computer-aided molecular modeling has the potential to reduce costs, increase the quality of data, and identify promising targets for drug development. Through the use of sophisticated methods, such as virtual screening, molecular docking, pharmacophore modeling, and quantitative structure-activity relationships, scientists can achieve higher levels of efficacy and safety for new drugs. Moreover, it can be used to understand the activity of known drugs and simplify the process of formulating, optimizing, and predicting the pharmacokinetics of new and existing drugs. In conclusion, computer-aided molecular modeling is an effective tool to rapidly progress drug discovery and design by predicting the interactions between molecules and anticipating the behavior of new drugs in the body.

Decision tree-based identification of important molecular fragments for protein-ligand binding

Fragment-based drug design is an emerging technology in pharmaceutical research and development. One of the key aspects of this technology is the identification and quantitative characterization of molecular fragments. This study presents a strategy for identifying important molecular fragments based on molecular fingerprints and decision tree algorithms and verifies its feasibility in predicting protein-ligand binding affinity. Specifically, the three-dimensional (3D) structures of protein-ligand complexes are encoded using extended-connectivity fingerprints (ECFP), and three decision tree models, namely Random Forest, XGBoost, and LightGBM, are used to quantitatively characterize the feature importance, thereby extracting important molecular fragments with high reliability. Few-shot learning reveals that the extracted molecular fragments contribute significantly and consistently to the binding affinity even with a small sample size. Despite the absence of location and distance information for molecular fragments in ECFP, 3D visualization, in combination with the reverse ECFP process, shows that the majority of the extracted fragments are located at the binding interface of the protein and the ligand. This alignment with the distance constraints critical for binding affinity further supports the reliability of the strategy for identifying important molecular fragments.

Mpox (formerly monkeypox): pathogenesis, prevention, and treatment

In 2022, a global outbreak of Mpox (formerly monkeypox) occurred in various countries across Europe and America and rapidly spread to more than 100 countries and regions. The World Health Organization declared the outbreak to be a public health emergency of international concern due to the rapid spread of the Mpox virus. Consequently, nations intensified their efforts to explore treatment strategies aimed at combating the infection and its dissemination. Nevertheless, the available therapeutic options for Mpox virus infection remain limited. So far, only a few numbers of antiviral compounds have been approved by regulatory authorities. Given the high mutability of the Mpox virus, certain mutant strains have shown resistance to existing pharmaceutical interventions. This highlights the urgent need to develop novel antiviral drugs that can combat both drug resistance and the potential threat of bioterrorism. Currently, there is a lack of comprehensive literature on the pathophysiology and treatment of Mpox. To address this issue, we conducted a review covering the physiological and pathological processes of Mpox infection, summarizing the latest progress of anti-Mpox drugs. Our analysis encompasses approved drugs currently employed in clinical settings, as well as newly identified small-molecule compounds and antibody drugs displaying potential antiviral efficacy against Mpox. Furthermore, we have gained valuable insights from the process of Mpox drug development, including strategies for repurposing drugs, the discovery of drug targets driven by artificial intelligence, and preclinical drug development. The purpose of this review is to provide readers with a comprehensive overview of the current knowledge on Mpox.

Computer-Aided Drug Design and Drug Discovery: A Prospective Analysis

In the dynamic landscape of drug discovery, Computer-Aided Drug Design (CADD) emerges as a transformative force, bridging the realms of biology and technology. This paper overviews CADDs historical evolution, categorization into structure-based and ligand-based approaches, and its crucial role in rationalizing and expediting drug discovery. As CADD advances, incorporating diverse biological data and ensuring data privacy become paramount. Challenges persist, demanding the optimization of algorithms and robust ethical frameworks. Integrating Machine Learning and Artificial Intelligence amplifies CADDs predictive capabilities, yet ethical considerations and scalability challenges linger. Collaborative efforts and global initiatives, exemplified by platforms like Open-Source Malaria, underscore the democratization of drug discovery. The convergence of CADD with personalized medicine offers tailored therapeutic solutions, though ethical dilemmas and accessibility concerns must be navigated. Emerging technologies like quantum computing, immersive technologies, and green chemistry promise to redefine the future of CADD. The trajectory of CADD, marked by rapid advancements, anticipates challenges in ensuring accuracy, addressing biases in AI, and incorporating sustainability metrics. This paper concludes by highlighting the need for proactive measures in navigating the ethical, technological, and educational frontiers of CADD to shape a healthier, brighter future in drug discovery.

Enhancing Opioid Bioactivity Predictions through Integration of Ligand-Based and Structure-Based Drug Discovery Strategies with Transfer and Deep Learning Techniques

The opioid epidemic has cast a shadow over public health, necessitating immediate action to address its devastating consequences. To effectively combat this crisis, it is crucial to discover better opioid drugs with reduced addiction potential. Artificial intelligence-based and other machine learning tools, particularly deep learning models, have garnered significant attention in recent years for their potential to advance drug discovery. However, using these tools poses challenges, especially when training samples are insufficient to achieve adequate prediction performance. In this study, we investigate the effectiveness of transfer learning in building robust deep learning models to enhance ligand bioactivity prediction for each individual opioid receptor (OR) subtype. This is achieved by leveraging knowledge obtained from pretraining a model using supervised learning on a larger data set of bioactivity data combined with ligand-based and structure-based molecular descriptors related to the entire OR subfamily. Our studies hold the potential to advance opioid research by enabling the rapid identification of novel chemical probes with specific bioactivities, which can aid in the study of receptor function and contribute to the future development of improved opioid therapeutics.

Artificial intelligence in small molecule drug discovery from 2018 to 2023: Does it really work?

Utilizing artificial intelligence (AI) in drug design represents an advanced approach for identifying targets and developing new drugs. Integrating AI techniques significantly reduces the workload involved in drug development and enhances the efficiency of early-stage drug discovery. This review aims to present a comprehensive overview of the utilization of AI methods in the field of small drug design, with a specific focus on four key areas: protein structure prediction, molecular virtual screening, molecular design, and absorption, distribution, metabolism, excretion, and toxicity (ADMET) prediction. Additionally, the role and limitations of AI in drug development are explored, and the impact of AI on decision-making processes is studied. It is important to note that while AI can bring numerous benefits to the early stage of drug development, the direction and quality of decision-making should still be emphasized, as AI should be considered as a tool rather than a decisive factor.

Design of new Mcl-1 inhibitors for cancer using fragments hybridization, molecular docking, and molecular dynamics studies

Apoptosis is a critical process that regulates cell survival and death and plays an essential role in cancer development. The Bcl-2 protein family, including myeloid leukemia 1 (Mcl-1), is a key regulator of the intrinsic apoptosis pathway, and its overexpression in many human cancers has prompted efforts to develop Mcl-1 inhibitors as potential anticancer agents. In this study, we aimed to design new Mcl-1 inhibitors using various computational techniques. First, we used the Mcl-1 receptor-ligand complex to build an e-pharmacophore hypothesis and screened a library of 567,000 fragments from the Enamine database. We obtained 410 fragments and used them to design 92,384 novel compounds, which we then docked into the Mcl-1 binding cavity using HTVS, SP, and XP docking modes of Glide. To assess their suitability as drug candidates, we conducted MM-GBSA calculations and ADME prediction, leading to the identification of 10 compounds with excellent binding affinity and favorable pharmacokinetic properties. To further investigate the interaction strength, we performed molecular dynamics simulations on the top three Mcl-1 receptor-ligand complexes to study their interaction stability. Overall, our findings suggest that these compounds have promising potential as anticancer agents, pending further experimental validation such as Mcl-1 apoptosis Assay. By combining experimental methods with various in silico approaches, these techniques prove to be invaluable for identifying novel drug candidates with distinct therapeutic applications using fragment-based drug design. This methodology has the potential to expedite the drug discovery process while also reducing its costs.Communicated by Ramaswamy H. Sarma.

Thermal titration molecular dynamics (TTMD): shedding light on the stability of RNA-small molecule complexes

Ribonucleic acids are gradually becoming relevant players among putative drug targets, thanks to the increasing amount of structural data exploitable for the rational design of selective and potent binders that can modulate their activity. Mainly, this information allows employing different computational techniques for predicting how well would a ribonucleic-targeting agent fit within the active site of its target macromolecule. Due to some intrinsic peculiarities of complexes involving nucleic acids, such as structural plasticity, surface charge distribution, and solvent-mediated interactions, the application of routinely adopted methodologies like molecular docking is challenged by scoring inaccuracies, while more physically rigorous methods such as molecular dynamics require long simulation times which hamper their conformational sampling capabilities. In the present work, we present the first application of Thermal Titration Molecular Dynamics (TTMD), a recently developed method for the qualitative estimation of unbinding kinetics, to characterize RNA-ligand complexes. In this article, we explored its applicability as a post-docking refinement tool on RNA in complex with small molecules, highlighting the capability of this method to identify the native binding mode among a set of decoys across various pharmaceutically relevant test cases.

Artificial Intelligence/Machine Learning-Driven Small Molecule Repurposing via Off-Target Prediction and Transcriptomics

The process of discovering small molecule drugs involves screening numerous compounds and optimizing the most promising ones, both in vitro and in vivo. However, approximately 90% of these optimized candidates fail during trials due to unexpected toxicity or insufficient efficacy. Current concepts with respect to drug-protein interactions suggest that each small molecule interacts with an average of 6-11 targets. This implies that approved drugs and even discontinued compounds could be repurposed by leveraging their interactions with unintended targets. Therefore, we developed a computational repurposing framework for small molecules, which combines artificial intelligence/machine learning (AI/ML)-based and chemical similarity-based target prediction methods with cross-species transcriptomics information. This repurposing methodology incorporates eight distinct target prediction methods, including three machine learning methods. By using multiple orthogonal methods for a "dataset" composed of 2766 FDA-approved drugs targeting multiple therapeutic target classes, we identified 27,371 off-target interactions involving 2013 protein targets (i.e., an average of around 10 interactions per drug). Relative to the drugs in the dataset, we identified 150,620 structurally similar compounds. The highest number of predicted interactions were for drugs targeting G protein-coupled receptors (GPCRs), enzymes, and kinases with 10,648, 4081, and 3678 interactions, respectively. Notably, 17,283 (63%) of the off-target interactions have been confirmed in vitro. Approximately 4000 interactions had an IC50 of <100 nM for 1105 FDA-approved drugs and 1661 interactions had an IC50 of <10 nM for 696 FDA-approved drugs. Together, the confirmation of numerous predicted interactions and the exploration of tissue-specific expression patterns in human and animal tissues offer insights into potential drug repurposing for new therapeutic applications.

Enhancing Opioid Bioactivity Predictions through Integration of Ligand-Based and Structure-Based Drug Discovery Strategies with Transfer and Deep Learning Techniques

The opioid epidemic has cast a shadow over public health, necessitating immediate action to address its devastating consequences. To effectively combat this crisis, it is crucial to discover better opioid drugs with reduced addiction potential. Artificial intelligence-based and other machine learning tools, particularly deep learning models, have garnered significant attention in recent years for their potential to advance drug discovery. However, utilizing these tools poses challenges, especially when training samples are insufficient to achieve adequate prediction performance. In this study, we investigate the effectiveness of transfer learning using combined ligand-based and structure-based molecular descriptors from the entire opioid receptor (OR) subfamily in building robust deep learning models for enhanced bioactivity prediction of opioid ligands at each individual OR subtype. Our studies hold the potential to greatly advance opioid research by enabling the rapid identification of novel chemical probes with specific bioactivities, which can aid in the study of receptor function and contribute to the future development of improved opioid therapeutics.

Molecular Modeling Unveils the Effective Interaction of B-RAF Inhibitors with Rare B-RAF Insertion Variants

The Food and Drug Administration (FDA) has approved MAPK inhibitors as a treatment for melanoma patients carrying a mutation in codon V600 of the BRAF gene exclusively. However, BRAF mutations outside the V600 codon may occur in a small percentage of melanomas. Although these rare variants may cause B-RAF activation, their predictive response to B-RAF inhibitor treatments is still poorly understood. We exploited an integrated approach for mutation detection, tumor evolution tracking, and assessment of response to treatment in a metastatic melanoma patient carrying the rare p.T599dup B-RAF mutation. He was addressed to Dabrafenib/Trametinib targeted therapy, showing an initial dramatic response. In parallel, in-silico ligand-based homology modeling was set up and performed on this and an additional B-RAF rare variant (p.A598_T599insV) to unveil and justify the success of the B-RAF inhibitory activity of Dabrafenib, showing that it could adeptly bind both these variants in a similar manner to how it binds and inhibits the V600E mutant. These findings open up the possibility of broadening the spectrum of BRAF inhibitor-sensitive mutations beyond mutations at codon V600, suggesting that B-RAF V600 WT melanomas should undergo more specific investigations before ruling out the possibility of targeted therapy.

Correction: Bassani, D.; Moro, S. Past, Present, and Future Perspectives on Computer-Aided Drug Design Methodologies. Molecules 2023, 28, 3906

In the published publication [...].