OpenFF Fragmenter
The main purpose of openff-fragmenter
is to fragment molecules for quantum chemical (QC) torsion drives.
Warning
openff-fragmenter
is still pre-alpha. It is not fully tested and the API is still in flux.
Currently two fragmentation schemes are supported:
a Wiberg Bond Order (WBO) sensitive fragmentation scheme [1] (recommended)
the fragmentation schema detailed by Rai et al [2] reffered to in this package as the ‘Pfizer’ scheme.
WBO Sensitive Fragmentation
The assumption when fragmenting molecules is that the chemistry is localized and that removing or changing remote
substituents (defined as substituents more than 2 bonds away from the central bond that is being driven in the torsion
drive) will not change the torsion potential around the bond of interest. However, that is not always the case.
openff-fragmenter
uses the Wiberg Bond Order (WBO) as a surrogate signal to determine if the chemistry around the bond of
interest was destroyed during fragmentation relative to the bond in the parent molecule.
The WBO is a measure of electronic population overlap between two atoms in a bond. It can be quickly calculated from an empirical QC calculation and is given by:
Where \(A\) and \(B\) are atoms \(A\) and \(B\) in a bond, \(D\) is the density matrix and \(\mu\) and \(\nu\) are occupied orbitals on atoms \(A\) and \(B\) respectively.
openff-fragmenter
calculates the WBO of the parent molecules, then fragments according to a set of rules and then
recalculates the WBO of the fragments. If the WBO for the fragment of the bond of interest changes more than a user’s
specified threshold, openff-fragmenter
will add more substituents until the WBO of the bond of interest is within the user
specified threshold.
Contributors
Acknowledgment
CDS is funded by a fellowship from The Molecular Sciences Software Institute
References
- 1
Chaya D Stern, Christopher I Bayly, Daniel G A Smith, Josh Fass, Lee-Ping Wang, David L Mobley, and John D Chodera. Capturing non-local through-bond effects when fragmenting molecules for quantum chemical torsion scans. bioRxiv, 2020. URL: https://www.biorxiv.org/content/early/2020/08/28/2020.08.27.270934, arXiv:https://www.biorxiv.org/content/early/2020/08/28/2020.08.27.270934.full.pdf, doi:10.1101/2020.08.27.270934.
- 2
Brajesh K. Rai, Vishnu Sresht, Qingyi Yang, Ray Unwalla, Meihua Tu, Alan M. Mathiowetz, and Gregory A. Bakken. Comprehensive assessment of torsional strain in crystal structures of small molecules and protein–ligand complexes using ab initio calculations. Journal of Chemical Information and Modeling, 59(10):4195–4208, 2019. PMID: 31573196. URL: https://doi.org/10.1021/acs.jcim.9b00373, arXiv:https://doi.org/10.1021/acs.jcim.9b00373, doi:10.1021/acs.jcim.9b00373.