Development of an In Silico Prediction Model for P-glycoprotein Efflux Potential in Brain Capillary Endothelial Cells toward the Prediction of Brain Penetration

A Practical In Silico Method for Predicting Compound Brain Concentration-Time Profiles: Combination of PK Modeling and Machine Learning

Given the aging populations in advanced countries globally, many pharmaceutical companies have focused on developing central nervous system (CNS) drugs. However, due to the blood-brain barrier, drugs do not easily reach the target area in the brain. Although conventional screening methods for drug discovery involve the measurement of (unbound fraction of drug) brain-to-plasma partition coefficients, it is difficult to consider nonequilibrium between plasma and brain compound concentration-time profiles. To truly understand the pharmacokinetics/pharmacodynamics of CNS drugs, compound concentration-time profiles in the brain are necessary; however, such analyses are costly and time-consuming and require a significant number of animals. Therefore, in this study, we attempted to develop an in silico prediction method that does not require a large amount of experimental data by combining modeling and simulation (M&S) with machine learning (ML). First, we constructed a hybrid model linking plasma concentration-time profile to the brain compartment that takes into account the transit time and brain distribution of each compound. Using mouse plasma and brain time experimental values for 103 compounds, we determined the brain kinetic parameters of the hybrid model for each compound; this case was defined as scenario I (a positive control experiment) and included the full brain concentration-time profile data. Next, we built an ML model using chemical structure descriptors as explanatory variables and rate parameters as the target variable, and we then input the predicted values from 5-fold cross-validation (CV) into the hybrid model; this case was defined as scenario II, in which no brain compound concentration-time profile data exist. Finally, for scenario III, assuming that the brain concentration is obtained at only one time point, we used the brain kinetic parameters from the result of the 5-fold CV in scenario II as the initial values for the hybrid model and performed parameter refitting against the observed brain concentration at that time point. As a result, the RMSE/R2-values of the brain compound concentration-time profiles over time were 0.445/0.517 in scenario II and 0.246/0.805 in scenario III, indicating the method provides high accuracy and suggesting that it is a practical method for predicting brain compound concentration-time profiles.

Artificial intelligence assisted identification of potential tau aggregation inhibitors: ligand- and structure-based virtual screening, in silico ADME, and molecular dynamics study

Alzheimer's disease (AD) is a severe, growing, multifactorial disorder affecting millions of people worldwide characterized by cognitive decline and neurodegeneration. The accumulation of tau protein into paired helical filaments is one of the major pathological hallmarks of AD and has gained the interest of researchers as a potential drug target to treat AD. Lately, Artificial Intelligence (AI) has revolutionized the drug discovery process by speeding it up and reducing the overall cost. As a part of our continuous effort to identify potential tau aggregation inhibitors, and leveraging the power of AI, in this study, we used a fully automated AI-assisted ligand-based virtual screening tool, PyRMD to screen a library of 12 million compounds from the ZINC database to identify potential tau aggregation inhibitors. The preliminary hits from virtual screening were filtered for similar compounds and pan-assay interference compounds (the compounds containing reactive functional groups which can interfere with the assays) using RDKit. Further, the selected compounds were prioritized based on their molecular docking score with the binding pocket of tau where the binding pockets were identified using replica exchange molecular dynamics simulation. Thirty-three compounds showing good docking scores for all the tau clusters were selected and were further subjected to in silico pharmacokinetic prediction. Finally, top 10 compounds were selected for molecular dynamics simulation and MMPBSA binding free energy calculations resulting in the identification of UNK_175, UNK_1027, UNK_1172, UNK_1173, UNK_1237, UNK_1518, and UNK_2181 as potential tau aggregation inhibitors.

Integrated use of ligand and structure-based virtual screening, molecular dynamics, free energy calculation and ADME prediction for the identification of potential PTP1B inhibitors

