Vitamin D supplementation affects bone marrow-derived mesenchymal stem cells differentiation into insulin-producing cells.

Background this study was conducted to assess the effects of vitamin D on differentiation of bone marrow- derived mesenchymal stem cells (BM-MSCs) into insulin producing cells (IPCs). Method BM-MSCs were isolated from femur and tibia of rats and incubated in low (LG) or high glucose (HG) (5mM or 25mM), or high glucose DMEM media supplemented with vitamin D (0.2nM) (HGD) for 14 days. Cells viability was analysis by MTT assay. Differentiation of SCs was confirmed using measuring genes expression level of pdx1 and insulin, and insulin secretion, glucose stimulated insulin secretion, and insulin content by ELISA method. Results Cell viability was significantly higher in HGD than LG (p < 0.05) in day 3, also, in HG and HGD than LG (p < 0.001), and HGD vs. HG (p < 0.001) in day 7. Pdx1 and insulin level was markedly higher in HGD than LG (p < 0.05 and p < 0.01). pdx1 expression was markedly higher in HGD (p < 0.05) than LG, also insulin expression the HG (p < 0.05), and HGD (p < 0.01) groups compared to the LG group. Insulin release at 5mM glucose was notably higher in the HGD group compared to LG (p < 0.05), and at 25mM glucose, both HG and HGD showed significant increases vs. LG (p < 0.05 and p < 0.01, respectively). Insulin content was significantly higher in both 5mM and 25mM glucose for HG and HGD vs. LG (p < 0.01 and p < 0.001, respectively). In conclusion, treatment BM-MSCs with vitamin D could increase their differentiation into IPCs and it can be considered as a potential supplementary agent in enhancing differentiation SCs into insulin generating cells.

Prediction of the First-Pass Metabolism of a Drug After Oral Intake Based on Structural Parameters and Physicochemical Properties.

BACKGROUND AND OBJECTIVE: The oral first-pass metabolism is a crucial factor that plays a key role in a drug's pharmacokinetic profile. Prediction of the oral first-pass metabolism based on chemical structural parameters can be useful in the drug-design process. Developing an orally administered drug with an acceptable pharmacokinetic profile is necessary to reduce the cost and time associated with evaluating the extent of the first-pass metabolism of a candidate compound in preclinical studies. The aim of this study is to estimate the first-pass metabolism of an orally administered drug. METHODS: A set of compounds with reported first-pass metabolism data were collected. Moreover, human intestinal absorption percentage and oral bioavailability data were extracted from the literature to propose a classification system that split the drugs up based on their first-pass metabolism extents. Various structural parameters were calculated for each compound. The relations of the structural and physicochemical values of each compound to the class the compound belongs to were obtained using logistic regression. RESULTS: Initial analysis showed that compounds with logD7.4 > 1 or a rugosity factor of > 1.5 are more likely to have high first-pass metabolism. Four different models that can predict the oral first-pass metabolism with acceptable error were introduced. The overall accuracies of the models were in the range of 72% (for models with simple descriptors) to 78% (for models with complex descriptors). Although the models with simple descriptors have lower accuracies compared to complex models, they are more interpretable and easier for researchers to utilize. CONCLUSION: A novel classification of drugs based on the extent of the oral first-pass metabolism was introduced, and mechanistic models were developed to assign candidate compounds to the appropriate proposed classes.

Construction of a bacteriophage-derived vector with potential applications in targeted drug delivery and cell imaging.

There is a strong relationship between the dysregulation of epidermal growth factor receptor (EGFR) and the development of epithelial-derived cancers. Therefore, EGFR has usually been considered the desired target for gene therapy. Here, we propose an approach for targeting EGFR-expressing cells by phage particles capable of displaying EGF and GFP as tumor-targeting and reporting elements, respectively. For this purpose, the superfolder GFP-EGF (sfGFP-EGF) coding sequence was inserted at the N-terminus of the pIII gene in the pIT2 phagemid. The capability of the constructed phage to recognize EGFR-overexpressing cells was monitored by fluorescence microscopy, fluorescence-activated cell sorting (FACS), and cell-based ELISA experiments. FACS analysis showed a significant shift in the mean fluorescence intensity (MFI) of the cells treated with phage displaying sfGFP-EGF compared to phage displaying only sfGFP. The binding of phage displaying sfGFP-EGF to A-431 cells, monitored by fluorescence microscopy, indicated the formation of the sfGFP-EGF-EGFR complex on the surface of the treated cells. Cell-based ELISA experiments showed that phages displaying either EGF or sfGFP-EGF can specifically bind EGFR-expressing cells. The vector constructed in the current study has the potential to be engineered for gene delivery purposes as well as cell-based imaging for tumor detection.

