Next Article in Journal
A Safe GDNF and GDNF/BDNF Controlled Delivery System Improves Migration in Human Retinal Pigment Epithelial Cells and Survival in Retinal Ganglion Cells: Potential Usefulness in Degenerative Retinal Pathologies
Next Article in Special Issue
Hemorphins Targeting G Protein-Coupled Receptors
Previous Article in Journal
Effectiveness of Riboflavin and Rose Bengal Photosensitizer Modified Adhesive Resin for Orthodontic Bonding
Previous Article in Special Issue
Synthesis, Biological, and Computational Evaluation of Antagonistic, Chiral Hydrobenzoin Esters of Arecaidine Targeting mAChR M1
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Pharmaceuticals
Volume 14
Issue 1
10.3390/ph14010049
pharmaceuticals-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Modifications on the Tetrahydroquinoline Scaffold Targeting a Phenylalanine Cluster on GPER as Antiproliferative Compounds against Renal, Liver and Pancreatic Cancer Cells

by
David Méndez-Luna
1,2,†,
Loreley Araceli Morelos-Garnica
1,†,
Juan Benjamín García-Vázquez
1,*,
Martiniano Bello
1,
Itzia Irene Padilla-Martínez
3,
Manuel Jonathan Fragoso-Vázquez
4,
Alfonso Dueñas González
5,
Nuria De Pedro
6,
José Antonio Gómez-Vidal
7,
Humberto Lubriel Mendoza-Figueroa
1 and
José Correa-Basurto
1,*
1
Laboratorio de Diseño y Desarrollo de Nuevos Fármacos e Innovación Biotecnológica (Laboratory for the Design and Development of New Drugs and Biotechnological Innovation), Escuela Superior de Medicina, Instituto Politécnico Nacional, Plan de San Luis y Díaz Mirón, s/n, Col. Casco de Santo Tomas, Ciudad de México 11340, Mexico
2
Departamento de Fisiología, Escuela Nacional de Ciencias Biológicas, Instituto Politécnico Nacional Zacatenco, Av. Wilfrido Massieu 399, Col. Nueva Industrial Vallejo, Alcaldía Gustavo A. Madero, Ciudad de México 07738, Mexico
3
Laboratorio de Química Supramolecular y Nanociencias, Unidad Profesional Interdisciplinaria de Biotecnología, Instituto Politécnico Nacional, Av. Acueducto s/n., Barrio La Laguna Ticomán, Ciudad de México 07340, Mexico
4
Departamento de Química Orgánica, Escuela Nacional de Ciencias Biológicas, Instituto Politécnico Nacional Prolongación de Carpio y Plan de Ayala S/N. Col. Casco de Santo Tomas, Ciudad de México 11340, Mexico
5
Genomic Medicine and Environmental Toxicology, Biomedical Research Institute, UNAM, National Cancer Institute, Av San Fernando 22, Tlalpan, Mexico City 14080, Mexico
6
Fundación MEDINA, Parque Tecnológico Ciencias de la Salud, Avenida del Conocimiento 34, 18016 Granada, Spain
7
Facultad de Farmacia, Departamento de Química Farmacéutica y Orgánica, Universidad de Granada, 18071 Granada, Spain
*
Authors to whom correspondence should be addressed.
†
These authors contributed equally to this work.
Pharmaceuticals 2021, 14(1), 49; https://doi.org/10.3390/ph14010049
Submission received: 14 October 2020 / Revised: 30 October 2020 / Accepted: 10 November 2020 / Published: 10 January 2021
(This article belongs to the Special Issue GPCRs: Ligands and beyond 2022)
Browse Figures
Review Reports Versions Notes

Abstract

:
The implementation of chemo- and bioinformatics tools is a crucial step in the design of structure-based drugs, enabling the identification of more specific and effective molecules against cancer without side effects. In this study, three new compounds were designed and synthesized with suitable absorption, distribution, metabolism, excretion and toxicity (ADME-tox) properties and high affinity for the G protein-coupled estrogen receptor (GPER) binding site by in silico methods, which correlated with the growth inhibitory activity tested in a cluster of cancer cell lines. Docking and molecular dynamics (MD) simulations accompanied by a molecular mechanics/generalized Born surface area (MMGBSA) approach yielded the binding modes and energetic features of the proposed compounds on GPER. These in silico studies showed that the compounds reached the GPER binding site, establishing interactions with a phenylalanine cluster (F206, F208 and F278) required for GPER molecular recognition of its agonist and antagonist ligands. Finally, a 3-(4,5-dimethylthiazol-2-yl)2,5-diphenyltetrazolium bromide (MTT) assay showed growth inhibitory activity of compounds 4, 5 and 7 in three different cancer cell lines—MIA Paca-2, RCC4-VA and Hep G2—at micromolar concentrations. These new molecules with specific chemical modifications of the GPER pharmacophore open up the possibility of generating new compounds capable of reaching the GPER binding site with potential growth inhibitory activities against nonconventional GPER cell models.

1. Introduction

It is widely known that steroid hormones play a crucial role in carcinogenesis and cancer progression [1]. Among this group of hormones, estrogen regulates several physiological and disease processes in humans at a transcriptional level, known as the genomic mechanism [2]. Additionally, there is a nongenomic level that involves rapid intracellular signaling, mediated through second messengers that are initiated by the activation of membrane receptors, specifically G protein-coupled receptors (GPCRs) [3]. Recently, a new estrogen membrane receptor was revealed as a pharmacological target in breast cancer, encouraging researchers to find new ligands against breast cancer [3,4]. Originally, this receptor was named as GPR30 but is now referred to as G protein-coupled estrogen receptor (GPER) by the International Union of Basic and Clinical Pharmacology (IUPHAR; www.iuphar.org/). GPER belongs to the superfamily of seven-transmembrane-spanning GPCRs, which are the most important group of pharmacological targets in the pharmaceutical market [5,6]. Specifically, GPER belongs to the class A rhodopsin-like GPCR family; its primary transduction mechanism is led through coupling with Gαi/o proteins and its secondary through the Gαs protein. GPER activates a network of transduction pathways, such as formation of the epidermal growth factor receptor (EGFR), the mitogen-activated protein kinases (MAPK), increment of the intracellular cyclic AMP (cAMP) and calcium mobilization [7,8,9]. Dihydroquinoline derivative G1, the first synthetic agonist ligand targeting GPER [10] that does not show affinity for estrogen receptors (ERα and ERβ), contains a structural scaffold homologous to other GPER antagonists (G15 and G36), and they differ only by the chemical moieties that could be responsible for their biological responses as either agonists or antagonists [10,11,12,13].
As a promising agent for use in cancer therapy, GPER could be a biological target in a wide range of cancer cells, producing positive results as quantified by their growth inhibitory activity, including renal [14], liver [15] and pancreatic [16] cancers. It is important to note that the biological activity of GPER is tissue-dependent as well as whether the cells are healthy or cancerous and the ability of tissues to express GPER [17,18]. Recent structural insights regarding the agonist/antagonist activities of GPER through in silico methods to explain how and why the activation of this receptor occurs have been applied to find new and more efficient ligands that target GPER [19,20]. In a previous study, docking and molecular dynamics (MD) simulations were combined with a molecular mechanics/generalized Born surface area (MMGBSA) approach to reveal the structural and energetic basis of the molecular recognition between agonists or antagonists on the GPER binding site. This study elucidated key residues of GPER involved in molecular recognition, including a phenylalanine cluster (F206, F208 and F278) and a polar residue (N310). The latter is considered to be a highly important residue involved in the activation or inactivation mechanism of GPER [21].
In accordance with the knowledge of key amino acid residues involved in the molecular recognition that dictates the binding mode of agonists and antagonists in the active site of GPER, we focused our efforts on proposing strategic chemical modifications on the pharmacophore that would improve molecular recognition through additional interactions and, in turn, explore additional cavities that would increase selectivity and affinity at the GPER binding site. First, we substituted the bromine atom with a phenyl ring, which was in turn substituted with electron withdrawing (m-NO2) (compound 4) or electron donating (m-OCH3) (compound 5) groups. The second chemical modification was to obtain a molecule with a tert-butyl group off the piperidine ring (compound 7), seeking to increase the hydrophobic environment that can establish interactions with the aforementioned cluster of aromatic residues. Thus, we started with a G1 analog previously prepared by our research group [22], which is subsequently referred to as G1-PABA. The chemical synthesis of G1-PABA involved modification of the p-amino acetophenone of G1 using p-aminobenzoic acid, which has shown inhibitory activity against breast cancer cells [22,23]. Once G1-PABA was synthetized, compounds 4, 5 and 7 were synthesized and tested in renal, liver and pancreatic cancer cells. The results showed growth inhibitory activities (half maximal inhibitory concentration (IC50) values) less than 50 µM. In addition, the structural and energetic properties of 4, 5 and 7 in the GPER binding site were explored by combining docking and MD simulations with a MMGBSA approach, enabling us to correlate the experimental and theoretical results.

2. Results and Discussion

2.1. Docking Calculations