Protein tyrosine phosphatases (PTPs) are the group of enzymes that control both cellular activity and the dephosphorylation of tyrosine (Tyr)-phosphorylated proteins. Dysregulation of PTP1B has contributed to numerous diseases including Diabetes Mellitus, Alzheimer's disease, and obesity rendering PTP1B as a legitimate target for therapeutic applications. It is highly challenging to target this enzyme because of its highly conserved and positively charged active-site pocket motivating researchers to find novel lead compounds against it. The present work makes use of an integrated approach combining ligand-based and structure-based virtual screening to find hit compounds targeting PTP1B. Initially, pharmacophore modeling was performed to find common features like two hydrogen bond acceptors, an aromatic ring and one hydrogen bond donor from the potent PTP1B inhibitors. The dataset of compounds matching with the common pharmacophoric features was filtered to remove Pan-Assay Interference substructure and to match the Lipinski criteria. Then, compounds were further prioritized using molecular docking and top fifty compounds with good binding affinity were selected for absorption, distribution, metabolism, and excretion (ADME) predictions. The top five compounds with high solubility, absorption and permeability holding score of - 10 to - 9.3 kcal/mol along with Ertiprotafib were submitted to all-atom molecular dynamic (MD) studies. The MD studies and binding free energy calculations showed that compound M4, M5 and M8 were having better binding affinity for PTP1B enzyme with ∆Gtotal score of - 24.25, - 31.47 and - 33.81 kcal/mol respectively than other compounds indicating that compound M8 could be a suitable lead compound as PTP1B inhibitor.

Accurate prediction of Kp,uu,brain based on experimental measurement of Kp,brain and computed physicochemical properties of candidate compounds in CNS drug discovery

A mathematical equation model was developed by building the relationship between the fu,b/fu,p ratio and the computed physicochemical properties of candidate compounds, thereby predicting Kp,uu,brain based on a single experimentally measured Kp,brain value. A total of 256 compounds and 36 marketed published drugs including acidic, basic, neutral, zwitterionic, CNS-penetrant, and non-CNS penetrant compounds with diverse structures and physicochemical properties were involved in this study. A strong correlation was demonstrated between the fu,b/fu,p ratio and physicochemical parameters (CLogP and ionized fraction). The model showed good performance in both internal and external validations. The percentages of compounds with Kp,uu,brain predictions within 2-fold variability were 80.0 %-83.3 %, and more than 90 % were within a 3-fold variability. Meanwhile, "black box" QSAR models constructed by machine learning approaches for predicting fu,b/fu,p ratio based on the chemical descriptors are also presented, and the ANN model displayed the highest accuracy with an RMSE value of 0.27 and 86.7 % of the test set drugs fell within a 2-fold window of linear regression. These models demonstrated strong predictive power and could be helpful tools for evaluating the Kp,uu,brain by a single measurement parameter of Kp,brain during lead optimization for CNS penetration evaluation and ranking CNS drug candidate molecules in the early stages of CNS drug discovery.

Effects of Aqueous Boundary Layers and Paracellular Transport on the Efflux Ratio as a Measure of Active Transport Across Cell Layers

The efflux ratio (ER), determined by Caco-2/MDCK assays, is the standard in vitro metric to establish qualitatively whether a compound is a substrate of an efflux transporter. However, others have also enabled the utilisation of this metric quantitatively by deriving a relationship that expresses the ER as a function of the intrinsic membrane permeability of the membrane (P0) as well as the permeability of carrier-mediated efflux (Ppgp). As of yet, Ppgp cannot be measured directly from transport experiments or otherwise, but the ER relationship provides easy access to this value if P0 is known. However, previous derivations of this relationship failed to consider the influence of additional transport resistances such as the aqueous boundary layers (ABLs) and the filter on which the monolayer is grown. Since single fluxes in either direction can be heavily affected by these experimental artefacts, it is crucial to consider the potential impact on the ER. We present a model that includes these factors and show both mathematically and experimentally that this simple ER relationship also holds for the more realistic scenario that does not neglect the ABLs/filter. Furthermore, we also show mathematically how paracellular transport affects the ER, and we experimentally confirm that paracellular dominance reduces the ER to unity and can mask potential efflux.