The Effect of Mesenchymal Stem Cells Derived-Conditioned Media in Combination with Oral Anti-Androgenic Drugs on Male Pattern Baldness: An Animal Study.

OBJECTIVE: Androgenetic alopecia (AGA) is a prevalent form of hair loss, mainly caused by follicular sensitivity to androgens. Despite developing different anti-androgen treatment options, the success rate of these treatments has been limited. Using animal models, this study evaluated the therapeutic effects of umbilical cord (UC) stem cell conditioned media (CM) combined with oral anti-androgens for hair regeneration. MATERIALS AND METHODS: In this experimental study, Poloxamer 407 (P407) was used as a drug carrier for subcutaneous testosterone injection. AGA models were treated with oral finasteride, oral flutamide, and CM injections. Samples were thoroughly evaluated and compared using histological, stereological, and molecular analyses. RESULTS: Injecting CM-loaded hydrogel alone or combined with oral intake of anti-androgens improved hair regeneration. These treatments could promote hair growth by inducing hair follicles in the anagen stage and shortening the telogen and catagen phases. Furthermore, the combination treatment led to an upregulation of hair induction gene expression with a downregulation of inflammation genes. CONCLUSION: Through a reduction in inflammation, injection of CM-loaded hydrogel alone or combined with oral intake of anti-androgens induces the hair cell cycle with regeneration in damaged follicles. Hence, this could be a promising therapeutic method for AGA patients.

Expression, Purification and Characterization of Functional Teduglutide Using GST Fusion System in Prokaryotic Cells.

Purpose: Teduglutide is the first and only FDA-approved drug for long-term treatment of short bowel syndrome (SBS). The current study aimed to present an approach for production of teduglutide using recombinant DNA technology. Methods: The coding gene for teduglutide was cloned into pGEX-2T vector, where coding sequence for factor Xa cleavage site was added between GST and teduglutide coding genes. The GST-teduglutide protein was overexpressed in E. coli BL21 (DE3) strain and affinity purified using glutathione sepharose affinity column. Results: On-column proteolytic activity of factor Xa followed by size exclusion chromatography resulted in the pure teduglutide. Circular dichroism (CD) spectropolarimetry showed that the produced teduglutide folds into mainly α-helical structure (>50%), as expected. In mass spectroscopy analysis, the fragments of teduglutide resulted by cyanogen bromide cleavage as well as those expected theoretically due to mass fragmentation were identified. The functionality of the produced peptide was evaluated by measuring its proliferative effect on Caco2 intestinal epithelial cells, and the results indicated that produced teduglutide induces cell proliferation by 19±0.30 and 33±7.82 % at 1.21 and 3.64 µM concentrations, respectively, compared to untreated cells. Conclusion: Teduglutide was successfully expressed and purified and its functionality and structural integrity were confirmed by in vitro experiments. We believe that the experimental-scale method presented in the current study can be useful for pilot-scale and also industrial-scale production of teduglutide.

Advantages of Material Biofunctionalization Using Nucleic Acid Aptamers in Tissue Engineering and Regenerative Medicine.

Material engineering is a fundamental issue in the applications of materials in the medical field. One of the aspects of material engineering is incorporating recognition sites on the surface of biomaterials, which plays an essential role in increasing the efficiency of tissue engineering scaffolds in various aspects. The application of peptides and antibodies to establish the recognition and adhesion sites has limitations, such as fragility and instability under physical and chemical processes. Therefore, synthetic ligands such as nucleic acid aptamers have received much attention for easy synthesis, minimal immunogenicity, high specificity, and stability under processing. Due to the effective role of these ligands in increasing the efficiency of engineered constructs in this study, the advantages of nucleic acid aptamers in tissue engineering will be reviewed. Aptamer-functionalized biomaterials can attract endogenous stem cells to wounded areas and organize their actions to facilitate tissue regeneration. This approach harnesses the body's inherent regeneration potential to treat many diseases. Also, increased efficacy in controlled release, slow and targeted drug delivery are important issues in drug delivery for tissue engineering approaches which can be achieved by incorporating aptamers in drug delivery systems. Aptamer-functionalized scaffolds have very applications, such as diagnosis of cancer, hematological infections, narcotics, heavy metals, toxins, controlled release from the scaffolds, and in vivo cell tracing. Aptasensors, as a result of many advantages over other traditional assay methods, can replace older methods. Furthermore, their unique targeting mechanism also targets compounds with no particular receptors. Targeting cell homing, local and targeted drug delivery, cell adhesion efficacy, cytocompatibility and bioactivity of scaffolds, aptamer-based biosensor, and aptamer-functionalized scaffolds are the topics that will be examined in this review study.

Drug Repurposing for Identification of S1P1 Agonists with Potential Application in Multiple Sclerosis Using In Silico Drug Design Approaches.

Purpose: Drug repurposing is an approach successfully used for discovery of new therapeutic applications for the existing drugs. The current study was aimed to use the combination of in silico methods to identify FDA-approved drugs with possible S1P1 agonistic activity useful in multiple sclerosis (MS). Methods: For this, a 3D-QSAR model for the known 21 S1P1 agonists were generated based on 3D-QSAR approach and used to predict the possible S1P1 agonistic activity of FDA-approved drugs. Then, the selected compounds were screened by docking into S1P1 and S1P3 receptors to select the S1P1 potent and selective compounds. Further evaluation was carried out by molecular dynamics (MD) simulation studies where the S1P1 binding energies of selected compounds were calculated. Results: The analyses resulted in identification of cobicistat, benzonatate and brigatinib as the selective and potent S1P1 agonists with the binding energies of -85.93, -69.77 and -67.44 kcal. mol-1, calculated using MM-GBSA algorithm based on 50 ns MD simulation trajectories. These values are better than that of siponimod (-59.35 kcal mol-1), an FDA approved S1P1 agonist indicated for MS treatment. Furthermore, similarity network analysis revealed that cobicistat and brigatinib are the most structurally favorable compounds to interact with S1P1. Conclusion: The findings in this study revealed that cobicistat and brigatinib can be evaluated in experimental studies as potential S1P1 agonist candidates useful in the treatment of MS.

S1PR1 modulators in multiple sclerosis: Efficacy, safety, comparison, and chemical structure insights.

Multiple sclerosis (MS) is a neurological disease that leads to severe physical and cognitive disabilities. Drugs used in the treatment of MS vary from small synthetic molecules to large macromolecules such as antibodies. Sphingosine 1-phosphate receptor modulators are frequently used for the treatment of MS. These medicines prevent the egress of lymphocytes from secondary lymphoid organs leading to immune system suppression. Currently, four S1PR modulators are on the market and several potential drug candidates are in clinical trials for the treatment of MS. These compounds differ in chemical structure, adverse effects, and efficacy points of view. The current article reviews the latest studies on S1PR1 modulators and compares them with other MS drugs in terms of efficacy, tolerability, and safety. A special focus was dedicated to discussing the structure-activity relationships of these compounds and performing a three-dimensional quantitative structure-activity relationship (3D-QSAR) analysis to gain better insight into the ligand-receptor interaction mode.

Synthesis of Mn-ofloxacin complex, experimental and in silico DNA binding evaluation, biological activity assessment.

Metal-fluoroquinolones have more antibacterial and cytotoxic effects compared to free fluoroquinolones. In this work, a bidentated Mn (II) complex with ofloxacin (MOC) was synthesized and its cytotoxicity activity, oxidative stress and DNA binding were studied. Anti- proliferative and cytotoxic tests revealed that MOC exhibits better anti proliferative and cytotoxic activities compared to OFL which was attributed to the more interaction of MOC with DNA. Therefore, the interaction of MOC with DNA was investigated by using voltammetry, UV-Vis, fluorescence, and in silico methods. The results revealed that MOC interacts with DNA via electrostatic and outside hydrogen binding via minor groove. The proposed DNA binding modes may support the greater in-vitro cytotoxicity of MOC compared to OFL alone.

Engineered cells along with smart scaffolds: critical factors for improving tissue engineering approaches.

