How Well Will This Compound Dissolve in Water?
From descriptors computed from its structure, what is this compound's solubility in water? Poorly soluble drug candidates fail late and expensively, so chemists screen for it early.
The table
9,982 compounds with 2D molecular descriptors computed by RDKit and their measured aqueous solubility (log S). It's in Welt's research catalogue as AqSolDB Aqueous Solubility Dataset, 9,982 rows by 25 columns, licensed CC0: Public Domain. The catalogue measures its signal at 0.74. At 9,982 rows and 24 features the table is larger than the model limits allow, so the walkthrough evaluates and fits on a seeded sample of 4,000 rows. TabICLv2 and Kumo Tabular Medium compete.
Walkthrough
Load it from the catalogue
Find the table, read its recorded target and print the column names.
import pandas as pd from welt import Client with Client() as client: item = client.research_datasets(q="AqSolDB aqueous solubility", limit=1)["items"][0] info = client.research_dataset(item["id"]) data = client.load_research_dataset(item["id"]) target = info["targets"][0] target = target["name"] if isinstance(target, dict) else target print([column["name"] for column in info["schema"]], target)Pick the rows
Set aside 20 rows to predict at the end. The models learn from a seeded sample of 4,000 of the rest, the most this table's width allows.
known = data.dropna(subset=[target]) new_rows = known.sample(20, random_state=1) table = known.drop(new_rows.index).sample(4000, random_state=7)Compare the models and predict
The models compete on the same folds against a simple guess, then the top-ranked one is fitted on the whole table. RMSE and MAE are average errors in the target's own units, so lower is better. The last table puts each prediction next to the real answer.
with Client() as client: evaluation = client.evaluate(table, target, task="regression", seed=7, models=["tabicl-v2", "kumo-medium"]) model = evaluation.fit() print(evaluation.summary()[["name", "rank", "rmse_mean", "mae_mean"]]) features = new_rows.drop(columns=target) print(pd.DataFrame({"actual": new_rows[target].values, "predicted": model.predict(features)}))
What you get
Try it
Use Open with to run the research datasets notebook in Colab, Kaggle or Jupyter and search for "AqSolDB aqueous solubility", or run the code above after pip install "welt-client[parquet]".
Open with
- Type
- Use case
- Area
- Science
- Question
- How Well Will This Compound Dissolve in Water?
- Model
- TabICLv2, Kumo Tabular Medium
- Data
- AqSolDB Aqueous Solubility Dataset, 9,982 rows
- For
- Medicinal chemists and materials scientists
- Licence
- CC0: Public Domain