The Trends and Future Prospective of In Silico Models from the Viewpoint of ADME Evaluation in Drug Discovery

Drug discovery and development are aimed at identifying new chemical molecular entities (NCEs) with desirable pharmacokinetic profiles for high therapeutic efficacy. The plasma concentrations of NCEs are a biomarker of their efficacy and are governed by pharmacokinetic processes such as absorption, distribution, metabolism, and excretion (ADME). Poor ADME properties of NCEs are a major cause of attrition in drug development. ADME screening is used to identify and optimize lead compounds in the drug discovery process. Computational models predicting ADME properties have been developed with evolving model-building technologies from a simplified relationship between ADME endpoints and physicochemical properties to machine learning, including support vector machines, random forests, and convolution neural networks. Recently, in the field of in silico ADME research, there has been a shift toward evaluating the in vivo parameters or plasma concentrations of NCEs instead of using predictive results to guide chemical structure design. Another research hotspot is the establishment of a computational prediction platform to strengthen academic drug discovery. Bioinformatics projects have produced a series of in silico ADME models using free software and open-access databases. In this review, we introduce prediction models for various ADME parameters and discuss the currently available academic drug discovery platforms.

Study Models of Drug-Drug Interactions Involving P-Glycoprotein: The Potential Benefit of P-Glycoprotein Modulation at the Kidney and Intestinal Levels

P-glycoprotein (P-gp) is a crucial membrane transporter situated on the cell's apical surface, being responsible for eliminating xenobiotics and endobiotics. P-gp modulators are compounds that can directly or indirectly affect this protein, leading to changes in its expression and function. These modulators can act as inhibitors, inducers, or activators, potentially causing drug-drug interactions (DDIs). This comprehensive review explores diverse models and techniques used to assess drug-induced P-gp modulation. We cover several approaches, including in silico, in vitro, ex vivo, and in vivo methods, with their respective strengths and limitations. Additionally, we explore the therapeutic implications of DDIs involving P-gp, with a special focus on the renal and intestinal elimination of P-gp substrates. This involves enhancing the removal of toxic substances from proximal tubular epithelial cells into the urine or increasing the transport of compounds from enterocytes into the intestinal lumen, thereby facilitating their excretion in the feces. A better understanding of these interactions, and of the distinct techniques applied for their study, will be of utmost importance for optimizing drug therapy, consequently minimizing drug-induced adverse and toxic effects.

Predictive Models Based on Molecular Images and Molecular Descriptors for Drug Screening

Various toxicity and pharmacokinetic evaluations as screening experiments are needed at the drug discovery stage. Currently, to reduce the use of animal experiments and developmental expenses, the development of high-performance predictive models based on quantitative structure-activity relationship analysis is desired. From these evaluation targets, we selected 50% lethal dose (LD50), blood-brain barrier penetration (BBBP), and the clearance (CL) pathway for this investigation and constructed predictive models for each target using 636-11,886 compounds. First, we constructed predictive models using the DeepSnap-deep learning (DL) method and images of compounds as features. The calculated area under the curve (AUC) and balanced accuracy (BAC) were, respectively, 0.887 and 0.818 for LD50, 0.893 and 0.824 for BBBP, and 0.883 and 0.763 for the CL pathway. Next, molecular descriptors (MDs) of compounds were calculated using Molecular Operating Environment, alvaDesc, and ADMET Predictor to construct predictive models using the MD-based method. Using these MDs, we constructed predictive models using DataRobot. The calculated AUC and BAC were, respectively, 0.931 and 0.805 for LD50, 0.919 and 0.849 for BBBP, and 0.900 and 0.807 for the CL pathway. In this investigation, we constructed predictive models combining the DeepSnap-DL and MD-based methods. In ensemble models using the mean predictive probability of the DeepSnap-DL and MD-based methods, the calculated AUC and BAC were, respectively, 0.942 and 0.842 for LD50, 0.936 and 0.853 for BBBP, and 0.908 and 0.832 for the CL pathway, with improved predictive performance observed for all variables compared with either single method alone. Moreover, in consensus models that adopted only compounds for which the results of the two methods agreed, the calculated BAC for LD50, BBBP, and the CL pathway were 0.916, 0.918, and 0.847, respectively, indicating higher predictive performance than the ensemble models for all three variables. The predictive models combining the DeepSnap-DL and MD-based methods displayed high predictive performance for LD50, BBBP, and the CL pathway. Therefore, the application of this approach to prediction targets in various drug discovery screenings is expected to accelerate drug discovery.