In this review, gene delivery and its applications are discussed in tissue engineering (TE); also, new techniques such as the CRISPR-Cas9 system, synthetics biology and molecular dynamics simulation to improve the efficiency of the scaffolds have been studied. CRISPR-Cas9 is expected to make significant advances in TE in the future. The fundamentals of synthetic biology have developed powerful and flexible methods for programming cells via artificial genetic circuits. The combination of regenerative medicine and artificial biology allows the engineering of cells and organisms for use in TE, biomaterials, bioprocessing and scaffold development. The dynamics of protein adsorption at the scaffold surface at the atomic level can provide valuable guidelines for the future design of TE scaffolds /implants.

Development of Bio-artificial Esophageal Tissue Engineering Utilization for Circumferential Lesion Transplantation: A Narrative Review.

The esophagus is the gastrointestinal tract's primary organ that transfers bolus into the stomach with peristaltic motion. Therefore, its lesions cause a significant disturbance in the nutrition and digestive system. Esophageal disease treatment sometimes requires surgical procedures that involve removal and circumferential full-thickness replacement. Unlike other organs, the esophagus has a limited regeneration ability and cannot be transplanted from donors. There are various methods of restoring the esophageal continuity; however, they are associated with certain flaws that lead to a non-functional recovery. As an exponentially growing science, tissue engineering has become a leading technique for the development of tissue replacement to repair damaged esophageal segments. Scaffold plays a significant role in the process of tissue engineering, as it acts as a template for the regeneration of growing tissue. A variety of scaffolds have been studied to replace the esophagus. Due to the many tissue quality challenges, the results are still inadequate and need to be improved. The success of esophageal tissue regeneration will finally depend on the scaffold's capability to mimic natural tissue properties and provide a qualified environment for regeneration. Thereby, scaffold fabrication techniques are fundamental. This article reviews the recent developments in esophageal tissue engineering for the treatment of circumferential lesions based on scaffold biomaterial engineering approaches.

In silico models to predict tubular secretion or reabsorption clearance pathway using physicochemical properties and structural characteristics.

Renal clearance is one of the main pathways for a drug to be cleared from plasma. The aim of this study is to develop in-silico models to find out the relationship between the type of renal clearance, and structural parameters.Literature data were used to categorise the drugs into those that undergo tubular secretion and those that undergo reabsorption. Different structural descriptors (VolSurf descriptors, Abraham solvation parameters, data warrior descriptors, logarithm of distribution coefficient at pH = 7.4 (logD7.4)) were applied to develop a mechanistic model for estimating renal clearance class whether its secretion or reabsorption.The results of this study show that logD7.4 and the number of hydrogen bond donors, as well as available uncharged species (AUS7.4), are the most effective descriptors to establish mechanistic models for predicting renal clearance class. The classification models were established with a level of accuracy of more than 75%.Developed models of this study can be helpful to predict renal clearance class for new drug candidates with an acceptable error. Hydrophilicity and hydrogen bond formation ability of drugs are among the main descriptors.

A mechanistic perspective, clinical applications, and phage-display-assisted discovery of TNFα inhibitors.

TNFα participates in a variety of physiological processes, but at supra-physiological concentrations it has been implicated in the pathology of inflammatory and autoimmune diseases. Therefore, much attention has been devoted to the development of strategies that overcome the effects of aberrant TNFα concentration. Promising strategies include drugs that destabilize the active (trimeric) form of TNFα and antagonists of TNFα receptor type I. Underpinning these strategies is the successful application of phage-display technology to identify anti-TNFα peptides and antibodies. Here, we review the development of inhibitors of the TNFα-TNF receptor system, with particular focus on the phage-display-assisted identification of molecules that interfere with this system by acting as inhibitors of TNFα or by sequestering TNFα away from its receptor.

Identification of Novel Mutations in Arabidopsis thaliana DOF 4.2 Coding Gene.

