Understanding kinase‐inhibitor selectivity continues to be a major objective in kinase drug discovery. We probe the molecular basis of selectivity of an allosteric inhibitor (MSC1609119A‐1) of the insulin‐like growth factor‐I receptor kinase (IGF1RK), which has been shown to be ineffective for the homologous insulin receptor kinase (IRK). Specifically, we investigated the structural and energetic basis of the allosteric binding of this inhibitor to each kinase by combining molecular modeling, molecular dynamics (MD) simulations, and thermodynamic calculations. We predict the inhibitor conformation in the binding pocket of IRK and highlight that the charged residues in the histidine‐arginine‐aspartic acid (HRD) and aspartic acid‐phenylalanine‐glycine (DFG) motifs and the nonpolar residues in the binding pocket govern inhibitor interactions in the allosteric pocket of each kinase. We suggest that the conformational changes in the IGF1RK residues M1054 and M1079, movement of the ⍺C‐helix, and the conformational stabilization of the DFG motif favor the selectivity of the inhibitor toward IGF1RK. Our thermodynamic calculations reveal that the observed selectivity can be rationalized through differences observed in the electrostatic interaction energy of the inhibitor in each inhibitor/kinase complex and the hydrogen bonding interactions of the inhibitor with the residue V1063 in IGF1RK that are not attained with the corresponding residue V1060 in IRK. Overall, our study provides a rationale for the molecular basis of recognition of this allosteric inhibitor by IGF1RK and IRK, which is potentially useful in developing novel inhibitors with improved affinity and selectivity.
- Award ID(s):
- 2137558
- NSF-PAR ID:
- 10349233
- Date Published:
- Journal Name:
- Cells
- Volume:
- 11
- Issue:
- 15
- ISSN:
- 2073-4409
- Page Range / eLocation ID:
- 2451
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
Molecular basis for differential recognition of an allosteric inhibitor by receptor tyrosine kinases
Abstract -
Abstract Phosphodiesterase‐5 (PDE5) is responsible for regulating the concentration of the second messenger molecule cGMP by hydrolyzing it into 5′‐GMP. PDE5 is implicated in erectile dysfunction and cardiovascular diseases. The substrate binding site in the catalytic domain of PDE5 is surrounded by several dynamic structural motifs (including the α14 helix, M‐loop, and H‐loop) that are known to switch between inactive and active conformational states via currently unresolved structural intermediates. We evaluated the conformational dynamics of these structural motifs in the apo state and upon binding of an allosteric inhibitor (
evodiamine ) oravanafil , a competitive inhibitor. We employed enhanced sampling‐based replica exchange solute scaling (REST2) method, principal component analysis (PCA), time‐lagged independent component analysis (tICA), molecular dynamics (MD) simulations, and well‐tempered metadynamics simulations to probe the conformational changes in these structural motifs. Our results support a regulatory mechanism for PDE5, where the α14 helix alternates between an inward (lower activity) conformation and an outward (higher activity) conformation that is accompanied by the folding/unfolding of the α8′ and α8″ helices of the H‐loop. When the allosteric inhibitor evodiamine is bound to PDE5, the inward (inactive) state of the α14 helix is preferred, thus preventing substrate access to the catalytic site. In contrast, competitive inhibitors of PDE5 block catalysis by occupying the active site accompanied by stabilization of the outward conformation of the α14 helix. Defining the conformational dynamics underlying regulation of PDE5 activation will be helpful in rational design of next‐generation small molecules modulators of PDE5 activity. -
Abstract Previous studies of the N-terminal PDZ tandem from PSD-95 produced divergent models and failed to identify interdomain contacts stabilizing the structure. We used ensemble and single-molecule FRET along with replica-exchange molecular dynamics to fully characterize the energy landscape. Simulations and experiments identified two conformations: an open-like conformation with a small contact interface stabilized by salt bridges, and a closed-like conformation with a larger contact interface stabilized by surface-exposed hydrophobic residues. Both interfaces were confirmed experimentally. Proximity of interdomain contacts to the binding pockets may explain the observed coupling between conformation and binding. The low-energy barrier between conformations allows submillisecond dynamics, which were time-averaged in previous NMR and FRET studies. Moreover, the small contact interfaces were likely overridden by lattice contacts as crystal structures were rarely sampled in simulations. Our hybrid approach can identify transient interdomain interactions, which are abundant in multidomain proteins yet often obscured by dynamic averaging.
-
Gram-positive pathogenic bacteria Staphylococcus express and secret staphylococcal peroxidase inhibitor (SPIN) proteins to help evade neutrophil-mediated immunity by inhibiting the activity of the main oxidative-defense player myeloperoxidase (MPO) enzyme. SPIN contains a structured 3-helix bundle C-terminal domain, which can specifically bind to MPO with high affinity, and an intrinsically disordered N-terminal domain (NTD), which folds into a structured β-hairpin and inserts itself into the active site of MPO for inhibition. Mechanistic insights of the coupled folding and binding process are needed in order to better understand how residual structures and/or conformational flexibility of NTD contribute to the different strengths of inhibition of SPIN homologs. In this work, we applied atomistic molecular dynamics simulations on two SPIN homologs, from S. aureus and S. delphini , respectively, which share high sequence identity and similarity, to explore the possible mechanistic basis for their different inhibition efficacies on human MPO. Direct simulations of the unfolding and unbinding processes at 450 K reveal that these two SPIN/MPO complexes systems follow surprisingly different mechanisms of coupled binding and folding. While coupled binding and folding of SPIN- aureus NTD is highly cooperative, SPIN- delphini NTD appears to mainly utilize a conformational selection-like mechanism. These observations are in contrast to an overwhelming prevalence of induced folding-like mechanisms for intrinsically disordered proteins that fold into helical structures upon binding. Further simulations of unbound SPIN NTDs at room temperature reveal that SPIN- delphini NTD has a much stronger propensity of forming β-hairpin like structures, consistent with its preference to fold and then bind. These may help explain why the inhibition strength is not well correlated with binding affinity for different SPIN homologs. Altogether, our work establishes the relationship between the residual conformational stability of SPIN-NTD and their inhibitory function, which can help us develop new strategies towards treating Staphylococcal infections.more » « less
-
Dynamic allosterism allows the propagation of signal throughout a protein. The PDZ (PSD-95/Dlg1/ZO-1) family has been named as a classic example of dynamic allostery in small modular domains. While the PDZ family consists of more than 200 domains, previous efforts have primarily focused on a few well-studied PDZ domains, including PTP-BL PDZ2, PSD-95 PDZ3, and Par6 PDZ. Taken together, experimental and computational studies have identified regions of these domains that are dynamically coupled to ligand binding. These regions include the αA helix, the αB lower-loop, and the αC helix. In this review, we summarize the specific residues on the αA helix, the αB lower-loop, and the αC helix of PTP-BL PDZ2, PSD-95 PDZ3, and Par6 PDZ that have been identified as participants in dynamic allostery by either experimental or computational approaches. This review can serve as an index for researchers to look back on the previously identified allostery in the PDZ family. Interestingly, our summary of previous work reveals clear consistencies between the domains. While the PDZ family has a low sequence identity, we show that some of the most consistently identified allosteric residues within PTP-BL PDZ2 and PSD-95 PDZ3 domains are evolutionarily conserved. These residues include A46/A347, V61/V362, and L66/L367 on PTP-BL PDZ2 and PSD-95 PDZ3, respectively. Finally, we expose a need for future work to explore dynamic allostery within (1) PDZ domains with multiple binding partners and (2) multidomain constructs containing a PDZ domain.more » « less