Activation of AMPD2 drives metabolic dysregulation and liver disease in mice with hereditary fructose intolerance.

Hereditary fructose intolerance (HFI) is a painful and potentially lethal genetic disease caused by a mutation in aldolase B resulting in accumulation of fructose-1-phosphate (F1P). No cure exists for HFI and treatment is limited to avoid exposure to fructose and sugar. Using aldolase B deficient mice, here we identify a yet unrecognized metabolic event activated in HFI and associated with the progression of the disease. Besides the accumulation of F1P, here we show that the activation of the purine degradation pathway is a common feature in aldolase B deficient mice exposed to fructose. The purine degradation pathway is a metabolic route initiated by adenosine monophosphate deaminase 2 (AMPD2) that regulates overall energy balance. We demonstrate that very low amounts of fructose are sufficient to activate AMPD2 in these mice via a phosphate trap. While blocking AMPD2 do not impact F1P accumulation and the risk of hypoglycemia, its deletion in hepatocytes markedly improves the metabolic dysregulation induced by fructose and corrects fat and glycogen storage while significantly increasing the voluntary tolerance of these mice to fructose. In summary, we provide evidence for a critical pathway activated in HFI that could be targeted to improve the metabolic consequences associated with fructose consumption.

Baldwin and Whitehead's Manzamine Alkaloids Biosynthesis Hypothesis Involves a Finely Tuned Reactivity of Acrolein: Automated Extraction of Reactivity Patterns from LC-MS2 Data.

Inspired by the multicomponent reaction-type scenario involving fatty dialdehydes, a nitrogen source, and acrolein, as a key C3 unit, put forward by Baldwin and Whitehead to explain the formation of manzamine-type alkaloids, 96 multicomponent reactions were designed, and their analytical readouts were deconvoluted using a herein-provided chemoinformatic workflow. This strategy pinpointed relevant conditions tuning the reactivity of acrolein to fulfill Baldwin and Whitehead's manzamine alkaloids biosynthetic hypothesis. This strategy can become part of a general method for the high-content analysis of multicomponent reactions applied to a natural product biosynthetic scenario.

Chemodiverse monoterpene indole alkaloids from Kopsia teoi, inhibitory potential against α-amylase, and their molecular docking studies.

Diabetes mellitus stands as a metabolic ailment marked by heightened blood glucose levels due to inadequate insulin secretion. The primary aims of this investigative inquiry encompassed the isolation of phytochemical components from the bark of Kopsia teoi, followed by the assessment of their α-amylase inhibition. The phytochemical composition of the K. teoi culminated in the discovery of a pair of new indole alkaloids; which are 16-epi-deacetylakuammiline N(4)-methylene chloride (akuammiline) (1), and N(1)-methoxycarbonyl-11-methoxy-12-hydroxy-Δ14-17-kopsinine (aspidofractinine) (2), together with five known compounds i.e. kopsiloscine G (aspidofractinine) (3), akuammidine (sarpagine) (4), leuconolam (aspidosperma) (5), N-methoxycarbonyl-12-methoxy-Δ16, 17-kopsinine (aspidofractinine) (6), and kopsininate (aspidofractinine) (7). All compounds were determined via spectroscopic analyses. The in vitro evaluation against α-amylase showed good inhibitory activities for compounds 5-7 with the inhibitory concentration (IC50) values of 21.7 ± 1.2, 34.1 ± 0.1, and 30.0 ± 0.8 µM, respectively compared with the reference acarbose (IC50 = 34.4 ± 0.1 µM). The molecular docking outputs underscored the binding interactions of compounds 5-7 ranging from -8.1 to -8.8 kcal/mol with the binding sites of α-amylase. Consequently, the outcomes highlighted the anti-hyperglycemic attributes of isolates from K. teoi.

Structure-based molecular networking, molecular docking, dynamics simulation and pharmacokinetic studies of Olax subscorpioidea for identification of potential inhibitors against selected cancer targets.