Purpose: DOF (DNA-binding with One Finger) proteins are plant-specific transcription factors which mediate numerous biological processes. The purpose of the current study is to report new naturally occurring mutations in the gene encoding for one of the members of DOF proteins named DOF 4.2. Methods: The expression of zinc finger domain of DOF 4.2 (DOF 4.2-ZF) was investigated by first synthesis of cDNA library using different parts of Arabidopsis thaliana plant. Then the coding sequence for zinc finger domain of DOF 4.2 protein was prepared using nested PCR experiment and cloned into pGEX-6P-1 expression vector. Finally, the prepared construct was used for protein expression. Furthermore, molecular dynamics (MD) simulation was carried out to predict DNA binding affinity of DOF 4.2-ZF using AMBER package. Results: For the first time a new variant of DOF 4.2-ZF protein with three mutations was detected. One of the mutations is silent while the other two mutations lead to amino acid replacement (S18G) as well as introduction of a stop codon ultimately resulting in a truncated protein production. In order to investigate whether the truncated form is able to recognize DNA binding motif, MD simulations were carried out and the results showed that the chance of binding of DOF 4.2-ZF protein to cognate DNA in its truncated form is very low. Conclusion: The findings demonstrated that the observed mutations adversely affect the DNA binding ability of the truncated form of DOF4.2 if it is expressed in the mutant variant of A. thaliana used in this study.

Expression, purification and molecular dynamics simulation of extracellular domain of glucagon-like peptide-2 receptor linked to teduglutide.

Teduglutide is the only drug approved for long-term treatment of short bowel syndrome (SBS). This drug exerts its pharmacological effects via binding to the GLP-2 receptors (ECD-GLP2R) located in intestinal tissue. The three dimensional (3D) structure of ECD-GLP2R hasn't been determined yet and hence its mode of interaction with agonists/antagonists is not clear. Therefore, it would be of great importance to develop a structural scaffold for investigation of ECD-GLP2R interactions with its binders. For this, the current study aimed to produce fusion protein of ECD-GLP2R-teduglutide. The ECD-GLP2R-teduglutide protein was expressed in bacterial expression system and purified using affinity and size exclusion chromatography techniques. Using circular dichroism the secondary structure content of purified protein was determined which was comparable to that of theoretical calculations. The low structural stability of purified protein (ΔG = 3.64 kJ.mol-1) was elucidated by monitoring its fluorescence emission at the presence of various concentrations of GdnHCl as a denaturant. Finally, a 3D model for ECD-GLP2R-teduglutide protein was generated and validated using molecular dynamics simulation whose information alongside the experimental studies can be useful for providing new insight into the mode of interaction of ECD-GLP2R with its specific ligands in order to design potent and specific GLP2R agonists.

A quantitative approach to predicting lung deposition profiles of pharmaceutical powder aerosols.

Dry powder inhalers (DPI) are widely used systems for pulmonary delivery of therapeutics. The inhalation performance of DPIs is influenced by formulation features, inhaler device and inhalation pattern. The current review presents the affecting factors with great focus on powder characteristics which include particle size, shape, surface, density, hygroscopicity and crystallinity. The properties of a formulation are greatly influenced by a number of physicochemical factors of drug and added excipients. Since available particle engineering techniques result in particles with a set of modifications, it is difficult to distinguish the effect of an individual feature on powder deposition behavior. This necessitates developing a predictive model capable of describing all influential factors on dry powder inhaler delivery. Therefore, in the current study, a model was constructed to correlate the inhaler device properties, inhalation flow rate, particle characteristics and drug/excipient physicochemical properties with the resultant fine particle fraction. The r2 value of established correlation was 0.74 indicating 86% variability in FPF values is explained by the model with the mean absolute errors of 0.22 for the predicted values. The authors believe that this model is capable of predicting the lung deposition pattern of a formulation with an acceptable precision when the type of inhaler device, inhalation flow rate, physicochemical behavior of active and inactive ingredients and the particle characteristics of DPI formulations are considered.

Synthesis and Biological Evaluation of 1,3,5-Trisubstituted 2-Pyrazolines as Novel Cyclooxygenase-2 Inhibitors with Antiproliferative Activity.

A new series of 1,3,5-trisubstituted 2-pyrazolines for the inhibition of cyclooxygenase-2 (COX-2) were synthesized. The designed structures include a COX-2 pharmacophore SO2 CH3 at the para-position of the phenyl ring located at C-5 of a pyrazoline scaffold. The synthesized compounds were tested for in vitro COX-1/COX-2 inhibition and cell toxicity against human colorectal adenocarcinoma cell lines HT-29. The lead compound (4-chlorophenyl){5-[4-(methanesulfonyl)phenyl]-3-phenyl-4,5-dihydro-1H-pyrazol-1-yl}methanone (16) showed significant COX-2 inhibition (IC50 =0.05±0.01 µM), and antiproliferative activity (IC50 =5.46±4.71 µM). Molecular docking studies showed that new pyrazoline-based compounds interact via multiple hydrophobic and hydrogen-bond interactions with key binding site residues of the COX-2 enzyme.

