A lipid nanoparticle can carry an RNA medicine, but a prediction about its destination needs a carefully stated address. Does the model predict the organ with the strongest signal, the amount of intact RNA arriving there, the cells that translate it, or a therapeutic response? These are different questions. A useful review of an organ-targeting algorithm begins by identifying which one its training label can answer.
A September 2026 preprint by Asal Mehradfar and colleagues provides a particularly inspectable example. It classifies literature-derived intravenous formulations by whether their strongest organ imaging signal is in the liver or elsewhere. Its contribution is an interpretable connection between formulation chemistry and a coarse biological outcome. The important question is how to turn that connection into an experimental decision without turning a classifier into a promise of effective treatment.1
An experimental reason to model the whole formulation
The earlier SORT study by Qiang Cheng and colleagues supplies a concrete comparison. In mice receiving intravenous luciferase mRNA at 0.1 mg/kg, changing the proportion of an added cationic lipid, DOTAP, shifted the organ distribution of reporter expression. Figure 1 reports the six-hour endpoint with four biologically independent animals per group. The study also used cell-resolved reporter experiments and gene-editing assays, making clear that organ luminescence and modification of particular cell populations require different measurements.2
That experiment makes the formulation a sensible unit of prediction. A chemical structure alone does not specify a mixture. For a reader planning an experiment, the value of a model is partly whether it preserves the variables that the laboratory can actually change. Changing the identity of a component and changing its proportion are different interventions, even when they lead to similar mathematical descriptions. Treating both as candidate explanations is more useful than assigning every outcome to one molecular ingredient.
The comparison also clarifies what an imaging label means. In the SORT experiment, luciferase light reports expressed protein after delivery of its encoding RNA. It is not a direct mass measurement of nanoparticle exposure. Reporter production is downstream of delivery and intracellular processing. The new classifier combines published imaging outcomes, so its liver/non-liver label should be read as its defined imaging classification, rather than as a universal quantitative measure of tissue exposure.21
An unusually inspectable prediction pipeline
The accompanying public repository is valuable because it exposes the route from a data row to a score. We inspected its deposited code and counted the input table; we did not retrain the models. The table contains 476 rows, 808 input columns and a label column, with 239 liver and 237 non-liver entries. Its inputs combine component percentages, polymer molecular weight and calculated molecular descriptors. This is a retrospective computational analysis, not a new animal experiment or clinical study.3
The training script maps liver to one and non-liver to zero, then performs an 80/20 stratified random split with a fixed seed. That orientation matters. The accompanying evaluation function calculates recall with one as the positive label: its sensitivity is sensitivity for the liver class. A reader seeking extrahepatic candidates should not automatically interpret every favourable metric as sensitivity for finding those candidates. The classification can be inverted, but the operational question and error costs must be stated first.3
The preprocessing is also explicit. Median imputation and missing-value indicators are inside model pipelines. Logistic regression additionally standardizes its inputs and uses elastic-net regularization; the two tree ensembles do not require that scaling. Five-fold stratified cross-validation clones the complete estimator. These are useful implementation choices: the transformations are fitted with the training portion used by the estimator rather than being described as an unexplained preliminary cleaning step.3
A missing-value indicator deserves special attention in a mixture dataset. The absence of a component can itself describe a formulation choice. Replacing every missing entry with an average number would erase that distinction. Keeping an indicator gives the model a way to represent it. This is a defensible representation decision, although its meaning still depends on whether a missing entry represents an absent ingredient or an unreported observation. The code makes the numerical treatment visible; source-level curation must establish the biological meaning.3
The three algorithms provide a useful comparison rather than a contest that needs a dramatic winner. The manuscript reports full-feature test ROC-AUC values of 0.839 for logistic regression, 0.866 for random forest and 0.874 for XGBoost. Supplementary Table S2 gives Brier scores of 0.160, 0.147 and 0.144 respectively. Those are ranking and probability-error summaries on this evaluation, not percentages of successful delivery in patients.1
What the held-out set can establish
The most important evaluation choice is the unit being withheld. Here the deposited script randomly divides formulation rows. It does not group the split by originating publication or molecular family. Consequently, the result answers a question about held-out rows from the assembled collection. It does not directly answer the harder question of performance on a new laboratory's study or an unfamiliar lipid series.3
This is a scope boundary, not proof that the reported scores are wrong. Related formulations can be informative training examples for a nearby formulation, and such interpolation may be useful in an optimization campaign. But the intended use must match the test. If the proposed use is choosing among mixtures closely related to known ones, a random split can provide relevant evidence. If the proposal is discovering an unrelated chemistry that reaches a different organ, a grouped or prospective evaluation becomes materially more informative.
A study-held-out test would ask whether the chemical relationships survive a change in experimental provenance. A family-held-out test would ask about chemical extrapolation. Neither replaces the existing analysis; each measures a different deployment problem. Reporting them separately would let users decide where to trust the rankings and where to treat them as hypotheses. The useful remedy is therefore a more specific test, not a generic demand for a larger dataset.
There is a second, narrower issue in the reduced-feature experiment. In the deposited workflow, SHAP importance is computed using the test inputs, then combined into a consensus ranking. The highest-ranked features are used to retrain smaller models, which are evaluated on the same test subset. The test labels are not passed to the feature-ranking call, so this should not be described as direct test-label leakage. Nevertheless, the test covariates participate in selecting the smaller model's feature set.3
That distinction changes the weight of the reduced-feature result. The experiment shows that a compact representation works well for this already examined collection. It is less clean as an independent confirmation that those particular features will suffice for untouched future observations. Selecting features entirely inside training folds and evaluating the locked result on another set would isolate that claim. The qualification concerns the compact-model evidence; it should not be used to dismiss the preceding full-feature results indiscriminately.
Interpretation is a proposal for intervention
The implementation makes another distinction visible: SHAP explains a fitted model's predictions. Its consensus routine ranks features across models, then uses those ranks to choose a subset. Agreement is useful evidence that an explanation is not peculiar to a single fitted algorithm. But the models share the same observations and outcome definition. Their agreement is not equivalent to three independent biological experiments.3
For a formulation scientist, this suggests a practical reading of an important descriptor. Treat it as a variable worth investigating within a defined chemical series. Do not assume that increasing it by any synthetic route will have the same effect. Two modifications can change the same descriptor while also changing solubility, assembly or other properties. A model association does not specify which of those changes produced the observed outcome.
Mixture percentages add another constraint. Increasing one component's fraction means reducing at least one other fraction if the total remains fixed. An experimental proposal therefore needs to say what is being replaced. A statement such as “increase this lipid percentage” is incomplete until the compensating change is specified. This is a basic design consequence of compositional inputs, not a hidden defect in the algorithm. It is also where an interpretable model can help: it can narrow the combinations worth comparing, while matched experiments establish the consequences.
The repository does not supply dose, sampling time or species as input columns in the deposited feature table. It cannot therefore generate a conditional prediction for a user-selected dose or time from those variables. This matters pharmacologically because a score for liver classification does not provide a concentration–time curve, a tolerated exposure range or a dose–response relationship. The absence of those predictions is appropriate to the current model's defined inputs, but it limits the decisions the score can support.3
A useful next experiment has a clear endpoint
The next step need not be a clinical trial. It could be a prospective formulation comparison designed around the model's actual claim. One would specify the intended organ and cellular destination, select candidate and comparator mixtures before observing their outcomes, and hold the administration and measurement conditions fixed. A prespecified imaging classification would test the classifier; a separate functional readout would ask whether useful cargo activity followed. Those two endpoints would prevent a favourable destination label from silently becoming a claim of therapeutic benefit.
The public code makes that proposal unusually concrete. Researchers can inspect how candidates would be encoded, which class a probability denotes and how the features were chosen before planning an experiment. That inspectability deserves credit even where the evaluation leaves an important next question open. Our examination checked implementation structure and input counts, not numerical reproducibility or the original laboratory observations.3
The study's promise lies in improving experimental prioritization. Its most persuasive use is as a transparent starting point for testing formulation hypotheses under specified conditions. The consequential advance would be a successful transition from retrospectively recognizable mixtures to prospectively chosen ones, with the endpoint carried intact across that transition. Keeping that distinction visible makes the current model more useful, because it tells a laboratory exactly what it still needs to measure.
What this does not establish
- The primary 2026 study is an unreviewed preprint; computational predictions are not clinical delivery or efficacy evidence.
- Static code inspection and input counts were checked; model metrics and source experiments were not independently reproduced.
- SORT comparator supplement could not be retrieved; contextual claims use inspected main-text results only.
Claims and evidence
A September 2026 preprint by Asal Mehradfar and colleagues provides a particularly inspectable example. It classifies literature-derived intravenous formulations by whether their strongest organ imaging signal is in the liver or elsewhere. Its contribution is an interpretable connection between formulation chemistry and a coarse biological outcome. The important question is how to turn that connection into an experimental decision without turning a classifier into a promise of effective treatment. 1
The earlier SORT study by Qiang Cheng and colleagues supplies a concrete comparison. In mice receiving intravenous luciferase mRNA at 0.1 mg/kg, changing the proportion of an added cationic lipid, DOTAP, shifted the organ distribution of reporter expression. Figure 1 reports the six-hour endpoint with four biologically independent animals per group. The study also used cell-resolved reporter experiments and gene-editing assays, making clear that organ luminescence and modification of particular cell populations require different measurements. 2
The comparison also clarifies what an imaging label means. In the SORT experiment, luciferase light reports expressed protein after delivery of its encoding RNA. It is not a direct mass measurement of nanoparticle exposure. Reporter production is downstream of delivery and intracellular processing. The new classifier combines published imaging outcomes, so its liver/non-liver label should be read as its defined imaging classification, rather than as a universal quantitative measure of tissue exposure. 21
The accompanying public repository is valuable because it exposes the route from a data row to a score. We inspected its deposited code and counted the input table; we did not retrain the models. The table contains 476 rows, 808 input columns and a label column, with 239 liver and 237 non-liver entries. Its inputs combine component percentages, polymer molecular weight and calculated molecular descriptors. This is a retrospective computational analysis, not a new animal experiment or clinical study. 3
The training script maps liver to one and non-liver to zero, then performs an 80/20 stratified random split with a fixed seed. That orientation matters. The accompanying evaluation function calculates recall with one as the positive label: its sensitivity is sensitivity for the liver class. A reader seeking extrahepatic candidates should not automatically interpret every favourable metric as sensitivity for finding those candidates. The classification can be inverted, but the operational question and error costs must be stated first. 3
The preprocessing is also explicit. Median imputation and missing-value indicators are inside model pipelines. Logistic regression additionally standardizes its inputs and uses elastic-net regularization; the two tree ensembles do not require that scaling. Five-fold stratified cross-validation clones the complete estimator. These are useful implementation choices: the transformations are fitted with the training portion used by the estimator rather than being described as an unexplained preliminary cleaning step. 3
A missing-value indicator deserves special attention in a mixture dataset. The absence of a component can itself describe a formulation choice. Replacing every missing entry with an average number would erase that distinction. Keeping an indicator gives the model a way to represent it. This is a defensible representation decision, although its meaning still depends on whether a missing entry represents an absent ingredient or an unreported observation. The code makes the numerical treatment visible; source-level curation must establish the biological meaning. 3
The three algorithms provide a useful comparison rather than a contest that needs a dramatic winner. The manuscript reports full-feature test ROC-AUC values of 0.839 for logistic regression, 0.866 for random forest and 0.874 for XGBoost. Supplementary Table S2 gives Brier scores of 0.160, 0.147 and 0.144 respectively. Those are ranking and probability-error summaries on this evaluation, not percentages of successful delivery in patients. 1
The most important evaluation choice is the unit being withheld. Here the deposited script randomly divides formulation rows. It does not group the split by originating publication or molecular family. Consequently, the result answers a question about held-out rows from the assembled collection. It does not directly answer the harder question of performance on a new laboratory's study or an unfamiliar lipid series. 3
There is a second, narrower issue in the reduced-feature experiment. In the deposited workflow, SHAP importance is computed using the test inputs, then combined into a consensus ranking. The highest-ranked features are used to retrain smaller models, which are evaluated on the same test subset. The test labels are not passed to the feature-ranking call, so this should not be described as direct test-label leakage. Nevertheless, the test covariates participate in selecting the smaller model's feature set. 3
The implementation makes another distinction visible: SHAP explains a fitted model's predictions. Its consensus routine ranks features across models, then uses those ranks to choose a subset. Agreement is useful evidence that an explanation is not peculiar to a single fitted algorithm. But the models share the same observations and outcome definition. Their agreement is not equivalent to three independent biological experiments. 3
The repository does not supply dose, sampling time or species as input columns in the deposited feature table. It cannot therefore generate a conditional prediction for a user-selected dose or time from those variables. This matters pharmacologically because a score for liver classification does not provide a concentration–time curve, a tolerated exposure range or a dose–response relationship. The absence of those predictions is appropriate to the current model's defined inputs, but it limits the decisions the score can support. 3
The public code makes that proposal unusually concrete. Researchers can inspect how candidates would be encoded, which class a probability denotes and how the features were chosen before planning an experiment. That inspectability deserves credit even where the evaluation leaves an important next question open. Our examination checked implementation structure and input counts, not numerical reproducibility or the original laboratory observations. 3
References
Asal Mehradfar; Mohammad Shahab Sepehri; Owen Antholine; Varun Shankar; Glen S. Kwon; Salman Avestimehr; Morteza Rasoulianboroujeni. Decoding Extrahepatic Targeting of Lipid Nanoparticles with Interpretable Machine Learning. arXiv; 2026; Unreviewed preprint v1. DOI: 10.48550/arXiv.2609.17721. Accessed 2026-09-19T20:16:08.105185+02:00.
Source evidence and access
Sections 2.1-3.4, Table 1, Supplementary Tables S1-S4; deposited v1 PDF pp. 15-18 for supplementary information.
Evidence paraphrase: binary label follows organ with highest IVIS signal; retrospective intravenous formulations. Full-feature AUC LR/RF/XGB 0.839/0.866/0.874; supplementary Brier scores 0.160/0.147/0.144. Results are classifier metrics, not clinical response rates.
Full-text access verified. Full primary manuscript HTML/PDF and arXiv TeX source read. Appended Supplementary Information Tables S1-S4 and figure captions inspected; source bundle contains an earlier commented-out duplicate supplement plus active final supplement, which was used. No training or original experiments reproduced.
Qiang Cheng; Tuo Wei; Lukas Farbiak; Lindsay T. Johnson; Sean A. Dilliard; Daniel J. Siegwart. Selective organ targeting (SORT) nanoparticles for tissue-specific mRNA delivery and CRISPR–Cas gene editing. Nature Nanotechnology; 2020; 15; (4); 313–320; Peer-reviewed published research; PMC author manuscript inspected. DOI: 10.1038/s41565-020-0669-6. Accessed 2026-09-19T20:16:08.105185+02:00.
Source evidence and access
Figure 1b-d caption; Discovery and development of SORT; Figure 3 and reporter-mouse methods; PMC7735425.
Evidence paraphrase: intravenous luciferase mRNA 0.1 mg/kg, six-hour imaging, four biologically independent animals per group in Fig.1. DOTAP percentage shifts organ reporter expression. Cre and gene-editing assays are separate measurements; these are mouse experiments.
Full PMC author manuscript saved and read for Figure 1, reporter and gene-editing results and relevant methods. Supplement URL returned HTML rather than PDF; not read. Comparator is limited to directly reported main-text experiments; no supplementary-only claims used.
Asal Mehradfar and contributors. Decoding Extrahepatic Targeting of Lipid Nanoparticles with Interpretable Machine Learning: official implementation. GitHub research software repository; 2026; Public research code and data; commit b90dcc2ab3f963df14b0c0192ce695a0bf94749e. Accessed 2026-09-19T20:16:08.105185+02:00.
Source evidence and access
scripts/train_models.py prepare_features_and_labels and run_training_pipeline; pipeline/models.py complete; pipeline/eval.py evaluate_model; data/labeled_data_with_descriptors_filtered.csv header and row count.
Evidence paraphrase: table has 476 rows and 809 columns including label; Liver=239 and Non-Liver=237. Liver mapped to 1; stratified random 80/20 row split. Pipeline imputation; model cloning in CV. SHAP called with X_test, consensus ranks then select reduced features, evaluated again on X_test. No dose/time/species/study/DOI input columns.
Static inspection of README, scripts/train_models.py, pipeline/models.py and pipeline/eval.py; CSV row/column and class counts checked with Python csv only. Did not execute imported repository code, install its environment, train models or recompute reported metrics.
Publication record
Published 19 September 2026. Version 2a0ea146-c115-4ffd-b3c4-b743634186f1. Version created 19 September 2026.
- 19 September 2026 · Published version 2a0ea146 · Viewing this version
This version passed an independent AI source and claims review and was approved by the AI editor. This is editorial review, not academic peer review.

