New Insights into the Binding Features of F508del CFTR Potentiators: A Molecular Docking, Pharmacophore Mapping and QSAR Analysis Approach.
Giada RighettiMonica CasaleMichele TonelliNara LiessiPaola FossaNicoletta PedemonteEnrico MilloElena CicheroPublished in: Pharmaceuticals (Basel, Switzerland) (2020)
Cystic fibrosis (CF) is the autosomal recessive disorder most recurrent in Caucasian populations. To combat this disease, many life-prolonging therapies are required and deeply investigated, including the development of the so-called cystic fibrosis transmembrane conductance regulator (CFTR) modulators, such as correctors and potentiators. Combination therapy with the two series of drugs led to the approval of several multi-drug effective treatments, such as Orkambi, and to the recent promising evaluation of the triple-combination Elexacaftor-Tezacaftor-Ivacaftor. This scenario enlightened the effectiveness of the multi-drug approach to pave the way for the discovery of novel therapeutic agents to contrast CF. The recent X-crystallographic data about the human CFTR in complex with the well-known potentiator Ivacaftor (VX-770) opened the possibility to apply a computational study aimed to explore the key features involved in the potentiator binding. Herein, we discussed molecular docking studies performed onto the chemotypes so far discussed in the literature as CFTR potentiator, reporting the most relevant interactions responsible for their mechanism of action, involving Van der Waals interactions and π-π stacking with F236, Y304, F305 and F312, as well as H-bonding F931, Y304, S308 and R933. This kind of positioning will stabilize the effective potentiator at the CFTR channel. These data have been accompanied by pharmacophore analyses, which promoted the design of novel derivatives endowed with a main (hetero)aromatic core connected to proper substituents, featuring H-bonding moieties. A highly predictive quantitative-structure activity relationship (QSAR) model has been developed, giving a cross-validated r2 (r2cv) = 0.74, a non-cross validated r2 (r2ncv) = 0.90, root mean square error (RMSE) = 0.347, and a test set r2 (r2pred) = 0.86. On the whole, the results are expected to gain useful information to guide the further development and optimization of new CFTR potentiators.
Keyphrases
- cystic fibrosis
- molecular docking
- pseudomonas aeruginosa
- lung function
- molecular dynamics simulations
- combination therapy
- structure activity relationship
- systematic review
- small molecule
- high resolution
- endothelial cells
- electronic health record
- magnetic resonance
- adverse drug
- randomized controlled trial
- big data
- healthcare
- magnetic resonance imaging
- emergency department
- high throughput
- molecular dynamics
- transcription factor
- autism spectrum disorder
- binding protein
- health information
- muscular dystrophy
- amino acid