The rationale at the basis of targeted approach in oncology is radically shifting-from development of highly specific agents aiming at a single target towards molecules interfering with multiple targets. This study was performed to isolate and characterize bioactive molecules from Olax subscorpioidea stem and investigate their potentials as multi-targeted inhibitors against selected non-small cell lung cancer, breast cancer and chronic myelogenous leukemia oncogenic targets. Three compounds: ß-sitosterol (1), α-amyrin (2) and stigmasterol (3) were isolated. The structures of 1 - 3 were elucidated by analysis of their spectroscopic data (NMR, MS and IR). To the best of our knowledge, this is the first time these compounds were isolated from O. subscorpioidea stems. Furthermore, integrated analysis of MS/MS data using the Global Natural Products Social Molecular Networking (GNPS) workflow enabled dereplication and identification of 26 compounds, including alkaloids (remerine, boldine), terpenoids (3-hydroxy-11-ursen-28,13-olide, oleanolic acid), flavonoids (kaempferitrin, olax chalcone A) and saponins in O. subscorpioidea stem. Molecular docking studies revealed that some of the compounds, including olax chalcone A (-9.2 to -10.9 kcal/mol), 3-Hydroxy-11-ursen-28,13-olide (-6.6 to -10.2 kcal/mol), α-amyrin (-6.6 to -10.2 kcal/mol), stigmasterol (-7.7 to -10.1 kcal/mol), ß-Sitosterol (-7 to -9.9 kcal/mol) and kaempferitrin (-7.7 to -9 kcal/mol) possessed good inhibitory potentials against selected cancer targets, when compared with reference inhibitors (-8.4 to -13.7 kcal/mol). A few of these compounds were shown to have considerable to favorable pharmacokinetic and drug-likeness properties. This study provides some rationale for the use of O. subscorpioidea in ethnomedicinal management of cancer and identifies some potential anticancer agents.Communicated by Ramaswamy H. Sarma.

Combination of Machine Learning and Empirical Computation for the Structural Validation of Trirosaline, a Natural Trimeric Monoterpene Indole Alkaloid from Catharanthus roseus.

Chemical investigation of the emblematic Catharanthus roseus led to the discovery of trirosaline (1), the first example of a tris-ajmalicine-type monoterpene indole alkaloid and the first natural trimeric MIA ever reported from this deeply dug plant species. Its structure was primarily elucidated based on NMR and HRESIMS analyses, and the nature of its unique intermonomeric linkages was firmly confirmed based on a combination of empirical computation and ML-J-DP4 study. Its absolute configuration was mitigated by comparison of experimental and TDDFT-simulated electronic circular dichroism (ECD) spectra. A possible biosynthetic pathway for trirosaline (1) was postulated.

Reproducible and fully automated testing of nocifensive behavior in mice.

Pain in rodents is often inferred from their withdrawal from noxious stimulation. Threshold stimulus intensity or response latency is used to quantify pain sensitivity. This usually involves applying stimuli by hand and measuring responses by eye, which limits reproducibility and throughput. We describe a device that standardizes and automates pain testing by providing computer-controlled aiming, stimulation, and response measurement. Optogenetic and thermal stimuli are applied using blue and infrared light, respectively. Precise mechanical stimulation is also demonstrated. Reflectance of red light is used to measure paw withdrawal with millisecond precision. We show that consistent stimulus delivery is crucial for resolving stimulus-dependent variations in withdrawal and for testing with sustained stimuli. Moreover, substage video reveals "spontaneous" behaviors for consideration alongside withdrawal metrics to better assess the pain experience. The entire process was automated using machine learning. RAMalgo (reproducible automated multimodal algometry) improves the standardization, comprehensiveness, and throughput of preclinical pain testing.

Artificial intelligence for natural product drug discovery.

Developments in computational omics technologies have provided new means to access the hidden diversity of natural products, unearthing new potential for drug discovery. In parallel, artificial intelligence approaches such as machine learning have led to exciting developments in the computational drug design field, facilitating biological activity prediction and de novo drug design for molecular targets of interest. Here, we describe current and future synergies between these developments to effectively identify drug candidates from the plethora of molecules produced by nature. We also discuss how to address key challenges in realizing the potential of these synergies, such as the need for high-quality datasets to train deep learning algorithms and appropriate strategies for algorithm validation.

Flindissone, a Limonoid Isolated from Trichilia prieuriana, Is an LXR Agonist.

In this study, the ability of six limonoids from Trichilia prieuriana (Meliaceae) to activate the liver X receptor (LXR) was assessed. One of these limonoids, flindissone, was shown to activate LXR by reporter-gene assays. Flindissone is a ring-intact limonoid, structurally similar to sterol-like LXR ligands. In endogenous cellular settings, flindissone showed an activity profile that is characteristic of LXR agonists. It induced cholesterol efflux in THP-1 macrophages by increasing the cholesterol transporter ABCA1 and ABCG1 gene expression. In HepG2 cells, flindissone induced the expression of IDOL, an LXR-target gene that is associated with the downregulation of the LDL receptor. However, unlike synthetic and similarly to sterol-based LXR agonists, flindissone did not induce the expression of the SREBP1c gene, a major transcription factor regulating de novo lipogenesis. Additionally, flindissone also appeared to be able to inhibit post-translational activation of SREBP1c. The results presented here reveal a natural product as a new LXR agonist and point to an additional property of T. prieuriana and other plant extracts containing flindissone.

