Next Article in Journal
Transcriptional Control of Seed Life: New Insights into the Role of the NAC Family
Next Article in Special Issue
Structure-Based De Novo Design for the Discovery of Miniprotein Inhibitors Targeting Oncogenic Mutant BRAF
Previous Article in Journal
Apolipoprotein-CIII O-Glycosylation Is Associated with Micro- and Macrovascular Complications of Type 2 Diabetes
Previous Article in Special Issue
Gaining Insights into Key Structural Hotspots within the Allosteric Binding Pockets of Protein Kinases
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 25
Issue 10
10.3390/ijms25105368
ijms-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

Structural Models for a Series of Allosteric Inhibitors of IGF1R Kinase

by
Jyoti Verma
1 and
Harish Vashisth
1,2,3,4,*
1
Department of Chemical Engineering and Bioengineering, University of New Hampshire, Durham, NH 03824, USA
2
Department of Chemistry, University of New Hampshire, Durham, NH 03824, USA
3
Integrated Applied Mathematics Program, University of New Hampshire, Durham, NH 03824, USA
4
Molecular and Cellular Biotechnology Program, University of New Hampshire, Durham, NH 03824, USA
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2024, 25(10), 5368; https://doi.org/10.3390/ijms25105368
Submission received: 7 April 2024 / Revised: 9 May 2024 / Accepted: 11 May 2024 / Published: 14 May 2024
(This article belongs to the Special Issue Protein Kinase in Disease)
Browse Figures
Versions Notes

Abstract

:
The allosteric inhibition of insulin-like growth factor receptor 1 kinase (IGF1RK) is a potential strategy to overcome selectivity barriers for targeting receptor tyrosine kinases. We constructed structural models of a series of 12 indole-butyl-amine derivatives that have been reported as allosteric inhibitors of IGF1RK. We further studied the dynamics and interactions of each inhibitor in the allosteric pocket via all-atom explicit-solvent molecular dynamics (MD) simulations. We discovered that a bulky carbonyl substitution at the R1 indole ring is structurally unfavorable for inhibitor binding in the IGF1RK allosteric pocket. Moreover, we found that the most potent derivative (termed C11) acquires a distinct conformation: forming an allosteric pocket channel with better shape complementarity and interactions with the receptor. In addition to a hydrogen-bonding interaction with V1063, the cyano derivative C11 forms a stable hydrogen bond with M1156, which is responsible for its unique binding conformation in the allosteric pocket. Our findings show that the positioning of chemical substituents with different pharmacophore features at the R1 indole ring influences molecular interactions and binding conformations of indole-butyl-amine derivatives and, hence, dramatically affects their potencies. Our results provide a structural framework for the design of allosteric inhibitors with improved affinities and specificities against IGF1RK.

1. Introduction

The type 1 insulin-like growth factor receptor (IGF1R) is a receptor tyrosine kinase (RTK) widely expressed in most somatic cells and tissues of all vertebrates [1]. IGF1R plays a critical role in cell growth, development, and differentiation [2]. IGF1R is a homodimer consisting of two α and two β chains connected by disulfide bonds [3]. It is synthesized as a precursor peptide of 180 kDa that undergoes posttranslational modifications to yield the mature α 2 β 2 receptor [4,5]. The extracellular domain is comprised of the α chain and part of the β chain. The remaining β chain makes up the single-pass transmembrane domain and a cytoplasmic tyrosine kinase domain. Further, six different protein subdomains are identified within the extracellular region: three type III fibronectin domains designated FnIII-1, FnIII-2, and FnIII-3; an N-terminal leucine-rich domain (L1); a cysteine-rich domain (CR); and a second leucine-rich domain (L2) [5]. Unlike many other RTKs, IGF1R is activated by the extracellular binding of secreted growth factor peptides IGF-1 and IGF-2, resulting in autophosphorylation of the intracellular kinase domain and the formation of docking sites for signaling molecules [6].
As a growth mediator, IGF1R is ubiquitously expressed and promotes normal tissue development [7]. However, it has also been implicated in mitogenesis and suppression of apoptosis in tumor cells [8,9,10]. Various cancer types are associated with elevated levels of IGF1R and circulating IGF-1/IGF-2 [11,12]. Different mechanisms can influence IGF1R activity to produce cancer phenotypes, including (i) IGF1R overexpression caused by IGF1R gene aberrations, (ii) regulation of IGF1R expression by modifier genes, (iii) increased levels of ligands, and (iv) transactivation by other membrane receptors [13]. Hence, various studies have demonstrated that IGF1R inhibition and downregulation cause substantial in vivo tumor cell death, impede tumorigenesis, and trigger a host response that eliminates surviving tumor cells [14,15,16,17]. Experimental investigations also suggest that inhibition of IGF1R kinase activity will initiate the disruption of tumor signaling, leading to inhibition of tumor growth [17]. Therefore, IGF1R is an attractive target for chemotherapeutic interventions. Moreover, targeting the RTK family of proteins is efficacious and well-tolerated, as demonstrated by the successful development of anti-tumor drugs against VEGFR, HER-2, and EGFR [18].
The therapeutics against IGF1R include antibodies to the extracellular domain that block receptor signaling and small-molecule inhibitors of the tyrosine kinase domain. A recent analysis reported that about 16 IGF1R drugs (both monoclonal antibodies and small-molecule inhibitors) entered oncology clinical trials between 2003 and 2021 [19]. However, none were approved for cancer treatment [19]. To date, several monoclonal antibodies against IGF1R, such as ganitumab, figitumumab, dalotuzumab, dusigitumab, and xentuzumab, have reached pre-clinical and clinical trials [20,21,22,23,24]. Small-molecule IGF1RK inhibitors have also shown considerable efficacy in pre-clinical trials [25]. The ATP-competitive inhibitors BMS-754807 and OSI906 (Linsitinib) are well-studied, although numerous trials indicated unsatisfactory outcomes [26,27]. AXL-1717 (picropodophyllotoxin), a non-ATP-competitive IGF1R kinase inhibitor, has shown some potential in treating patients with relapsed malignant astrocytomas; nevertheless, its precise mode of action is yet unknown [28]. Many small-molecule chemotypes still lack structural information and the mechanism of inhibition.
Despite persistent efforts to develop therapeutically useful IGF1R blockers with diverse molecular natures and mechanisms of action, the results of clinical trials have not been encouraging thus far. The most often cited explanations for the failure of IGF1R-targeted treatments are crosstalk with other receptors, IR-IGF1R hybrid receptors, a high degree of homology among the IGF1R and IR catalytic domains, and IR signaling compensation for IGF1R inhibition [25,29,30,31,32]. Hence, the specificity and selectivity of inhibitors continue to be further explored. Furthermore, allosteric inhibition has been considered as a possible strategy for increasing specificity and selectivity for kinase inhibitors [33,34]. However, experimental evidence of allosteric interaction has been demonstrated only for a small number of allosteric kinase inhibitors [35]. Henrich et al. described the design of a series of indole-butyl-amine derivatives as allosteric IGF1RK inhibitors and measured their biochemical activity against IGF1RK and IRK [36]. The allosteric binding site in the IGF1RK was characterized by X-ray crystallographic studies with one of the inhibitors (C10) from the series. However, the binding modes of other inhibitors in the series and their interactions with the IGF1RK remain unresolved. The inhibitor activity is invariably the consequence of various forces, including hydrophobic interactions, polar interactions, and shape and charge complementarity of the binding pocket, and hence, it is challenging to interpret the differences in potencies of closely related compounds without a thorough structural analysis. Here, we report the structural model of each indole-butyl-amine derivative in complex with IGF1RK and provide detailed structural analyses to provide a rationale for the observed differences in the activities of these allosteric inhibitors. We perform a structural and energetic comparison of the interactions of each allosteric inhibitor using molecular docking, all-atom MD simulations, and free-energy calculations. In addition, we report the unique binding conformation of the most potent derivative, C11, and unravel the chemical basis for the design of potent allosteric inhibitors of IGF1RK.

