Strain Calculation

Python API · stjames models · API example

How it works

Rowan's strain calculation estimates the energy cost of a molecular pose, such as a docked ligand conformation. It gently relaxes the submitted geometry under positional restraints, then compares its energy with the lowest-energy conformer found in a conformer search. Both are scored using the same final method and solvent settings. The reported strain is an energy difference in kcal/mol; the restraint penalty itself is excluded.

Settings

  • Conformer search: open "Conformer search settings" and choose a "Conformer generator": openconf (the default), ETKDG, iMTD, or MCMM. Adjust search effort, energy window, and "Max conformers" where available. A broader search can find a lower reference energy, increasing the estimated strain.
  • Energy method: "Optimization method" controls geometry relaxation and final energy scoring. The default uses GFN2-xTB with ALPB water for optimization, followed by g-xTB with CPCM-X water for scoring. AIMNet2 is another preset; "Custom" lets you edit optimization stages and the final single-point calculation. Retain at least one optimization stage.
  • Solvent: "Solvent (all steps)" applies an implicit solvent to steps that support it. Solvent choice can change the relative stability of conformers and therefore the strain.
  • Restraint strength: "Harmonic constraint spring constant" defaults to 5 kcal/(mol·Å²). Larger values keep atoms closer to the starting pose; smaller values allow more relaxation and may change the conformation substantially.
  • Hydrogens: "Constrain hydrogens?" is off by default. Heavy atoms are restrained, while hydrogens can relax freely. Enable it when hydrogen positions must also stay close to their starting coordinates.

Notes

Positional restraints let bond lengths and angles relax toward the chosen method's preferences while discouraging major changes to the pose. This reduces artificially high energies caused by small geometric differences between, for example, a docking force field and the energy method. Inspect the relaxed structure to check that it still represents the pose you intended to evaluate.

Results include the strain energy, the restrained geometry, and the reference conformers. The reference is the best conformer found, rather than a guaranteed global minimum. If the relaxed pose lies below the search minimum, an additional free optimization helps update the reference.

Strain describes the molecule's internal energetic cost under the selected model. It is not a binding affinity or a conformational free energy. Compare results using consistent protonation states, methods, solvent settings, and restraint strengths.

Further reading