Abstract
Using a combination of state-of-the-art free energy computational approaches together with in vitro IC50 assays, we evaluated the polypharmacological potential of a scaffold currently in clinical trials for inhibition of human neutrophil elastase (HNE), targeting cardiopulmonary disease, for efficacious inhibition of Proteinase 3 (PR3), a related neutrophil serine proteinase. Based on benchmarking of a series of ligands with known inhibitory activity against HNE using multisite λ-dynamics (MSλD) and multiple force field models, we demonstrated the best agreement between measured IC50 values and computed binding free energies using a hybrid force field of CHARMM c36m for the protein and ligand parameters from OPLS LigParGen that utilizes a CM1A-LBCC charge model. Using this computational model, the ligands were assessed as potential inhibitors of PR3 by anchoring the computed relative binding affinities to a single measured IC50. This led to suggested potential additional inhibitors for PR3 from those tested and characterized to be good binders to HNE. We synthesized and assayed a number of these compounds to further test our predicted binding affinities and discovered that both prediction and experiment demonstrated that the compounds bound less well to PR3, we observed a range of IC50s from 7.7-190 nM in PR3 to 0.5-1.3 nM for the same inhibitors for HNE. However, the affinities we observe suggest that the scaffold can serve as a suitable starting point for the establishment of polypharmacologically targeting both enzymes and enhancing the potential for treatments addressing diseases like chronic obstructive pulmonary disease.
Supplementary materials
Title
Supporting Information text
Description
Supporting information: computational and experimental details including IC50 assays and compounds synthesis. Detailed results and tables reporting predicted relative binding free energy for HNE (Table S2) and PR3 (Tables S4) with the four combinations of force fields, the associated predicted absolute binding free energies (HNE: Tables S3, S6, PR3: Tables S5) and their comparison against experimental data (HNE: Tables S3, S9 , PR3: Tables S8, S10) using different schemes to estimate the shift for minimization of systematic errors as experimental data became available. The predicted IC50 values for PR3 are reported in Table S7.
Simulation files (https://github.com/reuter-group/bayer_compounds_msld): Ligands mol2 files along with original FF parameters from CGenFF, GAFF2 and OPLS; MSλD prep containing hybrid ligands, ALF biases, and CHARMM production run files used for relative binding free energy calculations; CHARMM input files for bookending corrections.
Actions