scholarly journals Novel Pyrazolo[3,4-c]pyridine Antagonists with Nanomolar affinity for A1 / A3 Adenosine Receptors: Binding Kinetics and Exploration of their Binding Profile Using Mutagenesis Experiments, MD simulations and TI/MD calculations

Author(s):  
Margarita Stampelou ◽  
Anna Suchankova ◽  
Efpraxia Tzortzini ◽  
Lakshiv Dhingra ◽  
Kerry Barkan ◽  
...  

Drugs targeting the four adenosine receptor (AR) subtypes can provide “soft" treatment of various significant diseases. Even for the two experimentally resolved AR subtypes the description of the orthosteric binding area and structure-activity relationships of ligands remains a demanding task due to the high similar amino acids sequence but also the broadness and flexibility of the ARs binding area. The identification of new pharmacophoric moieties and nanomolar leads and the exploration of their binding area with mutagenesis and state-of-the-art computational methods useful also for drug design purposes remains a challenging aim for all ARs. Here, we identified several low nanomolar ligands and potent competitive antagonists against A1R / A3R, containing the novel pyrazolo[3,4-c]pyridine pharmacophore for ARs, from a screen of an in-house library of only 52 compounds, originally designed for anti-proliferative activity. We identified L2-L10, A15, A17 with 3-aryl, 7-anilino and a electronegative group at 5-position as low micromolar to low nanomolar A1R / A3R antagonists. A17 has for A1R Kd = 5.62 nM and a residence time (RT) 41.33 min and for A3R Kd = 13.5 nM, RT = 47.23 min. The kinetic data showed that compared to the not potent or mediocre congeners the active compounds have similar association, for example at A1R Kon = 13.97 x106 M-1 (A17) vs Kon = 3.36 x106 M-1 (A26) but much lower dissociation rate Koff = 0.024 min-1 (A17) vs 0.134 min-1 (A26). Using molecular dynamics (MD) simulations and mutagenesis experiments we investigated the binding site of A17 showing that it can interact with an array of residues in transmembrane helix 5 (TM5), TM6, TM7 of A1R or A3R including residues E5.30, E5.28, T7.35 in A1R instead of Q5.28, V5.30 , L7.35 in A3R. A striking observation for drug design purposes is that for L2506.51A the binding affinity of A17 significantly increased at A1R. A17 provides a lead representative of a promising series and by means of the Thermodynamics Integration coupled with MD simulations (TI/MD) method, first applied here on whole GPCR- membrane system and showing a very good agreement between calculated and experimental relative binding free energies for A1R and A3R (spearman rank correlation p = 0.82 and 0.84, respectively), and kinetic experiments can lead to ligands with improved profile against ARs.

2021 ◽  
Author(s):  
Margarita Stampelou ◽  
Anna Suchankova ◽  
Eva Tzortzini ◽  
Lakshiv Dhingra ◽  
Kerry Barkan ◽  
...  

Drugs targeting the four adenosine receptor (AR) subtypes can provide “soft" treatment of various significant diseases. Even for the two experimentally resolved AR subtypes the description of the orthosteric binding area and structure-activity relationships of ligands remains a demanding task due to the high similar amino acids sequence but also the broadness and flexibility of the ARs binding area. The identification of new pharmacophoric moieties and nanomolar leads and the exploration of their binding area with mutagenesis and state-of-the-art computational methods useful also for drug design purposes remains a challenging aim for all ARs. Here, we identified several low nanomolar ligands and potent competitive antagonists against A1R / A3R, containing the novel pyrazolo[3,4-c]pyridine pharmacophore for ARs, from a screen of an in-house library of only 52 compounds, originally designed for anti-proliferative activity. We identified L2-L10, A15, A17 with 3-aryl, 7-anilino and a electronegative group at 5-position as low micromolar to low nanomolar A1R / A3R antagonists. A17 has for A1R Kd = 5.62 nM and a residence time (RT) 41.33 min and for A3R Kd = 13.5 nM, RT = 47.23 min. The kinetic data showed that compared to the not potent or mediocre congeners the active compounds have similar association, for example at A1R Kon = 13.97 x106 M-1 (A17) vs Kon = 3.36 x106 M-1 (A26) but much lower dissociation rate Koff = 0.024 min-1 (A17) vs 0.134 min-1 (A26). Using molecular dynamics (MD) simulations and mutagenesis experiments we investigated the binding site of A17 showing that it can interact with an array of residues in transmembrane helix 5 (TM5), TM6, TM7 of A1R or A3R including residues E5.30, E5.28, T7.35 in A1R instead of Q5.28, V5.30 , L7.35 in A3R. A striking observation for drug design purposes is that for L2506.51A the binding affinity of A17 significantly increased at A1R. A17 provides a lead representative of a promising series and by means of the Thermodynamics Integration coupled with MD simulations (TI/MD) method, first applied here on whole GPCR- membrane system and showing a very good agreement between calculated and experimental relative binding free energies for A1R and A3R (spearman rank correlation p = 0.82 and 0.84, respectively), and kinetic experiments can lead to ligands with improved profile against ARs.