To analyze the possible affinity and binding posture of the designed ligands in the GPER binding site, docking studies were performed, and it was found that the designed ligands were capable of reaching cavities in addition to those observed in previous results [20,21]. As we reported previously, the molecular recognition conducted by GPER involves essential residues, mainly aromatic in nature, comprising a phenylalanine cluster (F206, F208 and F278) and a polar residue (N310). Accordingly, in this study, in which great chemical modifications of the chemical scaffold were explored (hindrance properties, electron donation and electron withdrawing as well as hydrophobic effects), either agonist or antagonist GPER ligands examined in silico showed binding modes on GPER that were equivalent to our previous findings [20,21]. Indeed, ligand 4 (m-NO2 phenyl group) appeared to be in close contact with the phenylalanine cluster due to the ionic–π interaction exerted by the electroattractant NO2 group facing the π electron density of the phenyl groups in the phenylalanine cluster. This property was reflected in the better binding energy of 4 (−8.46 Kcal/mol) compared to 5 (−8.36 Kcal/mol) due to the latter having an electron-donating OCH3 group (Table 1). Ligand 7 showed a binding mode in which the Boc (di-tert-butyl dicarbonate) moiety was in close contact with N310 (less than 5 Å). The other important feature of 7 was that the orientation of its binding mode overlapped with the binding mode of 5, sharing almost the same contacts with the residues involved in the binding site (Table 1).
It is important to note that the pivotal residues (phenylalanine cluster and N310) were not found at the same time to be involved in nonbonding interactions with the ligands during the first docking recognition study, although it is well known that GPER contains more than one binding site formed by residues of the same chemical nature as an orthosteric binding site [20]. The resulting binding pose of the ligands obtained from the docking on GPER (for an atomistic depiction please see the Figure S1 in Supplementary Material) were employed as input coordinates for MD simulations combined with energetic calculations using the MMGBSA approach. In addition, the ligands showed acceptable physicochemical and toxicological–biological properties (absorption, distribution, metabolism, excretion and toxicity (ADME-tox)) according to Lipinski’s rules, which are the parameters required for new drug design (see Table 2). Toxicological–biological evaluations showed no effect on any property analyzed.

2.2. Molecular Dynamics Simulations

2.2.1. Evaluation of the Equilibrium of Membrane-Embedded MD GPER–Ligand Complexes

Docking studies demonstrated that the compounds studied established a different map of contacts to those previously identified due to hindrance effects [20,21] involving interactions with the phenylalanine cluster (F206, F208 and F278) and a polar residue (N310) (data not shown). However, to provide more consistent information about the prevalence of the interactions observed through docking calculations as well as to provide conformational information, the GPER–ligand complexes were membrane-embedded and submitted through 100 ns long MD simulations, accompanied by binding free energy calculations using the MMGBSA approach to obtain a system with better physiological conditions than rigid strategies.
Before performing any analysis of the structural and energetic values of the simulated membrane-embedded GPER–ligand complexes, we evaluated the MD simulation time required to perform this analysis once the equilibration stages were reached (Figure 1). Therefore, three geometrical values were evaluated to determine whether the systems reached equilibrium: area per lipid of the 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphorylcholine (POPC) membrane; root mean square deviation (RMSD) of heavy atoms; and the radius of gyration (Rg) of the GPER with respect to the starting conformation. Figure 1A shows the MD simulation time evolution of the area per lipid for the three GPER–ligand complexes. It was observed that for the first 30 ns of each MD simulation, the area per lipid values were high but then decreased and converged at 60 ns to values of 61.03 ± 0.26, 61.34 ± 0.30 and 60.76 ± 0.32 Å for GPER-4, GPER-5 and GPER-7, respectively. RMSD analysis showed that the three GPER–ligand complexes reached equilibrium at similar MD simulation times (approximately 50 ns) with RMSD values of 2.6 ± 0.12, 2.2 ± 0.17 and 2.5 ± 0.11 Å for GPER-4, GPER-5 and GPER-7, respectively (see Figure 1B). Similarly, the Rg analysis illustrated that the three systems reached equilibrium at 60 ns with Rg values of 26.5 ± 0.1, 26.4 ± 0.13 and 26.5 ± 0.12 Å for GPER-4, GPER-5 and GPER-7, respectively (see Figure 1C). Based on these observations, it was evident that the three systems reached equilibrium at 60 ns, from which root mean square fluctuation (RMSF) calculations and subsequent analyses were performed. The RMSF analysis over the equilibrated simulations (the last 40 ns of the MD production runs) revealed that the molecular flexibility of the three systems were similar (see Figure 1D). The loop (residues 20–50) that flanked the receptor binding site was the region that exhibited the highest conformational mobility for the three systems, whereas loop 5 (residues 245–251) exhibited the highest flexibility for GPER-7 only (see Figure 1D).

2.2.2. Clustering Analysis

Clustering analysis over the equilibrated time was carried out for the GPER-4, GPER-5 and GPER-7 complexes to detect all conformational states of each complex. Clustering analysis using a cut-off of 2.0 Å demonstrated that GPER-4, GPER-5 and GPER-7 complexes were characterized by one, three and five clusters of conformers, respectively, with 100%, 73% and 56% of the clusters of each GPER-ligand complex being represented by the first cluster, respectively. This analysis allowed us to obtain the most populated conformation, which was useful in evaluating the map of interactions between the ligands and GPER.

2.2.3. Principal Component Analysis (PCA)

The total fluctuation of the GPER-4, GPER-5 and GPER-7 complexes was evaluated using Cartesian principal component analysis (cPCA) (see Methods). This analysis revealed that the total motion of the GPER-4, GPER-5 and GPER-7 complexes was dispersed over 3375 eigenvectors; however, the total MD simulations were described by 15 eigenvectors. Figure 2 shows the first 15 collective modes (Figure 2A) for the three systems and their individual percentages (Figure 2B). The collective modes captured 69.5%, 69.4% and 75.3% of the total motion for the GPER-4, GPER-5 and GPER-7 complexes, respectively. The projection of GPER-4, GPER-5 and GPER-7 onto eigenvectors 1 and 2 (PC1 vs. PC2) showed that the fluctuation of the three systems was confined within these two eigenvectors. However, GPER-7 covered a larger region of phase space along the PC1 and PC2 plane than GPER-4 and GPER-5, and this latter complex covered a larger region of phase space along the PC1 plane than GPER-4. The clusters of stable states were more defined in GPER-4 than GPER-7 and GPER-5. Quantification of the total flexibility among the different protein–ligand complexes, shown in Figure 2C, was obtained by analyzing the trace of the diagonalized covariance matrix of the backbone atomic positional fluctuations. Based on this parameter, the values for GPER-7 (12.71 nm2), GPER-4 (9.35 nm2) and GPER-5 (9.94 nm2) revealed differences in conformational flexibility in the order of GPER-7 > GPER-5 > GPER-4. This order was in line with the conformational complexity identified through the cluster analysis.

2.2.4. Structural Analysis of the Empty and Bound GPER Systems

Structural details of a map of the interactions for the GPER–ligand complexes were analyzed considering the most populated conformation obtained through cluster analysis (Figure 3). A comparison of the map of interactions of the three complexes showed that 7 (Figure 3B) and 4 (Figure 3C) shared contact with five of the same residues: Gln53, Leu59, Arg122, Met133 and Leu137. Ligands 7 and 5 also shared contacts with Gln53, Leu59, Glu275, Arg122 and Gly58, and 5 and 4 shared contacts with Gln53, Leu59, Leu119, His120 and Arg122. However, only three residues (Gln53, Leu59 and Arg122) were shared by all three compounds, indicating that small modifications in the pharmacophore group can modify the map of interactions.
Overall, this analysis revealed that the three compounds were mostly stabilized by hydrophobic residues. Residues belonging to the phenylalanine cluster (F206, F208 and F278) were identified as important stabilizers of the complexes with GPER [20,21]. Among the three complexes, only 4 (F206 and F208) and 5 (F278) were able to establish contact with some members of this cluster. With respect to polar interactions, only 7 and 4 were able to establish hydrogen bonds and salt bridges through their polar moieties at the GPER binding site. Ligand 7 formed a hydrogen bond with Gln53 and a salt bridge with Arg122 through its carboxyl group (Figure 3B), whereas 4 was stabilized by two hydrogen bonds through its carboxyl group with His120 and its NO2 group with Cys207 (Figure 3C). The latter polar interaction has been reported to be responsible for inducing the activated state [20]. Interestingly, the lack of a bromine (Br) atom in these compounds means that the polar interactions observed in G1 and G15 through Asn307 or Asn310, respectively [21], were not present. Instead, these two residues appeared in GPER-7 and GPER-5, forming polar interactions through backbone atoms.

2.2.5. Absolute Binding Free Energy Calculations

The absolute binding free energy (ΔGbind) was calculated to measure the energetic contributions in terms of noncovalent energies at the binding site of the GPER-4, GPER-5 and GPER-7 complexes (see Methods). This analysis revealed that the ΔGbind values for the three GPER–ligand systems, considering the entropic component, were thermodynamically favorable for all complexes (Table 3). The main energetic contribution to ΔGbind was guided by the nonpolar contributions (ΔEnonpolar = ΔEvwd + ΔGnpol,sol), whereas the polar contributions (ΔEpolar = ΔEele + ΔGpol,sol) displayed a thermodynamically unfavorable performance that contrasted with the binding process for the three complexes. Based on the entropic values (Table 3), it was observed that the three complexes exhibited a significant conformational reduction upon complex formation, notably impacting the ΔGbind value, suggesting that this molecular recognition was guided through an unfavorable entropic component. The GPER-4 system exhibited a more energetically favorable ΔGbind than the GPER-7 and GPER-5 complexes, suggesting that the NO2-containing compound established a more energetic map of interactions into the GPER binding site than the other two compounds.

2.3. Antiproliferative Assays

