Prerpints.org logo
Preprint
Article

Evaluation of Cx43 gap junction inhibitors using Quantitative Structure-Activity Relationship model

Altmetrics

Downloads

168

Views

59

Comments

0

Submitted:

03 April 2023

Posted:

05 April 2023

You are already at the latest version

Alerts
Abstract
Gap junctions (GJs) made of connexin-43 (Cx43) are important for the conduction of electrical impulses in the heart. Modulation of Cx43 activity may be useful in the treatment of cardiac arrhythmias and other dysfunctions. Search of novel GJ modulating agents using molecular docking approach allows to predict the binding affinities that often significantly differ from experimentally estimated their potencies. The objective of this study was to demonstrate that Quantitative Structure-Activity Relationship (QSAR) model could be used for more precise identification of potent Cx43 GJ inhibitors. Using the QSAR model and molecular docking, we evaluated the known Cx43 GJ inhibitors, suggested a new one, and tested it experimentally. Our QSAR model predicted the concentrations required to produce 50% of the maximal effect (IC50) for each of these compounds and estimated the correlation between predicted and experimental IC50ies (pIC50 and eIC50). This led to suggestion of monocyclic monoterpene d-limonene as putative Cx43 inhibitor with pIC50 = 1.07. In turn, dual whole-cell patch-clamp measurements provided eIC50=1.41. The pIC50ies of d-limonene and other Cx43 GJ inhibitors examined by our QSAR model well correlated with their eIC50ies (R = 0.88) in contrast to pIC50s obtained from molecular docking (R = 0.42). However, the molecular docking suggested that inhibitor potency may depend on their docking sites on Cx43.
Keywords: 
Subject: Medicine and Pharmacology  -   Medicine and Pharmacology

1. Introduction

Gap junctions (GJs) are intercellular channels indispensable for electrical interaction between cardiac myocytes and synchronized cardiac contraction [1]. Connexin (Cx)-based GJ channels are formed of two opposing hemichannels in the contiguous cells (reviewed in [2]). Six Cx subunits compose a hemichannel with ion-selective pore. Each Cx protein has four alpha helical transmembrane domains, intracellular N- and C-termini, two extracellular loops, and a cytoplasmic loop. Cxs are expressed in all tissues except differentiated skeletal muscle, erythrocytes, and mature sperm cells, and the family of Cx genes consists of 21 members in the human genome. The prevailing connexin isoform in human cardiac tissue is Cx43 [3]. Changes in the function, expression or localization of the Cx43 are associated with a higher frequency and severity of arrhythmias and sudden death in patients with some cardiovascular diseases [3]. On the other hand, modulation of GJ function is onerous due to the shortage of specific and Cx type-selective GJ inhibitors. Most GJ inhibitors are non-specific compounds, such as antimalarial drugs, polyamines, glycyrrhetinic acid, volatile anaesthetics, arachidonic acid, cyclodextrins, anti-cancer drugs cisplatin and oxaliplatin, fatty acid amides, terpenes [4,5,6], and other drugs. In addition, most of these compounds have been shown to inhibit not only intercellular communication through GJs but also affect the membranous Na+, K+, and Ca2+ channels that are crucial for the generation and spread of action potential [7,8,9].
Molecular docking approaches can relatively well predict the binding affinities of compounds; however, receptor occupancy often is not directly translated to its biological effect that also depends on the intrinsic activity and intrinsic efficacy of the compound [10]. The quantitative structure-activity relationship (QSAR) is a computational modeling method used for identifying the relationship between structural characteristics of chemical compounds and their binding affinities and biological activities [11]. QSAR analysis is particularly useful in the pharmacy industry as this method helps to select substances according to their desirable biological activity, thus significantly reducing the number of substances that need to be tested in vitro and in vivo [11]. QSAR was used to develop simple, accurate, and robust models useful to identify new inhibitors [12]. For instance, in the absence of a reliable 3D receptor-bound structure more complex 3D modeling methods might face limitations [13]. Modulation of Cx43 GJ activity (activation or inhibition) is expected to be beneficial in management of the cardiac pathologies [14], such as ischemic heart disease, heart failure, hypertrophic cardiomyopathy, arrhythmogenic right ventricular cardiomyopathy, and others [15]. Therefore, the aim of this study was to evaluate known and putative Cx43 inhibitors using the molecular docking and QSAR models and compare predicted IC50 (pIC50) values with experimental ones (eIC50) obtained by dual whole-cell patch-clamp experiments.

2. Materials and Methods

2.1. Bioinformatic analysis

The open access databases were searched for studies on the inhibitory effects of various substances on Cx43 GJs. The following substances with their half maximal inhibitory concentrations (IC50) were found [14,16,17,18,19,20,21]: 2-aminoethoxydiphenyl borate (2-APB), 3,3′-methylene-bis-4-hydroxycoumarin, α-copaene, α-pinene, sabinene, carbenoxolone, digoxin, dihydrogambogic acid (DGBA), heptanol, linoleic acid, meclofenamic acid, mefloquine, quinidine, warfarin. In such a way, the dataset for Cx43 inhibition QSAR was made of 14 compounds belonging to different classes (Table 1). These classes include terpenes (α-pinene, α-copaene, sabinene), benzenes (meclofenamic acid), quinolines (mefloquine), Cinchona alkaloids (quinidine), fatty acids (linoleic acid), coumarins and derivatives (dicumarol, warfarin), triterpenoids (carbenoxolone), fatty alcohols (heptanol), cardenolide glycosides (digoxin), organoboron compound (2-APB), pyranoxanthones (DGBA). For sake of clarity, PubChem CID of all inhibitors is also provided.
Given our previous research of terpenes as Cx43 inhibitors [19] widely mentioned in literature monoterpene d-limonene [22,23,24] gained our attention as a possible Cx43 inhibitor. We decided to test this hypothesis with molecular modeling and further with dual whole-cell patch-clamp experiments.