2021 ◽  
Author(s):  
Margarita Stampelou ◽  
Anna Suchankova ◽  
Efpraxia Tzortzini ◽  
Lakshiv Dhingra ◽  
Kerry Barkan ◽  
...  

Drugs targeting the four adenosine receptor (AR) subtypes can provide “soft" treatment of various significant diseases. Even for the two experimentally resolved AR subtypes the description of the orthosteric binding area and structure-activity relationships of ligands remains a demanding task due to the high similar amino acids sequence but also the broadness and flexibility of the ARs binding area. The identification of new pharmacophoric moieties and nanomolar leads and the exploration of their binding area with mutagenesis and state-of-the-art computational methods useful also for drug design purposes remains a challenging aim for all ARs. Here, we identified several low nanomolar ligands and potent competitive antagonists against A1R / A3R, containing the novel pyrazolo[3,4-c]pyridine pharmacophore for ARs, from a screen of an in-house library of only 52 compounds, originally designed for anti-proliferative activity. We identified L2-L10, A15, A17 with 3-aryl, 7-anilino and a electronegative group at 5-position as low micromolar to low nanomolar A1R / A3R antagonists. A17 has for A1R Kd = 5.62 nM and a residence time (RT) 41.33 min and for A3R Kd = 13.5 nM, RT = 47.23 min. The kinetic data showed that compared to the not potent or mediocre congeners the active compounds have similar association, for example at A1R Kon = 13.97 x106 M-1 (A17) vs Kon = 3.36 x106 M-1 (A26) but much lower dissociation rate Koff = 0.024 min-1 (A17) vs 0.134 min-1 (A26). Using molecular dynamics (MD) simulations and mutagenesis experiments we investigated the binding site of A17 showing that it can interact with an array of residues in transmembrane helix 5 (TM5), TM6, TM7 of A1R or A3R including residues E5.30, E5.28, T7.35 in A1R instead of Q5.28, V5.30 , L7.35 in A3R. A striking observation for drug design purposes is that for L2506.51A the binding affinity of A17 significantly increased at A1R. A17 provides a lead representative of a promising series and by means of the Thermodynamics Integration coupled with MD simulations (TI/MD) method, first applied here on whole GPCR- membrane system and showing a very good agreement between calculated and experimental relative binding free energies for A1R and A3R (spearman rank correlation p = 0.82 and 0.84, respectively), and kinetic experiments can lead to ligands with improved profile against ARs.


2021 ◽  
Author(s):  
Margarita Stampelou ◽  
Anna Suchankova ◽  
Eva Tzortzini ◽  
Lakshiv Dhingra ◽  
Kerry Barkan ◽  
...  

Drugs targeting the four adenosine receptor (AR) subtypes can provide “soft" treatment of various significant diseases. Even for the two experimentally resolved AR subtypes the description of the orthosteric binding area and structure-activity relationships of ligands remains a demanding task due to the high similar amino acids sequence but also the broadness and flexibility of the ARs binding area. The identification of new pharmacophoric moieties and nanomolar leads and the exploration of their binding area with mutagenesis and state-of-the-art computational methods useful also for drug design purposes remains a challenging aim for all ARs. Here, we identified several low nanomolar ligands and potent competitive antagonists against A1R / A3R, containing the novel pyrazolo[3,4-c]pyridine pharmacophore for ARs, from a screen of an in-house library of only 52 compounds, originally designed for anti-proliferative activity. We identified L2-L10, A15, A17 with 3-aryl, 7-anilino and a electronegative group at 5-position as low micromolar to low nanomolar A1R / A3R antagonists. A17 has for A1R Kd = 5.62 nM and a residence time (RT) 41.33 min and for A3R Kd = 13.5 nM, RT = 47.23 min. The kinetic data showed that compared to the not potent or mediocre congeners the active compounds have similar association, for example at A1R Kon = 13.97 x106 M-1 (A17) vs Kon = 3.36 x106 M-1 (A26) but much lower dissociation rate Koff = 0.024 min-1 (A17) vs 0.134 min-1 (A26). Using molecular dynamics (MD) simulations and mutagenesis experiments we investigated the binding site of A17 showing that it can interact with an array of residues in transmembrane helix 5 (TM5), TM6, TM7 of A1R or A3R including residues E5.30, E5.28, T7.35 in A1R instead of Q5.28, V5.30 , L7.35 in A3R. A striking observation for drug design purposes is that for L2506.51A the binding affinity of A17 significantly increased at A1R. A17 provides a lead representative of a promising series and by means of the Thermodynamics Integration coupled with MD simulations (TI/MD) method, first applied here on whole GPCR- membrane system and showing a very good agreement between calculated and experimental relative binding free energies for A1R and A3R (spearman rank correlation p = 0.82 and 0.84, respectively), and kinetic experiments can lead to ligands with improved profile against ARs.