Behavioral factors and SARS-CoV-2 transmission heterogeneity within a household cohort in Costa Rica.

INTRODUCTION: Variability in household secondary attack rates and transmission risks factors of SARS-CoV-2 remain poorly understood. METHODS: We conducted a household transmission study of SARS-CoV-2 in Costa Rica, with SARS-CoV-2 index cases selected from a larger prospective cohort study and their household contacts were enrolled. A total of 719 household contacts of 304 household index cases were enrolled from November 21, 2020, through July 31, 2021. Blood specimens were collected from contacts within 30-60 days of index case diagnosis; and serum was tested for presence of spike and nucleocapsid SARS-CoV-2 IgG antibodies. Evidence of SARS-CoV-2 prior infections among household contacts was defined based on the presence of both spike and nucleocapsid antibodies. We fitted a chain binomial model to the serologic data, to account for exogenous community infection risk and potential multi-generational transmissions within the household. RESULTS: Overall seroprevalence was 53% (95% confidence interval (CI) 48-58%) among household contacts. The estimated household secondary attack rate is 34% (95% CI 5-75%). Mask wearing by the index case is associated with the household transmission risk reduction by 67% (adjusted odds ratio = 0.33 with 95% CI: 0.09-0.75) and not sharing bedroom with the index case is associated with the risk reduction of household transmission by 78% (adjusted odds ratio = 0.22 with 95% CI 0.10-0.41). The estimated distribution of household secondary attack rates is highly heterogeneous across index cases, with 30% of index cases being the source for 80% of secondary cases. CONCLUSIONS: Modeling analysis suggests that behavioral factors are important drivers of the observed SARS-CoV-2 transmission heterogeneity within the household.

When living in the same house with known SARS-CoV-2 cases, household members may change their behavior and adopt preventive measures to reduce the spread of SARS-CoV-2. To understand how behavioral factors affect SARS-CoV-2 spreading in household settings, we focused on household members of individuals with laboratory-confirmed SARS-CoV-2 infections and followed the way SARS-CoV-2 spread within the household, by looking at who had antibodies against the virus, which means they were infected. We also asked participants detailed questions about their behavior and applied mathematical modeling to evaluate its impact on SARS-CoV-2 transmission. We found that mask-wearing by the SARS-CoV-2 cases, and avoiding sharing a bedroom with the infected individuals, reduces SARS-CoV-2 transmission. However, caring for SARS-CoV-2 cases, and prolonged interaction with infected individuals facilitate SARS-CoV-2 spreading. Our study helps inform what behaviors can help reduce SARS-CoV-2 transmission within a household.

Macrophage Xanthine Oxidoreductase Links LPS Induced Lung Inflammatory Injury to NLRP3 Inflammasome Expression and Mitochondrial Respiration.

Acute lung injury (ALI) and the acute respiratory distress syndrome (ARDS) remain poorly treated inflammatory lung disorders. Both reactive oxygen species (ROS) and macrophages are involved in the pathogenesis of ALI/ARDS. Xanthine oxidoreductase (XOR) is an ROS generator that plays a central role in the inflammation that contributes to ALI. To elucidate the role of macrophage-specific XOR in endotoxin induced ALI, we developed a conditional myeloid specific XOR knockout in mice. Myeloid specific ablation of XOR in LPS insufflated mice markedly attenuated lung injury demonstrating the essential role of XOR in this response. Macrophages from myeloid specific XOR knockout exhibited loss of inflammatory activation and increased expression of anti-inflammatory genes/proteins. Transcriptional profiling of whole lung tissue of LPS insufflated XOR fl/fl//LysM-Cre mice demonstrated an important role for XOR in expression and activation of the NLRP3 inflammasome and acquisition of a glycolytic phenotype by inflammatory macrophages. These results identify XOR as an unexpected link between macrophage redox status, mitochondrial respiration and inflammatory activation.

Discovery of anti-SARS-CoV-2 secondary metabolites from the heartwood of Pterocarpus santalinus using multi-informative molecular networking.

