@kent-tokyo/chematic 0.2.10 → 0.4.19

package/README.md

@@ -116,6 +116,12 @@ console.log(ifg); // [{"atoms":[1,2,3],"types":"OC=O"}, ...]
 - True MHFP (structural fragment hashing)
 - True ERG (Ertl 2017 functional group detection)
 
+## Bundle Size
+
+~500 KB gzip / ~1.3 MB raw (reduced from ~819 KB gzip in v0.4.17, -38.5%).
+
+PNG rasterization (`tiny_skia`) is excluded from the WASM build — use SVG output instead. All SVG depiction APIs remain fully available.
+
 ## Building from source
 
 ```sh

package/chematic_wasm.d.ts

@@ -169,6 +169,15 @@ export class MolHandle {
      * Number of heavy atoms (explicit atoms in the graph; does not count implicit H).
      */
     atom_count(): number;
+    /**
+     * Returns true when TPSA < 90 Ų, MW < 400, HBD ≤ 3.
+     */
+    bbb_passes(): boolean;
+    /**
+     * Clark (2000) blood-brain barrier logBB score.
+     * logBB > −1.0 = likely CNS penetrant.
+     */
+    bbb_score(): number;
     /**
      * Bertz complexity index (BertzCT).
      */
@@ -177,6 +186,11 @@ export class MolHandle {
      * Number of bonds.
      */
     bond_count(): number;
+    /**
+     * Palm (1997) Caco-2 intestinal permeability (logPCaco2).
+     * > −5.5 = high permeability.
+     */
+    caco2_permeability(): number;
     /**
      * Canonical SMILES string.
      */
@@ -222,8 +236,12 @@ export class MolHandle {
      */
     chi4v(): number;
     /**
-     * 2D PNG depiction (rasterized from SVG).
-     * Returns PNG data as base64-encoded string for embedding in HTML/JS.
+     * CYP3A4 metabolic inhibition risk score (0.0–1.0).
+     */
+    cyp3a4_inhibition_risk(): number;
+    /**
+     * 2D PNG depiction — not available in the WASM build (PNG stack disabled to reduce bundle size).
+     * Use `depict_svg()` in browser contexts; rasterize client-side if needed.
      */
     depict_png(): Uint8Array;
     /**
@@ -272,6 +290,10 @@ export class MolHandle {
      * Number of non-hydrogen heavy atoms.
      */
     heavy_atom_count(): number;
+    /**
+     * hERG cardiac toxicity risk score (0.0–1.0).
+     */
+    herg_risk_score(): number;
     /**
      * Isotope distribution as JSON.
      *
@@ -390,6 +412,14 @@ export class MolHandle {
      * Returns `true` if the molecule has no PAINS structural alerts.
      */
     pains_passes(): boolean;
+    /**
+     * Most acidic pKa in the molecule, or NaN if no acidic site.
+     */
+    pka_acid_value(): number;
+    /**
+     * Most basic pKa in the molecule, or NaN if no basic site.
+     */
+    pka_base_value(): number;
     /**
      * Quantitative Estimate of Drug-likeness (QED); range [0, 1].
      */
@@ -448,6 +478,18 @@ export class MolHandle {
  */
 export function add_hydrogens(mol: MolHandle): MolHandle;
 
+/**
+ * Compute a full ADMET property profile for a molecule.
+ *
+ * Returns a JSON object with fields:
+ * `bbb_score`, `bbb_passes`, `caco2`, `herg_risk`, `cyp3a4_risk`,
+ * `pka_acid` (null if absent), `pka_base` (null if absent),
+ * `esol`, `logd74`, `mw`, `logp`, `tpsa`, `hbd`, `hba`, `rotatable_bonds`
+ *
+ * Returns `{"error":"..."}` on parse failure.
+ */
+export function admet_profile_json(smiles: string): string;
+
 /**
  * AtomPair fingerprint as a bit-packed byte vector (256 bytes = 2048 bits).
  */
@@ -472,6 +514,14 @@ export function autocorr_3d_json(mol: MolHandle): string;
  */
 export function balance_check_json(reaction_smiles: string): string;
 
+/**
+ * Predict GI absorption and BBB penetration using the BOILED-Egg method
+ * (Daina & Zoete 2016).
+ *
+ * Returns JSON: `{"gi_absorbed":bool,"bbb_penetrant":bool,"logp":f64,"tpsa":f64}`
+ */
+export function boiled_egg_json(smiles: string): string;
+
 /**
  * Number of BRICS fragments produced by fragmenting the molecule.
  *
@@ -566,21 +616,10 @@ export function compare_molecules_batch_json(smiles_batch: string, delimiter: st
 export function compare_molecules_json(smiles1: string, smiles2: string): string;
 
 /**
- * Compute direct Coulomb energy for a molecule with Gasteiger partial charges.
- *
- * Returns JSON object: `{ "coulomb_energy": E, "unit": "kcal/mol" }`
- *
- * # Arguments
- * * `mol` - Molecule to evaluate
- *
- * # Example (JavaScript)
- * ```js
- * const mol = parse_smiles("CCO");
- * const result = coulomb_energy_json(mol);
- * // { "coulomb_energy": -12.34, "unit": "kcal/mol" }
- * ```
+ * Generate multiple conformers with RMSD-based pruning.
+ * Returns JSON: `{"conformers": [[[x,y,z],...], ...], "count": int}`.
  */
-export function coulomb_energy_json(mol: MolHandle): string;
+export function conformer_ensemble_json(mol: MolHandle, n: number, rmsd_threshold: number): string;
 
 /**
  * Return the CPK color (CSS hex string) for the given element symbol.
@@ -660,21 +699,6 @@ export function depict_svg_grid_highlighted(smiles_block: string, cols: number,
  */
 export function detect_functional_groups(mol: MolHandle): string;
 
-/**
- * Infer bond connectivity and bond orders from an XYZ-format string.
- *
- * Explicit hydrogen atoms must be present in the XYZ for reliable bond-order
- * assignment (without H, carbonyl C=O cannot be distinguished from C-O).
- *
- * Returns JSON on success: `{"smiles":"CCO","atom_count":3,"bond_count":2}`.
- * `atom_count` and `bond_count` refer to the heavy-atom skeleton (H removed).
- *
- * Returns JSON on error: `{"error":"molecule has 450 atoms; maximum is 300"}`.
- *
- * Safe: never freezes. All internal loops are O(n²). Capped at 300 atoms.
- */
-export function determine_bonds_from_xyz_json(xyz_str: string): string;
-
 /**
  * Dice similarity between `a` and `b` using ECFP4 fingerprints.
  */
@@ -805,6 +829,19 @@ export function find_reaction_center_json(reaction_smiles: string): string;
  */
 export function gasteiger_charges_json(mol: MolHandle): string;
 
+/**
+ * Generate 3D coordinates as raw JSON array [[x,y,z], ...].
+ *
+ * Unlike `generate_3d_pdb`, this returns coordinates that can be passed
+ * to descriptor functions like `whim_descriptors_json` or `shape_descriptors_json`.
+ */
+export function generate_3d_coords_json(mol: MolHandle): string;
+
+/**
+ * Generate 3D coordinates using ETKDG as raw JSON array [[x,y,z], ...].
+ */
+export function generate_3d_etkdg_coords_json(mol: MolHandle): string;
+
 /**
  * Generate 3D coordinates using ETKDG and minimize with DREIDING force field.
  */
@@ -941,10 +978,15 @@ export function get_descriptors_json(mol: MolHandle): string;
 export function get_dihedral_json(smiles: string, a: number, b: number, c: number, d: number): any;
 
 /**
- * Compute GETAWAY descriptors (GEometric, Topologic And wAveleT descriptors) from 3D coordinates.
- * Returns JSON array of 9 values: [G1, G2, G3, D1, D2, D3, T, V, A]
- * where G* = geometric autocorrelations (lag-1,2,3), D* = topologic distances,
- * T = total pairwise distance, V = bounding-box volume, A = anisotropy ratio.
+ * Compute GETAWAY descriptors (GEometry, Topology and Atom-Weights AssemblY) from 3D coords.
+ *
+ * Returns a JSON array of **19** values:
+ * - `[0..7]`  H[1..8]  — leverage autocorrelation at topological lags 1–8
+ * - `[8..15]` R[1..8]  — H[k] normalised by pair count W_k
+ * - `[16]` Hmax, `[17]` Hmean, `[18]` Htot — per-atom leverage statistics
+ *
+ * Note: requires 3D coordinates (non-planar); for flat/2D structures the hat matrix
+ * is degenerate and descriptors reflect squared centroid distances, not true leverage.
  */
 export function getaway_descriptors_json(mol: MolHandle): string;
 
@@ -1036,6 +1078,16 @@ export function maxmin_picks_ecfp4_json(smiles_json: string, n: number): string;
  */
 export function mcs_smiles_json(smiles_json: string): string;
 
+/**
+ * MCS with ring-awareness constraints.
+ *
+ * `smiles_json` — JSON array of at least 2 SMILES strings.
+ * `ring_matches_ring_only` — ring atoms may only match ring atoms.
+ * `complete_rings_only` — partial ring inclusion is removed from the result.
+ * Returns the MCS SMILES, or `"null"` when no common substructure was found.
+ */
+export function mcs_smiles_json_with_ring_config(smiles_json: string, ring_matches_ring_only: boolean, complete_rings_only: boolean): string;
+
 /**
  * MinHash fingerprint (128 hashes) as JSON.
  *
@@ -1071,6 +1123,17 @@ export function minimize_mmff94_json(mol: MolHandle, max_iter: number): string;
  */
 export function minimize_mmff94_lbfgs_json(mol: MolHandle, max_iter: number): string;
 
+/**
+ * Minimise a molecule's geometry using the Universal Force Field (UFF).
+ *
+ * `coords_json` — JSON array of `[x,y,z]` arrays (Å), one per atom.
+ * `max_iter` — maximum iterations (0 = default 500).
+ *
+ * Returns JSON: `{"coords":[[x,y,z],...], "energy":float, "iterations":int, "converged":bool}`
+ * or `{"error":"<msg>"}` on failure.
+ */
+export function minimize_uff_json(smiles: string, coords_json: string, max_iter: number): string;
+
 /**
  * MMFF94 partial charges (BCI table, ±0.1e accuracy) as a JSON array of f64.
  *
@@ -1357,6 +1420,15 @@ export function pharmacophore_fp_2d_summary(mol: MolHandle): string;
  */
 export function pharmacophore_fp_3d_summary(mol: MolHandle): string;
 
+/**
+ * Predict pKa for all ionizable sites in a molecule.
+ *
+ * Returns a JSON array: `[{"atom_idx":8,"pka":4.0,"type":"acid","group":"carboxylic_acid"},...]`
+ *
+ * Returns `[]` if no ionizable sites are found, or `{"error":"..."}` on parse failure.
+ */
+export function predict_pka_json(smiles: string): string;
+
 /**
  * Generate `count` random SMILES from a SMILES string using the given seed.
  * Atoms are permuted based on xorshift64 RNG. Each variant should parse back
@@ -1409,15 +1481,6 @@ export function rgroup_decompose_json(smiles_json: string, core_smarts: string):
  */
 export function ring_families_json(mol: MolHandle): string;
 
-/**
- * Run molecular dynamics simulation and return trajectory as JSON.
- *
- * Returns JSON object with trajectory frames: `{ "frames": [{ "step": N, "potential": E, "kinetic": K, "temp": T }, …] }`
- * Uses NVT ensemble (Berendsen thermostat) at 300 K by default.
- * Note: Limited to molecules with ~50 atoms or fewer for practical WASM performance.
- */
-export function run_md_json(mol: MolHandle, steps: number, temp_k: number): string;
-
 /**
  * Apply a SMIRKS reaction template and return product SMILES as a JSON string.
  *
@@ -1551,6 +1614,17 @@ export function smiles_array_to_sdf(smiles_json: string): string;
  */
 export function smiles_to_mol2(smiles: string): string;
 
+/**
+ * Write a molecule to AutoDock PDBQT format.
+ *
+ * `coords_json` — JSON array of `[x,y,z]` arrays (Å). Pass `"[]"` for zero coords.
+ * `charges_json` — JSON array of partial charges. Pass `"[]"` to write zeros.
+ * `name` — ligand name for the REMARK header.
+ *
+ * Returns the PDBQT string, or `"error:<msg>"` on failure.
+ */
+export function smiles_to_pdbqt(smiles: string, coords_json: string, charges_json: string, name: string): string;
+
 /**
  * Render a highlighted SVG from a SMILES string in one call.
  *
@@ -1688,12 +1762,6 @@ export function to_xyz(mol: MolHandle): string;
  */
 export function torsion_bitvec(mol: MolHandle): Uint8Array;
 
-/**
- * Scan a torsion dihedral i-j-k-l from 0° to 360° in `steps` increments.
- * Returns JSON array: [{"angle":0.0,"energy":E},...] or {"error":"..."}.
- */
-export function torsion_scan_json(mol: MolHandle, i: number, j: number, k: number, l: number, steps: number): string;
-
 /**
  * Virtual screen a query SMILES against a database of SMILES using ECFP4 Tanimoto.
  *
@@ -1707,14 +1775,15 @@ export function virtual_screen_ecfp4_json(query_smi: string, db_smiles_json: str
 
 /**
  * Compute WHIM descriptors (Weighted Holistic Invariant Molecular) from 3D coordinates.
- * Returns JSON array of 10 values: [L1, L2, L3, P1, P2, P3, ALPHA, BETA, GAMMA, DELTA]
- * where L* = inertia tensor eigenvalues, P* = principal moments, ALPHA = sum of moments,
- * BETA = average pairwise interaction, GAMMA = geometric mean, DELTA = anisotropy.
+ * Returns JSON array of 22 values: 11 unit-weight descriptors followed by 11 mass-weight
+ * descriptors. Each 11-element block is [λ₁, λ₂, λ₃, ν₁, ν₂, ν₃, T, A, V, K, D].
  */
 export function whim_descriptors_json(mol: MolHandle): string;
 
 /**
- * Compute combined WHIM + GETAWAY descriptors (19 values total) as JSON array.
+ * Compute combined WHIM + GETAWAY descriptors (**41** values total) as JSON array.
+ *
+ * Returns WHIM[0..21] (22 values) followed by GETAWAY[0..18] (19 values) = 41 total.
  * Useful for ML pipelines requiring both shape and topologic features.
  */
 export function whim_getaway_combined_json(mol: MolHandle): string;
@@ -1728,6 +1797,18 @@ export function whim_getaway_combined_json(mol: MolHandle): string;
  */
 export function write_smiles(mol: MolHandle): string;
 
+/**
+ * XLogP3 partition coefficient (alternative to Crippen LogP).
+ * Returns JSON: `{"xlogp3": float}`.
+ */
+export function xlogp3_json(mol: MolHandle): string;
+
+/**
+ * Per-atom XLogP3 contributions.
+ * Returns JSON array of floats (one per heavy atom).
+ */
+export function xlogp3_per_atom_json(mol: MolHandle): string;
+
 export type InitInput = RequestInfo | URL | Response | BufferSource | WebAssembly.Module;
 
 export interface InitOutput {
@@ -1737,10 +1818,12 @@ export interface InitOutput {
     readonly __wbg_mhfplshhandle_free: (a: number, b: number) => void;
     readonly __wbg_molhandle_free: (a: number, b: number) => void;
     readonly add_hydrogens: (a: number) => number;
+    readonly admet_profile_json: (a: number, b: number) => [number, number];
     readonly atom_pair_bitvec: (a: number) => [number, number];
     readonly autocorr_2d_json: (a: number) => [number, number];
     readonly autocorr_3d_json: (a: number) => [number, number];
     readonly balance_check_json: (a: number, b: number) => [number, number];
+    readonly boiled_egg_json: (a: number, b: number) => [number, number];
     readonly brics_fragment_count: (a: number) => number;
     readonly brics_fragments_json: (a: number) => [number, number];
     readonly butina_cluster_ecfp4_json: (a: number, b: number, c: number) => [number, number, number, number];
@@ -1748,6 +1831,9 @@ export interface InitOutput {
     readonly canonical_tautomer_with_blocked_atoms_json: (a: number, b: number, c: number) => [number, number];
     readonly cdxml_to_smiles_json: (a: number, b: number) => [number, number, number, number];
     readonly cip_assignments_json: (a: number) => [number, number];
+    readonly compare_molecules_batch_json: (a: number, b: number, c: number, d: number) => [number, number, number, number];
+    readonly compare_molecules_json: (a: number, b: number, c: number, d: number) => [number, number, number, number];
+    readonly conformer_ensemble_json: (a: number, b: number, c: number) => [number, number];
     readonly conformerhandle_add_generated_conformer: (a: number) => number;
     readonly conformerhandle_add_minimized_conformer: (a: number) => number;
     readonly conformerhandle_cluster_conformers_json: (a: number, b: number) => [number, number];
@@ -1758,7 +1844,6 @@ export interface InitOutput {
     readonly conformerhandle_mol: (a: number) => number;
     readonly conformerhandle_new: (a: number, b: number) => [number, number, number];
     readonly conformerhandle_remove_conformer: (a: number, b: number) => number;
-    readonly coulomb_energy_json: (a: number) => [number, number];
     readonly cpk_color: (a: number, b: number) => [number, number];
     readonly depict_data_json: (a: number) => [number, number];
     readonly depict_data_with_coords_json: (a: number, b: number, c: number) => [number, number];
@@ -1779,7 +1864,6 @@ export interface InitOutput {
     readonly depictoptions_set_show_atom_indices: (a: number, b: number) => void;
     readonly depictoptions_set_width: (a: number, b: number) => void;
     readonly detect_functional_groups: (a: number) => [number, number];
-    readonly determine_bonds_from_xyz_json: (a: number, b: number) => [number, number];
     readonly dice_ecfp4: (a: number, b: number) => number;
     readonly dice_ecfp6: (a: number, b: number) => number;
     readonly dice_maccs: (a: number, b: number) => number;
@@ -1797,9 +1881,13 @@ export interface InitOutput {
     readonly fcfp6_bitvec: (a: number) => [number, number];
     readonly find_reaction_center_json: (a: number, b: number) => [number, number];
     readonly gasteiger_charges_json: (a: number) => [number, number];
+    readonly generate_3d_coords_json: (a: number) => [number, number];
+    readonly generate_3d_etkdg_coords_json: (a: number) => [number, number];
     readonly generate_3d_etkdg_minimized_pdb: (a: number) => [number, number];
     readonly generate_3d_etkdg_pdb: (a: number) => [number, number];
+    readonly generate_3d_from_smiles: (a: number, b: number) => [number, number, number, number];
     readonly generate_3d_minimized_pdb: (a: number) => [number, number];
+    readonly generate_3d_optimized_pdb: (a: number, b: number) => [number, number, number, number];
     readonly generate_3d_pdb: (a: number) => [number, number];
     readonly generic_murcko_scaffold: (a: number) => number;
     readonly get_atom_info: (a: number, b: number) => [number, number];
@@ -1821,6 +1909,7 @@ export interface InitOutput {
     readonly match_smarts_smiles: (a: number, b: number, c: number, d: number) => [number, number, number, number];
     readonly maxmin_picks_ecfp4_json: (a: number, b: number, c: number) => [number, number, number, number];
     readonly mcs_smiles_json: (a: number, b: number) => [number, number, number, number];
+    readonly mcs_smiles_json_with_ring_config: (a: number, b: number, c: number, d: number) => [number, number, number, number];
     readonly mhfp_hashes_json: (a: number) => [number, number];
     readonly mhfplshhandle_add_smiles: (a: number, b: number, c: number) => [number, number, number];
     readonly mhfplshhandle_is_empty: (a: number) => number;
@@ -1830,6 +1919,7 @@ export interface InitOutput {
     readonly minimize_dreiding_json: (a: number) => [number, number];
     readonly minimize_mmff94_json: (a: number, b: number) => [number, number];
     readonly minimize_mmff94_lbfgs_json: (a: number, b: number) => [number, number];
+    readonly minimize_uff_json: (a: number, b: number, c: number, d: number, e: number) => [number, number];
     readonly mmff94_charges_json: (a: number) => [number, number];
     readonly mmff94_charges_typed_json: (a: number) => [number, number];
     readonly mmff94_energy_breakdown_json: (a: number) => [number, number];
@@ -1851,10 +1941,15 @@ export interface InitOutput {
     readonly mol_with_atom_removed: (a: number, b: number) => [number, number, number];
     readonly mol_with_bond_added: (a: number, b: number, c: number, d: number) => [number, number, number];
     readonly mol_with_bond_removed: (a: number, b: number) => [number, number, number];
+    readonly molecule_report_json: (a: number, b: number) => [number, number, number, number];
     readonly molhandle_aromatic_ring_count: (a: number) => number;
     readonly molhandle_assign_cip_json: (a: number) => [number, number];
+    readonly molhandle_atom_count: (a: number) => number;
+    readonly molhandle_bbb_passes: (a: number) => number;
+    readonly molhandle_bbb_score: (a: number) => number;
     readonly molhandle_bertz_ct: (a: number) => number;
     readonly molhandle_bond_count: (a: number) => number;
+    readonly molhandle_caco2_permeability: (a: number) => number;
     readonly molhandle_canonical_smiles: (a: number) => [number, number];
     readonly molhandle_chi0: (a: number) => number;
     readonly molhandle_chi0v: (a: number) => number;
@@ -1866,6 +1961,7 @@ export interface InitOutput {
     readonly molhandle_chi3v: (a: number) => number;
     readonly molhandle_chi4: (a: number) => number;
     readonly molhandle_chi4v: (a: number) => number;
+    readonly molhandle_cyp3a4_inhibition_risk: (a: number) => number;
     readonly molhandle_depict_png: (a: number) => [number, number];
     readonly molhandle_depict_svg: (a: number) => [number, number];
     readonly molhandle_depict_svg_opts: (a: number, b: number) => [number, number];
@@ -1878,6 +1974,7 @@ export interface InitOutput {
     readonly molhandle_hba_count: (a: number) => number;
     readonly molhandle_hbd_count: (a: number) => number;
     readonly molhandle_heavy_atom_count: (a: number) => number;
+    readonly molhandle_herg_risk_score: (a: number) => number;
     readonly molhandle_isotope_distribution_json: (a: number, b: number) => [number, number];
     readonly molhandle_iupac_name: (a: number) => [number, number];
     readonly molhandle_kappa1: (a: number) => number;
@@ -1904,6 +2001,8 @@ export interface InitOutput {
     readonly molhandle_num_stereocenters: (a: number) => number;
     readonly molhandle_num_unspecified_stereocenters: (a: number) => number;
     readonly molhandle_pains_passes: (a: number) => number;
+    readonly molhandle_pka_acid_value: (a: number) => number;
+    readonly molhandle_pka_base_value: (a: number) => number;
     readonly molhandle_qed: (a: number) => number;
     readonly molhandle_randic_index: (a: number) => number;
     readonly molhandle_reos_passes: (a: number) => number;
@@ -1931,13 +2030,14 @@ export interface InitOutput {
     readonly pharmacophore_features_json: (a: number) => [number, number];
     readonly pharmacophore_fp_2d_summary: (a: number) => [number, number];
     readonly pharmacophore_fp_3d_summary: (a: number) => [number, number];
+    readonly predict_pka_json: (a: number, b: number) => [number, number];
     readonly random_smiles_json: (a: number, b: number, c: number, d: bigint) => [number, number, number, number];
     readonly remove_hydrogens: (a: number) => number;
     readonly rgroup_decompose_json: (a: number, b: number, c: number, d: number) => [number, number, number, number];
     readonly ring_families_json: (a: number) => [number, number, number, number];
-    readonly run_md_json: (a: number, b: number, c: number) => [number, number];
     readonly run_reactants: (a: number, b: number, c: number, d: number) => [number, number, number, number];
     readonly sa_score: (a: number) => number;
+    readonly screen_smiles_json: (a: number, b: number, c: number, d: number) => [number, number];
     readonly sdf_from_records_json: (a: number, b: number, c: number, d: number, e: number, f: number) => [number, number, number, number];
     readonly sdf_to_records_json: (a: number, b: number) => [number, number];
     readonly sdf_to_smiles_json: (a: number, b: number) => [number, number];
@@ -1948,11 +2048,13 @@ export interface InitOutput {
     readonly smarts_match_atoms_with_chirality: (a: number, b: number, c: number, d: number) => [number, number, number, number];
     readonly smiles_array_to_sdf: (a: number, b: number) => [number, number, number, number];
     readonly smiles_to_mol2: (a: number, b: number) => [number, number];
+    readonly smiles_to_pdbqt: (a: number, b: number, c: number, d: number, e: number, f: number, g: number, h: number) => [number, number];
     readonly smiles_to_svg_highlighted: (a: number, b: number, c: number, d: number, e: number, f: number, g: number, h: number) => [number, number, number, number];
     readonly smr_vsa_json: (a: number) => [number, number];
     readonly sssr_rings_json: (a: number) => [number, number];
     readonly standardize_smiles: (a: number, b: number) => [number, number];
     readonly standardize_smiles_report_json: (a: number, b: number, c: number, d: number, e: number, f: number) => [number, number];
+    readonly start: () => void;
     readonly tanimoto_atom_pair: (a: number, b: number) => number;
     readonly tanimoto_ecfp4: (a: number, b: number) => number;
     readonly tanimoto_ecfp6: (a: number, b: number) => number;
@@ -1969,25 +2071,16 @@ export interface InitOutput {
     readonly to_mol_v3000_block: (a: number) => [number, number];
     readonly to_xyz: (a: number) => [number, number];
     readonly torsion_bitvec: (a: number) => [number, number];
-    readonly torsion_scan_json: (a: number, b: number, c: number, d: number, e: number, f: number) => [number, number];
     readonly virtual_screen_ecfp4_json: (a: number, b: number, c: number, d: number, e: number) => [number, number];
     readonly whim_descriptors_json: (a: number) => [number, number];
     readonly whim_getaway_combined_json: (a: number) => [number, number];
     readonly write_smiles: (a: number) => [number, number];
-    readonly start: () => void;
-    readonly molhandle_atom_count: (a: number) => number;
-    readonly compare_molecules_batch_json: (a: number, b: number, c: number, d: number) => [number, number, number, number];
-    readonly compare_molecules_json: (a: number, b: number, c: number, d: number) => [number, number, number, number];
-    readonly generate_3d_from_smiles: (a: number, b: number) => [number, number, number, number];
-    readonly generate_3d_optimized_pdb: (a: number, b: number) => [number, number, number, number];
-    readonly molecule_report_json: (a: number, b: number) => [number, number, number, number];
-    readonly screen_smiles_json: (a: number, b: number, c: number, d: number) => [number, number];
+    readonly xlogp3_json: (a: number) => [number, number];
+    readonly xlogp3_per_atom_json: (a: number) => [number, number];
     readonly __wbindgen_malloc: (a: number, b: number) => number;
     readonly __wbindgen_realloc: (a: number, b: number, c: number, d: number) => number;
-    readonly __wbindgen_free: (a: number, b: number, c: number) => void;
-    readonly __wbindgen_exn_store: (a: number) => void;
-    readonly __externref_table_alloc: () => number;
     readonly __wbindgen_externrefs: WebAssembly.Table;
+    readonly __wbindgen_free: (a: number, b: number, c: number) => void;
     readonly __externref_table_dealloc: (a: number) => void;
     readonly __wbindgen_start: () => void;
 }

package/chematic_wasm.js

@@ -415,6 +415,23 @@ export class MolHandle {
         const ret = wasm.molhandle_atom_count(this.__wbg_ptr);
         return ret >>> 0;
     }
+    /**
+     * Returns true when TPSA < 90 Ų, MW < 400, HBD ≤ 3.
+     * @returns {boolean}
+     */
+    bbb_passes() {
+        const ret = wasm.molhandle_bbb_passes(this.__wbg_ptr);
+        return ret !== 0;
+    }
+    /**
+     * Clark (2000) blood-brain barrier logBB score.
+     * logBB > −1.0 = likely CNS penetrant.
+     * @returns {number}
+     */
+    bbb_score() {
+        const ret = wasm.molhandle_bbb_score(this.__wbg_ptr);
+        return ret;
+    }
     /**
      * Bertz complexity index (BertzCT).
      * @returns {number}
@@ -431,6 +448,15 @@ export class MolHandle {
         const ret = wasm.molhandle_bond_count(this.__wbg_ptr);
         return ret >>> 0;
     }
+    /**
+     * Palm (1997) Caco-2 intestinal permeability (logPCaco2).
+     * > −5.5 = high permeability.
+     * @returns {number}
+     */
+    caco2_permeability() {
+        const ret = wasm.molhandle_caco2_permeability(this.__wbg_ptr);
+        return ret;
+    }
     /**
      * Canonical SMILES string.
      * @returns {string}
@@ -528,8 +554,16 @@ export class MolHandle {
         return ret;
     }
     /**
-     * 2D PNG depiction (rasterized from SVG).
-     * Returns PNG data as base64-encoded string for embedding in HTML/JS.
+     * CYP3A4 metabolic inhibition risk score (0.0–1.0).
+     * @returns {number}
+     */
+    cyp3a4_inhibition_risk() {
+        const ret = wasm.molhandle_cyp3a4_inhibition_risk(this.__wbg_ptr);
+        return ret;
+    }
+    /**
+     * 2D PNG depiction — not available in the WASM build (PNG stack disabled to reduce bundle size).
+     * Use `depict_svg()` in browser contexts; rasterize client-side if needed.
      * @returns {Uint8Array}
      */
     depict_png() {
@@ -654,6 +688,14 @@ export class MolHandle {
         const ret = wasm.molhandle_heavy_atom_count(this.__wbg_ptr);
         return ret >>> 0;
     }
+    /**
+     * hERG cardiac toxicity risk score (0.0–1.0).
+     * @returns {number}
+     */
+    herg_risk_score() {
+        const ret = wasm.molhandle_herg_risk_score(this.__wbg_ptr);
+        return ret;
+    }
     /**
      * Isotope distribution as JSON.
      *
@@ -914,6 +956,22 @@ export class MolHandle {
         const ret = wasm.molhandle_pains_passes(this.__wbg_ptr);
         return ret !== 0;
     }
+    /**
+     * Most acidic pKa in the molecule, or NaN if no acidic site.
+     * @returns {number}
+     */
+    pka_acid_value() {
+        const ret = wasm.molhandle_pka_acid_value(this.__wbg_ptr);
+        return ret;
+    }
+    /**
+     * Most basic pKa in the molecule, or NaN if no basic site.
+     * @returns {number}
+     */
+    pka_base_value() {
+        const ret = wasm.molhandle_pka_base_value(this.__wbg_ptr);
+        return ret;
+    }
     /**
      * Quantitative Estimate of Drug-likeness (QED); range [0, 1].
      * @returns {number}
@@ -1043,6 +1101,33 @@ export function add_hydrogens(mol) {
     return MolHandle.__wrap(ret);
 }
 
+/**
+ * Compute a full ADMET property profile for a molecule.
+ *
+ * Returns a JSON object with fields:
+ * `bbb_score`, `bbb_passes`, `caco2`, `herg_risk`, `cyp3a4_risk`,
+ * `pka_acid` (null if absent), `pka_base` (null if absent),
+ * `esol`, `logd74`, `mw`, `logp`, `tpsa`, `hbd`, `hba`, `rotatable_bonds`
+ *
+ * Returns `{"error":"..."}` on parse failure.
+ * @param {string} smiles
+ * @returns {string}
+ */
+export function admet_profile_json(smiles) {
+    let deferred2_0;
+    let deferred2_1;
+    try {
+        const ptr0 = passStringToWasm0(smiles, wasm.__wbindgen_malloc, wasm.__wbindgen_realloc);
+        const len0 = WASM_VECTOR_LEN;
+        const ret = wasm.admet_profile_json(ptr0, len0);
+        deferred2_0 = ret[0];
+        deferred2_1 = ret[1];
+        return getStringFromWasm0(ret[0], ret[1]);
+    } finally {
+        wasm.__wbindgen_free(deferred2_0, deferred2_1, 1);
+    }
+}
+
 /**
  * AtomPair fingerprint as a bit-packed byte vector (256 bytes = 2048 bits).
  * @param {MolHandle} mol
@@ -1118,6 +1203,29 @@ export function balance_check_json(reaction_smiles) {
     }
 }
 
+/**
+ * Predict GI absorption and BBB penetration using the BOILED-Egg method
+ * (Daina & Zoete 2016).
+ *
+ * Returns JSON: `{"gi_absorbed":bool,"bbb_penetrant":bool,"logp":f64,"tpsa":f64}`
+ * @param {string} smiles
+ * @returns {string}
+ */
+export function boiled_egg_json(smiles) {
+    let deferred2_0;
+    let deferred2_1;
+    try {
+        const ptr0 = passStringToWasm0(smiles, wasm.__wbindgen_malloc, wasm.__wbindgen_realloc);
+        const len0 = WASM_VECTOR_LEN;
+        const ret = wasm.boiled_egg_json(ptr0, len0);
+        deferred2_0 = ret[0];
+        deferred2_1 = ret[1];
+        return getStringFromWasm0(ret[0], ret[1]);
+    } finally {
+        wasm.__wbindgen_free(deferred2_0, deferred2_1, 1);
+    }
+}
+
 /**
  * Number of BRICS fragments produced by fragmenting the molecule.
  *
@@ -1360,28 +1468,19 @@ export function compare_molecules_json(smiles1, smiles2) {
 }
 
 /**
- * Compute direct Coulomb energy for a molecule with Gasteiger partial charges.
- *
- * Returns JSON object: `{ "coulomb_energy": E, "unit": "kcal/mol" }`
- *
- * # Arguments
- * * `mol` - Molecule to evaluate
- *
- * # Example (JavaScript)
- * ```js
- * const mol = parse_smiles("CCO");
- * const result = coulomb_energy_json(mol);
- * // { "coulomb_energy": -12.34, "unit": "kcal/mol" }
- * ```
+ * Generate multiple conformers with RMSD-based pruning.
+ * Returns JSON: `{"conformers": [[[x,y,z],...], ...], "count": int}`.
  * @param {MolHandle} mol
+ * @param {number} n
+ * @param {number} rmsd_threshold
  * @returns {string}
  */
-export function coulomb_energy_json(mol) {
+export function conformer_ensemble_json(mol, n, rmsd_threshold) {
     let deferred1_0;
     let deferred1_1;
     try {
         _assertClass(mol, MolHandle);
-        const ret = wasm.coulomb_energy_json(mol.__wbg_ptr);
+        const ret = wasm.conformer_ensemble_json(mol.__wbg_ptr, n, rmsd_threshold);
         deferred1_0 = ret[0];
         deferred1_1 = ret[1];
         return getStringFromWasm0(ret[0], ret[1]);
@@ -1584,36 +1683,6 @@ export function detect_functional_groups(mol) {
     }
 }
 
-/**
- * Infer bond connectivity and bond orders from an XYZ-format string.
- *
- * Explicit hydrogen atoms must be present in the XYZ for reliable bond-order
- * assignment (without H, carbonyl C=O cannot be distinguished from C-O).
- *
- * Returns JSON on success: `{"smiles":"CCO","atom_count":3,"bond_count":2}`.
- * `atom_count` and `bond_count` refer to the heavy-atom skeleton (H removed).
- *
- * Returns JSON on error: `{"error":"molecule has 450 atoms; maximum is 300"}`.
- *
- * Safe: never freezes. All internal loops are O(n²). Capped at 300 atoms.
- * @param {string} xyz_str
- * @returns {string}
- */
-export function determine_bonds_from_xyz_json(xyz_str) {
-    let deferred2_0;
-    let deferred2_1;
-    try {
-        const ptr0 = passStringToWasm0(xyz_str, wasm.__wbindgen_malloc, wasm.__wbindgen_realloc);
-        const len0 = WASM_VECTOR_LEN;
-        const ret = wasm.determine_bonds_from_xyz_json(ptr0, len0);
-        deferred2_0 = ret[0];
-        deferred2_1 = ret[1];
-        return getStringFromWasm0(ret[0], ret[1]);
-    } finally {
-        wasm.__wbindgen_free(deferred2_0, deferred2_1, 1);
-    }
-}
-
 /**
  * Dice similarity between `a` and `b` using ECFP4 fingerprints.
  * @param {MolHandle} a
@@ -1947,6 +2016,47 @@ export function gasteiger_charges_json(mol) {
     }
 }
 
+/**
+ * Generate 3D coordinates as raw JSON array [[x,y,z], ...].
+ *
+ * Unlike `generate_3d_pdb`, this returns coordinates that can be passed
+ * to descriptor functions like `whim_descriptors_json` or `shape_descriptors_json`.
+ * @param {MolHandle} mol
+ * @returns {string}
+ */
+export function generate_3d_coords_json(mol) {
+    let deferred1_0;
+    let deferred1_1;
+    try {
+        _assertClass(mol, MolHandle);
+        const ret = wasm.generate_3d_coords_json(mol.__wbg_ptr);
+        deferred1_0 = ret[0];
+        deferred1_1 = ret[1];
+        return getStringFromWasm0(ret[0], ret[1]);
+    } finally {
+        wasm.__wbindgen_free(deferred1_0, deferred1_1, 1);
+    }
+}
+
+/**
+ * Generate 3D coordinates using ETKDG as raw JSON array [[x,y,z], ...].
+ * @param {MolHandle} mol
+ * @returns {string}
+ */
+export function generate_3d_etkdg_coords_json(mol) {
+    let deferred1_0;
+    let deferred1_1;
+    try {
+        _assertClass(mol, MolHandle);
+        const ret = wasm.generate_3d_etkdg_coords_json(mol.__wbg_ptr);
+        deferred1_0 = ret[0];
+        deferred1_1 = ret[1];
+        return getStringFromWasm0(ret[0], ret[1]);
+    } finally {
+        wasm.__wbindgen_free(deferred1_0, deferred1_1, 1);
+    }
+}
+
 /**
  * Generate 3D coordinates using ETKDG and minimize with DREIDING force field.
  * @param {MolHandle} mol
@@ -2267,10 +2377,15 @@ export function get_dihedral_json(smiles, a, b, c, d) {
 }
 
 /**
- * Compute GETAWAY descriptors (GEometric, Topologic And wAveleT descriptors) from 3D coordinates.
- * Returns JSON array of 9 values: [G1, G2, G3, D1, D2, D3, T, V, A]
- * where G* = geometric autocorrelations (lag-1,2,3), D* = topologic distances,
- * T = total pairwise distance, V = bounding-box volume, A = anisotropy ratio.
+ * Compute GETAWAY descriptors (GEometry, Topology and Atom-Weights AssemblY) from 3D coords.
+ *
+ * Returns a JSON array of **19** values:
+ * - `[0..7]`  H[1..8]  — leverage autocorrelation at topological lags 1–8
+ * - `[8..15]` R[1..8]  — H[k] normalised by pair count W_k
+ * - `[16]` Hmax, `[17]` Hmean, `[18]` Htot — per-atom leverage statistics
+ *
+ * Note: requires 3D coordinates (non-planar); for flat/2D structures the hat matrix
+ * is degenerate and descriptors reflect squared centroid distances, not true leverage.
  * @param {MolHandle} mol
  * @returns {string}
  */
@@ -2546,6 +2661,39 @@ export function mcs_smiles_json(smiles_json) {
     }
 }
 
+/**
+ * MCS with ring-awareness constraints.
+ *
+ * `smiles_json` — JSON array of at least 2 SMILES strings.
+ * `ring_matches_ring_only` — ring atoms may only match ring atoms.
+ * `complete_rings_only` — partial ring inclusion is removed from the result.
+ * Returns the MCS SMILES, or `"null"` when no common substructure was found.
+ * @param {string} smiles_json
+ * @param {boolean} ring_matches_ring_only
+ * @param {boolean} complete_rings_only
+ * @returns {string}
+ */
+export function mcs_smiles_json_with_ring_config(smiles_json, ring_matches_ring_only, complete_rings_only) {
+    let deferred3_0;
+    let deferred3_1;
+    try {
+        const ptr0 = passStringToWasm0(smiles_json, wasm.__wbindgen_malloc, wasm.__wbindgen_realloc);
+        const len0 = WASM_VECTOR_LEN;
+        const ret = wasm.mcs_smiles_json_with_ring_config(ptr0, len0, ring_matches_ring_only, complete_rings_only);
+        var ptr2 = ret[0];
+        var len2 = ret[1];
+        if (ret[3]) {
+            ptr2 = 0; len2 = 0;
+            throw takeFromExternrefTable0(ret[2]);
+        }
+        deferred3_0 = ptr2;
+        deferred3_1 = len2;
+        return getStringFromWasm0(ptr2, len2);
+    } finally {
+        wasm.__wbindgen_free(deferred3_0, deferred3_1, 1);
+    }
+}
+
 /**
  * MinHash fingerprint (128 hashes) as JSON.
  *
@@ -2639,6 +2787,36 @@ export function minimize_mmff94_lbfgs_json(mol, max_iter) {
     }
 }
 
+/**
+ * Minimise a molecule's geometry using the Universal Force Field (UFF).
+ *
+ * `coords_json` — JSON array of `[x,y,z]` arrays (Å), one per atom.
+ * `max_iter` — maximum iterations (0 = default 500).
+ *
+ * Returns JSON: `{"coords":[[x,y,z],...], "energy":float, "iterations":int, "converged":bool}`
+ * or `{"error":"<msg>"}` on failure.
+ * @param {string} smiles
+ * @param {string} coords_json
+ * @param {number} max_iter
+ * @returns {string}
+ */
+export function minimize_uff_json(smiles, coords_json, max_iter) {
+    let deferred3_0;
+    let deferred3_1;
+    try {
+        const ptr0 = passStringToWasm0(smiles, wasm.__wbindgen_malloc, wasm.__wbindgen_realloc);
+        const len0 = WASM_VECTOR_LEN;
+        const ptr1 = passStringToWasm0(coords_json, wasm.__wbindgen_malloc, wasm.__wbindgen_realloc);
+        const len1 = WASM_VECTOR_LEN;
+        const ret = wasm.minimize_uff_json(ptr0, len0, ptr1, len1, max_iter);
+        deferred3_0 = ret[0];
+        deferred3_1 = ret[1];
+        return getStringFromWasm0(ret[0], ret[1]);
+    } finally {
+        wasm.__wbindgen_free(deferred3_0, deferred3_1, 1);
+    }
+}
+
 /**
  * MMFF94 partial charges (BCI table, ±0.1e accuracy) as a JSON array of f64.
  *
@@ -3434,6 +3612,30 @@ export function pharmacophore_fp_3d_summary(mol) {
     }
 }
 
+/**
+ * Predict pKa for all ionizable sites in a molecule.
+ *
+ * Returns a JSON array: `[{"atom_idx":8,"pka":4.0,"type":"acid","group":"carboxylic_acid"},...]`
+ *
+ * Returns `[]` if no ionizable sites are found, or `{"error":"..."}` on parse failure.
+ * @param {string} smiles
+ * @returns {string}
+ */
+export function predict_pka_json(smiles) {
+    let deferred2_0;
+    let deferred2_1;
+    try {
+        const ptr0 = passStringToWasm0(smiles, wasm.__wbindgen_malloc, wasm.__wbindgen_realloc);
+        const len0 = WASM_VECTOR_LEN;
+        const ret = wasm.predict_pka_json(ptr0, len0);
+        deferred2_0 = ret[0];
+        deferred2_1 = ret[1];
+        return getStringFromWasm0(ret[0], ret[1]);
+    } finally {
+        wasm.__wbindgen_free(deferred2_0, deferred2_1, 1);
+    }
+}
+
 /**
  * Generate `count` random SMILES from a SMILES string using the given seed.
  * Atoms are permuted based on xorshift64 RNG. Each variant should parse back
@@ -3559,31 +3761,6 @@ export function ring_families_json(mol) {
     }
 }
 
-/**
- * Run molecular dynamics simulation and return trajectory as JSON.
- *
- * Returns JSON object with trajectory frames: `{ "frames": [{ "step": N, "potential": E, "kinetic": K, "temp": T }, …] }`
- * Uses NVT ensemble (Berendsen thermostat) at 300 K by default.
- * Note: Limited to molecules with ~50 atoms or fewer for practical WASM performance.
- * @param {MolHandle} mol
- * @param {number} steps
- * @param {number} temp_k
- * @returns {string}
- */
-export function run_md_json(mol, steps, temp_k) {
-    let deferred1_0;
-    let deferred1_1;
-    try {
-        _assertClass(mol, MolHandle);
-        const ret = wasm.run_md_json(mol.__wbg_ptr, steps, temp_k);
-        deferred1_0 = ret[0];
-        deferred1_1 = ret[1];
-        return getStringFromWasm0(ret[0], ret[1]);
-    } finally {
-        wasm.__wbindgen_free(deferred1_0, deferred1_1, 1);
-    }
-}
-
 /**
  * Apply a SMIRKS reaction template and return product SMILES as a JSON string.
  *
@@ -3959,6 +4136,41 @@ export function smiles_to_mol2(smiles) {
     }
 }
 
+/**
+ * Write a molecule to AutoDock PDBQT format.
+ *
+ * `coords_json` — JSON array of `[x,y,z]` arrays (Å). Pass `"[]"` for zero coords.
+ * `charges_json` — JSON array of partial charges. Pass `"[]"` to write zeros.
+ * `name` — ligand name for the REMARK header.
+ *
+ * Returns the PDBQT string, or `"error:<msg>"` on failure.
+ * @param {string} smiles
+ * @param {string} coords_json
+ * @param {string} charges_json
+ * @param {string} name
+ * @returns {string}
+ */
+export function smiles_to_pdbqt(smiles, coords_json, charges_json, name) {
+    let deferred5_0;
+    let deferred5_1;
+    try {
+        const ptr0 = passStringToWasm0(smiles, wasm.__wbindgen_malloc, wasm.__wbindgen_realloc);
+        const len0 = WASM_VECTOR_LEN;
+        const ptr1 = passStringToWasm0(coords_json, wasm.__wbindgen_malloc, wasm.__wbindgen_realloc);
+        const len1 = WASM_VECTOR_LEN;
+        const ptr2 = passStringToWasm0(charges_json, wasm.__wbindgen_malloc, wasm.__wbindgen_realloc);
+        const len2 = WASM_VECTOR_LEN;
+        const ptr3 = passStringToWasm0(name, wasm.__wbindgen_malloc, wasm.__wbindgen_realloc);
+        const len3 = WASM_VECTOR_LEN;
+        const ret = wasm.smiles_to_pdbqt(ptr0, len0, ptr1, len1, ptr2, len2, ptr3, len3);
+        deferred5_0 = ret[0];
+        deferred5_1 = ret[1];
+        return getStringFromWasm0(ret[0], ret[1]);
+    } finally {
+        wasm.__wbindgen_free(deferred5_0, deferred5_1, 1);
+    }
+}
+
 /**
  * Render a highlighted SVG from a SMILES string in one call.
  *
@@ -4362,31 +4574,6 @@ export function torsion_bitvec(mol) {
     return v1;
 }
 
-/**
- * Scan a torsion dihedral i-j-k-l from 0° to 360° in `steps` increments.
- * Returns JSON array: [{"angle":0.0,"energy":E},...] or {"error":"..."}.
- * @param {MolHandle} mol
- * @param {number} i
- * @param {number} j
- * @param {number} k
- * @param {number} l
- * @param {number} steps
- * @returns {string}
- */
-export function torsion_scan_json(mol, i, j, k, l, steps) {
-    let deferred1_0;
-    let deferred1_1;
-    try {
-        _assertClass(mol, MolHandle);
-        const ret = wasm.torsion_scan_json(mol.__wbg_ptr, i, j, k, l, steps);
-        deferred1_0 = ret[0];
-        deferred1_1 = ret[1];
-        return getStringFromWasm0(ret[0], ret[1]);
-    } finally {
-        wasm.__wbindgen_free(deferred1_0, deferred1_1, 1);
-    }
-}
-
 /**
  * Virtual screen a query SMILES against a database of SMILES using ECFP4 Tanimoto.
  *
@@ -4419,9 +4606,8 @@ export function virtual_screen_ecfp4_json(query_smi, db_smiles_json, k) {
 
 /**
  * Compute WHIM descriptors (Weighted Holistic Invariant Molecular) from 3D coordinates.
- * Returns JSON array of 10 values: [L1, L2, L3, P1, P2, P3, ALPHA, BETA, GAMMA, DELTA]
- * where L* = inertia tensor eigenvalues, P* = principal moments, ALPHA = sum of moments,
- * BETA = average pairwise interaction, GAMMA = geometric mean, DELTA = anisotropy.
+ * Returns JSON array of 22 values: 11 unit-weight descriptors followed by 11 mass-weight
+ * descriptors. Each 11-element block is [λ₁, λ₂, λ₃, ν₁, ν₂, ν₃, T, A, V, K, D].
  * @param {MolHandle} mol
  * @returns {string}
  */
@@ -4440,7 +4626,9 @@ export function whim_descriptors_json(mol) {
 }
 
 /**
- * Compute combined WHIM + GETAWAY descriptors (19 values total) as JSON array.
+ * Compute combined WHIM + GETAWAY descriptors (**41** values total) as JSON array.
+ *
+ * Returns WHIM[0..21] (22 values) followed by GETAWAY[0..18] (19 values) = 41 total.
  * Useful for ML pipelines requiring both shape and topologic features.
  * @param {MolHandle} mol
  * @returns {string}
@@ -4481,10 +4669,50 @@ export function write_smiles(mol) {
         wasm.__wbindgen_free(deferred1_0, deferred1_1, 1);
     }
 }
+
+/**
+ * XLogP3 partition coefficient (alternative to Crippen LogP).
+ * Returns JSON: `{"xlogp3": float}`.
+ * @param {MolHandle} mol
+ * @returns {string}
+ */
+export function xlogp3_json(mol) {
+    let deferred1_0;
+    let deferred1_1;
+    try {
+        _assertClass(mol, MolHandle);
+        const ret = wasm.xlogp3_json(mol.__wbg_ptr);
+        deferred1_0 = ret[0];
+        deferred1_1 = ret[1];
+        return getStringFromWasm0(ret[0], ret[1]);
+    } finally {
+        wasm.__wbindgen_free(deferred1_0, deferred1_1, 1);
+    }
+}
+
+/**
+ * Per-atom XLogP3 contributions.
+ * Returns JSON array of floats (one per heavy atom).
+ * @param {MolHandle} mol
+ * @returns {string}
+ */
+export function xlogp3_per_atom_json(mol) {
+    let deferred1_0;
+    let deferred1_1;
+    try {
+        _assertClass(mol, MolHandle);
+        const ret = wasm.xlogp3_per_atom_json(mol.__wbg_ptr);
+        deferred1_0 = ret[0];
+        deferred1_1 = ret[1];
+        return getStringFromWasm0(ret[0], ret[1]);
+    } finally {
+        wasm.__wbindgen_free(deferred1_0, deferred1_1, 1);
+    }
+}
 function __wbg_get_imports() {
     const import0 = {
         __proto__: null,
-        __wbg___wbindgen_string_get_72bdf95d3ae505b1: function(arg0, arg1) {
+        __wbg___wbindgen_string_get_71bb4348194e31f0: function(arg0, arg1) {
             const obj = arg1;
             const ret = typeof(obj) === 'string' ? obj : undefined;
             var ptr1 = isLikeNone(ret) ? 0 : passStringToWasm0(ret, wasm.__wbindgen_malloc, wasm.__wbindgen_realloc);
@@ -4492,33 +4720,11 @@ function __wbg_get_imports() {
             getDataViewMemory0().setInt32(arg0 + 4 * 1, len1, true);
             getDataViewMemory0().setInt32(arg0 + 4 * 0, ptr1, true);
         },
-        __wbg___wbindgen_throw_1506f2235d1bdba0: function(arg0, arg1) {
+        __wbg___wbindgen_throw_ea4887a5f8f9a9db: function(arg0, arg1) {
             throw new Error(getStringFromWasm0(arg0, arg1));
         },
-        __wbg_error_a6fa202b58aa1cd3: function(arg0, arg1) {
-            let deferred0_0;
-            let deferred0_1;
-            try {
-                deferred0_0 = arg0;
-                deferred0_1 = arg1;
-                console.error(getStringFromWasm0(arg0, arg1));
-            } finally {
-                wasm.__wbindgen_free(deferred0_0, deferred0_1, 1);
-            }
-        },
-        __wbg_getRandomValues_3f44b700395062e5: function() { return handleError(function (arg0, arg1) {
-            globalThis.crypto.getRandomValues(getArrayU8FromWasm0(arg0, arg1));
-        }, arguments); },
-        __wbg_new_227d7c05414eb861: function() {
-            const ret = new Error();
-            return ret;
-        },
-        __wbg_stack_3b0d974bbf31e44f: function(arg0, arg1) {
-            const ret = arg1.stack;
-            const ptr1 = passStringToWasm0(ret, wasm.__wbindgen_malloc, wasm.__wbindgen_realloc);
-            const len1 = WASM_VECTOR_LEN;
-            getDataViewMemory0().setInt32(arg0 + 4 * 1, len1, true);
-            getDataViewMemory0().setInt32(arg0 + 4 * 0, ptr1, true);
+        __wbg_error_933f449d72fef598: function(arg0) {
+            console.error(arg0);
         },
         __wbindgen_cast_0000000000000001: function(arg0) {
             // Cast intrinsic for `F64 -> Externref`.
@@ -4559,12 +4765,6 @@ const MolHandleFinalization = (typeof FinalizationRegistry === 'undefined')
     ? { register: () => {}, unregister: () => {} }
     : new FinalizationRegistry(ptr => wasm.__wbg_molhandle_free(ptr, 1));
 
-function addToExternrefTable0(obj) {
-    const idx = wasm.__externref_table_alloc();
-    wasm.__wbindgen_externrefs.set(idx, obj);
-    return idx;
-}
-
 function _assertClass(instance, klass) {
     if (!(instance instanceof klass)) {
         throw new Error(`expected instance of ${klass.name}`);
@@ -4604,15 +4804,6 @@ function getUint8ArrayMemory0() {
     return cachedUint8ArrayMemory0;
 }
 
-function handleError(f, args) {
-    try {
-        return f.apply(this, args);
-    } catch (e) {
-        const idx = addToExternrefTable0(e);
-        wasm.__wbindgen_exn_store(idx);
-    }
-}
-
 function isLikeNone(x) {
     return x === undefined || x === null;
 }

package/package.json

@@ -5,7 +5,7 @@
     "kent-tokyo <kent-tokyo@users.noreply.github.com>"
   ],
   "description": "WebAssembly bindings for chematic — use chematic from JavaScript/TypeScript",
-  "version": "0.2.10",
+  "version": "0.4.19",
   "license": "MIT OR Apache-2.0",
   "repository": {
     "type": "git",
@@ -29,4 +29,4 @@
     "smiles",
     "chemistry"
   ]
-}
+}