Human Brain Penetration Prediction Using Scaling Approach from Animal Machine Learning Models

Machine learning (ML) approaches have been applied to predicting drug pharmacokinetic properties. Previously, we predicted rat unbound brain-to-plasma ratio (Kpuu,brain) by ML models. In this study, we aimed to predict human Kpuu,brain through animal ML models. First, we re-evaluated ML models for rat Kpuu,brain prediction by using trendy open-source packages. We then developed ML models for monkey Kpuu,brain prediction. Leave-one-out cross validation was utilized to rationally build models using a relatively small dataset. After establishing the monkey and rat ML models, human Kpuu,brain prediction was achieved by implementing the animal models considering appropriate scaling methods. Mechanistic NeuroPK models for the identical monkey and human dataset were treated as the criteria for comparison. Results showed that rat Kpuu,brain predictivity was successfully replicated. The optimal ML model for monkey Kpuu,brain prediction was superior to the NeuroPK model, where accuracy within 2-fold error was 78% (R2 = 0.76). For human Kpuu,brain prediction, rat model using relative expression factor (REF), scaled transporter efflux ratios (ERs), and monkey model using in vitro ERs can provide comparable predictivity to the NeuroPK model, where accuracy within 2-fold error was 71% and 64% (R2 = 0.30 and 0.52), respectively. We demonstrated that ML models can deliver promising Kpuu,brain prediction with several advantages: (1) predict reasonable animal Kpuu,brain; (2) prospectively predict human Kpuu,brain from animal models; and (3) can skip expensive monkey studies for human prediction by using the rat model. As a result, ML models can be a powerful tool for drug Kpuu,brain prediction in the discovery stage.

DruMAP: A Novel Drug Metabolism and Pharmacokinetics Analysis Platform

We developed a novel drug metabolism and pharmacokinetics (DMPK) analysis platform named DruMAP. This platform consists of a database for DMPK parameters and programs that can predict many DMPK parameters based on the chemical structure of a compound. The DruMAP database includes curated DMPK parameters from public sources and in-house experimental data obtained under standardized conditions; it also stores predicted DMPK parameters produced by our prediction programs. Users can predict several DMPK parameters simultaneously for novel compounds not found in the database. Furthermore, the highly flexible search system enables users to search for compounds as they desire. The current version of DruMAP comprises more than 30,000 chemical compounds, about 40,000 activity values (collected from public databases and in-house data), and about 600,000 predicted values. Our platform provides a simple tool for searching and predicting DMPK parameters and is expected to contribute to the acceleration of new drug development. DruMAP can be freely accessed at: https://drumap.nibiohn.go.jp/.

A Computational Physics-based Approach to Predict Unbound Brain-to-Plasma Partition Coefficient, Kp,uu

The blood-brain barrier (BBB) plays a critical role in preventing harmful endogenous and exogenous substances from penetrating the brain. Optimal brain penetration of small-molecule central nervous system (CNS) drugs is characterized by a high unbound brain/plasma ratio (K_p,uu). While various medicinal chemistry strategies and in silico models have been reported to improve BBB penetration, they have limited application in predicting K_p,uu directly. We describe a physics-based computational approach, a quantum mechanics (QM)-based energy of solvation (E-sol), to predict K_p,uu. Prospective application of this method in internal CNS drug discovery programs highlights the utility and accuracy of this new method, which showed a categorical accuracy of 79% and an R² of 0.61 from a linear regression model.