2. Results and Discussion

2.1. Pharmacophore Modeling and Molecular Docking

In the current work, we elucidate the binding conformations of a series of 12 allosteric IGF1RK inhibitors (Table 1) while determining the rationale for differences in their potencies [36]. The series of indole-butyl-amine derivatives (C1–C12) have a large common substructure (1-[4(5-cyano-1H-indol-3-yl)-butyl]-piperidin-4-yl)-amide) with modifications in the indole ring located at one of the termini (Figure 1 and Table 1). We categorized these allosteric inhibitors based on the chemical moieties substituted at the R group (Figure 1). Group A constitutes the indole amine derivatives, C1, C2, and C3, which differ in the amide bond attachment on the indole ring. The compounds C4 and C5 are alkyl derivatives in Group B with methyl and ethyl substituents added to the indolyl nitrogen. Group C includes the compounds C6, C7, C8, and C9, which have more bulkier aldehyde and ketone substituents. The cyano derivatives, C10, C11, and C12, are in Group D and have a cyano group attached to the indole ring.
Next, we identified the pharmacophore features of all compounds in the series: C1 to C12. This analysis revealed the following features: an H-bond donor (D), an H-bond acceptor (A), a hydrophobic site (H), an aromatic ring (R), and a positive ionic (P) group (Figure 2). Since the compounds have a large common substructure, the pharmacophore features defining the parent structure were identical among all inhibitors. The unique features of all compounds are highlighted in Figure 2. The presence of two aromatic rings was also consistent for all the R group substituents because of the common indole ring. The compounds from Group A (C1–C3), have an H-bond donor group along with two aromatic rings. The indole ring of compounds C1, C2, and C3 has an amide group attached to it, which acts as an electron-donating substituent. On the other hand, the Group B compounds have a hydrophobic site with the presence of methyl and ethyl substituents on 4-indoyl nitrogen for C4 and C5, respectively. However, the compounds belonging to Group C have distinct pharmacophore features. For example, C6 and C7 have an H-bond donor as well as an acceptor group attached to their indole ring. Compound C8 has a bulkier substituent, trifluoromethyl ketone, where the trifluoromethyl group is the hydrophobic site, the carbonyl oxygen is the H-bond acceptor, and the amide at the indole ring acts as the H-bond donor. Furthermore, C9 also has the carbonyl oxygen atom as the acceptor and the 7-indolyl as the donor as well as two additional H-bond donor groups from the carboxylic amide substituent. The cyano derivatives (Group D) C10 and C12 are isomers that have both H-bond donor and acceptor groups. Besides a donor and an acceptor group, C11 has an additional hydrophobic site with a fluorine atom attached to the indole ring at the fifth position. Following the chemical characterization, we further compared these allosteric inhibitors structurally and energetically in the IGF1RK allosteric pocket.
Using the molecular docking approach, we initially modeled and predicted the conformation of each allosteric inhibitor in the IGF1RK allosteric site. We used the X-ray crystal structure of the IGF1RK domain in complex with the allosteric inhibitor C10 (PDB ID 3LW0) as the basis for the initial modeling of the remaining 11 inhibitors. For each inhibitor, the conformation with the lowest docking score was selected. We further compared the docked conformations of all 11 derivatives with the crystallized conformation of C10 in the IGF1RK allosteric pocket (Figure 3). Overall, the conformations of the common substructure of these derivatives attained a similar configuration. All compounds showed an identical configuration of the 5-cyano indole ring (R2), while the conformational differences were observed along the R substituent at the other end (R1) (Figure 3B). The docking results reflect that hydrogen bonding with the carbonyl oxygen atom of V1063 could be a significant factor for inhibitor binding to the IGF1RK allosteric pocket. Hence, an H-bond donor group at the 5-cyano indole ring (R2) is a key feature of an allosteric inhibitor candidate for IGF1RK. The electron-withdrawing cyano substituent at the R2 indole ring potentially also has a critical role in interaction with residues in the α C-helix. The presence of these chemical features at the R2 indole ring (5-cyano indole) resulted in its identical conformation in all inhibitors at the IGF1RK pocket. To corroborate these observations and further investigate the significance of chemical modifications at the R1 indole with regard to inhibitor potency, we studied the conformational dynamics of each complex using all-atom MD simulations.

2.2. Conformational Dynamics of Complexes of IGF1RK/Inhibitors

We performed MD simulations for each complex to understand the conformational changes in the IGF1RK structure upon inhibitor binding in the allosteric pocket. In parallel, the dynamics of each allosteric inhibitor were studied to comprehend differences in their binding affinities. We first computed the RMSD values of the IGF1RK backbone atoms and the non-hydrogen atoms of each inhibitor. For Group A inhibitors, C1, C2, and C3, the IGF1RK backbone RMSD distribution curves have single peaks between 2.0 Å and 3.5 Å (Figure S1A). The protein backbone RMSD for the C2 complex has a narrow distribution curve that shifts further to the left, indicating a lower RMSD value. The Group B inhibitor C4 has an average protein backbone RMSD value of 2.5 Å, whereas C5 has a broader distribution curve ranging between 2 Å and 4.5 Å, reflecting higher deviations in the protein backbone atoms (Figure S1A). The Group C inhibitors C6, C7, C8, and C9 have relatively higher IGF1RK backbone RMSD values with wider distribution curves along the RMSD axis (Figure S1A). Further, in Group D, the cyano derivative C10 has its IGF1RK backbone RMSD distributed between 1.5 Å and 4.5 Å. In contrast, its isomer C12 has the lowest protein backbone RMSD of about 1.9 Å. Also, the inhibitor C11, on the other hand, has an IGF1RK backbone RMSD value of 2.2 Å. Overall, the IGF1RK backbone was more flexible with Group C inhibitors and was relatively less flexible with Group A and Group D inhibitors.
We further report the mean RMSD values derived for the non-hydrogen atoms of each inhibitor (Figure 4). Among Group A inhibitors, C2 has the lowest mean RMSD value of 0.8 Å, followed by C3 and C1 with 1.8 Å and 2.3 Å, respectively (Figure 4). Also, the RMSD distribution curve for C2 spans between 0.5 Å and 1.5 Å (Figure S1B). However, the inhibitors C1 and C3 have their RMSD values distributed between 1 Å and 3.5 Å and from 1 Å to 3 Å, respectively (Figure 4 and Figure S1B). These observations suggest that the inhibitor C2 (with a difference in the amide position at the indole ring R1 compared to C1) has lower deviations, leading to a relatively more stable conformation (Figure 4 and Figure S2). The mean RMSD for the C4 inhibitor is 0.86 Å, whereas for C5, it is 1.67 Å (Figure 4). Although C4 shows a single peak and a narrow distribution, C5 shows a distribution with two peaks spanning RMSD values between 1 Å and 2.5 Å (Figure S1B). Hence, in Group B, the inhibitor C4 (methyl substituted) exhibits overall low deviations in the complex. For Group C inhibitors, C6, C7 and C8, the mean RMSD values are 1.6 Å, 0.8 Å, and 1.6 Å, respectively. The restricted movement/flexibility of the substituted ketone and aldehyde groups of these inhibitors in the allosteric pocket may account for their lower RMSD values (Figure 4). However, the inhibitor C9 has an RMSD value of 2.3 Å, reflecting relatively more flexibility of the inhibitor in the allosteric pocket. Among the cyano derivatives (Group D), C10 (previously co-crystalized with IGF1RK) has a mean RMSD value of 1.7 Å. In contrast, its isomer C12 has the highest inhibitor RMSD value of 2.5 Å. The inhibitor C12 is highly flexible in the allosteric pocket and deviates from the allosteric region during simulations (Figure 4 and Figure S2). The inhibitor C11, on the other hand, has an RMSD value of 1.5 Å, which is the lowest among Group D compounds. These findings altogether imply that the position of the amide, as well as the cyano group or any electron-withdrawing group on the indole ring, has a critical role in inhibitor stability.
We also evaluated the RMSF values for the IGF1RK residues to characterize their dynamics. The overall trend of residue fluctuations in all complexes was relatively similar. There was an evident fluctuation peak corresponding to the loop region (residues 1093 to 1115) with RMSF values ranging between 4 Å and 6 Å in all complexes (Figure S1C). The fluctuation amplitudes corresponding to the RMSF values of the catalytic loop (residues 1133 to 1142) and the activation loop (residues 1152–1174) residues are comparable for all inhibitor complexes. Residues from these loops also specifically interact with the substituted R group of the inhibitors. The time evolution trace of the hydrogen bonds formed between the allosteric inhibitors and IGF1RK is shown in Figure S1D. The maximum number of hydrogen bonds that may form is 4–5; however, throughout the trajectory, all inhibitors consistently only display a small number (1–2) of hydrogen bonds. We further elucidated the stability of these hydrogen-bonding interactions in each complex. We also measured the radius of gyration ( R g ) for each complex to indicate the compactness of the protein structure. We observed that the mean Rg value for all complexes was in the range between 19.7 Å and 20.3 Å (Figure S3). The consistent value of Rg for all complexes indicates that the IGF1RK structure remained compact upon binding to inhibitors. Although a visual inspection of simulation trajectories provided insight into the binding conformations of the allosteric inhibitors in the IGF1RK pocket (Figure S2), we further investigated the predominant binding conformation for each inhibitor.

