skip to main content


Title: Coupling of conformational dynamics and inhibitor binding in the phosphodiesterase‐5 family
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.

 
more » « less
Award ID(s):
1757371
NSF-PAR ID:
10436197
Author(s) / Creator(s):
 ;  ;  
Publisher / Repository:
Wiley Blackwell (John Wiley & Sons)
Date Published:
Journal Name:
Protein Science
Volume:
32
Issue:
8
ISSN:
0961-8368
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    Structure and functions of S100 proteins are regulated by two distinct calcium binding EF hand motifs. In this work, we used solution‐state NMR spectroscopy to investigate the cooperativity between the two calcium binding sites and map the allosteric changes at the target binding site. To parse the contribution of the individual calcium binding events, variants of S100A12 were designed to selectively bind calcium to either the EF‐I (N63A) or EF‐II (E31A) loop, respectively. Detailed analysis of the backbone chemical shifts for wildtype protein and its mutants indicates that calcium binding to the canonical EF‐II loop is the principal trigger for the conformational switch between ‘closed’ apo to the ‘open’ Ca2+‐bound conformation of the protein. Elimination of binding in S100‐specific EF‐I loop has limited impact on the calcium binding affinity of the EF‐II loop and the concomitant structural rearrangement. In contrast, deletion of binding in the EF‐II loop significantly attenuates calcium affinity in the EF‐I loop and the structure adopts a ‘closed’ apo‐like conformation. Analysis of experimental amide nitrogen (15N) relaxation rates (R1,R2, and15N–{1H} NOE) and molecular dynamics (MD) simulations demonstrate that the calcium bound state is relatively floppy with pico–nanosecond motions induced in functionally relevant domains responsible for target recognition such as the hinge domain and the C‐terminal residues. Experimental relaxation studies combined with MD simulations show that while calcium binding in the EF‐I loop alone does not induce significant motions in the polypeptide chain, EF‐I regulates fluctuations in the polypeptide in the presence of bound calcium in the EF‐II loop. These results offer novel insights into the dynamic regulation of target recognition by calcium binding and unravels the role of cooperativity between the two calcium binding events in S100A12.

     
    more » « less
  2. The β2adrenergic receptor (β2AR) is an archetypal G protein coupled receptor (GPCR). One structural signature of GPCR activation is a large-scale movement (ca. 6 to 14 Å) of transmembrane helix 6 (TM6) to a conformation which binds and activates a cognate G protein. The β2AR exhibits a low level of agonist-independent G protein activation. The structural origin of this basal activity and its suppression by inverse agonists is unknown but could involve a unique receptor conformation that promotes G protein activation. Alternatively, a conformational selection model proposes that a minor population of the canonical active receptor conformation exists in equilibrium with inactive forms, thus giving rise to basal activity of the ligand-free receptor. Previous spin-labeling and fluorescence resonance energy transfer experiments designed to monitor the positional distribution of TM6 did not detect the presence of the active conformation of ligand-free β2AR. Here we employ spin-labeling and pressure-resolved double electron–electron resonance spectroscopy to reveal the presence of a minor population of unliganded receptor, with the signature outward TM6 displacement, in equilibrium with inactive conformations. Binding of inverse agonists suppresses this population. These results provide direct structural evidence in favor of a conformational selection model for basal activity in β2AR and provide a mechanism for inverse agonism. In addition, they emphasize 1) the importance of minor populations in GPCR catalytic function; 2) the use of spin-labeling and variable-pressure electron paramagnetic resonance to reveal them in a membrane protein; and 3) the quantitative evaluation of their thermodynamic properties relative to the inactive forms, including free energy, partial molar volume, and compressibility.

     
    more » « less
  3. Abstract

    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.

     
    more » « less
  4. Mechanistic understanding of DNA recombination in the Cre-loxPsystem has largely been guided by crystallographic structures of tetrameric synaptic complexes. Those studies have suggested a role for protein conformational dynamics that has not been well characterized at the atomic level. We used solution nuclear magnetic resonance (NMR) spectroscopy to discover the link between intrinsic flexibility and function in Cre recombinase. Transverse relaxation-optimized spectroscopy (TROSY) NMR spectra show the N-terminal and C-terminal catalytic domains (CreNTDand CreCat) to be structurally independent. Amide15N relaxation measurements of the CreCatdomain reveal fast-timescale dynamics in most regions that exhibit conformational differences in active and inactive Cre protomers in crystallographic tetramers. However, the C-terminal helix αN, implicated in assembly of synaptic complexes and regulation of DNA cleavage activity viatransprotein–protein interactions, is unexpectedly rigid in free Cre. Chemical shift perturbations and intra- and intermolecular paramagnetic relaxation enhancement (PRE) NMR data reveal an alternative autoinhibitory conformation for the αN region of free Cre, wherein it packsin cisover the protein DNA binding surface and active site. Moreover, binding toloxPDNA induces a conformational change that dislodges the C terminus, resulting in acis-to-transswitch that is likely to enable protein–protein interactions required for assembly of recombinogenic Cre intasomes. These findings necessitate a reexamination of the mechanisms by which this widely utilized gene-editing tool selects target sites, avoids spurious DNA cleavage activity, and controls DNA recombination efficiency.

     
    more » « less
  5. Abstract

    Human acyl protein thioesterases (APTs) catalyze the depalmitoylation ofS‐acylated proteins attached to the plasma membrane, facilitating reversible cycles of membrane anchoring and detachment. We previously showed that a bacterial APT homologue, FTT258 from the gram‐negative pathogenFrancisella tularensis, exists in equilibrium between a closed and open state based on the structural dynamics of a flexible loop overlapping its active site. Although the structural dynamics of this loop are not conserved in human APTs, the amino acid sequence of this loop is highly conserved, indicating essential but divergent functions for this loop in human APTs. Herein, we investigated the role of this loop in regulating the catalytic activity, ligand binding, and protein folding of human APT1, a depalmitoylase connected with cancer, immune, and neurological signaling. Using a combination of substitutional analysis with kinetic, structural, and biophysical characterization, we show that even in its divergent structural location in human APT1 that this loop still regulates the catalytic activity of APT1 through contributions to ligand binding and substrate positioning. We confirmed previously known roles for multiple residues (Phe72 and Ile74) in substrate binding and catalysis while adding new roles in substrate selectivity (Pro69), in catalytic stabilization (Asp73 and Ile75), and in transitioning between the membrane binding β‐tongue and substrate‐binding loops (Trp71). Even conservative substitution of this tryptophan (Trp71) fulcrum led to complete loss of catalytic activity, a 13°C decrease in total protein stability, and drastic drops in ligand affinity, indicating that the combination of the size, shape, and aromaticity of Trp71 are essential to the proper structure of APT1. Mixing buried hydrophobic surface area with contributions to an exposed secondary surface pocket, Trp71 represents a previously unidentified class of essential tryptophans within α/β hydrolase structure and a potential allosteric binding site within human APTs.

     
    more » « less