Hydrogen Sulfide Attenuates Lipopolysaccharide-Induced Inflammation via the P-glycoprotein and NF-κB Pathway in Astrocytes

Astrocyte activation is key in neurodegenerative diseases. Hydrogen sulfide (H2S) exhibits neuroprotective effects on astrocytes, although the underlying molecular mechanism remains unclear. Here, we explored the effects of H2S on lipopolysaccharide (LPS)-induced astrocyte activation and astrocyte-mediated neuroinflammation. After inducing primary astrocytes via LPS exposure, H2S levels were altered. The generation and secretion of inflammatory mediators by astrocytes and their interrelation with P-glycoprotein (P-gp), an important transporter belonging to the ABC transporter family, were assessed. Activated astrocytes showed upregulated glial fibrillary acidic protein (GFAP) mRNA expression, and significantly increased proinflammatory factor mRNA/protein expression and release. The secretory capacity of astrocytes was reduced, with significantly decreased proinflammatory factor levels in culture supernatant after P-gp inhibitor verapamil pretreatment. The increase in the intracellular H2S level inhibited LPS-induced GFAP expression and P65 nuclear entry in astrocytes. mRNA expression and release of proinflammatory factors were reduced significantly, with no significant changes in cytoplasmic protein expression. S-sulfhydration levels increased significantly with the increased concentration of sodium hydrosulfide or S-adenosyl-l-methionine addition, with only moderate changes in astrocyte P-gp expression. H2S regulates NF-κB activation, leads to S-sulfhydration of P-gp, and inhibits the biosynthesis and secretion of proinflammatory factors by astrocytes. The regulatory effects of H2S on astrocytes may have clinical value for exploring new therapeutic strategies against neurodegenerative diseases.

Recent Advances on P-Glycoprotein (ABCB1) Transporter Modelling with In Silico Methods

ABC transporters play a critical role in both drug bioavailability and toxicity, and with the discovery of the P-glycoprotein (P-gp), this became even more evident, as it plays an important role in preventing intracellular accumulation of toxic compounds. Over the past 30 years, intensive studies have been conducted to find new therapeutic molecules to reverse the phenomenon of multidrug resistance (MDR) ), that research has found is often associated with overexpression of P-gp, the most extensively studied drug efflux transporter; in MDR, therapeutic drugs are prevented from reaching their targets due to active efflux from the cell. The development of P-gp inhibitors is recognized as a good way to reverse this type of MDR, which has been the subject of extensive studies over the past few decades. Despite the progress made, no effective P-gp inhibitors to reverse multidrug resistance are yet on the market, mainly because of their toxic effects. Computational studies can accelerate this process, and in silico models such as QSAR models that predict the activity of compounds associated with P-gp (or analogous transporters) are of great value in the early stages of drug development, along with molecular modelling methods, which provide a way to explain how these molecules interact with the ABC transporter. This review highlights recent advances in computational P-gp research, spanning the last five years to 2022. Particular attention is given to the use of machine-learning approaches, drug-transporter interactions, and recent discoveries of potential P-gp inhibitors that could act as modulators of multidrug resistance.

Strategies for targeting the P2Y12 receptor in the central nervous system

The purinergic 2Y type 12 receptor (P2Y₁₂R) is a well-known biological target for anti-thrombotic drugs due to its role in platelet aggregation and blood clotting. While the importance of the P2Y₁₂R in the periphery has been known for decades, much less is known about its expression and roles in the central nervous system (CNS), where it is expressed exclusively on microglia - the first responders to brain insults and neurodegeneration. Several seminal studies have shown that P2Y₁₂ is a robust, translatable biomarker for anti-inflammatory and neuroprotective microglial phenotypes in models of degenerative diseases such as multiple sclerosis and Alzheimer's disease. An enduring problem for studying this receptor in vivo, however, is the lack of selective, high-affinity small molecule ligands that can bypass the blood-brain barrier and accumulate in the CNS. In this Digest, we discuss previous attempts by researchers to target the P2Y₁₂R in the CNS and opine on strategies that may be employed to design and assess the suitability of novel P2Y₁₂ ligands for this purpose going forward.

Comparison of the Blood-Brain Barrier Transport and Vulnerability to P-Glycoprotein-Mediated Drug-Drug Interaction of Domperidone versus Metoclopramide Assessed Using In Vitro Assay and PET Imaging

Domperidone and metoclopramide are widely prescribed antiemetic drugs with distinct neurological side effects. The impact of P-glycoprotein (P-gp)-mediated efflux at the blood−brain barrier (BBB) on brain exposure and BBB permeation was compared in vitro and in vivo using positron emission tomography (PET) imaging in rats with the radiolabeled analogs [11C]domperidone and [11C]metoclopramide. In P-gp-overexpressing cells, the IC50 of tariquidar, a potent P-gp inhibitor, was drastically different using [11C]domperidone (221 nM [198−248 nM]) or [11C]metoclopramide (4 nM [2−8 nM]) as the substrate. Complete P-gp inhibition led to a 1.8-fold higher increase in the cellular uptake of [11C]domperidone compared with [11C]metoclopramide (p < 0.0001). Brain PET imaging revealed that the baseline brain exposure (AUCbrain) of [11C]metoclopramide was 2.4-fold higher compared with [11C]domperidone (p < 0.001), consistent with a 1.8-fold higher BBB penetration (AUCbrain/AUCplasma). The maximal increase in the brain exposure (2.9-fold, p < 0.0001) and BBB penetration (2.9-fold, p < 0.0001) of [11C]metoclopramide was achieved using 8 mg/kg of tariquidar. In comparison, neither 8 nor 15 mg/kg of tariquidar increased the brain exposure of [11C]domperidone (p > 0.05). Domperidone is an avid P-gp substrate that was in vitro compared with metoclopramide. Domperidone benefits from a lower brain exposure and a limited risk for P-gp-mediated drug−drug interaction involving P-gp inhibition at the BBB.

In Silico Approaches for Addressing Challenges in CNS Radiopharmaceutical Design

Positron emission tomography (PET) is a highly sensitive and versatile molecular imaging modality that leverages radiolabeled molecules, known as radiotracers, to interrogate biochemical processes such as metabolism, enzymatic activity, and receptor expression. The ability to probe specific molecular and cellular events longitudinally in a noninvasive manner makes PET imaging a particularly powerful technique for studying the central nervous system (CNS) in both health and disease. Unfortunately, developing and translating a single CNS PET tracer for clinical use is typically an extremely resource-intensive endeavor, often requiring synthesis and evaluation of numerous candidate molecules. While existing in vitro methods are beginning to address the challenge of derisking molecules prior to costly in vivo PET studies, most require a significant investment of resources and possess substantial limitations. In the context of CNS drug development, significant time and resources have been invested into the development and optimization of computational methods, particularly involving machine learning, to streamline the design of better CNS therapeutics. However, analogous efforts developed and validated for CNS radiotracer design are conspicuously limited. In this Perspective, we overview the requirements and challenges of CNS PET tracer design, survey the most promising computational methods for in silico CNS drug design, and bridge these two areas by discussing the potential applications and impact of computational design tools in CNS radiotracer design.

The Application of in silico Methods for Prediction of Blood-Brain Barrier Permeability of Small Molecule PET Tracers