2.3. Dominant Conformation of Each Inhibitor in the Allosteric Pocket

We performed principal component analysis of all MD trajectories to understand the key conformational motions and extract the structural information for each IGF1RK-allosteric inhibitor complex. The first few principal components define the collective motions of the localized atomic fluctuations. The two-dimensional FES plots of each complex for the first two principal components (PC1 and PC2) are represented in Figure 5. The conformational phase space and the energy basins for each complex were unique and different. From the free-energy maps, we identified a global minimum conformation for each inhibitor that is associated with the lowest free energy. The allosteric inhibitors C1 and C2 from Group A and C4 and C5 from Group B have two well-separated minima with higher energy barriers (Figure 5). Hence, conformations corresponding to both minima were analyzed for these inhibitors. Despite a widespread conformational phase space, a single global free-energy minimum was observed for Group C and Group D inhibitors. The inhibitors C6 and C7 have well-separated global minima, while C8 and C9 have multiple local minima and a single global minimum with lower energy barriers. The cyano derivatives (Group D), on the other hand, have a single global free-energy minimum (Figure 5). We further obtained and evaluated the conformations of all allosteric inhibitors with IGF1RK corresponding to these minima.
We assessed the predominant conformation of each of the allosteric inhibitors in the IGF1RK allosteric pocket. The structure superimposition of all 12 complexes is shown in Figure 6. An overall assessment of the binding modes of all allosteric inhibitors showed that the Group A (C1–C3), B (C4–C5), and C (C6–C9) derivatives occupy similar configurations clustered at the same region of the allosteric pocket. As shown in Figure 6A, the shared indole R2 moiety in all inhibitors has similar conformations. Furthermore, the local conformational changes of the A-loop and the α C-helix caused by inhibitor binding affect the IGF1RK allosteric pocket characteristics. As a consequence, Group D molecules (cyano derivatives, C11 and C12) protrude away from the proximal allosteric binding pocket region and attain different conformations with respect to their initial docked conformations (Figure 4 and Figure 6). A detailed assessment of the conformations of these derivatives is discussed in the following section.

2.4. Effect of Substitutions at the R1 Indole Ring and Unique Binding Mode of C11

The effects of chemical substitutions in indole-butyl-amine derivatives on their biochemical activities (IC50) against IGF1RK was investigated by Henrich et al. [36]. We further recapitulate the substitutions at the R1 indole ring with molecular and structural aspects. The comparison of the whole series of derivatives emphasizes the molecular basis for the activity differences among them. We observed a favorable H-bond interaction between the acceptor carbonyl oxygen atom of V1063 and the donor group NH9 of the R2 indole ring in all inhibitor complexes (Figures S4–S7). This interaction, however, was disrupted in the complexes of C1, C6, C8, and C9, but it was stable and consistent for all other inhibitors (Figures S4 and S6). Furthermore, the atoms of the R2 indole ring participate in interactions with the hydrophobic residues V1063, M1054, and M1079, and these interactions are conserved in all complexes (Figure S9). This suggests that the 5-cyano indole (R2) at one end of the inhibitor establishes a stable hydrogen bond and is buried in a hydrophobic cavity created by M1054 and M1079, making it suitable for the IGF1RK allosteric pocket. Hence, we show that an electron-donating group with two aromatic rings is responsible for crucial interactions with the 5-cyano indole (R2) in the binding pocket. The substitution at the R1 indole ring, however, accounts for the differences in their affinities for IGF1RK. Residues H1133, R1134, I1151, G1152, D1153, and L1174 are common residues involved in interactions with the R1 indole ring of each derivative. In particular, the sidechain of the residue R1134 interacts with the R1 indole ring, influencing its conformation. This cation– π interaction with R1134 is mostly conserved in all inhibitors. For C1, the position of the amide at the indole ring disfavors its initial interactions at the binding pocket, resulting in dissociation of the H-bond with V1063 while forming new H-bonds with K1171 and G1172 (Figure S4). Furthermore, the dissociation of the H-bond in C6, C8, and C9 could be due to the inability of bulky carbonyl groups like ketone, aldehyde, and amide to fit into the allosteric pocket (Figure S6). Further, the cyano derivative C12 interacts with the α C-helix by extending out of the allosteric region. This observation emphasizes that the position of the cyano group in C12 might cause steric hindrance and, hence, is difficult to accommodate in the allosteric pocket (Figure S7). Altogether, the position of each substitution at the R1 moiety is critical, as it is responsible for both electronic properties and steric effects in the structure. An R1 indole ring with an electron-donating group (amide group) at position 7 and an electron-acceptor group at position 3 likely favors inhibitor binding as well as its activity against IGF1RK.
We performed a comprehensive analysis of the binding free energy using the MMGBSA approach to describe the molecular interactions and energetics of the binding of allosteric inhibitors. The total binding free energy ( Δ G) of each allosteric inhibitor with IGF1RK is shown in Figure S8. The time evolution traces for the binding energy estimated at each frame indicate a relatively stable curve for the cyano derivatives, C10 and C11 (Figure S8). The Δ Gbind for individual residues involved in the interaction with each derivative are also given in Figure S9. The C10 and C11 compounds are the most potent derivatives reported against IGF1RK, with IC50 values of 0.4 μ M and 0.2 μ M, respectively (Table 1). The derivative C10 has been co-crystallized with the IGF1RK domain with the most favorable conformation of the 3-cyano indolyl carboxamide substituent (R1) in the X-ray structure. Additional simulations of C10 and C11 (Figure S10) showed that the RMSD of the protein backbone atoms and the RMSF of the C α atoms were comparable in independent MD simulations. Our structure was determined by analyzing ensemble conformations of C10 and has its 3-cyano indolyl carboxamide substituent (R1) stacked parallel to the C-loop, forming favorable cation– π interactions with the residue R1134 ( Δ Gbind = −2.32 kcal/mol) (Figure 7). Further, the IGF1RK residues M1054, V1063, and M1079 significantly contribute to the binding of C10 with Δ Gbind values of −2.45 kcal/mol, −1.67 kcal/mol, and −1.56 kcal/mol, respectively (Figure 7).
Our findings show that the most potent C11 derivative achieves a unique conformation in the IGF1RK pocket compared to all other derivatives. Figure 6B shows that C11 induces a conformational change in the A-loop that transforms the shape of the allosteric pocket in IGF1RK and offers a cleft-like formation at the allosteric binding region. A detailed inspection of these two compounds bound to IGF1RK shows that the binding of C11 is influenced by an additional hydrogen bond with the residue M1156, providing better shape complementarity and interactions with the receptor (Figure 6 and Figure 7). The indole ring R1 of C11 is stacked between the sidechains of residues K1171 and T1157 and makes favorable interactions with residues G1155 and G1172 (Figure 7). Hence, the R1 substituent of C11 is stabilized by the favorable energetic contributions of these residues, with G1155 ( Δ Gbind = −1.68 kcal/mol), M1156 ( Δ Gbind = −1.84 kcal/mol), and T1157 ( Δ Gbind = −2.00 kcal/mol) contributing most to the binding free energy of C11 in the IGF1RK pocket (Figure 7). The presence of a fluorine atom (hydrophobic site; H) at the indole ring provides improved activity by supporting and directing a stable conformation of the R1 indole ring in a cleft-like formation of the A-loop, while the carbonyl oxygen (acceptor; A) secures a stable hydrogen bond with the donor group at M1156. We thereby propose a structural rationale at the molecular level behind the higher potency of the cyano derivative C11.

