4101444

Automating bespoke torsion parameterization on Orion ®

Date
August 20, 2024

The accuracy of force fields is vital in several molecular modeling efforts, including studying protein-ligand interactions for drug discovery. Despite having far fewer parameters than traditional force fields, those released by the Open Force Field (OpenFF) Initiative have been shown to accurately reproduce geometries and energies obtained from quantum chemistry (QC) calculations. However, the sensitivity of torsion parameters to local chemical environments in molecules sometimes leads to the underperformance of force fields for molecules with novel chemistry. Bespoke fitting of torsion parameters has been discussed as a solution to this problem. It has been shown that bespoke fitting of torsion parameters can improve the accuracy of property prediction, such as binding free energies.

In this work, we have developed a protocol for bespoke fitting of the torsion parameters for small molecules based on QC calculations. The QC torsion scans are performed on fragments constructed around rotatable bonds that are minimal yet chemically representative, thus reducing the computational cost of generating the reference data. Our effort of bespoke fitting torsion parameters of OpenFF force fields differs from the existing solutions in that we use a hybrid HF-3C/DFT QC protocol and a newly developed efficient parameter optimization procedure. We apply our newly fitted parameters to perform binding free energy calculations of small molecules using Non-Equilibrium switching (NES). We show that bespoke fitting of torsion parameters improves the accuracy of binding free energy predictions obtained using NES for multiple datasets with known binding affinities. Applying bespoke fitting of torsion parameters to improve binding free energy predictions via NES poses a promising pathway for more successful drug discovery campaigns.

Related Products

Thumbnail for Improved ligand-based virtual screening with Flexible-ROCS
Improved ligand-based virtual screening with Flexible-ROCS
Rapid Overlay of Chemical Structures (ROCS) is a molecular alignment software that overlays ligands based on their shape and chemistry. ROCS achieves its remarkable speed by employing a Gaussian representation of molecular volume and performing rigid overlay…
Thumbnail for Constant pH molecular dynamics on the GPUs using enveloping distribution sampling based multistate scheme | Poster Board #445
Constant pH molecular dynamics on the GPUs using enveloping distribution sampling based multistate scheme | Poster Board #445
We introduce a new GPU-accelerated implementation of constant pH simulations, for the computational exploration of protonation states in explicit solvent environments…
Thumbnail for Trans histone pathway of Ubiquitination of nucleosomal H2A by BRCA1/BARD1-UbcH5c Complex | Poster Board #444
Trans histone pathway of Ubiquitination of nucleosomal H2A by BRCA1/BARD1-UbcH5c Complex | Poster Board #444
The Breast Cancer Associated Protein 1 and its binding partner BRCA1 Associated Ring Domain Protein (BARD1) are a heterodimeric protein that, upon mutation, possess a high risk of breast and ovarian cancer…
Thumbnail for ShapeFit: Pose prediction using flexible fitting with shape, color, and force field
ShapeFit: Pose prediction using flexible fitting with shape, color, and force field
ShapeFit is a ligand-centric pose prediction methodology that performs both molecular alignment and scoring, during pose prediction, against the bound-ligand of the input receptor…