2.2. Molecular docking

The rat Cx43 homology model was created in the Phyre2 web portal for protein modeling, prediction, and analysis [25]. Human Cx43 structure was used as a template (PDB ID 7F93) [26] for the construction of rat Cx43 homology model as a target for molecular docking. 3D molecular models for all ligands except 2-APB, digoxin, and DGBA were available in the PubChem database [27]. These 3 ligands either without clear 3D structures (digoxin and DGBA) or organoboron compound (2-APB) (no modelling parameters for B atom) were excluded from docking. Gnina version 1.0 was used for molecular docking calculations, which utilizes an ensemble of convolutional neural networks (CNNs) as a scoring function [28]. The docking mode was set to whole protein docking. All ligand output conformations except top scoring ones were excluded from further analysis automatically with program settings leaving only the most reliable one for each ligand. 3D docking image was generated with ChimeraX 1.4 software [29] and 2D docking plots were generated with LigPlot+ 2.2 software [30]. Log(pIC50) from molecular docking was calculated using the equation of linear regression:
log ⁡ p I C 50 = m · D M A + n  
where m and n are fitting coefficients and DMA is Docking Minimized Affinity.

2.3. QSAR model development

QSAR model development requires molecular descriptors calculated from their chemical structure. PaDEL-Decriptor software [31] was used to calculate 1444 molecular descriptors for those 14 compounds. Some of these descriptors had identical values for all 14 compounds like C2SP1 – amount of carbon atoms in investigated compound chemical structure containing chemical bounds to two other nearby carbon atoms and exactly one more atom of a different type had identical values for all analysed compounds so were considered useless. In total 596 molecular descriptors remained excluding those unsuitable for further modeling.
Regression analysis methods are statistical instruments widely used for the determination of relationships between molecular descriptors and biological activities of compounds [32]. Multiple linear regression (MLR) [33] for QSAR was developed by using R software for statistical computing (https://www.r-project.org/) with a Leaps package. Stepwise regression using a Forward Selection Approach [34] and also other approaches available in that package were applied. MLR was chosen to obtain relatively easily understandable QSAR models as it describes predicted biological activities in the form:
log ⁡ p I C 50 = a 1 x 1 + a 2 x 2 + … + a n x n + b    
where pIC50 is predicted biological activity; xn, is molecular descriptor calculated by software; an is its fitting coefficient; b is the assumed log(pIC50) value when all molecular descriptors are equal to 0. Given the limited size of the data set, the number of chosen descriptors for the final model (n) was chosen to be just 2 to avoid overfitting [35]. All possible models (175806 in total) with those two descriptors were created and analysed searching for one providing the most significant correlation between pIC50 and eIC50 values. Promising QSAR models were transferred into Microsoft Excel software for final inspection and validation. Microsoft Excel Analysis ToolPak add-in was used to compute the final values of regression analysis: R - coefficient of correlation - measured both the strength and the direction of a linear relationship; R2 - coefficient of determination - provided percentage variation making it easier to compare between different models; adjusted R2 helped to spot issues of overfitting. SigmaPlot 12 was also used for validation and final QSAR plot preparation.

2.4. Cell lines and culture conditions

Experiments were performed on HeLa (human cervix carcinoma, ATCC CCL-2, Manassas, VA, USA) cells stably transfected with Cx43 tagged with a green fluorescent protein (Cx43-EGFP). A stable HeLa cell line expressing Cx43-EGFP was obtained in collaboration with Dr. F. Bukauskas (Albert Einstein College of Medicine, New York, USA). The construction protocol of the vector is described elsewhere [36]. Cell line expressing Cx43-EGFP was selected using 500 µg/ml G418/geneticin (Sigma-Aldrich Co.). Cells were grown in DMEM medium containing 10% fetal bovine serum (FBS), penicillin/streptomycin mix (100 U/ml penicillin and 100 μg/ml streptomycin; Gibco Laboratories) at 37 °C and 5 % CO2. Typically, the cells were analyzed on the second day after passage.

2.5. Electrophysiological Measurements

For electrophysiological recordings, the cells grown onto glass coverslips were transferred to an experimental chamber with constant flow-through perfusion mounted on the stage of the inverted microscope Olympus IX81 equipped with the Orca-R2 cooled digital camera. Junctional conductance gj between contiguous cells was measured using the dual whole-cell patch-clamp technique [37]. Cell-1 and Cell-2 of a cell pair (Figure 4C) were voltage clamped independently with the patch-clamp amplifier MultiClamp 700B (Molecular Devices, Inc., USA) at the same holding potential, V1 = V2. By applying repeated every 10ths voltage ramps (-10 mV, 20 ms) (Figure 4D) in the Cell-1 (ΔV1) and keeping the other constant, the junctional current was measured as the change in current in the Cell-2, Ij = ΔI2. Thus, gj was obtained from the ratio –Ij/ΔV1, where ΔV1 is equal to transjunctional voltage (Vj), and a negative sign indicates that the junctional current measured in the Cell-2 is oppositely oriented to the one measured in the Cell-1. Voltages and currents were digitized using the Digidata 1440A data acquisition system (Molecular Devices, Inc., USA) and acquired and analyzed using pClamp 10 software (Molecular Devices, Inc., USA). Patch pipettes were pulled from borosilicate glass capillary tubes with filaments. Experiments were performed at room temperature in a modified Krebs-Ringer solution (in mM): NaCl, 140; KCl, 4; CaCl2, 2; MgCl2, 1; glucose, 5; pyruvate, 2; HEPES, 5 (pH 7.4). Patch pipettes were filled with solution containing (in mM): KCl, 130; Na aspartate, 10; MgATP, 1; MgCl2, 1; CaCl2, 0.2; EGTA, 2; HEPES, 5 (pH 7.3). All chemical reagents were purchased from Sigma-Aldrich Corp. Stock solutions of d-limonene were prepared in DMSO at 100 mM concentration and later diluted with modified Krebs-Ringer solution into needed concentration (10, 30, 50 or 100 µM).

2.6. Statistical analysis

Dose–response curve obtained with different concentrations of d-limonene was fitted to a three-parameter Hill’s equation, and a concentration of the compound required to produce 50% of the maximal effect (IC50) was derived using SigmaPlot 12.0 software. Data are reported as means ± SEM.

3. Results

3.1. Molecular docking of Cx43 inhibitors

Pharmacological characterization of d-limonene and other compounds was performed in silico using molecular simulation software. Molecular docking showed a common docking site for the majority of examined compounds except for sabinene, α-pinene, mefloquine, and carbenoxolone. This common docking site is a large hydrophobic furrow on the Cx43 transmembrane domains (TMs 1, 2, and 3) between neighbouring subunits (Figure 1A and Supplementary Video). Sabinene and α-pinene docked on a different site of TM3 close to extracellular loops. This hydrophobic location is suitable for docking such small highly hydrophobic compounds, but not for larger ones. Carbenoxalone docked on the Cx43 TM4 at a unique location containing two contiguous phenylalanines well fitting for its more hydrophilic part. Mefloquine docked on the pore side of TM2 reaching the N terminal domain. The exact docking locations of investigated Cx43 inhibitors are shown in 2D plots (Figure 2). From those images, it can be concluded that compounds of a highly hydrophobic nature and appropriate size are promising as Cx43 inhibitors. The docking of d-limonene is shown in 3D model (Figure 1B) which shares the same common docking site with the majority of other compounds. A fair correlation (R = 0.42) between docking minimized affinities and eIC50ies was found presumably due to quite a large discrepancy between compound affinities and potencies (Figure 3A). From linear regression (Equation 1) m and n were determined equal to 0.39 and 3.99, respectively. Using these values, pIC50 was calculated from docking (Table 1). Docking minimized affinity for each compound is also presented in Table 1.

3.2. QSAR modeling of Cx43 inhibitors

Using the R leaps package-based approaches the molecular descriptor pair BCUTp-1l and SM1_Dzs was selected for Cx43 inhibition QSAR modeling. BCUTp-1l is the first lowest eigenvalue of the burden matrix weighted by polarization and SM1_Dzs is the spectral moment of order 1 from the Barysz matrix weighted by I-state [38,39]. Using them the following QSAR model was developed for Cx43 inhibition:
log ⁡ ( p I C 50 ) = − 1.02 × B C U T p 1 l − 0.41 × S M 1 _ D z s + 7.18
The developed QSAR model allowed to calculate pIC50 of each compound from two calculated molecular descriptor values and compare with their eIC50ies. The best QSAR model (Equation 3) achieved a very strong correlation between pIC50s and eIC50s: R = 0.88, R2 = 0.78, and R2adj = 0.76. Values of pIC50 calculated with this model are provided in Table 1 together with literature data for comparison. The same values were used in Figure 3B where eIC50 was plotted against pIC50 calculated by this model. The substance with the lowest values of selected molecular descriptors (BCUTp-1l and SM1_Dzs) was heptanol and it inhibited Cx43 with the highest IC50. Conversely, the compounds with the highest values of these descriptors - α-copaene and digoxin - inhibited Cx43 GJs with the lowest IC50 concentration.
Preprints 70490 g003
Figure 3. The plot of log(eIC50) values versus docking minimized affinities for the docked Cx43 inhibitors (A). The plot of log(eIC50) values versus log(pIC50) values from the QSAR model (Equation 3) for the Cx43 inhibitors (B).
Figure 3. The plot of log(eIC50) values versus docking minimized affinities for the docked Cx43 inhibitors (A). The plot of log(eIC50) values versus log(pIC50) values from the QSAR model (Equation 3) for the Cx43 inhibitors (B).
Preprints 70490 g003
Table 1. Values of selected molecular descriptors (BCUTp-1l and SM1_Dzs) and values of experimental log(eIC50), predicted log(pIC50), and docking minimized affinities of the examined substances.
Table 1. Values of selected molecular descriptors (BCUTp-1l and SM1_Dzs) and values of experimental log(eIC50), predicted log(pIC50), and docking minimized affinities of the examined substances.
PubChem CID Name Log(eIC50) BCUTp-1l SM1_Dzs Log(pIC50) Docking minimized affinity Log(pIC50) from docking
11240513 α-Pinene [19] 1.08 6.26 -0.33 0.92 -5.18 1.96
12303902 α-Copaene [19] 0.14 6.49 -0.49 0.75 -6.12 1.59
18818 Sabinene [19] 1.57 6.23 -0.64 1.08 -4.90 2.07
4037 Meclofenamic acid [18] 2.42 4.44 0.85 2.28 -6.41 1.47
4046 Mefloquine [17] 0.95 3.73 5.91 0.92 -6.66 1.38
441074 Quinidine [17] 2.60 4.91 0.92 1.77 -6.57 1.41
5280450 Linoleic acid [21] 2.00 3.72 1.19 2.88 -5.75 1.73
54676038 Dicumarol [14] 0.47 5.05 3.67 0.50 -7.54 1.03
54678486 Warfarin [14] 0.87 4.66 2.38 1.43 -6.73 1.35
636403 Carbenoxolone [18] 2.32 3.75 4.18 1.62 -7.09 1.21
8129 Heptanol [18] 3.34 3.93 0.66 2.89 -3.50 2.62
1598 2-APB [20] 1.71 4.64 0.41 2.27 - -
2724385 Digoxin [18] -0.86 4.26 8.45 -0.66 - -
6857793 DGBA [16] 1.11 4.30 4.20 1.05 - -
22311 d-Limonene 1.41 6.18 -0.50 1.07 -5.12 1.98

3.3. D-limonene dose-dependently inhibits Cx43 GJ conductance

In our earlier study, we have demonstrated that constituents of nutmeg essential oil - monoterpenes sabinene and α-pinene and sesquiterpene α-copaene - were potent and efficient Cx43 GJ inhibitors [19]. In the current study, we aimed to examine the effect of another constituent of nutmeg essential oil, suggested by QSAR modeling using the molecular descriptors, monocyclic monoterpene d-limonene, on the conductance of GJs composed of Cx43 (Figure 4 A,B) exogenously expressed in HeLa cells.
To determine the effect of d-limonene on Cx43 GJ conductance, we performed dual whole-cell patch-clamp experiments in pairs of HeLa cells expressing exogenous Cx43-EGFP (Figure 4C) applying -10 mV, 20 ms voltage ramps to Cell-1 and measuring junctional current in Cell-2 (Figure 4D). A threshold concentration of d-limonene for inhibition of gj was 10 µM (Figure 4E). Further, applying higher concentrations (30, 50 and 100 µM) (Figure 4F–H) we found that gj could be completely blocked by 100 µM of d-limonene. IC50 value of 30 µM was derived from the fit of the experimental points to the Hill’s equation (Figure 4I). Hill’s coefficient was 2.8 suggesting more than one binding site on Cx43 GJ channel, similar to that obtained for other terpenes in our earlier study [5].
Preprints 70490 g004
Figure 4. The effect of d-limonene on Cx43 GJ conductance. (A, B) Cx43 GJ homology model (side and top view, respectively). (C) Dual whole-cell patch-clamp measurement of Cx43-EGFP conductance in HeLa cells. (D) Gj was measured by applying repeated -10 mV, 20 ms Vj ramps which do not cause the voltage-dependent gating of Cx43 GJ channels. (E-H) Typical effects of d-limonene at indicated concentrations on Cx43 GJ conductance. (I) Dose-dependence of d-limonene effect on Cx43 GJ conductance (eIC50 = 30 µM; Hill’s coefficient = 2.8).
Figure 4. The effect of d-limonene on Cx43 GJ conductance. (A, B) Cx43 GJ homology model (side and top view, respectively). (C) Dual whole-cell patch-clamp measurement of Cx43-EGFP conductance in HeLa cells. (D) Gj was measured by applying repeated -10 mV, 20 ms Vj ramps which do not cause the voltage-dependent gating of Cx43 GJ channels. (E-H) Typical effects of d-limonene at indicated concentrations on Cx43 GJ conductance. (I) Dose-dependence of d-limonene effect on Cx43 GJ conductance (eIC50 = 30 µM; Hill’s coefficient = 2.8).
Preprints 70490 g004

4. Discussion

Cardiac remodeling, which involves structural and electrical changes in the heart, may be impacted by altered expression and localization of Cx43 GJs. Decreased expression of Cx43 proteins and heterogeneous arrangement of channels can impair cardiac conduction and lead to supraventricular or ventricular arrhythmias [40]. Cx channels are promising pharmacological targets because inhibitors of Cx channels could be useful for treating not only arrhythmias, but also other communication-dependent diseases affecting other body systems. The importance of Cx43 has been well established, particularly in the heart where a lack of Cx43 leads to abnormal cardiac development and death at birth [41]. Abnormalities in Cx43 organization and regulation have also been linked to myocardial ischemia [42]. Preservation of Cx43 in the intercalated discs joining cardiomyocytes may not only support conduction velocity but also impact the total normoxic Cx43 interactome [43]. So, Cx43 is a considerable drug target, especially during heart ischemia and reperfusion [14,44,45]. For example, Cx43 inhibitor digoxin (analyzed in QSAR here) isolated from Digitalis lanata is well known in cardiology [46]. It is used to treat both irregular heartbeat [47] and heart failure [48], but its side effects like gynecomastia also are significant which could be explained through its chemical similarity to estrogen [49]. Another Cx43 inhibitor quinidine also is a popular antiarrhythmic drug [50]. Another examined compound carbenoxolone is used for the treatment of peptic, esophageal and oral ulceration and inflammation; however, it also has been shown in humans to slow the myocardial conduction [51]. An anti-malarial drug mefloquine may lead to complete heart arrest [52] and number of neuropsychiatric effects including suicide [53]. More recently, Cx inhibiting peptides (antiarrhythmic peptides AAP10, ZP123; GAP-134; RXP-E; and the Cx43 mimetic peptides Gap 26 and Gap 27) were suggested for treatment of arrhythmias in patients with ischaemic cardiomyopathy [54]. On the other hand, peptides underperformed as drug candidates due to unfavorable characteristics mainly regarding their pharmacokinetic behavior, including plasma stability, membrane permeability, and circulation half-life [55]. The discovery of new modulators of GJ channel function lacking similarity to steroid hormones to avoid side effects is of interest for human health [14,44,45]. Also, a long-standing challenge in the study of GJs is the lack of specific, high-affinity activators and inhibitors of GJ channels [56,57]. Therefore, it is important to predict in silico which substances could effectively modulate GJ conductance. This also could serve as an innovative approach in repurposing of licensed drugs with predicted new Cx inhibitory properties for other communication-dependent illnesses.
The results obtained in molecular docking considering minimized affinity had much worse correlation with log(eIC50) (R was 0.42) when compared to QSAR modeling (R was 0.88) so it might be concluded that QSAR outperforms docking. Fair correlation in case of molecular docking can be explained by often happening discrepancy between affinity and potency of the compounds [58].
All top 3 Cx43 inhibitors docked into the same common docking site suggesting that this site should be explored further for more potent inhibitors. Thus, even if molecular docking accuracy is low it can be used in combination with QSAR. Using the developed QSAR model, the limonene log(pIC50) was calculated equal to 1.07 (Equation 3), and an experimental value of 1.41 was later determined by our patch-clamp experiments. Therefore, this new Cx43 inhibitor could be added the current GJ inhibitor nomenclature.

Supplementary Materials

The following supporting information can be downloaded at the website of this paper posted on Preprints.org, Video S1: 3D version of Figure 1.

Author Contributions

Conceptualization, V.R.; methodology, V.R., Ra.M. and E.I.; validation, V.R. and V.A.S.; formal analysis, V.R., Ra.M., E.I. and V.A.S.; investigation, V.R., Ra.M., E.I. and Ro.M.; resources, V.A.S.; data curation, Ro.M., V.R., Ra.M. and E.I.; writing—original draft preparation, Ra.M. and E.I.; writing—review and editing, V.R. and V.A.S.; visualization, V.R. and V.A.S.; supervision, V.R. and V.A.S.; project administration, V.A.S.; funding acquisition, V.A.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Research Council of Lithuania, grant number P-MIP-23-105 and the APC was funded by the same grant.

Data Availability Statement

All data are provided within the article or supplementary material.

Acknowledgments

We thank Ms. Kamilė Eitkevičiūtė for skillful technical assistance and preparation of the cells.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Del Ry, S.; Moscato, S.; Bianchi, F.; Morales, M.A.; Dolfi, A.; Burchielli, S.; Cabiati, M.; Mattii, L. Altered expression of connexin 43 and related molecular partners in a pig model of left ventricular dysfunction with and without dipyrydamole therapy. Pharmacological Research 2015, 95, 92–101. [Google Scholar] [CrossRef] [PubMed]
  2. Račkauskas, M.; Neverauskas, V.; Skeberdis, V.A. Diversity and properties of connexin gap junction channels. Medicina 2010, 46, 1. [Google Scholar] [CrossRef]
  3. Lucero, C.M.; Andrade, D.C.; Toledo, C.; Díaz, H.S.; Pereyra, K.V.; Diaz-Jara, E.; Schwarz, K.G.; Marcus, N.J.; Retamal, M.A.; Quintanilla, R.A. Cardiac remodeling and arrhythmogenesis are ameliorated by administration of Cx43 mimetic peptide Gap27 in heart failure rats. Scientific Reports 2020, 10, 1–12. [Google Scholar] [CrossRef] [PubMed]
  4. Srinivas, M. Pharmacology of connexin channels. In Connexins: a guide, Harris, A.L., Locke, D., Eds.; Humana Press: New York, 2009; pp. 207–224. [Google Scholar]
  5. Mickus, R.; Jančiukė, G.; Raškevičius, V.; Mikalayeva, V.; Matulytė, I.; Marksa, M.; Maciūnas, K.; Bernatonienė, J.; Skeberdis, V.A. The effect of nutmeg essential oil constituents on Novikoff hepatoma cell viability and communication through Cx43 gap junctions. Biomed Pharmacother 2021, 135, 111229. [Google Scholar] [CrossRef] [PubMed]
  6. Wang, Q.; You, T.; Yuan, D.; Han, X.; Hong, X.; He, B.; Wang, L.; Tong, X.; Tao, L.; Harris, A.L. Cisplatin and oxaliplatin inhibit gap junctional communication by direct action and by reduction of connexin expression, thereby counteracting cytotoxic efficacy. J Pharmacol Exp Ther 2010, 333, 903–911. [Google Scholar] [CrossRef]
  7. Du, Y.-m.; Xia, C.-k.; Zhao, N.; Dong, Q.; Lei, M.; Xia, J.-h. 18β-Glycyrrhetinic acid preferentially blocks late Na current generated by ΔKPQ Nav1.5 channels. Acta Pharmacologica Sinica 2012, 33, 752–760. [Google Scholar] [CrossRef]
  8. Han, J.; Su, G.-h.; Wang, Y.-h.; Lu, Y.-x.; Zhao, H.-l.; Shuai, X.-x. 18β-Glycyrrhetinic Acid Improves Cardiac Diastolic Function by Attenuating Intracellular Calcium Overload. Current Medical Science 2020, 40, 654–661. [Google Scholar] [CrossRef]
  9. Meves, H. Arachidonic acid and ion channels: an update. British journal of pharmacology 2008, 155, 4–16. [Google Scholar] [CrossRef]
  10. Kenakin, T. Pharmacologic Analysis of Drug-Receptor Interaction, 2 ed.; Raven Press, Ltd.: New York, NY 10036, 1993. [Google Scholar]
  11. Kwon, S.; Bae, H.; Jo, J.; Yoon, S. Comprehensive ensemble in QSAR prediction for drug discovery. BMC bioinformatics 2019, 20, 1–12. [Google Scholar] [CrossRef]
  12. Liang, Y.n.; Qin, D.; Zhang, Y.; Liu, W.; Liang, G. Comprehensive interactions of ACE inhibitors with their receptor by a Support Vector Machine model and molecular docking. Journal of the Chinese Chemical Society 2017, 64, 310–320. [Google Scholar] [CrossRef]
  13. Cramer, R.D. Topomer CoMFA: a design methodology for rapid lead optimization. Journal of medicinal chemistry 2003, 46, 374–388. [Google Scholar] [CrossRef] [PubMed]
  14. Salameh, A.; Dhein, S. Pharmacology of gap junctions. New pharmacological targets for treatment of arrhythmia, seizure and cancer? Biochimica et Biophysica Acta (BBA)-Biomembranes 2005, 1719, 36–58. [Google Scholar] [CrossRef] [PubMed]
  15. Marsh, S.R.; Williams, Z.J.; Pridham, K.J.; Gourdie, R.G. Peptidic connexin43 therapeutics in cardiac reparative medicine. Journal of cardiovascular development and disease 2021, 8, 52. [Google Scholar] [CrossRef]
  16. Choi, E.J.; Yeo, J.H.; Yoon, S.M.; Lee, J. Gambogic acid and its analogs inhibit gap junctional intercellular communication. Frontiers in pharmacology 2018, 9, 814. [Google Scholar] [CrossRef] [PubMed]
  17. Picoli, C.; Nouvel, V.; Aubry, F.; Reboul, M.; Duchêne, A.; Jeanson, T.; Thomasson, J.; Mouthon, F.; Charvériat, M. Human connexin channel specificity of classical and new gap junction inhibitors. Journal of biomolecular screening 2012, 17, 1339–1347. [Google Scholar] [CrossRef]
  18. Burnham, M.; Sharpe, P.; Garner, C.; Hughes, R.; Pollard, C.; Bowes, J. Investigation of connexin 43 uncoupling and prolongation of the cardiac QRS complex in preclinical and marketed drugs. British Journal of Pharmacology 2014, 171, 4808–4819. [Google Scholar] [CrossRef] [PubMed]
  19. Mickus, R.; Jančiukė, G.; Raškevičius, V.; Mikalayeva, V.; Matulytė, I.; Marksa, M.; Maciūnas, K.; Bernatonienė, J.; Skeberdis, V.A. The effect of nutmeg essential oil constituents on Novikoff hepatoma cell viability and communication through Cx43 gap junctions. Biomedicine & Pharmacotherapy 2021, 135, 111229. [Google Scholar]
  20. D'hondt, C.; Ponsaerts, R.; De Smedt, H.; Bultynck, G.; Himpens, B. Pannexins, distant relatives of the connexin family with specific cellular functions? Bioessays 2009, 31, 953–974. [Google Scholar] [CrossRef]
  21. Willebrords, J.; Maes, M.; Yanguas, S.C.; Vinken, M. Inhibitors of connexin and pannexin channels as potential therapeutics. Pharmacology & therapeutics 2017, 180, 144–160. [Google Scholar]
  22. Karr, L.L.; Coats, J.R. Insecticidal properties of d-limonene. Journal of Pesticide Science 1988, 13, 287–290. [Google Scholar] [CrossRef]
  23. Kim, Y.W.; Kim, M.J.; Chung, B.Y.; Bang, D.Y.; Lim, S.K.; Choi, S.M.; Lim, D.S.; Cho, M.C.; Yoon, K.; Kim, H.S. Safety evaluation and risk assessment of d-limonene. Journal of Toxicology and Environmental Health, Part B 2013, 16, 17–38. [Google Scholar] [CrossRef] [PubMed]
  24. Sun, J. D-Limonene: safety and clinical applications. Alternative medicine review 2007, 12, 259. [Google Scholar] [PubMed]
  25. Kelley, L.A.; Mezulis, S.; Yates, C.M.; Wass, M.N.; Sternberg, M.J. The Phyre2 web portal for protein modeling, prediction and analysis. Nature protocols 2015, 10, 845–858. [Google Scholar] [CrossRef] [PubMed]
  26. Woo, J.-S.; Lee, H.-J.; Cha, H.J.; Jeong, H.; Lee, S.-N.; Lee, C.-W.; Kim, M.; Yoo, J. Structural insights into the gating mechanism of human Cx43/GJA1 gap junction channel. 2021.
  27. Kim, S.; Chen, J.; Cheng, T.; Gindulyte, A.; He, J.; He, S.; Li, Q.; Shoemaker, B.A.; Thiessen, P.A.; Yu, B. PubChem 2019 update: improved access to chemical data. Nucleic acids research 2019, 47, D1102–D1109. [Google Scholar] [CrossRef] [PubMed]
  28. McNutt, A.T.; Francoeur, P.; Aggarwal, R.; Masuda, T.; Meli, R.; Ragoza, M.; Sunseri, J.; Koes, D.R. GNINA 1.0: molecular docking with deep learning. Journal of cheminformatics 2021, 13, 1–20. [Google Scholar] [CrossRef]
  29. Goddard, T.D.; Huang, C.C.; Meng, E.C.; Pettersen, E.F.; Couch, G.S.; Morris, J.H.; Ferrin, T.E. UCSF ChimeraX: Meeting modern challenges in visualization and analysis. Protein Science 2018, 27, 14–25. [Google Scholar] [CrossRef]
  30. Laskowski, R.A.; Swindells, M.B. LigPlot+: multiple ligand–protein interaction diagrams for drug discovery. 2011. [Google Scholar] [CrossRef]
  31. Yap, C.W. PaDEL-descriptor: An open source software to calculate molecular descriptors and fingerprints. Journal of computational chemistry 2011, 32, 1466–1474. [Google Scholar] [CrossRef]
  32. Draper, N.R.; Smith, H. Applied regression analysis; John Wiley & Sons, 1998; Volume 326. [Google Scholar]
  33. Eberly, L.E. Multiple linear regression. Topics in Biostatistics 2007, 165–187. [Google Scholar]
  34. Goodarzi, M.; Dejaegher, B.; Heyden, Y.V. Feature selection methods in QSAR studies. Journal of AOAC International 2012, 95, 636–651. [Google Scholar] [CrossRef] [PubMed]
  35. Scior, T.; Medina-Franco, J.; Do, Q.-T.; Martínez-Mayorga, K.; Yunes Rojas, J.; Bernard, P. How to recognize and workaround pitfalls in QSAR studies: a critical review. Current medicinal chemistry 2009, 16, 4297–4313. [Google Scholar] [CrossRef]
  36. Bukauskas, F.F.; Jordan, K.; Bukauskiene, A.; Bennett, M.V.; Lampe, P.D.; Laird, D.W.; Verselis, V.K. Clustering of connexin 43–enhanced green fluorescent protein gap junction channels and functional coupling in living cells. Proceedings of the National Academy of Sciences 2000, 97, 2556–2561. [Google Scholar] [CrossRef]
  37. Skeberdis, V.A.; Rimkute, L.; Skeberdyte, A.; Paulauskas, N.; Bukauskas, F.F. pH-dependent modulation of connexin-based gap junctional uncouplers. The Journal of Physiology 2011, 589, 3495–3506. [Google Scholar] [CrossRef] [PubMed]
  38. Consonni, V.; Todeschini, R. Molecular Descriptors for Chemoinformatics: Volume I: Alphabetical Listing/Volume II: Appendices, References; John Wiley & Sons, 2009. [Google Scholar]
  39. Randic, M. On molecular identification numbers. Journal of Chemical Information and Computer Sciences 1984, 24, 164–175. [Google Scholar] [CrossRef]
  40. Fontes, M.S.; van Veen, T.A.; de Bakker, J.M.; van Rijen, H.V. Functional consequences of abnormal Cx43 expression in the heart. Biochimica et Biophysica Acta (BBA)-Biomembranes 2012, 1818, 2020–2029. [Google Scholar] [CrossRef]
  41. Reaume, A.G.; de Sousa, P.A.; Kulkarni, S.; Langille, B.L.; Zhu, D.; Davies, T.C.; Juneja, S.C.; Kidder, G.M.; Rossant, J. Cardiac malformation in neonatal mice lacking connexin43. Science 1995, 267, 1831–1834. [Google Scholar] [CrossRef]
  42. Fontes, M.S.C.; van Veen, T.A.B.; de Bakker, J.M.T.; van Rijen, H.V.M. Functional consequences of abnormal Cx43 expression in the heart. Biochimica et Biophysica Acta (BBA)-Biomembranes 2012, 1818, 2020–2029. [Google Scholar] [CrossRef]
  43. Martins-Marques, T.; Anjo, S.I.; Pereira, P.; Manadas, B.; Girão, H. Interacting Network of the Gap Junction (GJ) Protein Connexin43 (Cx43) is Modulated by Ischemia and Reperfusion in the Heart*[S]. Molecular & Cellular Proteomics 2015, 14, 3040–3055. [Google Scholar] [CrossRef]
  44. Laird, D.W.; Lampe, P.D. Therapeutic strategies targeting connexins. Nature reviews Drug discovery 2018, 17, 905–921. [Google Scholar] [CrossRef]
  45. Picoli, C.; Soleilhac, E.; Journet, A.; Barette, C.; Comte, M.; Giaume, C.; Mouthon, F.; Fauvarque, M.-O.; Charvériat, M. High-content screening identifies new inhibitors of connexin 43 gap junctions. ASSAY and Drug Development Technologies 2019, 17, 240–248. [Google Scholar] [CrossRef]
  46. Hollman, A. Drugs for atrial fibrillation. Digoxin comes from Digitalis lanata. BMJ: British Medical Journal 1996, 312, 912. [Google Scholar] [CrossRef]
  47. Sticherling, C.; Oral, H.; Horrocks, J.; Chough, S.P.; Baker, R.L.; Kim, M.H.; Wasmer, K.; Pelosi, F.; Knight, B.P.; Michaud, G.F. Effects of Digoxin on acute, atrial fibrillation–induced changes in atrial refractoriness. Circulation 2000, 102, 2503–2508. [Google Scholar] [CrossRef]
  48. McDonagh, T.A.; Metra, M.; Adamo, M.; Gardner, R.S.; Baumbach, A.; Böhm, M.; Burri, H.; Butler, J.; Čelutkienė, J.; Chioncel, O.; et al. 2021 ESC Guidelines for the diagnosis and treatment of acute and chronic heart failure. Eur Heart J 2021, 42, 3599–3726. [Google Scholar] [CrossRef]
  49. Moscovitz, T.; Aldrighi, J.M.; Abrahanshon, P.A.; Zorn, T.M.; Logullo, A.F.; Gebara, O.C.; Rosano, G.G.; Ramires, J.F. Repercussions of digoxin, digitoxin and estradiol on the endometrial histomorphometry of oophorectomized mice. Gynecol Endocrinol 2005, 20, 213–220. [Google Scholar] [CrossRef]
  50. Grace, A.A.; Camm, A.J. Quinidine. New England Journal of Medicine 1998, 338, 35–45. [Google Scholar] [CrossRef]
  51. Kojodjojo, P.; Kanagaratnam, P.; Segal Oliver, R.; Hussain, W.; Peters Nicholas, S. The Effects of Carbenoxolone on Human Myocardial Conduction. Journal of the American College of Cardiology 2006, 48, 1242–1249. [Google Scholar] [CrossRef]
  52. Ladipo, G.O.A.; Essien, E.E.; Andy, J.J. Complete heart block in chronic chloroquine poisoning. International journal of cardiology 1983, 4, 198–200. [Google Scholar] [CrossRef]
  53. Nevin, R.L. Mefloquine blockade of connexin 36 and connexin 43 gap junctions and risk of suicide. Biological Psychiatry 2012, 71, e1–e2. [Google Scholar] [CrossRef]
  54. De Vuyst, E.; Boengler, K.; Antoons, G.; Sipido, K.R.; Schulz, R.; Leybaert, L. Pharmacological modulation of connexin-formed channels in cardiac pathophysiology. Br J Pharmacol 2011, 163, 469–483. [Google Scholar] [CrossRef]
  55. Lamers, C. Overcoming the shortcomings of peptide-based therapeutics. Future Drug Discovery 2022, 4, FDD75. [Google Scholar] [CrossRef]
  56. Harris, A.; Locke, D. Connexins. 2009. [Google Scholar]
  57. Evans, W.H.; Boitano, S. Connexin mimetic peptides: specific inhibitors of gap-junctional intercellular communication. Biochemical Society Transactions 2001, 29, 606–612. [Google Scholar] [CrossRef] [PubMed]
  58. Borgert, C.J.; Baker, S.P.; Matthews, J.C. Potency matters: thresholds govern endocrine activity. Regulatory Toxicology and Pharmacology 2013, 67, 83–88. [Google Scholar] [CrossRef]
Preprints 70490 g001
Figure 1. Docking of examined compounds to Cx43 GJ channel. (A) Cx43 hemichannel with highlighted in red and green two adjacent subunits docking the examined compounds that are indicated in different colors on the single subunit or between subunits (see Figure 2 for details). (B) Molecular docking conformation of d-limonene. Cx43 subunits are marked with numbers in parenthesis, different Cx43 transmembrane domains are marked as M1, M2 or M3 helixes, and amino acid residues interacting with d-limonene as balls.
Figure 1. Docking of examined compounds to Cx43 GJ channel. (A) Cx43 hemichannel with highlighted in red and green two adjacent subunits docking the examined compounds that are indicated in different colors on the single subunit or between subunits (see Figure 2 for details). (B) Molecular docking conformation of d-limonene. Cx43 subunits are marked with numbers in parenthesis, different Cx43 transmembrane domains are marked as M1, M2 or M3 helixes, and amino acid residues interacting with d-limonene as balls.
Preprints 70490 g001
Preprints 70490 g002
Figure 2. Specification of molecular docking of examined compounds to Cx43 GJ hemichannel. Interacting amino acids are presented with their number in the Cx43 sequence and Cx43 subunit number in parenthesis. Hydrophobic interactions are shown as brick-red spoked arcs.
Figure 2. Specification of molecular docking of examined compounds to Cx43 GJ hemichannel. Interacting amino acids are presented with their number in the Cx43 sequence and Cx43 subunit number in parenthesis. Hydrophobic interactions are shown as brick-red spoked arcs.
Preprints 70490 g002
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.
Copyright: This open access article is published under a Creative Commons CC BY 4.0 license, which permit the free download, distribution, and reuse, provided that the author and preprint are cited in any reuse.
Prerpints.org logo

Preprints.org is a free preprint server supported by MDPI in Basel, Switzerland.

facebook logotwitter logolinkedin logo
MDPI Initiatives
Important Links
Subscribe

Disclaimer

Terms of Use

Privacy Policy

© 2024 MDPI (Basel, Switzerland) unless otherwise stated