3. Materials and Methods

3.1. Generating 3D Chemical Structures and Pharmacophore Features of Allosteric Inhibitors

We first obtained the chemical structures of the series of indole-butyl amine derivatives (labeled 1–12) reported as a novel class of allosteric IGF1RK inhibitors [36] and generated their atomic coordinates using the ChemAxon MarvinSketch package, version 22.22 (https://www.chemaxon.com (12 October 2023)). The molecular structures were further categorized into four groups (Groups A–D) based on the ‘R’ group substituent (Figure 1). Furthermore, we generated the pharmacophore features for all inhibitors (1–12) using the Phase module of the Schrödinger software (version 2022.4) package [37,38]. The features include an H-bond donor (D), an H-bond acceptor (A), an aromatic ring (R), a hydrophobic site (H), a positive ionic (P) group, and a negative ionic (N) group.

3.2. Molecular Docking of Allosteric Inhibitors with IGF1RK

The crystal structure of IGF1RK in complex with the allosteric inhibitor C10 (PDB ID 3LW0) was obtained from the Protein Data Bank, which we used as a basis for our modeling of other inhibitor-docked structures [36]. We modeled the missing residues in the IGF1RK structure, A1097 through A1105 and G1169 through G1171, using the Chimera Interface to Modeller [39,40]. Eventually, we utilized the IGF1RK structure in complex with the inhibitor C10 for molecular docking of other inhibitors in the chemical series [36]. We optimized the structures of all inhibitors with correct molecular mechanics parameters before docking. The chemical structures were optimized using the OPLS3e force field with desalting and tautomer generation. The Epik module was used to generate all possible ionization states of the inhibitors at pH 7 ± 2, and at most 32 stereoisomers were generated for each inhibitor. All these ligand optimization steps were accomplished using the LigPrep module of the Schrödinger software package (version 2022.4) [41,42]. Next, the protein structure was optimized using the protein preparation wizard of Schrödinger [41]. The protein was first pre-processed by adding missing hydrogen atoms, removing water molecules, and adding any missing sidechains in residues. Next, H-bond optimization was performed to address any overlapping hydrogens added during pre-processing. Finally, a restrained minimization was carried out with the OPLS3e force field to remove steric clashes while allowing the displacement of heavy atoms up to 0.3 Å. Following minimization, the receptor grid generation panel was used to specify the protein as the receptor structure for docking and grid generation. A 20 Å3 receptor grid was generated around the IGF1RK allosteric pocket by selecting the bound inhibitor (C10) in the structure as the centroid of the grid. Finally, the docking of each inhibitor was performed in the extra-precision (XP) mode of the Glide module [43,44]. The binding affinity (docking score) and the energetically favorable binding conformations of each inhibitor were determined [45]. From the XP docking, we obtained two poses for each inhibitor in the binding pocket. Ultimately, for each inhibitor, the docked conformation with the lowest docking score (more negative) was selected for further MD simulations. The docked conformations were the initial models of interactions of inhibitors in the binding pocket, and further MD simulations were carried out to account for the flexibility of each inhibitor and the IGF1RK.

3.3. Molecular Dynamics (MD) Simulation Setup

We performed all-atom MD simulations of each IGF1RK-inhibitor complex using the GROMACS-2020.4 package [46]. The IGF1RK structure was parameterized using the AMBER99SB force field, while the inhibitors were parametrized with the general amber force field (GAFF) using the Antechamber module [47,48,49]. The protein and inhibitor coordinates were merged, and the inhibitor topology information was integrated into the system topology file. Subsequently, a cubic simulation domain was generated with the protein at the center and at least 10 Å from the edge of the simulation domain. The TIP3P model for water molecules was used to solvate each system, and then Na+ and Cl− ions were added to neutralize each system’s net charge [50]. Following that, each system was energy minimized for 50,000 steps using the steepest-descent algorithm [51]. The inhibitor molecules were subjected to a positional restraint prior to unrestrained system equilibration. Each system was equilibrated first in the NVT ensemble for 10 ns at 300 K using the velocity-rescale thermostat followed by equilibration in the NPT ensemble for 50 ns using the Berendsen barostat at 1 atm pressure and 300 K temperature [52,53]. Finally, restraints were removed, and a 200 ns production MD simulation was carried out for each IGF1RK/inhibitor complex. In all MD simulations, an integration time step of 2 fs was employed, and coordinates were saved every 100 ps.

3.4. Binding Free-Energy Calculations

We further computed the binding free energy of each inhibitor for the IGF1RK domain to provide a distinct comparison of their binding modes. We estimated the enthalpy change and solvation free energy for each complex using the Molecular Mechanics Generalized-Born Surface Area (MMGBSA) approach [54,55]. We used the gmxMMPBSA package, which integrates the AMBER tools MMPBSA.py script for performing free-energy calculations by utilizing MD trajectories [56,57]. We calculated the potential energy (bonded and non-bonded interactions), polar solvation energy, and non-polar solvation energy terms estimated over an average of 2000 frames in a 200 ns trajectory. The polar component was calculated utilizing the GB model with external and internal dielectric constants of 78.5 and 1.0, respectively. The non-polar solvation energy was estimated using the solvent-accessible surface area (SASA) approach, with the surface tension of the solvent set to 0.0072 and the surface tension offset set to 0.0. Furthermore, with per residue decomposition, the average contributions of residues to the binding free energy were also estimated.

3.5. Conformational and Interaction Analyses

The analyses of all MD trajectories were performed using the Gromacs built-in utilities from the Gromacs package (2020.4) and the Visual Molecular Dynamics (VMD) program [58,59]. We calculated the root mean squared deviation (RMSD) of the protein backbone atoms, taking the initial structure as a reference, and the root mean squared fluctuation (RMSF) was evaluated for the Cα atom of each residue from its mean position. The RMSD was calculated for the heavy atoms of the inhibitor, taking the initial docked conformation as a reference. The kernel density distribution function in MATLAB version 9.13.0.2193358 (R2022b) was used to evaluate the RMSD distribution curve [60].
We performed principal component analysis (PCA) to identify the essential conformational modes in the IGF1RK inhibitor complexes specifically at the allosteric pocket [61]. After building and diagonalizing a covariance matrix of the eigenvectors that represented the atomic variations in the protein and ligand atoms, a set of eigenvectors and their matching eigenvalues were obtained. The eigenvectors corresponding to essential movements in the protein–inhibitor complexes that have the greatest eigenvalues are called principal components, or PCs. We also constructed a free-energy surface (FES) projected along the top two eigenvectors corresponding to the first and second principal components, PC1 and PC2, respectively.
All hydrogen bonds (H-bonds) were computed with a donor–acceptor interatomic distance threshold of 3.5 Å in VMD [59]. The protein–inhibitor interactions were determined using the LigPlot+ program [62]. All the images were rendered using PyMOL software version 2.5.0 [63]. The MATLAB program was used to generate the plots [60].

4. Conclusions

We used pharmacophore modeling, molecular docking, and all-atom MD simulations to structurally rationalize the structure–activity relationships and binding affinities for a series of 12 allosteric inhibitors of IGF1RK. Our results provide structural explanations for the experimentally observed activities of the indole-butyl-amine derivatives against IGF1RK. We show that the key interaction with the V1063 residue was initially preserved in all complexes with a network of hydrophobic interactions involving the sulfur side chains of residues M1054 and M1079. Hence, the presence of an electron-donating group with two aromatic rings of the 5-cyano indole of the inhibitor is essential for these interactions. Our findings suggest that bulky carbonyl substitutions such as ketones, amides, or aldehydes at the R1 indole ring are less favorable for accommodating inhibitors in the IGF1RK allosteric pocket. The position of the electron-withdrawing (acceptor; A), electron-donating (donor; D), and hydrophobic (H) substituents at the R1 indole ring influence the binding conformations of allosteric inhibitors, impacting IGF1RK inhibition. Furthermore, we show that the binding of C11 induces the formation of a distinctive pocket with the A-loop and forms a unique H-bond with M1156. This allosteric pocket channel has a striking shape as well as chemical complementarity for the inhibitor C11. Thus, our results explain the reasoning underlying the higher potency of the cyano derivative C11. These structural investigations offer molecular explanations for variations in potencies of allosteric inhibitors and emphasize the conformational changes and molecular interactions that can be exploited for the design of further-improved allosteric inhibitors of IGF1RK and other kinases.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ijms25105368/s1.

Author Contributions

Conceptualization, J.V. and H.V.; methodology, J.V.; software, J.V.; validation, J.V.; formal analysis, J.V.; investigation, J.V.; data curation, J.V.; visualization, J.V.; writing—original draft preparation, J.V.; writing—review and editing, H.V.; supervision, H.V.; resources, H.V.; project administration, H.V.; funding acquisition, H.V. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Institutes of Health (NIH) through grant Nos. R35GM138217 and P20GM113131.

Data Availability Statement

The data that support the findings of this study are available in the supporting information of this article.

Acknowledgments

We acknowledge the financial support provided by the National Institutes of Health (NIH) through grants R35GM138217 and P20GM113131. The content is solely the responsibility of the authors and does not necessarily represent the official views of the NIH. We are grateful for computational support through Premise: a central shared HPC cluster at UNH supported by the Research Computing Center.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
CRCysteine-Rich Domain
EGFREpidermal Growth Factor Receptor
HER-2Human Epidermal Growth Factor Receptor 2
IGF1RInsulin-like Growth Factor 1 Receptor
IGF-1Insulin-like Growth Factor 1
IGF-2Insulin-like Growth Factor 2
IRInsulin Receptor
RTKReceptor Tyrosine Kinase
VEGFRVascular Endothelial Growth Factor Receptor

References

  1. Adams, T.; Epa, V.; Garrett, T.; Ward, C. Structure and function of the type 1 insulin-like growth factor receptor. Cell. Mol. Life Sci. 2000, 57, 1050–1093. [Google Scholar] [CrossRef] [PubMed]
  2. Chitnis, M.M.; Yuen, J.S.; Protheroe, A.S.; Pollak, M.; Macaulay, V.M. The type 1 insulin-like growth factor receptor pathway. Clin. Cancer Res. 2008, 14, 6364–6370. [Google Scholar] [CrossRef] [PubMed]
  3. Cabail, M.Z.; Li, S.; Lemmon, E.; Bowen, M.E.; Hubbard, S.R.; Miller, W.T. The insulin and IGF1 receptor kinase domains are functional dimers in the activated state. Nat. Commun. 2015, 6, 6406. [Google Scholar] [CrossRef] [PubMed]
  4. Garrett, T.P.; McKern, N.M.; Lou, M.; Frenkel, M.J.; Bentley, J.D.; Lovrecz, G.O.; Elleman, T.C.; Cosgrove, L.J.; Ward, C.W. Crystal structure of the first three domains of the type-1 insulin-like growth factor receptor. Nature 1998, 394, 395–399. [Google Scholar] [CrossRef] [PubMed]
  5. Li, J.; Choi, E.; Yu, H.; Bai, X.c. Structural basis of the activation of type 1 insulin-like growth factor receptor. Nat. Commun. 2019, 10, 4567. [Google Scholar] [CrossRef] [PubMed]
  6. Maruyama, I.N. Mechanisms of activation of receptor tyrosine kinases: Monomers or dimers. Cells 2014, 3, 304–330. [Google Scholar] [CrossRef] [PubMed]
  7. Larsson, O.; Girnita, A.; Girnita, L. Role of insulin-like growth factor 1 receptor signalling in cancer. Br. J. Cancer 2005, 92, 2097–2101. [Google Scholar] [CrossRef]
  8. Ekyalongo, R.C.; Yee, D. Revisiting the IGF-1R as a breast cancer target. NPJ Precis. Oncol. 2017, 1, 14. [Google Scholar] [CrossRef]
  9. Zhang, M.; Liu, J.; Li, M.; Zhang, S.; Lu, Y.; Liang, Y.; Zhao, K.; Li, Y. Insulin-like growth factor 1/insulin-like growth factor 1 receptor signaling protects against cell apoptosis through the PI3K/AKT pathway in glioblastoma cells. Exp. Ther. Med. 2018, 16, 1477–1482. [Google Scholar] [CrossRef]
  10. Sun, Y.; Sun, X.; Shen, B. Molecular imaging of IGF-1R in cancer. Mol. Imaging 2017, 16, 1536012117736648. [Google Scholar] [CrossRef]
  11. Dricu, A.; Kanter, L.; Wang, M.; Nilsson, G.; Hjertman, M.; Wejde, J.; Larsson, O. Expression of the insulin-like growth factor 1 receptor (IGF-1R) in breast cancer cells: Evidence for a regulatory role of dolichyl phosphate in the transition from an intracellular to an extracellular IGF-1 pathway. Glycobiology 1999, 9, 571–579. [Google Scholar] [CrossRef]
  12. Zhang, Y.; Gao, C.; Cao, F.; Wu, Y.; Chen, S.; Han, X.; Mo, J.; Qiu, Z.; Fan, W.; Zhou, P.; et al. Pan-cancer analysis of IGF-1 and IGF-1R as potential prognostic biomarkers and immunotherapy targets. Front. Oncol. 2021, 11, 755341. [Google Scholar] [CrossRef] [PubMed]
  13. Wang, P.; Mak, V.C.; Cheung, L.W. Drugging IGF-1R in cancer: New insights and emerging opportunities. Genes Dis. 2023, 10, 199–211. [Google Scholar] [CrossRef] [PubMed]
  14. Ofer, P.; Heidegger, I.; Eder, I.E.; Schöpf, B.; Neuwirt, H.; Geley, S.; Klocker, H.; Massoner, P. Both IGF1R and INSR knockdown exert antitumorigenic effects in prostate cancer in vitro and in vivo. Mol. Endocrinol. 2015, 29, 1694–1707. [Google Scholar] [CrossRef] [PubMed]
  15. Waraky, A.; Aleem, E.; Larsson, O. Downregulation of IGF-1 receptor occurs after hepatic linage commitment during hepatocyte differentiation from human embryonic stem cells. Biochem. Biophys. Res. Commun. 2016, 478, 1575–1581. [Google Scholar] [CrossRef] [PubMed]
  16. Khalil, A.; Jameson, M.J. Downregulation of IGF1R expression inhibits growth and enhances cisplatin sensitivity of head and neck squamous cell carcinoma cells in vitro. Horm. Cancer 2019, 10, 11–23. [Google Scholar] [CrossRef] [PubMed]
  17. Li, R.; Pourpak, A.; Morris, S.W. Inhibition of the insulin-like growth factor-1 receptor (IGF1R) tyrosine kinase as a novel cancer therapy approach. J. Med. Chem. 2009, 52, 4981–5004. [Google Scholar] [CrossRef]
  18. Press, M.F.; Lenz, H.J. EGFR, HER2 and VEGF pathways: Validated targets for cancer treatment. Drugs 2007, 67, 2045–2075. [Google Scholar] [CrossRef]
  19. Jentzsch, V.; Osipenko, L.; Scannell, J.W.; Hickman, J.A. Costs and causes of oncology drug attrition with the example of insulin-like growth factor-1 receptor inhibitors. JAMA Netw. Open 2023, 6, e2324977. [Google Scholar] [CrossRef]
  20. Langer, C.J.; Novello, S.; Park, K.; Krzakowski, M.; Karp, D.D.; Mok, T.; Benner, R.J.; Scranton, J.R.; Olszanski, A.J.; Jassem, J. Randomized, phase III trial of first-line figitumumab in combination with paclitaxel and carboplatin versus paclitaxel and carboplatin alone in patients with advanced non–small-cell lung cancer. J. Clin. Oncol. 2014, 32, 2059. [Google Scholar] [CrossRef]
  21. Sclafani, F.; Kim, T.Y.; Cunningham, D.; Kim, T.W.; Tabernero, J.; Schmoll, H.J.; Roh, J.K.; Kim, S.Y.; Park, Y.S.; Guren, T.K.; et al. A randomized phase II/III study of dalotuzumab in combination with cetuximab and irinotecan in chemorefractory, KRAS wild-type, metastatic colorectal cancer. J. Natl. Cancer Inst. 2015, 107, djv258. [Google Scholar] [CrossRef]
  22. DuBois, S.G.; Krailo, M.D.; Glade-Bender, J.; Buxton, A.; Laack, N.; Randall, R.L.; Chen, H.X.; Seibel, N.L.; Boron, M.; Terezakis, S.; et al. Randomized phase III trial of ganitumab with interval-compressed chemotherapy for patients with newly diagnosed metastatic Ewing sarcoma: A report from the Children’s Oncology Group. J. Clin. Oncol. 2023, 41, 2098–2107. [Google Scholar] [CrossRef]
  23. Fuchs, C.; Azevedo, S.; Okusaka, T.; Van Laethem, J.L.; Lipton, L.; Riess, H.; Szczylik, C.; Moore, M.; Peeters, M.; Bodoky, G.; et al. A phase 3 randomized, double-blind, placebo-controlled trial of ganitumab or placebo in combination with gemcitabine as first-line therapy for metastatic adenocarcinoma of the pancreas: The GAMMA trial. Ann. Oncol. 2015, 26, 921–927. [Google Scholar] [CrossRef]
  24. de Bono, J.; Lin, C.C.; Chen, L.T.; Corral, J.; Michalarea, V.; Rihawi, K.; Ong, M.; Lee, J.H.; Hsu, C.H.; Yang, J.C.H.; et al. Two first-in-human studies of xentuzumab, a humanised insulin-like growth factor (IGF)-neutralising antibody, in patients with advanced solid tumours. Br. J. Cancer 2020, 122, 1324–1332. [Google Scholar] [CrossRef]
  25. Chen, H.X.; Sharon, E. IGF-1R as an anti-cancer target—trials and tribulations. Chin. J. Cancer 2013, 32, 242. [Google Scholar] [CrossRef]
  26. Carboni, J.M.; Wittman, M.; Yang, Z.; Lee, F.; Greer, A.; Hurlburt, W.; Hillerman, S.; Cao, C.; Cantor, G.H.; Dell-John, J.; et al. BMS-754807, a small molecule inhibitor of insulin-like growth factor-1R/IR. Mol. Cancer Ther. 2009, 8, 3341–3349. [Google Scholar] [CrossRef]
  27. Gulbins, A.; Horstmann, M.; Flögel, U.; Oeverhaus, M.; Banga, J.P.; Krause, G.; Hammer, G.D.; Eckstein, A.; Görtz, G.E. Linsitinib, an IGF-1R inhibitor, attenuates disease development and progression in a model of thyroid eye disease. Front. Endocrinol. 2023, 14, 1211473. [Google Scholar] [CrossRef]
  28. Aiken, R.; Axelson, M.; Harmenberg, J.; Klockare, M.; Larsson, O.; Wassberg, C. Phase I clinical trial of AXL1717 for treatment of relapsed malignant astrocytomas: Analysis of dose and response. Oncotarget 2017, 8, 81501. [Google Scholar] [CrossRef]
  29. Buck, E.; Mulvihill, M. Small molecule inhibitors of the IGF-1R/IR axis for the treatment of cancer. Expert. Opin. Investig. Drugs 2011, 20, 605–621. [Google Scholar] [CrossRef]
  30. Baserga, R. The decline and fall of the IGF-I receptor. J. Cell. Physiol. 2013, 228, 675–679. [Google Scholar] [CrossRef]
  31. King, H.; Aleksic, T.; Haluska, P.; Macaulay, V.M. Can we unlock the potential of IGF-1R inhibition in cancer therapy? Cancer Treat. Rev. 2014, 40, 1096–1105. [Google Scholar] [CrossRef]
  32. Boone, D.N.; Lee, A.V. Targeting the insulin-like growth factor receptor: Developing biomarkers from gene expression profiling. Crit. Rev. Oncog. 2012, 17, 161–173. [Google Scholar] [CrossRef] [PubMed]
  33. Huggins, D.J.; Sherman, W.; Tidor, B. Rational approaches to improving selectivity in drug design. J. Med. Chem. 2012, 55, 1424–1444. [Google Scholar] [CrossRef]
  34. Wang, B.; Wu, H.; Hu, C.; Wang, H.; Liu, J.; Wang, W.; Liu, Q. An overview of kinase downregulators and recent advances in discovery approaches. Signal Transduct. Target. Ther. 2021, 6, 423. [Google Scholar] [CrossRef]
  35. Wu, P.; Clausen, M.H.; Nielsen, T.E. Allosteric small-molecule kinase inhibitors. Pharmacol. Ther. 2015, 156, 59–68. [Google Scholar] [CrossRef] [PubMed]
  36. Heinrich, T.; Grädler, U.; Böttcher, H.; Blaukat, A.; Shutes, A. Allosteric IGF-1R inhibitors. ACS Med. Chem. Lett. 2010, 1, 199–203. [Google Scholar] [CrossRef]
  37. Dixon, S.L.; Smondyrev, A.M.; Knoll, E.H.; Rao, S.N.; Shaw, D.E.; Friesner, R.A. PHASE: A new engine for pharmacophore perception, 3D QSAR model development, and 3D database screening: 1. Methodology and preliminary results. J. Comput. Aided Mol. Des. 2006, 20, 647–671. [Google Scholar] [CrossRef]
  38. Dixon, S.L.; Smondyrev, A.M.; Rao, S.N. PHASE: A novel approach to pharmacophore modeling and 3D database searching. Chem. Biol. Drug Des. 2006, 67, 370–372. [Google Scholar] [CrossRef]
  39. Pettersen, E.F.; Goddard, T.D.; Huang, C.C.; Couch, G.S.; Greenblatt, D.M.; Meng, E.C.; Ferrin, T.E. UCSF Chimera—A visualization system for exploratory research and analysis. J. Comput. Chem. 2004, 25, 1605–1612. [Google Scholar] [CrossRef]
  40. Yang, Z.; Lasker, K.; Schneidman-Duhovny, D.; Webb, B.; Huang, C.C.; Pettersen, E.F.; Goddard, T.D.; Meng, E.C.; Sali, A.; Ferrin, T.E. UCSF Chimera, MODELLER, and IMP: An integrated modeling system. J. Struct. Biol. 2012, 179, 269–278. [Google Scholar] [CrossRef]
  41. Madhavi Sastry, G.; Adzhigirey, M.; Day, T.; Annabhimoju, R.; Sherman, W. Protein and ligand preparation: Parameters, protocols, and influence on virtual screening enrichments. J. Comput. Aided Mol. Des. 2013, 27, 221–234. [Google Scholar] [CrossRef] [PubMed]
  42. Shelley, J.C.; Cholleti, A.; Frye, L.L.; Greenwood, J.R.; Timlin, M.R.; Uchimaya, M. Epik: A software program for pK a prediction and protonation state generation for drug-like molecules. J. Comput. Aided Mol. Des. 2007, 21, 681–691. [Google Scholar] [CrossRef] [PubMed]
  43. Halgren, T.A.; Murphy, R.B.; Friesner, R.A.; Beard, H.S.; Frye, L.L.; Pollard, W.T.; Banks, J.L. Glide: A new approach for rapid, accurate docking and scoring. 2. Enrichment factors in database screening. J. Med. Chem. 2004, 47, 1750–1759. [Google Scholar] [CrossRef] [PubMed]
  44. Friesner, R.A.; Banks, J.L.; Murphy, R.B.; Halgren, T.A.; Klicic, J.J.; Mainz, D.T.; Repasky, M.P.; Knoll, E.H.; Shelley, M.; Perry, J.K.; et al. Glide: A new approach for rapid, accurate docking and scoring. 1. Method and assessment of docking accuracy. J. Med. Chem. 2004, 47, 1739–1749. [Google Scholar] [CrossRef] [PubMed]
  45. Friesner, R.A.; Murphy, R.B.; Repasky, M.P.; Frye, L.L.; Greenwood, J.R.; Halgren, T.A.; Sanschagrin, P.C.; Mainz, D.T. Extra precision glide: Docking and scoring incorporating a model of hydrophobic enclosure for protein- ligand complexes. J. Med. Chem. 2006, 49, 6177–6196. [Google Scholar] [CrossRef] [PubMed]
  46. Berendsen, H.J.; 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]
  47. Hornak, V.; Abel, R.; Okur, A.; Strockbine, B.; Roitberg, A.; Simmerling, C. Comparison of multiple Amber force fields and development of improved protein backbone parameters. Proteins 2006, 65, 712–725. [Google Scholar] [CrossRef]
  48. 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] [PubMed]
  49. Wang, J.; Wang, W.; Kollman, P.A.; Case, D.A. Automatic atom type and bond type perception in molecular mechanical calculations. J. Mol. Graph. Modell. 2006, 25, 247–260. [Google Scholar] [CrossRef]
  50. Jorgensen, W.L.; Chandrasekhar, J.; Madura, J.D.; Impey, R.W.; Klein, M.L. Comparison of simple potential functions for simulating liquid water. J. Chem. Phys. 1983, 79, 926–935. [Google Scholar] [CrossRef]
  51. Chester, C.; Friedman, B.; Ursell, F. An extension of the method of steepest descents. In Mathematical Proceedings of the Cambridge Philosophical Society; Cambridge University Press: Cambridge, UK, 1957; Volume 53, pp. 599–611. [Google Scholar]
  52. Berendsen, H.J.; Postma, J.V.; Van Gunsteren, W.F.; DiNola, A.; Haak, J.R. Molecular dynamics with coupling to an external bath. J. Chem. Phys. 1984, 81, 3684–3690. [Google Scholar] [CrossRef]
  53. Bussi, G.; Donadio, D.; Parrinello, M. Canonical sampling through velocity rescaling. J. Chem. Phys. 2007, 126, 014101. [Google Scholar] [CrossRef] [PubMed]
  54. Tsui, V.; Case, D.A. Theory and applications of the generalized Born solvation model in macromolecular simulations. Biopolymers 2000, 56, 275–291. [Google Scholar] [CrossRef] [PubMed]
  55. Genheden, S.; Ryde, U. The MM/PBSA and MM/GBSA methods to estimate ligand-binding affinities. Expert Opin. Drug Discov. 2015, 10, 449–461. [Google Scholar] [CrossRef] [PubMed]
  56. Valdés-Tresanco, M.S.; Valdés-Tresanco, M.E.; Valiente, P.A.; Moreno, E. gmx_MMPBSA: A new tool to perform end-state free energy calculations with GROMACS. J. Chem. Theory Comput. 2021, 17, 6281–6291. [Google Scholar] [CrossRef] [PubMed]
  57. Miller, B.R., III; McGee, T.D., Jr.; Swails, J.M.; Homeyer, N.; Gohlke, H.; Roitberg, A.E. MMPBSA. py: An efficient program for end-state free energy calculations. J. Chem. Theory Comput. 2012, 8, 3314–3321. [Google Scholar] [CrossRef] [PubMed]
  58. Van Der Spoel, D.; Lindahl, E.; Hess, B.; Groenhof, G.; Mark, A.E.; Berendsen, H.J. GROMACS: Fast, flexible, and free. J. Comput. Chem. 2005, 26, 1701–1718. [Google Scholar] [CrossRef] [PubMed]
  59. Humphrey, W.; Dalke, A.; Schulten, K. VMD: Visual molecular dynamics. J. Mol. Graph. 1996, 14, 33–38. [Google Scholar] [CrossRef]
  60. The MathWorks Inc. MATLAB, version: 9.13. 0 (R2022b); The MathWorks Inc.: Natick, MA, USA, 2022. [Google Scholar]
  61. David, C.C.; Jacobs, D.J. Principal Component Analysis: A Method for Determining the Essential Dynamics of Proteins. In Protein Dynamics: Methods and Protocols; Livesay, D.R., Ed.; Humana Press: Totowa, NJ, USA, 2014; pp. 193–226. [Google Scholar]
  62. Laskowski, R.A.; Swindells, M.B. LigPlot+: Multiple Ligand–Protein Interaction Diagrams for Drug Discovery. J. Chem. Inf. Model. 2011, 51, 2778–2786. [Google Scholar] [CrossRef]
  63. Schrodinger, L. The PyMOL Molecular Graphics System; Schrödinger, Inc.: New York, NY, USA, 2015. [Google Scholar]