2018 ◽  
Author(s):  
Alejandro Lara ◽  
Maximiliano Riquelme ◽  
Esteban Vöhringer-Martinez

<div> <div> <div> <p>Partition coefficients serve in various areas as pharmacology and environmental sciences to predict the hydrophobicity of different substances. Recently, they have been also used to address the accuracy of force fields for various organic compounds and specifically the methylated DNA bases. In this study atomic charges were derived by different partitioning methods (Hirshfeld and Minimal Basis Iterative Stockholder) directly from the electron density obtained by electronic structure calculations in vac- uum, with an implicit solvation model or with explicit solvation taking the dynamics of the solute and the solvent into account. To test the ability of these charges to describe electrostatic interactions in force fields for condensed phases the original atomic charges of the AMBER99 force field were replaced with the new atomic charges and combined with different solvent models to obtain the hydration and chloroform solvation free energies by molecular dynamics simulations. Chloroform-water partition coefficients derived from the obtained free energies were compared to experimental and previously reported values obtained with the GAFF or the AMBER-99 force field. The results show that good agreement with experimental data is obtained when the polarization of the electron density by the solvent has been taken into account deriving the atomic charges of polar DNA bases and when the energy needed to polarize the electron den- sity of the solute has been considered in the transfer free energy. These results were further confirmed by hydration free energies of polar and aromatic amino acid side chain analogues. Comparison of the two partitioning methods Hirsheld-I and Minimal Basis Iterative Stockholder (MBIS) revealed some deficiencies in the Hirshfeld-I method related to nonexistent isolated anionic nitrogen pro-atoms used in the method. Hydration free energies and partitioning coefficients obtained with atomic charges from the MBIS partitioning method accounting for polarization by the implicit solvation model are in good agreement with the experimental values. </p> </div> </div> </div>


Author(s):  
Smriti Sharma ◽  
Vinayak Bhatia

: Pyrazole and its derivatives are a pharmacologically significant active scaffold that have innumerable physiological and pharmacological activities. They can be very good targets for the discovery of novel anti-bacterial, anticancer, anti-inflammatory, anti-fungal, anti-tubercular, antiviral, antioxidant, antidepressant, anti-convulsant and neuroprotective drugs. This review focuses on the importance of in silico manipulations of pyrazole and its derivatives for medicinal chemistry. The authors have discussed currently available information on the use of computational techniques like molecular docking, structure-based virtual screening (SBVS), molecular dynamics (MD) simulations, quantitative structure activity relationship (QSAR), comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) to drug design using pyrazole moieties. Pyrazole based drug design is mainly dependent on the integration of experimental and computational approaches. The authors feel that more studies need to be done to fully explore the pharmacological potential of the pyrazole moiety and in silico method can be of great help.


2020 ◽  
Vol 17 (9) ◽  
pp. 1084-1101
Author(s):  
Tingjuan Wu ◽  
Xu Yao ◽  
Guan Wang ◽  
Xiaohe Liu ◽  
Hongfei Chen ◽  
...  