A pigment-depleted extract from the heartwood of Pterocarpus santalinus L. f. (PS-DE) showed promising anti-SARS-CoV-2 activity with an IC50 of 29.9 µg/mL in Caco-2-F03 cells. To determine the potential active constituents within the extract prior to isolation, multi-informative molecular network (MN) was applied. Therefore, the extract was separated by high-performance counter-current chromatography (HPCCC) into 11 fractions which were subsequently tested for anti-SARS-CoV-2 activity and analysed by UPLC-tandem mass spectrometry (MS2). The resulting MN combines the bioactivity data of the fractions with the MS2 data. The MN analysis led to the targeted isolation of seven compounds including one pterocarpan (7) reported for the first time as constituent of P. santalinus and four so far undescribed natural products (NPs) that belong to the compound classes of arylpropanes (9), isoflavanones (10) coumestans (16) and 3-arylcoumarins (17), respectively. In total, 15 constituents from the heartwood of P. santalinus and one synthetic isoflavonoid that is structurally related to the natural metabolites were tested for anti-SARS-CoV-2 activity. Thereby, the two pterocarpans (-)-homopterocarpin (5) and (-)-medicarpin (2), the stilbene (E)-pterostilbene (1) and the isoflavonoid 7-O-methylgenistein (11) showed a distinct antiviral activity with IC50 values of 17.2, 33.4, 34.7, and 37.9 µM, respectively, and no cytotoxic effects against Caco-2-F03 cells (CC50 > 100 µM). In addition, a structure-activity relationship (SAR) was proposed indicating structural requirements of pterocarpans for anti-SARS-CoV-2 activity. The herein presented results support the implementation of multi-informative molecular networks as powerful tool for dereplication and targeted isolation of bioactive NPs.

Unifying the configuration of historical alkaloids from Borreria capitata through an extensive spectroscopic reinvestigation.

Reinvestigation of the structure of borrecapine and borreline through extensive spectroscopic analysis of their authentic samples led to the assignment of their absolute configurations. Newly acquired spectroscopic data determined that the previously assigned relative configuration for borrecapine was incorrect and that the claimed absolute configuration of borreline should be revised to its enantiomer.

From the Spectroscopic Reassessment of Authentic Alkaloid Samples to the Molecular Networking-Guided Discovery of Criophylline-Related Analogues from Callichilia inaequalis.

The molecular network-guided exploration of the alkaloid extract of Callichilia inaequalis stems revealed a cluster attributed tentatively to dimeric monoterpene indole alkaloids of the rare criophylline subtype, initiating the dual study reported herein. A patrimonial-themed portion of this work was aimed at performing a spectroscopic reassessment of criophylline (1), a monoterpene bisindole alkaloid for which the nature of the inter-monomeric connectivity and configurational assignments have remained dubious. A targeted isolation of the entity annotated as criophylline (1) was undertaken to strengthen the available analytical evidence. An extensive set of spectroscopic data was acquired from the authentic sample of criophylline (1a) isolated earlier by Cavé and Bruneton. These spectroscopic studies proved the samples to be identical, and the complete structure of criophylline could be assigned, half a century after it was first isolated. The absolute configuration of andrangine (2) was also ascertained based on a TDDFT-ECD approach from the authentic sample. The forward-looking aspect of this investigation resulted in the characterization of two new criophylline derivatives from C. inaequalis stems, namely, 14'-hydroxycriophylline (3) and 14'-O-sulfocriophylline (4). Their structures, including absolute configurations, were elucidated by analysis of NMR and MS spectroscopic data and by ECD analysis. Notably, 14'-O-sulfocriophylline (4) is the first sulfated monoterpene indole alkaloid to have been reported. The antiplasmodial activity against the chloroquine-resistant strain of Plasmodium falciparum FcB1 was determined for criophylline and its two new analogues.

On the structures of the penduflorines from Tabernaemontana penduliflora.

The structures of the recently published monoterpene indole alkaloids penduflorines A and B (1a and 1b), isolated from Tabernaemontana penduliflora (Apocynaceae), have been revised. Rather than an inseparable mixture of two compounds, they appear to be the known alkaloid vobasine (2). Although we could not comprehensively revise the structures of penduflorines C-E due to lacking spectral data, since their structural elucidations were based on that of 1a and 1b, their structures should also be treated with caution.

Dependence of cortical neuronal strength-duration properties on TMS pulse shape.

OBJECTIVE: The aim of present study was to explore the effects of different combinations of transcranial magnetic stimulation (TMS) pulse width and pulse shape on cortical strength-duration time constant (SDTC) and rheobase measurements. METHODS: Resting motor thresholds (RMT) at pulse widths (PW) of 30, 45, 60, 90 and 120 µs and M-ratios of 0.2, 0.1 and 0.025 were determined using figure-of-eight coil with initial posterior-to-anterior induced current. The M-ratio indicates the relative phases of the induced current with lower values signifying a more unidirectional stimulus. Strength-duration time constant (SDTC) and rheobase were estimated for each M-ratio and various PW combinations. Simulations of biophysically realistic cortical neuron models assessed underlying neuronal populations and physiological mechanisms mediating pulse shape effects on strength-duration properties. RESULTS: The M-ratio exerted significant effect on SDTC (F(2,44) = 4.386, P = 0.021), which was longer for M-ratio of 0.2 (243.4 ± 61.2 µs) compared to 0.025 (186.7 ± 52.5 µs, P = 0.034). Rheobase was significantly smaller when assessed with M-ratio 0.2 compared to 0.025 (P = 0.026). SDTC and rheobase values were most consistent with pulse width sets of 30/45/60/90/120 µs, 30/60/90/120 µs, and 30/60/120 µs. Simulation studies indicated that isolated pyramidal neurons in layers 2/3, 5, and large basket-cells in layer 4 exhibited SDTCs comparable to experimental results. Further, simulation studies indicated that reducing transient Na+ channel conductance increased SDTC with larger increases for higher M-ratios. CONCLUSIONS: Cortical strength-duration curve properties vary with pulse shape, and the modulating effect of the hyperpolarising pulse phase on cortical axonal transient Na+ conductances could account for these changes, although a shift in the recruited neuronal populations may contribute as well. SIGNIFICANCE: The dependence of the cortical strength-duration curve properties on the TMS pulse shape and pulse width selection underscores the need for consistent measurement methods across studies and the potential to extract information about pathophysiological processes.

Cortical inhibition and facilitation are mediated by distinct physiological processes.

A complex interaction of inhibitory and facilitatory interneuronal processes may underlie development of cortical excitability in the human motor cortex. To determine whether distinct interneuronal processes mediated cortical excitability, threshold tracking transcranial magnetic stimulation was utilised to assess cortical excitability, with figure-of-eight coil oriented in posterior-anterior (PA), anterior-posterior (AP) and latero-medial (LM) directions. Motor evoked potential (MEP) responses were recorded over the contralateral abductor pollicis brevis. Resting motor threshold (RMT), short interval intracortical inhibition (SICI), short interval intracortical facilitation (SICF) and intracortical facilitation were recorded. Significant effects of coil orientation were evident on SICI (F = 8.560, P = 0.002) and SICF (F = 7.132, P = 0.003). SICI was greater when recorded with PA (9.7 ± 10.9%, P = 0.029) and AP (13.1 ± 7.0%, P = 0.003) compared to LM (5.2 ± 7.3%) directed currents. SICF was significantly greater with PA (-14.7 ± 8.1%, P = 0.016) and LM (-14.7 ± 8.8%, P = 0.005) compared to AP (-9.1 ± 7.2%) coil orientations. SICI recorded with LM and PA coil orientations were correlated (R = 0.7, P = 0.002), as was SICF recorded with AP vs LM (R = 0.60, P = 0.019) and LM vs PA (R = 0.69, P = 0.002) coil orientations. RMT was significantly smaller with PA compared to AP (P < 0.001) and LM (P = 0.018) stimulation. Recruitment of distinct interneuronal processes with variable cortical orientation and thresholds underlies short interval intracortical inhibition and facilitation.

Ominoxanthone-The First Xanthone Linearly Fused to a γ-Lactone from Cortinarius ominosus Bidaud Basidiomata. CASE- and DFT-Based Structure Elucidation.

The UHPLC-HRMS analysis of Cortinarius ominosus basidiomata extract revealed that this mushroom accumulated elevated yields of an unreported specialized metabolite. The molecular formula of this unknown compound, C17H10O8, indicated that a challenging structure elucidation lay ahead, owing to its critically low H/C atom ratio. The structure of this new isolate, namely ominoxanthone (1), could not be solved from the interpretation of the usual set of 1D/2D NMR data that conveyed too limited information to afford a single, unambiguous structure. To remedy this, a Computer-Assisted Structure Elucidation (CASE) workflow was used to rank the different possible structure candidates consistent with our scarce spectroscopic data. DFT-based chemical shift calculations on a limited set of top-ranked structures further ascertained the determined structure for ominoxanthone. Although the determined scaffold of ominoxanthone is unprecedented as a natural product, a plausible biosynthetic scenario involving a precursor known from cortinariaceous sources and classical biogenetic reactions could be proposed.

Structure elucidation of an aspidofractinine-type monoterpene indole alkaloid from Melodinus reticulatus.

The structure and complete NMR assignments of aspidoreticulofractine, an aspidofractinine N-oxide, are reported. Its structure was elucidated based on a combination of spectroscopic techniques including 1D and 2D NMR, high-resolution mass spectrometry, and electronic circular dichroism.

S-adenosylmethionine blocks tumorigenesis and with immune checkpoint inhibitor enhances anti-cancer efficacy against BRAF mutant and wildtype melanomas.

Despite marked success in treatment with immune checkpoint inhibitor (CPI), only a third of patients are responsive. Thus, melanoma still has one of the highest prevalence and mortality rates; which has led to a search for novel combination therapies that might complement CPI. Aberrant methylomes are one of the mechanisms of resistance to CPI therapy. S-adenosylmethionine (SAM), methyl donor of important epigenetic processes, has significant anti-cancer effects in several malignancies; however, SAM's effect has never been extensively investigated in melanoma. We demonstrate that SAM modulates phenotype switching of melanoma cells and directs the cells towards differentiation indicated by increased melanogenesis (melanin and melanosome synthesis), melanocyte-like morphology, elevated Mitf and Mitf activators' expression, increased antigen expression, reduced proliferation, and reduced stemness genes' expression. Consistently, providing SAM orally, reduced tumor growth and progression, and metastasis of syngeneic BRAF mutant and wild-type (WT) melanoma mouse models. Of note, SAM and anti-PD-1 antibody combination treatment had enhanced anti-cancer efficacy compared to monotherapies, showed significant reduction in tumor growth and progression, and increased survival. Furthermore, SAM and anti-PD-1 antibody combination triggered significantly higher immune cell infiltration, higher CD8+ T cells infiltration and effector functions, and polyfunctionality of CD8+ T cells in YUMMER1.7 tumors. Therefore, SAM combined with CPI provides a novel therapeutic strategy against BRAF mutant and WT melanomas and provides potential to be translated into clinic.

A new vaccine supply chain network under COVID-19 conditions considering system dynamic: Artificial intelligence algorithms.

With the discovery of the COVID-19 vaccine, what has always been worrying the decision-makers is related to the distribution management, the vaccination centers' location, and the inventory control of all types of vaccines. As the COVID-19 vaccine is highly demanded, planning for its fair distribution is a must. University is one of the most densely populated areas in a city, so it is critical to vaccinate university students so that the spread of this virus is curbed. As a result, in the present study, a new stochastic multi-objective, multi-period, and multi-commodity simulation-optimization model has been developed for the COVID-19 vaccine's production, distribution, location, allocation, and inventory control decisions. In this study, the proposed supply chain network includes four echelons of manufacturers, hospitals, vaccination centers, and volunteer vaccine students. Vaccine manufacturers send the vaccines to the vaccination centers and hospitals after production. The students with a history of special diseases such as heart disease, corticosteroids, blood clots, etc. are vaccinated in hospitals because of accessing more medical care, and the rest of the students are vaccinated in the vaccination centers. Then, a system dynamic structure of the prevalence of COVID -19 in universities is developed and the vaccine demand is estimated using simulation, in which the demand enters the mathematical model as a given stochastic parameter. Thus, the model pursues some goals, namely, to minimize supply chain costs, maximize student desirability for vaccination, and maximize justice in vaccine distribution. To solve the proposed model, Variable Neighborhood Search (VNS) and Whale Optimization Algorithm (WOA) algorithms are used. In terms of novelties, the most important novelties in the simulation model are considering the virtual education and exerted quarantine effect on estimating the number of the vaccines. In terms of the mathematical model, one of the remarkable contributions is paying attention to social distancing while receiving the injection and the possibility of the injection during working and non-working hours, and regarding the novelties in the solution methodology, a new heuristic method based on a meta-heuristic algorithm called Modified WOA with VNS (MVWOA) is developed. In terms of the performance metrics and the CPU time, the MOWOA is discovered with a superior performance than other given algorithms. Moreover, regarding the data, a case study related to the COVID-19 pandemic period in Tehran/Iran is provided to validate the proposed algorithm. The outcomes indicate that with the demand increase, the costs increase sharply while the vaccination desirability for students decreases with a slight slope.