Co-electrospraying technology as a novel approach for dry powder inhalation formulation of montelukast and budesonide for pulmonary co-delivery.

In the current study electrospraying methodology was used for particle engineering of montelukast and budesonide to prepare a combined inhalable dry powder formulation applicable as a smart regimen in asthma treatment. For this, electrospraying was carried out using different solvents and drug concentrations. No carrier was added for the formulation of montelukast-budesonide combination as montelukast played the role of both active ingredient and carrier. Scanning electron microscopy, particle size analysis, gas chromatography, powder X-ray diffraction, Fourier transform infrared spectroscopy, and differential scanning calorimetry were used to evaluate the physicochemical properties of the produced drug particles. In vitro drug deposition pattern was assessed using next generation impactor, and the dissolution profile of the selected formulations was characterized via modified diffusion franz cell method. The FPF value for the co-electrosprayed carrier free formulation of montelukast-budesonide was 38% with a significantly enhanced dissolution rate for budesonide compared to the budesonide alone formulations. The pharmacological effects of hypothesized combined formulation was assessed by measuring its power to inhibit the production of reactive oxygen species in human normal lung cells. The results showed that the combination of montelukast and budesonide can exert a synergistic effect. The findings in the current study emphasize that using montelukast as a carrier for budesonide not only has greatly improved the aerosolization behavior and dissolution rate of budesonide but also has resulted in synergistic pharmacological effects, indicating the suitability of this combination as an anti-asthmatic therapeutic.

Application of bioinformatics and molecular dynamics simulation approaches for identification of fibroblast growth factor 10 analogues with potentially improved thermostability.

Fibroblast growth factor 10 functions as a paracrine mesenchymal molecule to initiate signalling pathways regarding to cellular development and health. However, the low thermal stability restricts it's functionality in the human body and the shelf-life of FGF10-based formulations. The current study aimed to employ rational design and bioinformatics approaches to identify some point mutations which may improve the thermal stability of FGF10. Bioinformatics analyses resulted in N105D, C106F, K144R, K153M and I156R as the potential stability conferring mutations. The identified mutants were subjected to MD simulation indicating that all mutations are both structurally and energetically favoured. Finally, the effects of the identified mutations on receptor binding of FGF10 were predicted and the results showed that K144R and K153M mutations may increase the binding affinity relative to the wild type. The findings of the current study propose potentially improved FGF10 analogues for further experimental investigations.

Alignment independent 3D-QSAR studies and molecular dynamics simulations for the identification of potent and selective S1P1 receptor agonists.

Sphingosine 1-phosphate type 1 (S1P1) receptors are expressed on lymphocytes and regulate immune cells trafficking. Sphingosine 1-phosphate and its analogues cause internalization and degradation of S1P1 receptors, preventing the auto reactivity of immune cells in the target tissues. It has been shown that S1P1 receptor agonists such as fingolimod can be suitable candidates for treatment of autoimmune diseases. The current study aimed to generate GRIND-based 3D-QSAR predictive models for agonistic activities of 2-imino-thiazolidin-4-one derivatives on S1P1 to be used in virtual screening of chemical libraries. The developed model for the S1P1 receptor agonists showed appropriate power of predictivity in internal (r2acc 0.93 and SDEC 0.18) and external (r2 0.75 and MAE (95% data), 0.28) validations. The generated model revealed the importance of variables DRY-N1 and DRY-O in the potency and selectivity of these compounds towards S1P1 receptor. To propose potential chemical entities with S1P1 agonistic activity, PubChem chemicals database was searched and the selected compounds were virtually tested for S1P1 receptor agonistic activity using the generated models, which resulted in four potential compounds with high potency and selectivity towards S1P1 receptor. Moreover, the affinities of the identified compounds towards S1P1 receptor were evaluated using molecular dynamics simulations. The results indicated that the binding energies of the compounds were in the range of -39.31 to -46.18 and -3.20 to -9.75 kcal mol-1, calculated by MM-GBSA and MM-PBSA algorithms, respectively. The findings in the current work may be useful for the identification of potent and selective S1P1 receptor agonists with potential use in diseases such as multiple sclerosis.