Designing positron emission tomography (PET) tracers for targets in the central nervous system (CNS) is challenging. Besides showing high affinity and high selectivity for their intended target, these tracers have to be able to cross the blood-brain barrier (BBB). Since only a small fraction of small molecules is estimated to be able to cross the BBB, tools that can predict permeability at an early stage during the development are of great importance. One such tool is in silico models for predicting BBB-permeability. Thus far, such models have been built based on CNS drugs, with one exception. Herein, we sought to discuss and analyze if in silico predictions that have been built based on CNS drugs can be applied for CNS PET tracers as well, or if dedicated models are needed for the latter. Depending on what is taken into account in the prediction, i.e., passive diffusion or also active influx/efflux, there may be a need for a model build on CNS PET tracers. Following a brief introduction, an overview of a few selected in silico BBB-permeability predictions is provided along with a short historical background to the topic. In addition, a combination of previously reported CNS PET tracer datasets were assessed in a couple of selected models and guidelines for predicting BBB-permeability. The selected models were either predicting only passive diffusion or also the influence of ADME (absorption, distribution, metabolism and excretion) parameters. To conclude, we discuss the potential need of a prediction model dedicated for CNS PET tracers and present the key issues in respect to setting up a such a model.

Nanomedicine for brain cancer

CNS tumors remain among the deadliest forms of cancer, resisting conventional and new treatment approaches, with mortality rates staying practically unchanged over the past 30 years. One of the primary hurdles for treating these cancers is delivering drugs to the brain tumor site in therapeutic concentration, evading the blood-brain (tumor) barrier (BBB/BBTB). Supramolecular nanomedicines (NMs) are increasingly demonstrating noteworthy prospects for addressing these challenges utilizing their unique characteristics, such as improving the bioavailability of the payloadsviacontrolled pharmacokinetics and pharmacodynamics, BBB/BBTB crossing functions, superior distribution in the brain tumor site, and tumor-specific drug activation profiles. Here, we review NM-based brain tumor targeting approaches to demonstrate their applicability and translation potential from different perspectives. To this end, we provide a general overview of brain tumor and their treatments, the incidence of the BBB and BBTB, and their role on NM targeting, as well as the potential of NMs for promoting superior therapeutic effects. Additionally, we discuss critical issues of NMs and their clinical trials, aiming to bolster the potential clinical applications of NMs in treating these life-threatening diseases.

Computational strategies to identify new drug candidates against neuroinflammation

The even more increasing application of computational approaches in these last decades has deeply modified the process of discovery and commercialization of new therapeutic entities. This is especially true in the field of neuroinflammation, in which both the peculiar anatomical localization and the presence of the blood-brain barrier makeit mandatory to finely tune the candidates' physicochemical properties from the early stages of the discovery pipeline. The aim of this review is therefore to provide a general overview to the readers about the topic of neuroinflammation, together with the most common computational strategies that can be exploited to discover and design small molecules controlling neuroinflammation, especially those based on the knowledge of the three-dimensional structure of the biological targets of therapeutic interest. The techniques used to describe the molecular recognition mechanisms, such as molecular docking and molecular dynamics, will therefore be eviscerated, highlighting their advantages and their limitations. Finally, we report several case studies in which computational methods have been applied in drug discovery on neuroinflammation, focusing on the last decade's research.

Drug Design Targeting the Muscarinic Receptors and the Implications in Central Nervous System Disorders

There is substantial evidence that cholinergic system function impairment plays a significant role in many central nervous system (CNS) disorders. During the past three decades, muscarinic receptors (mAChRs) have been implicated in various pathologies and have been prominent targets of drug-design efforts. However, due to the high sequence homology of the orthosteric binding site, many drug candidates resulted in limited clinical success. Although several advances in treating peripheral pathologies have been achieved, targeting CNS pathologies remains challenging for researchers. Nevertheless, significant progress has been made in recent years to develop functionally selective orthosteric and allosteric ligands targeting the mAChRs with limited side effect profiles. This review highlights past efforts and focuses on recent advances in drug design targeting these receptors for Alzheimer’s disease (AD), schizophrenia (SZ), and depression.