Background: Oleanolic Acid (OA) is a ubiquitous product of triterpenoid compounds. Due to its inexpensive availability, unique bioactivities, pharmacological effects and non-toxic properties, OA has attracted tremendous interest in the field of drug design and synthesis. Furthermore, many OA derivatives have been developed for ameliorating the poor water solubility and bioavailability. Objective: Over the past few decades, various modifications of the OA framework structure have led to the observation of enhancement in bioactivity. Herein, we focused on the synthesis and medicinal performance of OA derivatives modified on A-ring. Moreover, we clarified the relationship between structures and activities of OA derivatives with different functional groups in A-ring. The future application of OA in the field of drug design and development also was discussed and inferred. Conclusion: This review concluded the novel achievements that could add paramount information to the further study of OA-based drugs.


2012 ◽  
Vol 11 (02) ◽  
pp. 283-295 ◽  
Author(s):  
BAHRAM GHALAMI-CHOOBAR ◽  
ALI GHIAMI-SHOMAMI ◽  
PARIA NIKPARSA

In this work, calculations of p K b values have been performed for aniline and its substituted derivatives and sulfonamide drugs by using Gaussian 98 software package. Gas-phase energies were calculated with HF /6-31 G ** and B3LYP /6-31 G ** levels of theory. Free energies of solvation have been computed using the polarizable continuum model (PCM), conductor-like polarizable continuum model (CPCM) and the integral equation formalism-polarizable continuum model (IEFPCM) at the same levels which have been used for geometry determination in the gas-phase. The results show that the calculated p K b values using the B3LYP /6-31 G ** are better than those using the corresponding HF /6-31 G **. At first, the correlation equation was found to determine the p K b values of the investigated anilines. Then, this correlation equation was used to calculate the p K b values of the sulfonamide drugs. The results obtained indicate that the PCM model is a suitable solvation model for calculating p K b values in comparison to the other solvation models. For the investigated compounds a good agreement between the experimental and the calculated p K b values was also observed.


2017 ◽  
Vol 10 (5) ◽  
pp. 1823-1830 ◽  
Author(s):  
Marin S. Robinson ◽  
Irena Grgić ◽  
Vid S. Šelih ◽  
Martin Šala ◽  
Marsha Bitsui ◽  
...  

Abstract. A widely used instrument for collecting size-segregated particles is the micro-orifice uniform deposit impactor (MOUDI). In this work, a 10-stage MOUDI (cut-point diameter of 10 µm to 56 nm) was used to collect samples in Ljubljana, Slovenia, and Martinska, Croatia. Filters, collected with and without rotation, were cut in half and analyzed for nine elements (As, Cu, Fe, Ni, Mn, Pb, Sb, V, Zn) using laser ablation ICP-MS. Elemental image maps (created with ImageJ) were converted to concentrations using NIST SRM 2783. Statistical analysis of the elemental maps indicated that for submicron particles (stages 6–10), ablating 10 % of the filter (0.5 cm2, 20 min ablation time) was sufficient to give values in good agreement (±10 %) to analysis of larger parts of the filter and with good precision (RSE < 1 %). Excellent sensitivity was also observed (e.g., 20 ± 0.2 pg m−3 V). The novel use of LA-ICP-MS, together with image mapping, provided a fast and sensitive method for elemental analysis of size-segregated MOUDI filters, particularly for submicron particles.


2017 ◽  
Vol 53 (2) ◽  
pp. 85-93 ◽  
Author(s):  
J. Zhou ◽  
L. Zhang ◽  
L. Chen ◽  
Y. Du ◽  
Z.K. Liu

A critical thermodynamic assessment of the metastable c-TiAlZrN coatings, which are reported to spinodally decompose into triple domains, i.e., c-TiN, c-AlN, and c-ZrN, was performed via the CALculation of PHAse Diagram (CALPHAD) technique based on the limited experimental data as well as the first-principles computed free energies. The metastable c-TiAlZrN coatings were modeled as a pseudo-ternary phase consisting of c-TiN, c-AlN and c-ZrN species, and described using the substitutional solution model. The thermodynamic descriptions for the three boundary binaries were directly taken from either the CALPHAD assessment or the first-principles results available in the literature except for a re-adjustment of the pseudo-binary c-AlN/c-ZrN system based on the experimental phase equilibria in the pseudo-ternary system. The good agreement between the calculated phase equilibria and the experimental data over the wide temperature range was obtained, validating the reliability of the presently obtained thermodynamic descriptions for the c-TiAlZrN system. Based on the present thermodynamic description, different phase diagrams and thermodynamic properties can be easily predicted. It is anticipated that the present thermodynamic description of the metastable c-TiAlZrN coatings can serve as the important input for the later quantitative description of the microstructure evolution during service life.


Sign in / Sign up

Export Citation Format

Share Document