Investigating Polypharmacology Through Targeting Known Human Neutrophil Elastase Inhibitors to Proteinase 3

07 December 2023, Version 1
This content is a preprint and has not undergone peer review at the time of posting.

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.

Keywords

proteinase 3
neutrophil elastase
lambda dynamics
free energy

Supplementary materials

Title
Description
Actions
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

Comments

Comments are not moderated before they are posted, but they can be removed by the site moderators if they are found to be in contravention of our Commenting Policy [opens in a new tab] - please read this policy before you post. Comments should be used for scholarly discussion of the content in question. You can find more information about how to use the commenting feature here [opens in a new tab] .
This site is protected by reCAPTCHA and the Google Privacy Policy [opens in a new tab] and Terms of Service [opens in a new tab] apply.