Stability analysis

from kallima import KallimaClient

client = KallimaClient(api_key)

job = client.stability_analyses.submit(variant_id=variant_id)
job.wait(timeout=300)

result = job.results

print(result["solubility_score"])       # float 0–1, higher = more soluble
print(result["aggregation_score"])      # float 0–1, lower = less aggregation-prone
print(result["thermostability_score"])  # float 0–1
print(result["melting_temp_predicted"]) # predicted Tm in °C, e.g. 72.4
print(result["developability_score"])   # composite 0–1

ddg = result["ddg_matrix"] or {}
for position, ddg_value in ddg.items():
    if ddg_value < -1.0:
        print(f"{position}: stabilizing ({ddg_value:.2f} kcal/mol)")

for mut in result["stabilizing_mutations"] or []:
    print(mut["position"], mut["wild_type"], mut["mutation"], mut["ddg"])

for job in client.stability_analyses.list(variant_id=variant_id):
    print(job.id, job.status, job["created_at"])

client.stability_analyses.cancel(job.id)