GPER plays a pivotal role in triggering cancer development, and most studies in the literature have focused on its overexpression in breast cancer. The current findings relate to the properties of GPER, particularly that its activity depends specifically on the expression tissue and the cognate ligands (either agonist or antagonist) that activate or inactivate it. Thus, new possibilities have been generated in the evaluation of new compounds targeting GPER in cancer cells, unlike classical models of activity evaluation [24,25,26]. However, there are other types of cancer in which GPER is present, and it could be an interesting biological target in some other cancers, such as renal [14], liver [15] and pancreatic [16]. The evaluation of growth inhibition of the synthetized compounds was conducted using a 3-(4,5-dimethylthiazol-2-yl)2,5-diphenyltetrazolium bromide (MTT) assay in selected cancer cell lines. RCC4-VA, RCC4-VHL, MIA Paca-2 and Hep G2 were used to investigate the possible antiproliferative effects of these three compounds as these are diverse cancer cells that all have GPER in common. As shown in Table 4, MIA Paca-2 was the most sensitive cancer cell line to the compounds tested. In this cell line, all tested compounds had an antiproliferative effect in a concentration-dependent manner. The most effective ligand in the MIA Paca-2 cell line was 7, which has a Boc protecting group on the amine of the piperidine ring of the pharmacophore of GPER (see Figure 4). In this cell line, the expression pattern of estrogen receptors is known, and antiestrogen therapy is used to block ER activity in order to treat the cancer [27,28].
The GPER signaling pathway was studied and its role in both healthy and cancerous processes were examined, particularly with regard to the very important role in the tumor microenvironment for renal, liver and pancreatic cancers, which could be affected by the synthetized compounds. Interestingly, only two compounds showed an effect in mediating growth inhibition in the RCC4-VA cell line (7 and 5), but with higher concentrations than those assayed in MIA Paca-2 cells (see Figure 5). In this cell line, it is important to note that the trend of 4, which showed a non-dose-dependent behavior (evident at 6.3 µM concentration), might be due to intrinsic properties of the cell to adapt to increased concentrations of the tested compounds. In this case, the cell could be increasing the expression of proteins involved into externalize the drugs (e.g., P-glycoprotein (PGP) ATPase) or another mechanism, which conferred to the cell the capability to not respond to the treatments assayed. This correlated with the result yielded by 7 assayed at the highest concentration, although the result was slightly lower than the previous concentration (25 µM). This pharmacological phenomenon is named as a biphasic pattern and has been observed for other cancer cells (https://www.ncbi.nlm.nih.gov/pmc/articles/PMC4560104/). In contrast, in the RCC4-VHL cell line, all compounds showed activity at the concentrations assayed (data not shown). This result is promising due to the molecular differences by which these cell lines are employed, noting that VHL/HIF deregulated signaling and the distinctive molecular feature of RCC4-VA cells responsible for the high rate of tumors characteristic of Von Hippel-Lindau syndrome [29,30]. Finally, in Hep G2 (human liver carcinoma) cells, only compound 5 showed activity (Figure 6) with almost the same concentration found for RCC4-VA. This result might unmask an axis that involves the participation of GPER/HIF-1α/vascular endothelial growth factor (VEGF) as has been reported elsewhere [31]. Additionally, all cell models tested have been implicated in processes in which the participation of molecular mediators under hypoxic conditions are initiated once GPER activation has occurred [32,33]. These findings suggest that GPER is a potential pharmacological target to develop new therapeutic strategies for cancer.

3. Materials and Methods

3.1. Docking, Physicochemical and ADME-Tox Properties Analysis

The potential interaction between the designed compounds and GPER was analyzed using molecular docking studies. First, the two-dimensional (2D) chemical structures of the ligands were drawn with the ChemBioDraw Ultra 12.0 program [34]. The Z-matrices of the structures were generated using the Gauss View 5.0 program [35] to assign the correct interatomic distance and hydrogen types. This step was undertaken through energetic and geometric minimization under an AM1 semiempirical method that employed the quantum chemistry package Gaussian 09 [36]. The GPER three-dimensional (3D) model previously generated and validated by our research group was used in the molecular docking simulations [20]. These simulations were carried out using the AutoDock 4.2.6 software [37]. Polar hydrogen atoms were added to the polar atoms of the protein, and Kollman charges were assigned. Gasteiger charges were assigned for ligands, and a grid-based procedure was used to generate the affinity maps delimiting the protein binding site with the following grid box parameters: XYZ of 60 Å3 over the Cα atom of the N310 residue (focused procedure) with a grid spacing of 0.375 Å3. Scoring sampling used the Lamarckian genetic algorithm considering the following as search parameters: a randomized initial population of 100 individuals and a maximum number of energy evaluations of 1 × 107 cycles. The ligand–receptor interactions were analyzed with PyMOL v0.99 graphical viewer [38]. The physicochemical and ADME-tox properties were analyzed using the Molinspiration server [39] and DataWarrior software [40], respectively.

3.2. Embedding of the GPER–Ligand Complexes into the Membrane

Orientation of the GPER-4, GPER-5 and GPER-7 complexes with respect to the POPC membrane was performed using the Orientations of Proteins in Membranes (OPM) server [41]. Rectangular pre-equilibrated POPC bilayers were built for each complex using the membrane builder tool CHARMM [42] with dimensions 110.437 × 110.437 × 122.485 Å (xyz). GPER–ligand complexes were embedded into the POPC bilayer using the replacement method in CHARMM [43,44]. POPC bilayers constituted approximately 314 POPC phospholipids with 154 and 160 on the outer and the inner leaflets, respectively, and they were solvated using the TIP3 water model and neutralized with an ionic strength of 0.15 M by NaCl using the solvation and autoionize module [43].

3.3. MD Simulations of the Membrane-Embedded GPER–Ligand Complexes

The MD simulations were performed with pmemd.cuda AMBER 12 executable [45]. Topologies for the GPER–ligand complexes were built using the Leap module, employing the ff99SB, Lipid11 and GAFF force fields [46,47]. These complexes were minimized and equilibrated with the sander module through the following steps. The complexes were energy minimized through 10,000 steps with position restraints on the receptor and lipids to allow relaxation of the water molecules. Complexes were slowly heated under constant temperature and constant volume (NVT) ensemble through two sequential runs from 0 to 300 K for 1 ns, maintaining the receptor atom and lipid restraints. Afterwards, a 1 ns simulation under the NPT ensemble at 300 K and 1 bar pressure with the heavy atoms restrained was carried out to equilibrate the system, followed by 1 ns of equilibration with the fully unrestricted system. The equilibration runs were followed by 100 ns MD simulations, which were run for each system without position restraints under periodic boundary conditions (PBCs) using a NPT ensemble at 300 K and 1 bar pressure. The electrostatic interactions were computed by the particle mesh Ewald (PME) method [48]. A 10 Å cut-off was used for the van der Waals interactions, and the SHAKE algorithm was used to constrict the bonds between the heavy atoms and hydrogen atoms [49]. The temperature was conserved using Langevin dynamics, and pressure was controlled using a semi-isotropic constant surface tension to preserve a specific area per lipid. The pressure was maintained at 1 bar, and the pressure coupling constant was set to 1 ps with a collision frequency of 1 ps−1. The time step was set to 2 fs, and the coordinates were saved for analysis every 1 ps.

3.4. Analysis of MD Simulations

The RMSD and Rg were calculated only for the α-carbon atoms and the receptor, respectively, after removing the overall translational and rotational motion. The area per lipid was calculated using the following equation: area per lipid = (box X length × box Y length)/number of phospholipids per layer. The most populated conformers, which are by definition the most stabilized systems present in a physiological context, were obtained through a RMSD conformational clustering analysis using the g_cluster command with the GROMOS algorithm and a cut-off of 0.25 nm [50] to evaluate the conformational features linked to the most stable GPER–ligand complexes. Images and structural representations were prepared using PyMOL v0.99 [38].

3.5. Principal Component Analysis

Principal component analysis, the essential dynamics (ED) method [51] and the g_covar and g_anaeig GROMACS 4.5.3 package utilities [52,53] were implemented to predict the collective motion of atoms. The symmetric matrix was diagonalized to form a covariance matrix, which was used to obtain a set of eigenvectors and eigenvalues. The essential subspace relevant to the GPER–ligand systems was constructed by projecting the PCA along the most significant eigenvectors.

3.6. Calculation of Absolute Binding Free Energies

The binding free energy was computed for the GPER–ligand complexes using the MMGBSA method and the single trajectory approach [54,55,56,57,58,59,60,61]. To this end, the water molecules and counterions were removed from the trajectory, and the energy computations were performed by selecting 400 snapshots at time intervals of 100 ps from the last 40 ns of the simulation run for each system. The binding free energy (ΔGbind) of each GPER–ligand complex was calculated as follows:
∆ G b i n d = G c o m p l e x − G r e c e p t o r − G l i g a n d
= ∆ E M M + ∆ G G B + ∆ G S A − T ∆ S
where ΔEMM corresponds to the gas-phase interaction energy for the protein–ligand complex, which includes the van der Waals (ΔEvdw) and the electrostatic (ΔEele) interaction energies; ΔGGB and ΔGSA are the electrostatic and nonpolar contributions to desolvation upon ligand binding, respectively; and −TΔS is the entropy contribution, which is a result of the structural changes in the degrees of freedom between the free and bound partners forming the protein–ligand complex. Because −TΔS was considered in the computation of the binding free energy calculations, the values we obtained using the MMGBSA approach could be named as absolute binding free energies, ΔGbind.

3.7. Calculation of Entropy Contributions

The entropy contribution was evaluated for each complex using the MMPBSA.py module implemented in Amber 12. Due to the large memory required for this calculation, only 20 snapshots were selected for each GPER–ligand complex at time intervals of 2 ns from the last 40 ns of each simulation run.

3.8. Chemical Synthesis

All reagent-grade chemicals were obtained from commercial suppliers and were used without further purification. Reactions were monitored by thin-layer chromatography (TLC) on aluminum-backed sheets with silica gel 60 GF254 (HX805651) and a fluorescent indicator (visualized with UV light of 254 nm). Flash chromatography was performed using silica gel 60 (230–400 mesh). Melting points (mp) were determined on an Electrothermal IA 91000 apparatus (Electrothermal, Bibby Scientific, Staffordshire, ST15 OSA, UK) and were uncorrected. 1H and 13C NMR spectra were obtained on a Varian Mercury 300 MHz spectrometer or a Bruker Avance III 750 MHz spectrometer using deuterated dimethyl sulfoxide (DMSO-d6) or deuterated chloroform (CDCl3) as the solvent and tetramethylsilane (TMS) as the internal standard. Chemical shifts (δ) are reported in ppm downfield from the internal standard, and coupling constants are reported in Hertz (Hz) (Supplementary Material). Electrospray ionization high-resolution mass spectrometry positive mode (ESI-HRMS) was performed with an Agilent 6545 QTOF LC/MS instrument (Agilent Technologies, Santa Clara, CA, USA).

3.8.1. Synthesis of (3aS, 4R, 9bR)-4-(6-Bromobenzo[d][1,3]dioxol-5-yl)-3a,4,5,9b-tetrahydro-3H-cyclopenta[c]quinoline-8-carboxylic acid, G1-PABA 1

G1-PABA (1) and the tert-butyl derivative (compound 7) were obtained following the previous method reported by our research group [22,23].

3.8.2. Synthesis of (3aS, 4R, 9bR)-4-(6-(3-Nitrophenyl)benzo[d][1,3]dioxol-5-yl)-3a,4,5,9b-tetrahydro-3H-cyclopenta[c]quinoline-8-carboxylic acid, 4

Compound 4 was obtained by following the previously reported method [62]. G1-PABA (0.020 g, 0.048 mmol), m-nitro phenylboronic acid (0.0096 g, 0.057 mmol), PEG2000 (0.060 g, 0.10 mmol) and Pd(AcO)2 (0.004 g, 0.017 mmol) were added to a solution of K2CO3 (0.0266 g, 0.19 mmol) and MeOH/H2O (7 mL) in a 1:1 mixture with constant stirring (see Scheme 1). The reaction mixture was heated to reflux at 70 °C for 5 h under N2 atmosphere. The mixture was then enabled to cool to room temperature, and the resulting mixture was extracted with EtOAc (3 × 25 mL). The organic phase was washed with distilled water and brine, dried over anhydrous sodium sulfate (Na2SO4) and concentrated under reduced pressure. The residue was purified by flash chromatography using a mixture of hexane/EtOAc (7:3) as the eluent, affording 55% yield of compound 4 as a clay solid, Rf = 0.42 (hexane/EtOAc, 7:3), mp = 184–185 °C. HPLC purity = 97.19 %. 1H NMR (750 MHz, CDCl3) δ 8.22 (d, 1H, 7.5 Hz, H-4′′), 8.14 (bs, 1H, H-2′′), 7.68 (d, 1H, J = 7.5 Hz, H-7), 7.67 (s, 1H, H-9), 7.60 (d, 1H, J = 7.5 Hz, H6′′), 7.57 (t, 1H, J = 7.5 Hz, H5′′), 7.26 (s, 1H, H-7′), 7.23 (s, 1H, H-4′), 6.68 (s, 1H, NH), 6.56 (d, 1H, J = 7.5 Hz, H-6), 6.06 and 6.05 (AB, 2H, H-2′), 5.88 (bm, 1H, H-1), 5.69 (bm, 1H, H-2), 4.29 (d, 1H, J = 3.7 Hz, H-4), 4.22 (d, 1H, J = 6.0 Hz, H-9b), 2.67 (m, 1H, H-3a), 1.71 (m, 1H, H-3 down), 1.44 (m, 1H, H-3 up). 13C NMR (187.5 MHz, CDCl3) δ 167.9 (CO2H), 150.4 (C-7a’), 148.3 (C-3′’), 148.2 (C-3′), 146.9 (C-5a), 142.5 (C-1′′), 135.6 (C-6′′), 133.9 (C-2), 133.1 (C-6′), 132.6 (C-9), 132.0 (C-1), 130.6 (C-4′), 129.6 (C-5′′), 129.5 (C-2′′), 124.8 (C-7), 124.4 (C-4′′), 122.6 (C-9a), 119.4 (C-8), 115.3 (C-6), 110.4 (C-4′), 107.1 (C-7′), 101.8 (C-2′), 53.6 (C-4), 45.6 (C-9b), 39.0 (C-3a), 38.9 (C-3). HRMS (ESI) calculated for ([C26H20N2O6] + H)+: 457.1400; found: 457.13824.

3.8.3. Synthesis of (3aS, 4R, 9bR)-4-(6-(3-Methoxyphenyl)benzo[d][1,3]dioxol-5-yl)-3a,4,5,9b-tetrahydro-3H-cyclopenta[c]quinoline-8-carboxylic acid, 5

Compound 5 was obtained through a Suzuki–Miyaura cross-coupling reaction [63]. A solution of K2CO3 (0.0266 g, 0.19 mmol) in dimethylformamide (DMF, 1.5 mL) and water (1.5 mL) was stirred for 15 min under N2 atmosphere. Subsequently, G1-PABA (0.020 g, 0.048 mmol), methoxy-nitro phenyl boronic acid (0.0096, 0.057 mmol), triphenylphosphine (0.0015 g, 0.0057 mmol) and Pd(AcO)2 (0.0004 g, 0.017 mmol) were added. The reaction mixture was heated to 100 °C for 24 h (see Scheme 2). Once the reaction was complete, it was allowed to cool to room temperature, and the resulting mixture was extracted with EtOAc (3 × 25 mL). The organic layer was washed with distilled water and brine, dried over anhydrous Na2SO4 and concentrated under reduced pressure. The obtained residue was purified by flash chromatography with a mobile phase consisting of hexane/EtOAc (8:2), affording a 13% yield of compound 5 as a gray solid, Rf = 0.5 (hexane/EtOAc 8:2), mp = 200 °C. HPLC purity = 87.31 %. 1H NMR (300 MHz, CDCl3) δ 7.69 (d, 1H, J = 7.1 Hz, H-7), 7.68 (s, 1H, H-9), 7.31 (t, 1H, J = 8.8 Hz, H-5′′), 7.18 (s, 1H, H-7′), 6.88 (dd, 1H, J = 7.9, 2.9 Hz, H-6′′), 6.82 (d, 1H, J = 7.6 Hz, H-4′′), 6.78 (d, 1H, J = 1.7 Hz, H-2′′), 6.71 (s, 1H, H-4′), 6.54 (d, 1H, J = 8.8 Hz, H-6), 6.02 (m, 2H, H-2′), 5.90 (m, 1H, H-1), 5.69 (m, 1H, H-2), 4.70 (d, 1H, J = 3 Hz, H-4), 3.81 (s, 3H, OMe), 3.81 (dd, 1H, J = 17.6, 8.8 Hz, H-9b), 2.80 (dt, 1H, J = 8.3, 3 Hz, H-3a), 2.64 (dd, J = 15.2, 8.3 Hz, H-3 down), 1.97 (m, 1H, H-3 up). 13C NMR (75 MHz, CDCl3) δ 171.6 (CO2H), 159.3 (C-1′′), 150.6 (C-5a), 147.4 (C-3′a), 142.0 (C-7′a), 146.2 (C-1′′), 135.1 (C-5′), 132.6 (C-9a), 124.7 (C-6′), 118.8 (C-8), 101.2 (C-2′), 131.8 (C-7), 129.3 (C-9), 128.0 (C-5′′), 106.5 (C-7′), 112.3 (C-6′′), 121.7 (C-4′′), 115.2 (C-2′′), 110.3 (C-4′), 114.9 (C-6), 133.7 (C-1), 130.6 (C-2), 53.2 (C-4), 55.2 (OMe), 45.4 (C-9b), 44.5 (C-3a), 31.9 (C-3). HRMS (ESI) calculated for ([C27H23NO5] + H)+: 442.1654; found: 442.1646.

3.8.4. Synthesis of (3aS, 4R, 9bR)-4-(6-Bromobenzo[d][1,3]dioxol-5-yl)-5-(tert-butoxycarbonyl)-3a,4,5,9b-tetrahydro-3H-cyclopenta[c]quinoline-8-carboxylic acid, 7

To a cold solution of G1-PABA in DMF (3.0 mL), triethylamine (336 μL, 2.413 mmol) was added dropwise and with continuous stirring for 30 min at the same temperature. Subsequently, di-tert-butyl dicarbonate (Boc2O, 0.5268 g, 2.413 mmol) was added, the reaction mixture was allowed to reach room temperature and continued overnight with vigorous stirring (see Scheme 3). The progress of the reaction was monitored by TLC, and the resulting mixture was extracted with EtOAc (3 × 25 mL). The organic phase was washed with water and brine, dried over anhydrous Na2SO4 and concentrated to dryness. The residue obtained was purified by flash chromatography using a mobile phase of hexane/EtOAc (9:1), affording a 54% yield of compound 7 as a pearl yellow solid, Rf = 0.63 (hexane/EtOAc 9:1), mp = 184.3 °C. HPLC purity = 95.50 %. 1H NMR (300 MHz, CDCl3) δ 7.88 (d, 1H, J = 1.5 Hz, H-9), 7.81 (dd, 1H, J = 8.4, 1.5 Hz, H-7), 7.61 (d, 1H, J = 8.8 Hz, H-6), 6.94 (s, 1H, H-7′), 6.29 (s, 1H, H-4′), 6.17 (m, 1H, H-1), 6.05 (d, 1H, J = 9.7 Hz, H-4), 5.85 (AB, 2H, H-2′), 5.64, (m, 1H, H-2), 3.85 (bd, 1H, J = 8.4 Hz, H-9b), 3.44 (q, 1H, J = 8.3 Hz, H-3a), 2.20 (dd, 1H, J = 16.3, 8.3 Hz, H-3 down), 1.79 (dd, 1H, J = 16.5, 6.8, Hz, H-3 up), 1.40 (s, 9H, OtBu). 13C NMR (75 MHz, CDCl3) δ 165.6 (CO2H), 153.3 (C-carbamate), 147.0 (C-7′a), 146.9 (C-3′a), 141.4 (C-5a), 134.2 (C-9a), 132.2 (C-5′), 114.4 (C-6′), 81.7 (C-tBu), 129.0 (C-9), 127.4 (C-7), 124.7 (C-6), 112.3 (C-7′), 107.9 (C-4′), 133.0 (C-1), 56.0 (C-4), 101.6 (C-2′), 131.8 (C-2), 43.5 (C-9b), 41.6 (C-3a), 34.9 (C-3), 28.2 (3Me-tBu). HRMS (ESI) calculated for [C25H24BrNO6Na] + [M + Na]+: 536.0679; found: 536.0669.

3.9. Purity Analysis by HPLC

All HPLC analyses were performed using an Agilent 1260 infinity series liquid chromatograph (Agilent Technologies, Palo Alto, CA, USA) equipped with a quaternary pump delivery system (G1311B), robotic autosampler (G1316A), column thermostat (G1316A) and multiwavelength UV detector (G1315C), and the results were analyzed by OpenLab CDS EZChrom. The mobile phase with isocratic elution consisted of a mixture of (A) 0.2% acetic acid in water (v/v) at pH 3.0 and (B) CAN (acetonitrile) in a proportion of 40% A and 60% B. The system was ready for the next injection without further equilibration. The mobile phase was prepared fresh daily using deionized water that was filtered (0.22 μm) and degassed prior to use. The flow rate was 0.5 mL/min. The proposed method utilized a Zorbax SB-C18 column (5 μm, 4.6 × 150 mm) at a temperature of 25 °C. Acceptable purity was delimited when all compounds were confirmed by ≥95%.

3.10. Cell Culture

Four cell lines were used for this study. The human renal cell carcinoma cell line RCC4 was stably transfected with an empty vector, pcDNA3 (ECACC N-03112702; called RCC4-VA), which confers neomycin resistance, or with pcDNA3-VHL (ECACC N-0312703; called RCC4-VHL), which confers neomycin resistance and encodes the VHL tumor suppressor gene, pVHL, which is associated with resistance to chemotherapy. RCC4 plus the empty vector served as a negative control cell line to study the effects of pVHL expression by pcDNA3-VHL. The other two cell lines used in this study were the pancreatic cancer cell line MIA Paca-2 and the liver cancer cell line Hep G2. Cell culture medium, fetal bovine serum (FBS), L-glutamine, sodium pyruvate, minimum essential medium (MEM) nonessential amino acids, penicillin-streptomycin, triple express, phosphate-buffered saline (PBS) and bovine insulin were purchased from Life Technologies (Gibco, Invitrogen Corporation, CA, USA). MTT and methyl methanesulfonate (MMS) were purchased from Sigma Chemical Co. (St. Louis, MO, USA). Cells were maintained and plated with a robotic cell culture system, SelecT (TAP Biosystems). Hep G2 and MIA Paca-2 cells were obtained from the American Type Culture Collection (ATCC, Manassas, VA, USA). Hep G2 cells (human liver carcinoma, CCL-8065) were grown in ATCC-formulated Eagle’s MEM with 10% qualified FBS, 2 mM L-glutamine, 1 mM sodium pyruvate and 100 μM MEM nonessential amino acids. The MIA Paca-2 cell line (ATCC CRL-1420), a fibroblast primary pancreatic adenocarcinoma cell line that contains mutated K-RAS, P16 and P53 genes, was grown in Dulbecco’s modified Eagle’s medium supplemented with 10% fetal bovine serum, 2.5% fetal horse serum, 1% L-glutamine and 1% penicillin–streptomycin. RCC4-VA and RCC4-VHL cell lines were maintained in Dulbecco’s modified Eagle’s medium, supplemented with 10% FBS, 0.01% L-glutamine 200 mM, 0.01% penicillin–streptomycin (GIBCO) and 0.001% Geneticin G418 0.5 mg/mL (GIBCO). Cell cultures were grown at 37 °C in a humidified atmosphere of 5% CO2 and 95% air.

3.11. Antiproliferative Assays

The rate of MTT reduction is an indicator of the functional integrity of mitochondria and hence of cellular viability [64,65]. For this assay, cells were treated with compounds for 72 h. These cells were then seeded into 96-well plates at a density of 103 cells/well for Hep G2, MIA Paca-2, RCC4-VA and RCC4-VHL. Three microliters of each pure compound were dispensed into 597 μL of fresh medium (see SM for purity statement). From this mixture, 200 μL was transferred to three different cell plates with the Biomek FX (Beckman Coulter) automated liquid handling system. MMS was used as a positive control and 0.5% DMSO was used as a negative control. The maximum concentration of DMSO was 0.5% to minimize any background solvent toxicity. Compounds were tested in triplicate with 10-point dilutions. After compounds and controls had been added, the plates were incubated at 37 °C in a 5% CO2 incubator for 72 h. After this time, the MTT solution was prepared at 5 mg/mL in PBS 1× and then diluted to 0.5 mg/mL in MEM without phenol red. The sample solution in each well was skimmed off, and 100 µL of MTT dye was added to each well. The plates were gently shaken and incubated for 3 h at 37 °C in a 5% CO2 incubator. The supernatant was removed, and 100 µL of 100% DMSO was added. The plates were gently shaken to solubilize the formed formazan. The absorbance was measured using a VictorTM multireader at a wavelength of 570 nm. One-way ANOVA test was performed to determine any significant difference between each compound treatment and the DMSO control; p > 0.05 was considered as statistically significant.

4. Conclusions

Bioinformatic tools, such as molecular docking and molecular dynamics simulations, provide relevant structural information, which allows the characterization of structural features that carry out molecular recognition. This offers a broad panorama in the design of selective drugs for the treatment of diseases. Under this context, we designed three compounds that showed notable inhibitory activity against GPER when tested in nonconventional cell models. According to the chemical data of the structural modifications of the GPER pharmacophore, we were able to identify the key entities (bromine atom and NH group of the piperidine ring) that could be functionalized to improve the activity of the compounds designed to target this receptor. We found a very good correlation of the in vitro experimental results with the in silico approximations on GPER. Our results suggest that GPER has the capacity to recognize different structural scaffolds of flat and aromatic character, specifically ligands with simple chemical modifications that improve the stereoelectronic properties of the pharmacophore. A comparison of the results obtained in this study with our previous structural insights demonstrates the great capability of GPER to bind ligands of different chemical structures, specifically ligands with simple chemical modifications over their pharmacophores. Nevertheless, additional assays, for example, binding assays and intracellular Ca2+ mobilization, will be crucial to corroborate our hypothesis. The appropriate use of rational in silico methods and tools could yield new compounds for the promising list of ligands to be used as potential drugs in cancer therapy.

Supplementary Materials

The explanation will be as follow: The Supplementary Material for this paper is available online at https://www.mdpi.com/1424-8247/14/1/49/s1, 1H and 13C NMR spectra data, HPLC data to corroborate the purity statement, Figure S1. Binding pose yielded by docking analysis for (A) Compound 4, (B) Compound 5 and (C) Compound 7. H bond interactions are depicted as red dashes, and Force Field parameters employed in MD and Docking simulations.

Author Contributions

It would be as follow: Conceptualization, J.B.G.-V., M.J.F.-V. and J.C.-B.; methodology, D.M.-L., L.A.M.-G., M.B., A.D.G., I.I.P.-M., N.D.P., J.A.G.-V., H.L.M.-F.; software, D.M.-L., M.B.; validation, J.B.G.-V., M.J.F.-V., D.M.-L. and N.D.P.; formal analysis, D.M.-L., J.B.G.-V., M.J.F.-V., H.L.M.-F. and J.C.-B.; investigation, D.M.-L. and L.A.M.-G.; resources, J.B.G.-V., J.C.-B.; data curation, N.D.P.; writing—original draft preparation, D.M.-L., M.B., J.C.-B.; writing—review and editing, A.D.G. and J.C.-B.; visualization, D.M.-L., M.B., M.J.F.-V.; supervision, J.B.G.-V., M.J.F.-V., H.L.M.-F., J.C.-B.; project administration, J.C.-B.; funding acquisition, J.C.-B. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by grants from the CONACYT (Grants: CB-254600, APN-782 and SEP-CONACYT-ANUIES-ECOS Francia: 296636), Instituto Politécnico Nacional (Grant: Proyectos Insignia IPN-2015) and COFAA-SIP/IPN.

Acknowledgments

We gratefully acknowledge to Centro de Nanociencias y Micro y Nanotecnologías del IPN and MC Alma Ortiz and MC Alan Estrada for the assistance in HPCL assays.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Allred, D.C. Issues and updates: Evaluating estrogen receptor-α, progesterone receptor, and HER2 in breast cancer. Mod. Pathol. 2010, 23, S52–S59. [Google Scholar] [CrossRef] [Green Version]
  2. DeRoo, B.J.; Korach, K.S. Estrogen receptors and human disease. J. Clin. Investig. 2006, 116, 561–570. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Prossnitz, E.R.; Arterburn, J.B.; Sklar, L.A. GPR30: A G protein-coupled receptor for estrogen. Mol. Cell. Endocrinol. 2007, 265, 138–142. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Olde, B.; Leeb-Lundberg, L.M.F. GPR30/GPER1: Searching for a role in estrogen physiology. Trends Endocrinol. Metab. 2009, 20, 409–416. [Google Scholar] [CrossRef] [PubMed]
  5. Pierce, K.L.; Premont, R.T.; Lefkowitz, R.J. Signalling: Seven-transmembrane receptors. Nat. Rev. Mol. Cell Biol. 2002, 3, 639–650. [Google Scholar] [CrossRef]
  6. Dorsam, R.T.; Gutkind, J.S. G-protein-coupled receptors and cancer. Nat. Rev. Cancer 2007, 7, 79–94. [Google Scholar] [CrossRef]
  7. Revankar, C.M.; Cimino, D.F.; Sklar, L.A.; Arterburn, J.B.; Prossnitz, E.R. A Transmembrane Intracellular Estrogen Receptor Mediates Rapid Cell Signaling. Science 2005, 307, 1625–1630. [Google Scholar] [CrossRef] [Green Version]
  8. Thomas, P.; Pang, Y.; Filardo, E.J.; Dong, J. Identity of an Estrogen Membrane Receptor Coupled to a G Protein in Human Breast Cancer Cells. Endocrinology 2005, 146, 624–632. [Google Scholar] [CrossRef] [Green Version]
  9. Filardo, E.J.; Quinn, J.A.; Bland, K.I.; Frackelton, A.R. Estrogen-Induced Activation of Erk-1 and Erk-2 Requires the G Protein-Coupled Receptor Homolog, GPR30, and Occurs via Trans-Activation of the Epidermal Growth Factor Receptor through Release of HB-EGF. Mol. Endocrinol. 2000, 14, 1649–1660. [Google Scholar] [CrossRef]
  10. Bologa, C.G.; Revankar, C.M.; Young, S.M.; Edwards, B.S.; Arterburn, J.B.; Kiselyov, A.S.; Parker, M.A.; Tkachenko, S.E.; Savchuck, N.P.; Sklar, L.A.; et al. Virtual and biomolecular screening converge on a selective agonist for GPR30. Nat. Chem. Biol. 2006, 2, 207–212. [Google Scholar] [CrossRef]
  11. Burai, R.; Ramesh, C.; Shorty, M.; Curpan, R.; Bologa, C.G.; Sklar, L.A.; Oprea, T.I.; Prossnitz, E.R.; Arterburn, J.B. Highly efficient synthesis and characterization of the GPR30-selective agonist G-1 and related tetrahydroquinoline analogs. Org. Biomol. Chem. 2010, 8, 2252–2259. [Google Scholar] [CrossRef]
  12. Dennis, M.K.; Burai, R.; Ramesh, C.; Petrie, W.K.; Alcon, S.N.; Nayak, T.K.; Bologa, C.G.; Leitao, A.; Brailoiu, E.; Deliu, E.; et al. In vivo effects of a GPR30 antagonist. Nat. Chem. Biol. 2009, 5, 421–427. [Google Scholar] [CrossRef] [Green Version]
  13. Dennis, M.K.; Field, A.S.; Burai, R.; Ramesh, C.; Petrie, W.K.; Bologa, C.G.; Oprea, T.I.; Yamaguchi, Y.; Hayashi, S.-I.; Sklar, L.A.; et al. Identification of a GPER/GPR30 antagonist with improved estrogen receptor counterselectivity. J. Steroid Biochem. Mol. Biol. 2011, 127, 358–366. [Google Scholar] [CrossRef] [Green Version]
  14. Guan, B.-Z.; Yan, R.-L.; Huang, J.-W.; Rui-Ling, Y.; Zhong, Y.-X.; Chen, Y.; Liu, F.-N.; Hu, B.; Huang, S.-B.; Yin, L.-H. Activation of G protein coupled estrogen receptor (GPER) promotes the migration of renal cell carcinoma via the PI3K/AKT/MMP-9 signals. Cell Adhes. Migr. 2018, 12, 109–117. [Google Scholar] [CrossRef] [Green Version]
  15. Chaturantabut, S.; Shwartz, A.; Evason, K.J.; Cox, A.G.; Labella, K.; Schepers, A.G.; Yang, S.; Acuña, M.; Houvras, Y.; Mancio-Silva, L.; et al. Estrogen Activation of G-Protein–Coupled Estrogen Receptor 1 Regulates Phosphoinositide 3-Kinase and mTOR Signaling to Promote Liver Growth in Zebrafish and Proliferation of Human Hepatocytes. Gastroenterology 2019, 156, 1788–1804. [Google Scholar] [CrossRef] [Green Version]
  16. Cortes, E.; Sarper, M.; Robinson, B.; Lachowski, D.; Chronopoulos, A.; Thorpe, S.D.; Lee, D.A.; Hernández, A.E.D.R. GPER is a mechanoregulator of pancreatic stellate cells and the tumor microenvironment. EMBO Rep. 2019, 20, 46556. [Google Scholar] [CrossRef]
  17. Wang, C.; Lv, X.; Jiang, C.; Davis, J.S. The putative G-protein coupled estrogen receptor agonist G-1 suppresses proliferation of ovarian and breast cancer cells in a GPER-independent manner. Am. J. Transl. Res. 2012, 4, 390–402. [Google Scholar]
  18. Wang, C.; Lv, X.Y.; He, C.; Hua, G.; Tsai, M.-Y.; Davis, J.S. The G-protein-coupled estrogen receptor agonist G-1 suppresses proliferation of ovarian cancer cells by blocking tubulin polymerization. Cell Death Dis. 2013, 4, 1–11. [Google Scholar] [CrossRef] [Green Version]
  19. Rosano, C.; Ponassi, M.; Santolla, M.F.; Pisano, A.; Felli, L.; Vivacqua, A.; Maggiolini, M.; Lappano, R. Macromolecular Modelling and Docking Simulations for the Discovery of Selective GPER Ligands. AAPS J. 2015, 18, 41–46. [Google Scholar] [CrossRef] [Green Version]
  20. Méndez-Luna, D.; Martínez-Archundia, M.; Maroun, R.C.; Ceballos-Reyes, G.; Fragoso-Vázquez, M.; González-Juárez, D.; Correa-Basurto, J. Deciphering the GPER/GPR30-agonist and antagonists interactions using molecular modeling studies, molecular dynamics, and docking simulations. J. Biomol. Struct. Dyn. 2015, 33, 2161–2172. [Google Scholar] [CrossRef]
  21. Méndez-Luna, D.; Bello, M.; Correa-Basurto, J. Understanding the molecular basis of agonist/antagonist mechanism of GPER1/GPR30 through structural and energetic analyses. J. Steroid Biochem. Mol. Biol. 2016, 158, 104–116. [Google Scholar] [CrossRef]
  22. Zacarias-Lara, O.J.; Mendez-Luna, D.; Martinez-Ruiz, G.; Garcia-Sanchez, J.R.; Fragoso-Vazquez, M.J.; Bello, M.; Becerra-Martinez, E.; Garcia-Vazquez, J.B.; Correa-Basurto, J. Synthesis and In Vitro Evaluation of Tetrahydroquinoline Derivatives as Antiproliferative Compounds of Breast Cancer via Targeting the GPER. Anti Cancer Agents Med. Chem. 2019, 19, 760–771. [Google Scholar] [CrossRef]
  23. Martínez-Muñoz, A.; Prestegui-Martel, B.; Mendez-Luna, D.; Fragoso-Vazquez, M.J.; García-Sánchez, J.R.; Bello, M.; Martinez-Archundia, M.; Chávez-Blanco, A.D.; Duenas-Gonzalez, A.; Mendoza-Lujambio, I.; et al. Selection of a GPER1 Ligand via Ligand-based Virtual Screening Coupled to Molecular Dynamics Simulations and Its Anti-proliferative Effects on Breast Cancer Cells. Anti Cancer Agents Med. Chem. 2019, 18, 1629–1638. [Google Scholar] [CrossRef]
  24. Moreno-Ulloa, A.; Mendez-Luna, D.; Beltrán-Partida, E.; Castillo, C.; Guevara, G.; Ramírez-Sánchez, I.; Correa-Basurto, J.; Ceballos, G.; Villarreal, F. The effects of (-)-epicatechin on endothelial cells involve the G protein-coupled estrogen receptor (GPER). Pharmacol. Res. 2015, 100, 309–320. [Google Scholar] [CrossRef] [Green Version]
  25. Qian, H.; Xuan, J.; Liu, Y.; Shi, G. Function of G-Protein-Coupled Estrogen Receptor-1 in Reproductive System Tumors. J. Immunol. Res. 2016, 2016, 1–6. [Google Scholar] [CrossRef] [Green Version]
  26. Jacenik, D.; Cygankiewicz, A.I.; Krajewska, W.M. The G protein-coupled estrogen receptor as a modulator of neoplastic transformation. Mol. Cell. Endocrinol. 2016, 429, 10–18. [Google Scholar] [CrossRef]
  27. Abe, M.; Yamashita, J.; Ogawa, M. Medroxyprogesterone acetate inhibits human pancreatic carcinoma cell growth by inducing apoptosis in association with Bcl-2 phosphorylation. Cancer 2000, 88, 2000–2009. [Google Scholar] [CrossRef]
  28. Guo, J.-M.; Xiao, B.-X.; Dai, D.-J.; Liu, Q.; Ma, H.-H. Effects of daidzein on estrogen-receptor-positive and negative pancreatic cancer cellsin vitro. World J. Gastroenterol. 2004, 10, 860–863. [Google Scholar] [CrossRef]
  29. Lonser, R.R.; Glenn, G.M.; Walther, M.; Chew, E.Y.; Libutti, S.K.; Linehan, W.M.; Oldfield, E.H. von Hippel-Lindau disease. Lancet 2003, 361, 2059–2067. [Google Scholar] [CrossRef]
  30. Cautain, B.; De Pedro, N.; De Escalona, M.M.; Tormo, J.R.; Genilloud, O.; Vicente, F. HCS strategy targeting dysregulation of the VHL/HIF pathway for drug discovery. Adv. Biosci. Biotechnol. 2013, 4, 398–405. [Google Scholar] [CrossRef] [Green Version]
  31. Rigiracciolo, D.C.; Scarpelli, A.; Lappano, R.; Pisano, A.; Santolla, M.F.; De Marco, P.; Cirillo, F.; Cappello, A.R.; Dolce, V.; Belfiore, A.; et al. Copper activates HIF-1α/GPER/VEGF signalling in cancer cells. Oncotarget 2015, 6, 34158–34177. [Google Scholar] [CrossRef] [Green Version]
  32. De Francesco, E.M.; Lappano, R.; Santolla, M.F.; Marsico, S.; Caruso, A.; Maggiolini, M. HIF-1α/GPER signaling mediates the expression of VEGF induced by hypoxia in breast cancer associated fibroblasts (CAFs). Breast Cancer Res. 2013, 15, 1–18. [Google Scholar] [CrossRef] [Green Version]
  33. De Francesco, E.M.; Pellegrino, M.; Santolla, M.F.; Lappano, R.; Ricchio, E.; Abonante, S.; Maggiolini, M. GPER Mediates Activation of HIF1/VEGF Signaling by Estrogens. Cancer Res. 2014, 74, 4053–4064. [Google Scholar] [CrossRef] [Green Version]
  34. ChemBioDraw Ultra 12.0. Available online: http://www.cambridgesoft.com/software/overview.aspx (accessed on 12 November 2020).
  35. Denningtion, R.; Roy, T.; Millam, J. GaussView, Version 5; Semichem Inc.: Shawnee Mission, KS, USA, 2009. [Google Scholar]
  36. Frisch, M.J.; Trucks, G.W.; Schlegel, H.B.; Scuseria, G.E.; Robb, M.A.; Cheeseman, J.R.; Scalmani, G.; Barone, V.; Petersson, G.A.; Nakatsuji, H.; et al. Gaussian 09, Revision A.02; Gaussian, Inc.: Wallingford, CT, USA, 2016. [Google Scholar]
  37. Morris, G.M.; Huey, R.; Lindstrom, W.; Sanner, M.F.; Belew, R.K.; Goodsell, D.S.; Olson, A. AutoDock4 and AutoDockTools4: Automated docking with selective receptor flexibility. J. Comput. Chem. 2009, 30, 2785–2791. [Google Scholar] [CrossRef] [Green Version]
  38. DeLano, W.L. The PyMOL Molecular Graphics System; DeLano Scientific: San Carlos, CA, USA, 2002; Available online: http://www.pymol.org (accessed on 12 November 2020).
  39. Molinspiration Cheminformatics, Bratislava, Slovak Republic. Available online: http://www.molinspiration.com (accessed on 12 November 2020).
  40. Sander, T.; Freyss, J.; Von Korff, M.; Rufener, C. DataWarrior: An Open-Source Program for Chemistry Aware Data Visualization and Analysis. J. Chem. Inf. Model. 2015, 55, 460–473. [Google Scholar] [CrossRef]
  41. Lomize, M.A.; Lomize, A.L.; Pogozheva, I.D.; Mosberg, H.I. OPM: Orientations of Proteins in Membranes database. Bioinformatics 2006, 22, 623–625. [Google Scholar] [CrossRef]
  42. Jo, S.; Kim, T.; Im, W. Automated Builder and Database of Protein/Membrane Complexes for Molecular Dynamics Simulations. PLoS ONE 2007, 2, e880. [Google Scholar] [CrossRef] [Green Version]
  43. Jo, S.; Lim, J.B.; Klauda, J.B.; Im, W. CHARMM-GUI Membrane Builder for Mixed Bilayers and Its Application to Yeast Membranes. Biophys. J. 2009, 97, 50–58. [Google Scholar] [CrossRef] [Green Version]
  44. Woolf, T.B.; Roux, B. Structure, energetics, and dynamics of lipid-protein interactions: A molecular dynamics study of the gramicidin A channel in a DMPC bilayer. Proteins 1996, 24, 92–114. [Google Scholar] [CrossRef]
  45. Case, D.A.; Cheatham, T.E., III; Darden, T.; Gohlke, H.; Luo, R.; Merz, K.M., Jr.; Onufriev, A.; Simmerling, C.; Wang, B.; Woods, R.J. The Amber biomolecular simulation programs. J. Comput. Chem. 2005, 26, 1668–1688. [Google Scholar] [CrossRef] [Green Version]
  46. Skjevik, Å.A.; Madej, B.D.; Walker, R.C.; Eigen, K.T. LIPID11: A Modular Framework for Lipid Simulations Using Amber. J. Phys. Chem. B 2012, 116, 11124–11136. [Google Scholar] [CrossRef] [Green Version]
  47. Wang, J.; Wolf, R.M.; Caldwell, J.W.; Kollman, P.A.; Case, D.A. Development and testing of a general amber force field. J. Comput. Chem. 2004, 25, 1157–1174. [Google Scholar] [CrossRef]
  48. Darden, T.A.; York, D.; Pedersen, L. Particle mesh Ewald: An N-log(N) method for Ewald sums in large systems. J. Chem. Phys. 1993, 98, 10089–10092. [Google Scholar] [CrossRef] [Green Version]
  49. Van Gunsteren, W.; Berendsen, H. Algorithms for macromolecular dynamics and constraint dynamics. Mol. Phys. 1977, 34, 1311–1327. [Google Scholar] [CrossRef]
  50. Lindahl, E.; Hess, B.; Van Der Spoel, D. GROMACS 3.0: A package for molecular simulation and trajectory analysis. J. Mol. Model. 2001, 7, 306–317. [Google Scholar] [CrossRef]
  51. Amadei, A.; Linssen, A.B.M.; Berendsen, H.J.C. Essential dynamics of proteins. Proteins Struct. Funct. Bioinform. 1993, 17, 412–425. [Google Scholar] [CrossRef]
  52. Berendsen, H.; Van Der Spoel, D.; Van Drunen, R. GROMACS: A message-passing parallel molecular dynamics implementation. Comput. Phys. Commun. 1995, 91, 43–56. [Google Scholar] [CrossRef]
  53. Van Der Spoel, D.; Lindahl, E.; Hess, B.; Groenhof, G.; Mark, A.E.; Berendsen, H.J.C. GROMACS: Fast, flexible, and free. J. Comput. Chem. 2005, 26, 1701–1718. [Google Scholar] [CrossRef]
  54. Gohlke, H.; Case, D.A. Converging free energy estimates: MM-PB(GB)SA studies on the protein-protein complex Ras-Raf. J. Comput. Chem. 2004, 25, 238–250. [Google Scholar] [CrossRef]
  55. Kollman, P.A.; Massova, I.; Reyes, C.; Kuhn, B.; Huo, S.; Chong, L.; Lee, M.; Lee, T.; Duan, Y.; Wang, W.; et al. Calculating Structures and Free Energies of Complex Molecules: Combining Molecular Mechanics and Continuum Models. Accounts Chem. Res. 2000, 33, 889–897. [Google Scholar] [CrossRef]
  56. Hou, T.; Wang, J.; Li, Y.; Wang, W. Assessing the Performance of the MM/PBSA and MM/GBSA Methods. 1. The Accuracy of Binding Free Energy Calculations Based on Molecular Dynamics Simulations. J. Chem. Inf. Model. 2011, 51, 69–82. [Google Scholar] [CrossRef]
  57. Hou, T.; Wang, J.; Li, Y.; Wang, W. Assessing the performance of the molecular mechanics/Poisson Boltzmann surface area and molecular mechanics/generalized Born surface area methods. II. The accuracy of ranking poses generated from docking. J. Comput. Chem. 2011, 32, 866–877. [Google Scholar] [CrossRef] [Green Version]
  58. Xu, L.; Sun, H.; Li, Y.; Wang, J.; Hou, T. Assessing the Performance of MM/PBSA and MM/GBSA Methods. 3. The Impact of Force Fields and Ligand Charge Models. J. Phys. Chem. B 2013, 117, 8408–8421. [Google Scholar] [CrossRef]
  59. Sun, H.; Li, Y.; Tian, S.; Xu, L.; Hou, T. Assessing the performance of MM/PBSA and MM/GBSA methods. 4. Accuracies of MM/PBSA and MM/GBSA methodologies evaluated by various simulation protocols using PDBbind data set. Phys. Chem. Chem. Phys. 2014, 16, 16719–16729. [Google Scholar] [CrossRef]
  60. Hou, T.; Li, N.; Li, Y.; Wang, W. Characterization of Domain–Peptide Interaction Interface: Prediction of SH3 Domain-Mediated Protein–Protein Interaction Network in Yeast by Generic Structure-Based Models. J. Proteome Res. 2012, 11, 2982–2995. [Google Scholar] [CrossRef]
  61. Sun, H.; Li, Y.; Shen, M.; Tian, S.; Xu, L.; Pan, P.; Guan, Y.; Hou, T. Assessing the performance of MM/PBSA and MM/GBSA methods. 5. Improved docking performance using high solute dielectric constant MM/GBSA and MM/PBSA rescoring. Phys. Chem. Chem. Phys. 2014, 16, 22035–22045. [Google Scholar] [CrossRef]
  62. Liu, L.; Zhang, Y.; Wang, Y. Phosphine-free Palladium acetate catalyzed Suzuki reaction in water. J. Org. Chem. 2005, 70, 6122–6125. [Google Scholar] [CrossRef]
  63. Miyaura, N.; Yamada, K.; Suzuki, A. A new stereospecific cross-coupling by the palladium-catalyzed reaction of 1-alkenylboranes with 1-alkenyl or 1-alkynyl halides. Tetrahedron Lett. 1979, 20, 3437–3440. [Google Scholar] [CrossRef] [Green Version]
  64. Verhulst, C.; Coiffard, C.; Coiffard, L.J.; Rivalland, P.; De Roeck-Holtzhauer, Y. In vitro correlation between two colorimetric assays and the pyruvic acid consumption by fibroblasts cultured to determine the sodium laurylsulfate cytotoxicity. J. Pharmacol. Toxicol. Methods 1998, 39, 143–146. [Google Scholar] [CrossRef]
  65. Liu, Y.; Nair, M.G. An Efficient and Economical MTT Assay for Determining the Antioxidant Activity of Plant Natural Product Extracts and Pure Compounds. J. Nat. Prod. 2010, 73, 1193–1195. [Google Scholar] [CrossRef]
Pharmaceuticals 14 00049 g001 550
Figure 1. Equilibrium properties of GPER–ligand systems through molecular dynamics (MD) simulations. (A) Area per lipid values; (B) root mean square deviation (RMSD); (C) radius of gyration (Rg); (D) root mean square fluctuation (RMSF) for GPER-4 (red), GPER-5 (blue) and GPER-7 (black).
Figure 1. Equilibrium properties of GPER–ligand systems through molecular dynamics (MD) simulations. (A) Area per lipid values; (B) root mean square deviation (RMSD); (C) radius of gyration (Rg); (D) root mean square fluctuation (RMSF) for GPER-4 (red), GPER-5 (blue) and GPER-7 (black).
Pharmaceuticals 14 00049 g001
Pharmaceuticals 14 00049 g002 550
Figure 2. Projection of the GPER–ligand complexes in phase space. (A) The first 15 eigenvectors of the covariance matrix for GPER-4 (red), GPER-5 (blue) and GPER-7 (black). (B) The percentage of each eigenvector from panel A. (C) Projection of the motion in the phase space along the first and second eigenvectors (principal component (PC) 1 vs. PC2) for GPER-4 (red), GPER-5 (blue) and GPER-7 (black).
Figure 2. Projection of the GPER–ligand complexes in phase space. (A) The first 15 eigenvectors of the covariance matrix for GPER-4 (red), GPER-5 (blue) and GPER-7 (black). (B) The percentage of each eigenvector from panel A. (C) Projection of the motion in the phase space along the first and second eigenvectors (principal component (PC) 1 vs. PC2) for GPER-4 (red), GPER-5 (blue) and GPER-7 (black).
Pharmaceuticals 14 00049 g002
Pharmaceuticals 14 00049 g003 550
Figure 3. (A) Overlapping representation of the close contact of all compounds over GPER. Map of interactions for the (B) GPER-4, (C) GPER-5 and (D) GPER-7 systems. The first and most populated complex derived from a clustering analysis was used to obtain the map of interactions. Hydrogen bonds are depicted as green dashes.
Figure 3. (A) Overlapping representation of the close contact of all compounds over GPER. Map of interactions for the (B) GPER-4, (C) GPER-5 and (D) GPER-7 systems. The first and most populated complex derived from a clustering analysis was used to obtain the map of interactions. Hydrogen bonds are depicted as green dashes.
Pharmaceuticals 14 00049 g003
Pharmaceuticals 14 00049 g004 550
Figure 4. Growth inhibitory activities of 4, 5 and 7 on the MIA Paca-2 cancer cell line. The results are expressed as the mean of three individual experiments; p < 0.05.
Figure 4. Growth inhibitory activities of 4, 5 and 7 on the MIA Paca-2 cancer cell line. The results are expressed as the mean of three individual experiments; p < 0.05.
Pharmaceuticals 14 00049 g004
Pharmaceuticals 14 00049 g005 550
Figure 5. Growth inhibitory activities of 4, 5 and 7 on the RCC4-VA cancer cell line. The results are expressed as the mean of three individual experiments; p < 0.05.
Figure 5. Growth inhibitory activities of 4, 5 and 7 on the RCC4-VA cancer cell line. The results are expressed as the mean of three individual experiments; p < 0.05.
Pharmaceuticals 14 00049 g005
Pharmaceuticals 14 00049 g006 550
Figure 6. Growth inhibitory activities of 4, 5 and 7 on the Hep G2 cancer cell line. The results are expressed as the mean of three individual experiments; p < 0.05.
Figure 6. Growth inhibitory activities of 4, 5 and 7 on the Hep G2 cancer cell line. The results are expressed as the mean of three individual experiments; p < 0.05.
Pharmaceuticals 14 00049 g006
Pharmaceuticals 14 00049 sch001 550
Scheme 1. Synthesis of 4 using PEG2000 for cross-coupling reactions under N2 atmosphere.
Scheme 1. Synthesis of 4 using PEG2000 for cross-coupling reactions under N2 atmosphere.
Pharmaceuticals 14 00049 sch001
Pharmaceuticals 14 00049 sch002 550
Scheme 2. Synthesis of 5 by Suzuki–Miyaura cross-coupling reactions.
Scheme 2. Synthesis of 5 by Suzuki–Miyaura cross-coupling reactions.
Pharmaceuticals 14 00049 sch002
Pharmaceuticals 14 00049 sch003 550
Scheme 3. Synthesis of 7 by protection with a Boc protecting group.
Scheme 3. Synthesis of 7 by protection with a Boc protecting group.
Pharmaceuticals 14 00049 sch003
Table
Table 1. The binding energies (−ΔGbinding in Kcal/mol) and interacting residues of G protein-coupled estrogen receptor (GPER) towards the synthetized ligands yielded by docking analysis.
Table 1. The binding energies (−ΔGbinding in Kcal/mol) and interacting residues of G protein-coupled estrogen receptor (GPER) towards the synthetized ligands yielded by docking analysis.
LigandBinding EnergyNonbonded Interactions
4−8.46N44, H52, Q54 *, L59, L119, R122 *, F206
5−8.36Q53, Q54 *, G58, L59, H120*, F278 *, I279, H302 *, P303, G306, H307, N310
7−8.21E51, Q53 *, Q54, G58, L59, H120 *, F278 *, I279, H282 *, H302 *, P303, N310
* Residue forming H bond.
Table
Table 2. Physicochemical and toxicological–biological properties of the ligands analyzed through online servers. MW = molecular weight, HBD = hydrogen bond donor, HBA = hydrogen bond acceptor, Mut = mutagenesis, Ter = teratogenicity, Irri = irritability, Rep = reproductive effects.
Table 2. Physicochemical and toxicological–biological properties of the ligands analyzed through online servers. MW = molecular weight, HBD = hydrogen bond donor, HBA = hydrogen bond acceptor, Mut = mutagenesis, Ter = teratogenicity, Irri = irritability, Rep = reproductive effects.
LigandMWLog PHBDHBAMutTerIrriRep
4456.455.6228NoneNoneNoneNone
5441.485.7126NoneNoneNoneNone
7514.375.9717NoneNoneNoneNone
Table
Table 3. Binding free energy components of bound GPER systems (in units of Kcal/mol): the polar (ΔEpolar = ΔEele + ΔGele,sol) and nonpolar (ΔEnonpolar = ΔEvwd + ΔGnpol,sol) contributions. All the energies are averaged over 400 snapshots at time intervals of 100 ps from the last 60 ns production runs and are in Kcal/mol (±standard deviation). Results were obtained by MD computations.
Table 3. Binding free energy components of bound GPER systems (in units of Kcal/mol): the polar (ΔEpolar = ΔEele + ΔGele,sol) and nonpolar (ΔEnonpolar = ΔEvwd + ΔGnpol,sol) contributions. All the energies are averaged over 400 snapshots at time intervals of 100 ps from the last 60 ns production runs and are in Kcal/mol (±standard deviation). Results were obtained by MD computations.
SystemΔEvdwΔEeleΔGele,solΔGnpol,solΔEpolarΔEnonpolarΔGmm/gbsaTDSΔGbind
GPER-4−51.66 (2.9)15.72 (4.2)3.50 (0.10)−6.44 (0.22)19.22−58.10−38.88 (4.0)−26.43 (1.4)−12.43
GPER-5−47.3 (2.78)19.15 (3.0)4.70 (0.45)−6.10 (0.20)23.85−53.40−29.55 (3.2)−22.34 (2.6)−7.21
GPER-7−41.73 (2.5)−3.90 (0.3)22.94 (0.3)−5.52 (0.10)19.04−47.25−28.21 (2.5)−24.14 (1.3)−4.07
Table
Table 4. Half maximal inhibitory concentration (IC50) values obtained for the ligands tested in the cancer cell lines. Values >50 are displayed for information only; values are in µM range.
Table 4. Half maximal inhibitory concentration (IC50) values obtained for the ligands tested in the cancer cell lines. Values >50 are displayed for information only; values are in µM range.
IC50
LigandRCC4-VARCC4-VHLMIA Paca-2Hep G2
724.46>5014.78>50
4>50>5018.18>50
539.29>5021.1531.58
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Méndez-Luna, D.; Morelos-Garnica, L.A.; García-Vázquez, J.B.; Bello, M.; Padilla-Martínez, I.I.; Fragoso-Vázquez, M.J.; Dueñas González, A.; De Pedro, N.; Gómez-Vidal, J.A.; Mendoza-Figueroa, H.L.; et al. Modifications on the Tetrahydroquinoline Scaffold Targeting a Phenylalanine Cluster on GPER as Antiproliferative Compounds against Renal, Liver and Pancreatic Cancer Cells. Pharmaceuticals 2021, 14, 49. https://doi.org/10.3390/ph14010049

AMA Style

Méndez-Luna D, Morelos-Garnica LA, García-Vázquez JB, Bello M, Padilla-Martínez II, Fragoso-Vázquez MJ, Dueñas González A, De Pedro N, Gómez-Vidal JA, Mendoza-Figueroa HL, et al. Modifications on the Tetrahydroquinoline Scaffold Targeting a Phenylalanine Cluster on GPER as Antiproliferative Compounds against Renal, Liver and Pancreatic Cancer Cells. Pharmaceuticals. 2021; 14(1):49. https://doi.org/10.3390/ph14010049

Chicago/Turabian Style

Méndez-Luna, David, Loreley Araceli Morelos-Garnica, Juan Benjamín García-Vázquez, Martiniano Bello, Itzia Irene Padilla-Martínez, Manuel Jonathan Fragoso-Vázquez, Alfonso Dueñas González, Nuria De Pedro, José Antonio Gómez-Vidal, Humberto Lubriel Mendoza-Figueroa, and et al. 2021. "Modifications on the Tetrahydroquinoline Scaffold Targeting a Phenylalanine Cluster on GPER as Antiproliferative Compounds against Renal, Liver and Pancreatic Cancer Cells" Pharmaceuticals 14, no. 1: 49. https://doi.org/10.3390/ph14010049

APA Style

Méndez-Luna, D., Morelos-Garnica, L. A., García-Vázquez, J. B., Bello, M., Padilla-Martínez, I. I., Fragoso-Vázquez, M. J., Dueñas González, A., De Pedro, N., Gómez-Vidal, J. A., Mendoza-Figueroa, H. L., & Correa-Basurto, J. (2021). Modifications on the Tetrahydroquinoline Scaffold Targeting a Phenylalanine Cluster on GPER as Antiproliferative Compounds against Renal, Liver and Pancreatic Cancer Cells. Pharmaceuticals, 14(1), 49. https://doi.org/10.3390/ph14010049

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop