Cannabinoid receptor 1 antagonist genistein attenuates marijuana-induced vascular inflammation.

Epidemiological studies reveal that marijuana increases the risk of cardiovascular disease (CVD); however, little is known about the mechanism. Δ9-tetrahydrocannabinol (Δ9-THC), the psychoactive component of marijuana, binds to cannabinoid receptor 1 (CB1/CNR1) in the vasculature and is implicated in CVD. A UK Biobank analysis found that cannabis was an risk factor for CVD. We found that marijuana smoking activated inflammatory cytokines implicated in CVD. In silico virtual screening identified genistein, a soybean isoflavone, as a putative CB1 antagonist. Human-induced pluripotent stem cell-derived endothelial cells were used to model Δ9-THC-induced inflammation and oxidative stress via NF-κB signaling. Knockdown of the CB1 receptor with siRNA, CRISPR interference, and genistein attenuated the effects of Δ9-THC. In mice, genistein blocked Δ9-THC-induced endothelial dysfunction in wire myograph, reduced atherosclerotic plaque, and had minimal penetration of the central nervous system. Genistein is a CB1 antagonist that attenuates Δ9-THC-induced atherosclerosis.

Mapping transcriptomic vector fields of single cells.

Single-cell (sc)RNA-seq, together with RNA velocity and metabolic labeling, reveals cellular states and transitions at unprecedented resolution. Fully exploiting these data, however, requires kinetic models capable of unveiling governing regulatory functions. Here, we introduce an analytical framework dynamo (https://github.com/aristoteleo/dynamo-release), which infers absolute RNA velocity, reconstructs continuous vector fields that predict cell fates, employs differential geometry to extract underlying regulations, and ultimately predicts optimal reprogramming paths and perturbation outcomes. We highlight dynamo's power to overcome fundamental limitations of conventional splicing-based RNA velocity analyses to enable accurate velocity estimations on a metabolically labeled human hematopoiesis scRNA-seq dataset. Furthermore, differential geometry analyses reveal mechanisms driving early megakaryocyte appearance and elucidate asymmetrical regulation within the PU.1-GATA1 circuit. Leveraging the least-action-path method, dynamo accurately predicts drivers of numerous hematopoietic transitions. Finally, in silico perturbations predict cell-fate diversions induced by gene perturbations. Dynamo, thus, represents an important step in advancing quantitative and predictive theories of cell-state transitions.

Meiotic recombination mirrors patterns of germline replication in mice and humans.

Genetic recombination generates novel trait combinations, and understanding how recombination is distributed across the genome is key to modern genetics. The PRDM9 protein defines recombination hotspots; however, megabase-scale recombination patterning is independent of PRDM9. The single round of DNA replication, which precedes recombination in meiosis, may establish these patterns; therefore, we devised an approach to study meiotic replication that includes robust and sensitive mapping of replication origins. We find that meiotic DNA replication is distinct; reduced origin firing slows replication in meiosis, and a distinctive replication pattern in human males underlies the subtelomeric increase in recombination. We detected a robust correlation between replication and both contemporary and historical recombination and found that replication origin density coupled with chromosome size determines the recombination potential of individual chromosomes. Our findings and methods have implications for understanding the mechanisms underlying DNA replication, genetic recombination, and the landscape of mammalian germline variation.

Metabolic modeling of single Th17 cells reveals regulators of autoimmunity.

Metabolism is a major regulator of immune cell function, but it remains difficult to study the metabolic status of individual cells. Here, we present Compass, an algorithm to characterize cellular metabolic states based on single-cell RNA sequencing and flux balance analysis. We applied Compass to associate metabolic states with T helper 17 (Th17) functional variability (pathogenic potential) and recovered a metabolic switch between glycolysis and fatty acid oxidation, akin to known Th17/regulatory T cell (Treg) differences, which we validated by metabolic assays. Compass also predicted that Th17 pathogenicity was associated with arginine and downstream polyamine metabolism. Indeed, polyamine-related enzyme expression was enhanced in pathogenic Th17 and suppressed in Treg cells. Chemical and genetic perturbation of polyamine metabolism inhibited Th17 cytokines, promoted Foxp3 expression, and remodeled the transcriptome and epigenome of Th17 cells toward a Treg-like state. In vivo perturbations of the polyamine pathway altered the phenotype of encephalitogenic T cells and attenuated tissue inflammation in CNS autoimmunity.

Systematic Discovery of Short Linear Motifs Decodes Calcineurin Phosphatase Signaling.

Short linear motifs (SLiMs) drive dynamic protein-protein interactions essential for signaling, but sequence degeneracy and low binding affinities make them difficult to identify. We harnessed unbiased systematic approaches for SLiM discovery to elucidate the regulatory network of calcineurin (CN)/PP2B, the Ca2+-activated phosphatase that recognizes LxVP and PxIxIT motifs. In vitro proteome-wide detection of CN-binding peptides, in vivo SLiM-dependent proximity labeling, and in silico modeling of motif determinants uncovered unanticipated CN interactors, including NOTCH1, which we establish as a CN substrate. Unexpectedly, CN shows SLiM-dependent proximity to centrosomal and nuclear pore complex (NPC) proteins-structures where Ca2+ signaling is largely uncharacterized. CN dephosphorylates human and yeast NPC proteins and promotes accumulation of a nuclear transport reporter, suggesting conserved NPC regulation by CN. The CN network assembled here provides a resource to investigate Ca2+ and CN signaling and demonstrates synergy between experimental and computational methods, establishing a blueprint for examining SLiM-based networks.

Viral T-cell epitopes - Identification, characterization and clinical application.

T-cell immunity, mediated by CD4+ and CD8+ T cells, represents a cornerstone in the control of viral infections. Virus-derived T-cell epitopes are represented by human leukocyte antigen (HLA)-presented viral peptides on the surface of virus-infected cells. They are the prerequisite for the recognition of infected cells by T cells. Knowledge of viral T-cell epitopes provides on the one hand a diagnostic tool to decipher protective T-cell immune responses in the human population and on the other hand various prophylactic and therapeutic options including vaccination approaches and the transfer of virus-specific T cells. Such approaches have already been proven to be effective against various viral infections, particularly in immunocompromised patients lacking sufficient humoral, antibody-based immune response. This review provides an overview on the state of the art as well as current studies regarding the identification and characterization of viral T-cell epitopes and approaches of clinical application. In the first chapter in silico prediction tools and direct, mass spectrometry-based identification of viral T-cell epitopes is compared. The second chapter provides an overview of commonly used assays for further characterization of T-cell responses and phenotypes. The final chapter presents an overview of clinical application of viral T-cell epitopes with a focus on human immunodeficiency virus (HIV), human cytomegalovirus (HCMV) and severe acute respiratory syndrome coronavirus type 2 (SARS-CoV-2), being representatives of relevant viruses.

Advances in phage-host interaction prediction: in silico method enhances the development of phage therapies.

Phages can specifically recognize and kill bacteria, which lead to important application value of bacteriophage in bacterial identification and typing, livestock aquaculture and treatment of human bacterial infection. Considering the variety of human-infected bacteria and the continuous discovery of numerous pathogenic bacteria, screening suitable therapeutic phages that are capable of infecting pathogens from massive phage databases has been a principal step in phage therapy design. Experimental methods to identify phage-host interaction (PHI) are time-consuming and expensive; high-throughput computational method to predict PHI is therefore a potential substitute. Here, we systemically review bioinformatic methods for predicting PHI, introduce reference databases and in silico models applied in these methods and highlight the strengths and challenges of current tools. Finally, we discuss the application scope and future research direction of computational prediction methods, which contribute to the performance improvement of prediction models and the development of personalized phage therapy.

Integrating network pharmacology and in silico analysis deciphers Withaferin-A's anti-breast cancer potential via hedgehog pathway and target network interplay.

This study examines the remarkable effectiveness of Withaferin-A (WA), a withanolide obtained from Withania somnifera (Ashwagandha), in encountering the mortiferous breast malignancy, a global peril. The predominant objective is to investigate WA's intrinsic target proteins and hedgehog (Hh) pathway proteins in breast cancer targeting through the application of in silico computational techniques and network pharmacology predictions. The databases and webtools like Swiss target prediction, GeneCards, DisGeNet and Online Mendelian Inheritance in Man were exploited to identify the common target proteins. The culmination of the WA network and protein-protein interaction network were devised using Stitch and String web tools, through which the drug-target network of 30 common proteins was constructed employing Cytoscape-version 3.9. Enrichment analysis was performed by incorporating Gprofiler, Metascape and Cytoscape plugins. David compounded the Gene Ontology and Kyoto Encyclopedia of Genes and Genomes, and enrichment was computed through bioinformatics tools. The 20 pivotal proteins were docked harnessing Glide, Schrodinger Suite 2023-2. The investigation was governed by docking scores and affinity. The shared target proteins underscored the precise Hh and WA network roles with the affirmation enrichment P-value of <0.025. The implications for hedgehog and cancer pathways were profound with enrichment (P < 0.01). Further, the ADMET and drug-likeness assessments assisted the claim. Robust interactions were noticed with docking studies, authenticated through molecular dynamics, molecular mechanics generalized born surface area scores and bonds. The computational investigation emphasized WA's credible anti-breast activity, specifically with Hh proteins, implying stem-cell-level checkpoint restraints. Rigorous testament is imperative through in vitro and in vivo studies.

BEERS2: RNA-Seq simulation through high fidelity in silico modeling.

Simulation of RNA-seq reads is critical in the assessment, comparison, benchmarking and development of bioinformatics tools. Yet the field of RNA-seq simulators has progressed little in the last decade. To address this need we have developed BEERS2, which combines a flexible and highly configurable design with detailed simulation of the entire library preparation and sequencing pipeline. BEERS2 takes input transcripts (typically fully length messenger RNA transcripts with polyA tails) from either customizable input or from CAMPAREE simulated RNA samples. It produces realistic reads of these transcripts as FASTQ, SAM or BAM formats with the SAM or BAM formats containing the true alignment to the reference genome. It also produces true transcript-level quantification values. BEERS2 combines a flexible and highly configurable design with detailed simulation of the entire library preparation and sequencing pipeline and is designed to include the effects of polyA selection and RiboZero for ribosomal depletion, hexamer priming sequence biases, GC-content biases in polymerase chain reaction (PCR) amplification, barcode read errors and errors during PCR amplification. These characteristics combine to make BEERS2 the most complete simulation of RNA-seq to date. Finally, we demonstrate the use of BEERS2 by measuring the effect of several settings on the popular Salmon pseudoalignment algorithm.

Synthesis of geometrically realistic and watertight neuronal ultrastructure manifolds for in silico modeling.

Understanding the intracellular dynamics of brain cells entails performing three-dimensional molecular simulations incorporating ultrastructural models that can capture cellular membrane geometries at nanometer scales. While there is an abundance of neuronal morphologies available online, e.g. from NeuroMorpho.Org, converting those fairly abstract point-and-diameter representations into geometrically realistic and simulation-ready, i.e. watertight, manifolds is challenging. Many neuronal mesh reconstruction methods have been proposed; however, their resulting meshes are either biologically unplausible or non-watertight. We present an effective and unconditionally robust method capable of generating geometrically realistic and watertight surface manifolds of spiny cortical neurons from their morphological descriptions. The robustness of our method is assessed based on a mixed dataset of cortical neurons with a wide variety of morphological classes. The implementation is seamlessly extended and applied to synthetic astrocytic morphologies that are also plausibly biological in detail. Resulting meshes are ultimately used to create volumetric meshes with tetrahedral domains to perform scalable in silico reaction-diffusion simulations for revealing cellular structure-function relationships. Availability and implementation: Our method is implemented in NeuroMorphoVis, a neuroscience-specific open source Blender add-on, making it freely accessible for neuroscience researchers.

Cortical cell size regulates root metabolic cost.

It has been hypothesized that vacuolar occupancy in mature root cortical parenchyma cells regulates root metabolic cost and thereby plant fitness under conditions of drought, suboptimal nutrient availability, and increased soil mechanical impedance. However, the mechanistic role of vacuoles in reducing root metabolic cost was unproven. Here we provide evidence to support this hypothesis. We first show that root cortical cell size is determined by both cortical cell diameter and cell length. Significant genotypic variation for both cortical cell diameter (~1.1- to 1.5-fold) and cortical cell length (~ 1.3- to 7-fold) was observed in maize and wheat. GWAS and QTL analyses indicate cortical cell diameter and length are heritable and under independent genetic control. We identify candidate genes for both phenes. Empirical results from isophenic lines contrasting for cortical cell diameter and length show that increased cell size, due to either diameter or length, is associated with reduced root respiration, nitrogen content, and phosphorus content. RootSlice, a functional-structural model of root anatomy, predicts that an increased vacuolar: cytoplasmic ratio per unit cortical volume causes reduced root respiration and nutrient content. Ultrastructural imaging of cortical parenchyma cells with varying cortical diameter and cortical cell length confirms the in silico predictions and shows that an increase in cell size is correlated with increased vacuolar volume and reduced cytoplasmic volume. Vacuolar occupancy and its relationship with cell size merits further investigation as a phene for improving crop adaptation to edaphic stress.

Calibration of computational tools for missense variant pathogenicity classification and ClinGen recommendations for PP3/BP4 criteria.

Recommendations from the American College of Medical Genetics and Genomics and the Association for Molecular Pathology (ACMG/AMP) for interpreting sequence variants specify the use of computational predictors as "supporting" level of evidence for pathogenicity or benignity using criteria PP3 and BP4, respectively. However, score intervals defined by tool developers, and ACMG/AMP recommendations that require the consensus of multiple predictors, lack quantitative support. Previously, we described a probabilistic framework that quantified the strengths of evidence (supporting, moderate, strong, very strong) within ACMG/AMP recommendations. We have extended this framework to computational predictors and introduce a new standard that converts a tool's scores to PP3 and BP4 evidence strengths. Our approach is based on estimating the local positive predictive value and can calibrate any computational tool or other continuous-scale evidence on any variant type. We estimate thresholds (score intervals) corresponding to each strength of evidence for pathogenicity and benignity for thirteen missense variant interpretation tools, using carefully assembled independent data sets. Most tools achieved supporting evidence level for both pathogenic and benign classification using newly established thresholds. Multiple tools reached score thresholds justifying moderate and several reached strong evidence levels. One tool reached very strong evidence level for benign classification on some variants. Based on these findings, we provide recommendations for evidence-based revisions of the PP3 and BP4 ACMG/AMP criteria using individual tools and future assessment of computational methods for clinical interpretation.

Comprehensive assessment of nine target prediction web services: which should we choose for target fishing?

Identification of potential targets for known bioactive compounds and novel synthetic analogs is of considerable significance. In silico target fishing (TF) has become an alternative strategy because of the expensive and laborious wet-lab experiments, explosive growth of bioactivity data and rapid development of high-throughput technologies. However, these TF methods are based on different algorithms, molecular representations and training datasets, which may lead to different results when predicting the same query molecules. This can be confusing for practitioners in practical applications. Therefore, this study systematically evaluated nine popular ligand-based TF methods based on target and ligand-target pair statistical strategies, which will help practitioners make choices among multiple TF methods. The evaluation results showed that SwissTargetPrediction was the best method to produce the most reliable predictions while enriching more targets. High-recall similarity ensemble approach (SEA) was able to find real targets for more compounds compared with other TF methods. Therefore, SwissTargetPrediction and SEA can be considered as primary selection methods in future studies. In addition, the results showed that k = 5 was the optimal number of experimental candidate targets. Finally, a novel ensemble TF method based on consensus voting is proposed to improve the prediction performance. The precision of the ensemble TF method outperforms the individual TF method, indicating that the ensemble TF method can more effectively identify real targets within a given top-k threshold. The results of this study can be used as a reference to guide practitioners in selecting the most effective methods in computational drug discovery.

Bioinformatics toolbox for exploring target mutation-induced drug resistance.

Drug resistance is increasingly among the main issues affecting human health and threatening agriculture and food security. In particular, developing approaches to overcome target mutation-induced drug resistance has long been an essential part of biological research. During the past decade, many bioinformatics tools have been developed to explore this type of drug resistance, and they have become popular for elucidating drug resistance mechanisms in a low cost, fast and effective way. However, these resources are scattered and underutilized, and their strengths and limitations have not been systematically analyzed and compared. Here, we systematically surveyed 59 freely available bioinformatics tools for exploring target mutation-induced drug resistance. We analyzed and summarized these resources based on their functionality, data volume, data source, operating principle, performance, etc. And we concisely discussed the strengths, limitations and application examples of these tools. Specifically, we tested some predictive tools and offered some thoughts from the clinician's perspective. Hopefully, this work will provide a useful toolbox for researchers working in the biomedical, pesticide, bioinformatics and pharmaceutical engineering fields, and a good platform for non-specialists to quickly understand drug resistance prediction.

Ultraliser: a framework for creating multiscale, high-fidelity and geometrically realistic 3D models for in silico neuroscience.

Ultraliser is a neuroscience-specific software framework capable of creating accurate and biologically realistic 3D models of complex neuroscientific structures at intracellular (e.g. mitochondria and endoplasmic reticula), cellular (e.g. neurons and glia) and even multicellular scales of resolution (e.g. cerebral vasculature and minicolumns). Resulting models are exported as triangulated surface meshes and annotated volumes for multiple applications in in silico neuroscience, allowing scalable supercomputer simulations that can unravel intricate cellular structure-function relationships. Ultraliser implements a high-performance and unconditionally robust voxelization engine adapted to create optimized watertight surface meshes and annotated voxel grids from arbitrary non-watertight triangular soups, digitized morphological skeletons or binary volumetric masks. The framework represents a major leap forward in simulation-based neuroscience, making it possible to employ high-resolution 3D structural models for quantification of surface areas and volumes, which are of the utmost importance for cellular and system simulations. The power of Ultraliser is demonstrated with several use cases in which hundreds of models are created for potential application in diverse types of simulations. Ultraliser is publicly released under the GNU GPL3 license on GitHub (BlueBrain/Ultraliser). SIGNIFICANCE: There is crystal clear evidence on the impact of cell shape on its signaling mechanisms. Structural models can therefore be insightful to realize the function; the more realistic the structure can be, the further we get insights into the function. Creating realistic structural models from existing ones is challenging, particularly when needed for detailed subcellular simulations. We present Ultraliser, a neuroscience-dedicated framework capable of building these structural models with realistic and detailed cellular geometries that can be used for simulations.

Fast, accurate ranking of engineered proteins by target-binding propensity using structure modeling.

Deep-learning-based methods for protein structure prediction have achieved unprecedented accuracy, yet their utility in the engineering of protein-based binders remains constrained due to a gap between the ability to predict the structures of candidate proteins and the ability toprioritize proteins by their potential to bind to a target. To bridge this gap, we introduce Automated Pairwise Peptide-Receptor Analysis for Screening Engineered proteins (APPRAISE), a method for predicting the target-binding propensity of engineered proteins. After generating structural models of engineered proteins competing for binding to a target using an established structure prediction tool such as AlphaFold-Multimer or ESMFold, APPRAISE performs a rapid (under 1 CPU second per model) scoring analysis that takes into account biophysical and geometrical constraints. As proof-of-concept cases, we demonstrate that APPRAISE can accurately classify receptor-dependent vs. receptor-independent adeno-associated viral vectors and diverse classes of engineered proteins such as miniproteins targeting the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) spike, nanobodies targeting a G-protein-coupled receptor, and peptides that specifically bind to transferrin receptor or programmed death-ligand 1 (PD-L1). APPRAISE is accessible through a web-based notebook interface using Google Colaboratory (https://tiny.cc/APPRAISE). With its accuracy, interpretability, and generalizability, APPRAISE promises to expand the utility of protein structure prediction and accelerate protein engineering for biomedical applications.

Interleukin-2 superkines by computational design.

Affinity maturation of protein­protein interactions is an important approach in the development of therapeutic proteins such as cytokines. Typical experimental strategies involve targeting the cytokine-receptor interface with combinatorial libraries and then selecting for higher-affinity variants. Mutations to the binding scaffold are usually not considered main drivers for improved affinity. Here we demonstrate that computational design can provide affinity-enhanced variants of interleukin-2 (IL-2) "out of the box" without any requirement for interface engineering. Using a strategy of global IL-2 structural stabilization targeting metastable regions of the three-dimensional structure, rather than the receptor binding interfaces, we computationally designed thermostable IL-2 variants with up to 40-fold higher affinity for IL-2Rß without any library-based optimization. These IL-2 analogs exhibited CD25-independent activities on T and natural killer (NK) cells both in vitro and in vivo, mimicking the properties of the IL-2 superkine "super-2" that was engineered through yeast surface display [A. M. Levin et al., Nature, 484, 529­533 (2012)]. Structure-guided stabilization of cytokines is a powerful approach to affinity maturation with applications to many cytokine and protein­protein interactions.

A quantitative systems pharmacology workflow toward optimal design and biomarker stratification of atopic dermatitis clinical trials.

BACKGROUND: The development of atopic dermatitis (AD) drugs is challenged by many disease phenotypes and trial design options, which are hard to explore experimentally. OBJECTIVE: We aimed to optimize AD trial design using simulations. METHODS: We constructed a quantitative systems pharmacology model of AD and standard of care (SoC) treatments and generated a phenotypically diverse virtual population whose parameter distribution was derived from known relationships between AD biomarkers and disease severity and calibrated using disease severity evolution under SoC regimens. RESULTS: We applied this workflow to the immunomodulator OM-85, currently being investigated for its potential use in AD, and calibrated the investigational treatment model with the efficacy profile of an existing trial (thereby enriching it with plausible marker levels and dynamics). We assessed the sensitivity of trial outcomes to trial protocol and found that for this particular example the choice of end point is more important than the choice of dosing regimen and patient selection by model-based responder enrichment could increase the expected effect size. A global sensitivity analysis revealed that only a limited subset of baseline biomarkers is needed to predict the drug response of the full virtual population. CONCLUSIONS: This AD quantitative systems pharmacology workflow built around knowledge of marker-severity relationships as well as SoC efficacy can be tailored to specific development cases to optimize several trial protocol parameters and biomarker stratification and therefore has promise to become a powerful model-informed AD drug development and personalized medicine tool.

Immunoinformatic Approach for Rational Identification of Immunogenic Peptides Against Host Entry and/or Exit Mpox Proteins and Potential Multiepitope Vaccine Construction.

COVID-19 has intensified humanity's concern about the emergence of new pandemics. Since 2018, epidemic outbreaks of the mpox virus have become worrisome. In June 2022, the World Health Organization declared the disease a global health emergency, with 14 500 cases reported by the Centers for Disease Control and Prevention in 60 countries. Therefore, the development of a vaccine based on the current virus genome is paramount in combating new cases. In view of this, we hypothesized the obtainment of rational immunogenic peptides predicted from proteins responsible for entry of the mpox virus into the host (A17L, A26L/A30L, A33R, H2R, L1R), exit (A27L, A35R, A36R, C19L), and both (B5R). To achieve this, we aligned the genome sequencing data of mpox virus isolated from an infected individual in the United States in June 2022 (ON674051.1) with the reference genome dated 2001 (NC_003310.1) for conservation analysis. The Immune Epitope Database server was used for the identification and characterization of the epitopes of each protein related to major histocompatibility complex I or II interaction and recognition by B-cell receptors, resulting in 138 epitopes for A17L, 233 for A28L, 48 for A33R, 77 for H2R, 77 for L1R, 270 for A27L, 72 for A35R, A36R, 148 for C19L, and 276 for B5R. These epitopes were tested in silico for antigenicity, physicochemical properties, and allergenicity, resulting in 51, 40, 10, 34, 38, 57, 25, 7, 47, and 53 epitopes, respectively. Additionally, to select an epitope with the highest promiscuity of binding to major histocompatibility complexes and B-cell receptor simultaneously, all epitopes of each protein were aligned, and the most repetitive and antigenic regions were identified. By classifying the results, we obtained 23 epitopes from the entry proteins, 16 from the exit proteins, and 7 from both. Subsequently, 1 epitope from each protein was selected, and all 3 were fused to construct a chimeric protein that has potential as a multiepitope vaccine. The constructed vaccine was then analyzed for its physicochemical, antigenic, and allergenic properties. Protein modeling, molecular dynamics, and molecular docking were performed on Toll-like receptors 2, 4, and 8, followed by in silico immune simulation of the vaccine. Finally, the results indicate an effective, stable, and safe vaccine that can be further tested, especially in vitro and in vivo, to validate the findings demonstrated in silico.

Enhancement of intestinal calcium transport by short-chain fatty acids: roles of Na+/H+ exchanger 3 and transient receptor potential vanilloid subfamily 6.

Small organic molecules in the intestinal lumen, particularly short-chain fatty acids (SCFAs) and glucose, have long been postulated to enhance calcium absorption. Here, we used 45Ca radioactive tracer to determine calcium fluxes across the rat intestine after exposure to glucose and SCFAs. Confirming previous reports, glucose was found to increase the apical-to-basolateral calcium flux in the cecum. Under apical glucose-free conditions, SCFAs (e.g., butyrate) stimulated the cecal calcium fluxes by approximately twofold, while having no effect on proximal colon. Since SCFAs could be absorbed into the circulation, we further determined whether basolateral SCFA exposure rendered some positive actions. It was found that exposure of duodenum and cecum on the basolateral side to acetate or butyrate increased calcium fluxes. Under butyrate-rich conditions, cecal calcium transport was partially diminished by Na+/H+ exchanger 3 (NHE3) inhibitor (tenapanor) and nonselective transient receptor potential vanilloid subfamily 6 (TRPV6) inhibitor (miconazole). To confirm the contribution of TRPV6 to SCFA-stimulated calcium transport, we synthesized another TRPV6 inhibitor that was demonstrated by in silico molecular docking and molecular dynamics to occlude TRPV6 pore and diminish the glucose- and butyrate-induced calcium fluxes. Therefore, besides corroborating the importance of luminal molecules in calcium absorption, our findings provided foundation for development of more effective calcium-rich nutraceuticals in combination with various absorptive enhancers, e.g., glucose and SCFAs.NEW & NOTEWORTHY Organic molecules in the intestinal lumen, e.g., glucose and short-chain fatty acids (SCFAs), the latter of which are normally produced by microfloral fermentation, can stimulate calcium absorption dependent on transient receptor potential vanilloid subfamily 6 (TRPV6) and Na+/H+ exchanger 3 (NHE3). A selective TRPV6 inhibitor synthesized and demonstrated by in silico docking and molecular dynamics to specifically bind to the pore domain of TRPV6 was used to confirm a significant contribution of this channel. Our findings corroborate physiological significance of nutrients and SCFAs in enhancing calcium absorption.