Ijms 25 05368 g001
Figure 1. Chemical structures of the indole-butyl-amine derivatives. The derivatives are categorized into four groups (Groups A–D) based on their R-group substitutions. Group A: indole-amine derivatives; Group B: alkyl derivatives; Group C: aldehyde/ketone derivatives; Group D: cyano derivatives.
Figure 1. Chemical structures of the indole-butyl-amine derivatives. The derivatives are categorized into four groups (Groups A–D) based on their R-group substitutions. Group A: indole-amine derivatives; Group B: alkyl derivatives; Group C: aldehyde/ketone derivatives; Group D: cyano derivatives.
Ijms 25 05368 g001
Ijms 25 05368 g002
Figure 2. Pharmacophore features in each of the derivatives. Each feature is represented by a unique color and letter. A: Acceptor; D: Donor; H: Hydrophobic; P: Positive ionic; R: Ring.
Figure 2. Pharmacophore features in each of the derivatives. Each feature is represented by a unique color and letter. A: Acceptor; D: Donor; H: Hydrophobic; P: Positive ionic; R: Ring.
Ijms 25 05368 g002
Ijms 25 05368 g003
Figure 3. Docked conformations of allosteric inhibitors in the IGF1RK binding pocket. (A) The X-ray crystal structure of the C10 inhibitor bound to the IGF1RK. (B) Structural superimposition of the docked conformations of 11 indole-butyl-amine derivatives on the X-ray crystal structure of IGF1RK bound to the C10 derivative. (C) The docked conformation of each inhibitor in the IGF1RK pocket.
Figure 3. Docked conformations of allosteric inhibitors in the IGF1RK binding pocket. (A) The X-ray crystal structure of the C10 inhibitor bound to the IGF1RK. (B) Structural superimposition of the docked conformations of 11 indole-butyl-amine derivatives on the X-ray crystal structure of IGF1RK bound to the C10 derivative. (C) The docked conformation of each inhibitor in the IGF1RK pocket.
Ijms 25 05368 g003
Ijms 25 05368 g004
Figure 4. RMSD data for the heavy atoms of each inhibitor. The snapshot for each inhibitor with its initial docked conformation (white) superimposed on the most populated conformation in MD simulations is shown.
Figure 4. RMSD data for the heavy atoms of each inhibitor. The snapshot for each inhibitor with its initial docked conformation (white) superimposed on the most populated conformation in MD simulations is shown.
Ijms 25 05368 g004
Ijms 25 05368 g005
Figure 5. The free-energy surface (FES) projections along two principal components (PC1 and PC2) for each inhibitor complex with IGF1RK. The minima corresponding to the lowest free energy in each plot are marked by red circle(s). The free energy is given in kcal/mol and is indicated by the color palette from lower (blue) to higher (yellow) values.
Figure 5. The free-energy surface (FES) projections along two principal components (PC1 and PC2) for each inhibitor complex with IGF1RK. The minima corresponding to the lowest free energy in each plot are marked by red circle(s). The free energy is given in kcal/mol and is indicated by the color palette from lower (blue) to higher (yellow) values.
Ijms 25 05368 g005
Ijms 25 05368 g006
Figure 6. Binding conformation of the indole-butyl-amine derivatives (C1–C12) in the allosteric pocket of IGF1RK. (A) Structure superposition of inhibitor conformations corresponding to the lowest free-energy minimum (see Figure 5) for each IGF1RK–inhibitor complex. (B) The unique binding mode of the most potent derivative C11 in the IGF1RK allosteric pocket. (C) A surface view of the binding pocket with the C11 derivative.
Figure 6. Binding conformation of the indole-butyl-amine derivatives (C1–C12) in the allosteric pocket of IGF1RK. (A) Structure superposition of inhibitor conformations corresponding to the lowest free-energy minimum (see Figure 5) for each IGF1RK–inhibitor complex. (B) The unique binding mode of the most potent derivative C11 in the IGF1RK allosteric pocket. (C) A surface view of the binding pocket with the C11 derivative.
Ijms 25 05368 g006
Ijms 25 05368 g007
Figure 7. The binding conformations and the interactions of the allosteric inhibitors C10 and C11 in the IGF1RK pocket. (A) The binding modes of C10 and C11 in the allosteric pocket, with the residues involved in the interactions. (B) The residue−level energy contributions ( Δ Gbind) for each complex. The residues with higher/maximum energy contributions are represented in darker shades. (C) The distance traces for the unique hydrogen bonds between the IGF1RK residues and inhibitor atoms.
Figure 7. The binding conformations and the interactions of the allosteric inhibitors C10 and C11 in the IGF1RK pocket. (A) The binding modes of C10 and C11 in the allosteric pocket, with the residues involved in the interactions. (B) The residue−level energy contributions ( Δ Gbind) for each complex. The residues with higher/maximum energy contributions are represented in darker shades. (C) The distance traces for the unique hydrogen bonds between the IGF1RK residues and inhibitor atoms.
Ijms 25 05368 g007
Table 1. IGF1RK allosteric inhibitors with their respective IC50 (µM) values, as reported by Henrich et al. [36]. We categorized these derivates into four groups (Groups A–D) based on their ‘R’ group substituent.
Table 1. IGF1RK allosteric inhibitors with their respective IC50 (µM) values, as reported by Henrich et al. [36]. We categorized these derivates into four groups (Groups A–D) based on their ‘R’ group substituent.
GroupCompoundNameIC50 (µM)
AC11H-Indole-6-carboxylic acid 1-[4-(5-cyano-1H-indol-3-yl)-butyl]-piperidin-4-yl-amide10 ± 0.09
C21H-Indole-7-carboxylic acid 1-[4-(5-cyano-1H-indol-3-yl) butyl]-piperidin-4-yl-amide2.5 ± 0.002
C31H-Indole-4-carboxylic acid 1-[4-(5-cyano-1H-indol-3-yl)-butyl]-piperidin-4-yl-amide6 ± 0.20
BC41-Methyl-1H-indole-4-carboxylic acid 1-[4-(5-cyano-1H-indol-3-yl)-butyl]-piperidin-4-yl-amide2.7 ± 1.10
C51-Ethyl-1H-indole-4-carboxylic acid 1-[4-(5-cyano-1H-indol-3-yl)-butyl]-piperidin-4-yl-amide3.6 ± 1.10
CC63-Formyl-1H-indole-7-carboxylic acid 1-[4-(5-cyano-1H-indol-3-yl)-butyl]-piperidin-4-yl-amide0.9 ± 0.01
C73-Acetyl-1H-indole-7-carboxylic acid 1-[4-(5-cyano-1H-indol-3-yl)-butyl]-piperidin-4-yl-amide4.6 ± 0.08
C83-(2,2,2-Trifluoro-acetyl)-1H-indole-7-carboxylic acid 1-[4-(5-cyano-1H-indol-3-yl)-butyl]-piperidin-4-yl-amide4.0 ± 0.90
C9N7-1-[4-(5-cyano-1H-indol-3-yl)butyl]piperidin-4-yl-1H-indole-3,7-dicarboxamide5.8 ± 0.40
DC103-Cyano-1H-indole-7-carboxylic acid 1-[4-(5-cyano-1H-indol-3-yl)-butyl]-piperidin-4-yl-amide0.4 ± 0.07
C113-Cyano-5-fluoro-1H-indole-7-carboxylic acid 1-[4-(5-cyano-1H-indol-3-yl)-butyl]-piperidin-4-yl-amide0.2 ± 0.07
C123-Cyano-1H-indole-4-carboxylic acid 1-[4-(5-cyano-1H-indol-3-yl)-butyl]-piperidin-4-yl-amide3.5 ± 2.50
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.

Share and Cite

MDPI and ACS Style

Verma, J.; Vashisth, H. Structural Models for a Series of Allosteric Inhibitors of IGF1R Kinase. Int. J. Mol. Sci. 2024, 25, 5368. https://doi.org/10.3390/ijms25105368

AMA Style

Verma J, Vashisth H. Structural Models for a Series of Allosteric Inhibitors of IGF1R Kinase. International Journal of Molecular Sciences. 2024; 25(10):5368. https://doi.org/10.3390/ijms25105368

Chicago/Turabian Style

Verma, Jyoti, and Harish Vashisth. 2024. "Structural Models for a Series of Allosteric Inhibitors of IGF1R Kinase" International Journal of Molecular Sciences 25, no. 10: 5368. https://doi.org/10.3390/ijms25105368

APA Style

Verma, J., & Vashisth, H. (2024). Structural Models for a Series of Allosteric Inhibitors of IGF1R Kinase. International Journal of Molecular Sciences, 25(10), 5368. https://doi.org/10.3390/ijms25105368

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