Flexibility and inhibitor binding in cdc25 phosphatases

Cdc25 phosphatases involved in cell cycle checkpoints are now active targets for the development of anti‐cancer therapies. Rational drug design would certainly benefit from detailed structural information for Cdc25s. However, only apo‐ or sulfate‐bound crystal structures of the Cdc25 catalytic domain have been described so far. Together with previously available crystalographic data, results from molecular dynamics simulations, bioinformatic analysis, and computer‐generated conformational ensembles shown here indicate that the last 30–40 residues in the C‐terminus of Cdc25B are partially unfolded or disordered in solution. The effect of C‐terminal flexibility upon binding of two potent small molecule inhibitors to Cdc25B is then analyzed by using three structural models with variable levels of flexibility, including an equilibrium distributed ensemble of Cdc25B backbone conformations. The three Cdc25B structural models are used in combination with flexible docking, clustering, and calculation of binding free energies by the linear interaction energy approximation to construct and validate Cdc25B‐inhibitor complexes. Two binding sites are identified on top and beside the Cdc25B active site. The diversity of interaction modes found increases with receptor flexibility. Backbone flexibility allows the formation of transient cavities or compact hydrophobic units on the surface of the stable, folded protein core that are unexposed or unavailable for ligand binding in rigid and densely packed crystal structures. The present results may help to speculate on the mechanisms of small molecule complexation to partially unfolded or locally disordered proteins. Proteins 2010. © 2010 Wiley‐Liss, Inc.     

Solution NMR studies reveal no global flexibility in the catalytic domain of CDC25B

The CDC25B phosphatase is a critical regulator of the cell cycle and has been validated as an important therapeutic target in cancer. Previous studies using molecular dynamics simulations have concluded that the catalytic domain of CDC25B may experience a significant degree of dynamics or be partially disordered in solution, a finding that has a pronounced impact on the structure‐based development of CDC25B inhibitors. We have probed the backbone dynamics of the CDC25B catalytic domain in solution using NMR relaxation experiments and found that the core of the protein is relatively rigid and does not experience any large‐scale dynamics over a broad range of time scales. Furthermore, based on residual dipolar coupling measurements we have concluded that the conformation in solution is very similar to that observed in the crystal form. Importantly, these findings rationalize the application of the CDC25B crystal structure in structure‐based drug development. Proteins 2014; 82:2889–2895. © 2014 Wiley Periodicals, Inc.     

Characterization of the flexibility of the peripheral stalk of prokaryotic rotary A‐ATPases by atomistic simulations

Rotary ATPases are involved in numerous physiological processes, with the three distinct types (F/A/V‐ATPases) sharing functional properties and structural features. The basic mechanism involves the counter rotation of two motors, a soluble ATP hydrolyzing/synthesizing domain and a membrane‐embedded ion pump connected through a central rotor axle and a stator complex. Within the A/V‐ATPase family conformational flexibility of the EG stators has been shown to accommodate catalytic cycling and is considered to be important to function. For the A‐ATPase three EG structures have been reported, thought to represent conformational states of the stator during different stages of rotary catalysis. Here we use long, detailed atomistic simulations to show that those structures are conformers explored through thermal fluctuations, but do not represent highly populated states of the EG stator in solution. We show that the coiled coil tail domain has a high persistence length (～100 nm), but retains the ability to adapt to different conformational states through the presence of two hinge regions. Moreover, the stator network of the related V‐ATPase has been suggested to adapt to subunit interactions in the collar region in addition to the nucleotide occupancy of the catalytic domain. The MD simulations reported here, reinforce this observation showing that the EG stators have enough flexibility to adapt to significantly different structural re‐arrangements and accommodate structural changes in the catalytic domain whilst resisting the large torque generated by catalytic cycling. These results are important to understand the role the stators play in the rotary‐ATPase mechanism. Proteins 2016; 84:1203–1212. © 2016 The Authors. Proteins: Structure, Function, and Bioinformatics Published by Wiley Periodicals, Inc.     

Segment swapping aided the evolution of enzyme function: The case of uroporphyrinogen III synthase

In an earlier study, we showed that two‐domain segment‐swapped proteins can evolve by domain swapping and fusion, resulting in a protein with two linkers connecting its domains. We proposed that a potential evolutionary advantage of this topology may be the restriction of interdomain motions, which may facilitate domain closure by a hinge‐like movement, crucial for the function of many enzymes. Here, we test this hypothesis computationally on uroporphyrinogen III synthase, a two‐domain segment‐swapped enzyme essential in porphyrin metabolism. To compare the interdomain flexibility between the wild‐type, segment‐swapped enzyme (having two interdomain linkers) and circular permutants of the same enzyme having only one interdomain linker, we performed geometric and molecular dynamics simulations for these species in their ligand‐free and ligand‐bound forms. We find that in the ligand‐free form, interdomain motions in the wild‐type enzyme are significantly more restricted than they would be with only one interdomain linker, while the flexibility difference is negligible in the ligand‐bound form. We also estimated the entropy costs of ligand binding associated with the interdomain motions, and find that the change in domain connectivity due to segment swapping results in a reduction of this entropy cost, corresponding to ～20% of the total ligand binding free energy. In addition, the restriction of interdomain motions may also help the functional domain‐closure motion required for catalysis. This suggests that the evolution of the segment‐swapped topology facilitated the evolution of enzyme function for this protein by influencing its dynamic properties. Proteins 2016; 85:46–53. © 2016 Wiley Periodicals, Inc.     

An unexpected role for PI4,5P2 in EGF receptor endosomal trafficking

In mammalian cells, three Cdc25 phosphatases A, B, C coordinate cell cycle progression through activating dephosphorylation of Cyclin-dependent kinases. Whereas Cdc25B is believed to trigger entry into mitosis, Cdc25C is thought to act at a later stage of mitosis and in the nucleus. We report that a fraction of Cdc25C localises to centrosomes in a cell cycle-dependent fashion, as of late S phase and throughout G2 and mitosis. Moreover, Cdc25C colocalises with Cyclin B1 at centrosomes in G2 and in prophase and Fluorescence Recovery after Photobleaching experiments reveal that they are both in dynamic exchange between the centrosome and the cytoplasm. The centrosomal localisation of Cdc25C is essentially mediated by its catalytic C-terminal domain, but does not require catalytic activity. In fact phosphatase-dead and substrate-binding hotspot mutants of Cdc25C accumulate at centrosomes together with phosphoTyr15-Cdk1 and behave as dominant negative forms that impair entry into mitosis. Taken together, our data suggest an unexpected function for Cdc25C at the G2/M transition, in dephosphorylation of Cdk1. We propose that Cdc25C may participate in amplification of Cdk1-Cyclin B1 activity following initial activation by Cdc25B, and that this process is initiated at the centrosome, then further propagated throughout the cytoplasm thanks to the dynamic behavior of both Cdc25C and Cyclin B1.     

Crystal structure of the catalytic subunit of Cdc25B required for G2/M phase transition of the cell cycle.

Cdc25B is a dual specificity phosphatase involved in the control of cyclin-dependent kinases and the progression of cells through the cell cycle. A series of minimal domain Cdc25B constructs maintaining catalytic activity have been expressed. The structure of a minimum domain construct binding sulfate was determined at 1.9 A resolution and a temperature of 100 K. Other forms of the same co?nstruct were determined at lower resolution and room temperature. The overall folding and structure of the domain is similar to that found for Cdc25A. An important difference between the two is that the Cdc25B domain binds oxyanions in the catalytic site while that of Cdc25A appears unable to bind oxyanions. There are also important conformational differences in the C-terminal region. In Cdc25B, both sulfate and tungstate anions are shown to bind in the catalytic site containing the signature motif (HCxxxxxR) in a conformation similar to that of other protein tyrosine phosphatases and dual specificity phosphatases, with the exception of the Cdc25A. The Cdc25B constructs, with various truncations of the C-terminal residues, are shown to have potent catalytic activity. When cut back to the site at which the Cdc25A structure begins to deviate from the Cdc25B structure, the activity is considerably less. There is a pocket extending from the catalytic site to an anion-binding site containing a chloride about 14 A away. The catalytic cysteine residue, Cys473, can be oxidized to form a disulfide linkage to Cys426. A readily modifiable cysteine residue, Cys484, resides in another pocket that binds a sulfate but not in the signature motif conformation. This region of the structure is highly conserved between the Cdc25 molecules and could serve some unknown function.     

A molecular dynamics study of the binary complexes of APP,JIP1, and the cargo binding domain of KLC

Mutations in the amyloid precursor protein (APP) are responsible for the formation of amyloid‐β peptides. These peptides play a role in Alzheimer's and other dementia‐related diseases. The cargo binding domain of the kinesin‐1 light chain motor protein (KLC1) may be responsible for transporting APP either directly or via interaction with C‐jun N‐terminal kinase‐interacting protein 1 (JIP1). However, to date there has been no direct experimental or computational assessment of such binding at the atomistic level. We used molecular dynamics and free energy estimations to gauge the affinity for the binary complexes of KLC1, APP, and JIP1. We find that all binary complexes (KLC1:APP, KLC1:JIP1, and APP:JIP1) contain conformations with favorable binding free energies. For KLC1:APP the inclusion of approximate entropies reduces the favorability. This is likely due to the flexibility of the 42‐residue APP protein. In all cases we analyze atomistic/residue driving forces for favorable interactions. Proteins 2017; 85:221–234. © 2016 Wiley Periodicals, Inc.     

The C-terminal tail of the dual-specificity Cdc25B phosphatase mediates modular substrate recognition.

Cdc25 is a dual-specificity phosphatase that catalyzes the activation of the cyclin-dependent kinases (Cdk/cyclins), thus triggering initiation and progression of successive phases of the cell cycle. In our efforts to elucidate the interaction between Cdc25B and the natural substrate, bis-phosphorylated Cdk2/CycA (Cdk2-pTpY/CycA), we have previously found that the 17 residues of the C-terminal tail mediate a factor of 10 in substrate recognition. In the studies reported here, we localize the majority of this interaction using site-directed mutagenesis to two arginine residues (Arg556 and Arg562) located within this C-terminal region. We also show that the catalytic domain of Cdc25C, which differs most significantly from Cdc25B in this tail region, has a 100-fold lower activity toward Cdk2-pTpY/CycA. We further demonstrate that the proper presentation of the C-terminal tail of Cdc25B can be achieved in a "gain-of-function" chimeric protein consisting of the C-terminal tail of Cdc25B fused onto the catalytic core of Cdc25C. The >10-fold increase in activity seen only in the chimeric protein containing the two critical arginine residues demonstrates that the modular C-terminal tail of Cdc25B is the basis for most of the catalytic advantage of Cdc25B versus Cdc25C toward the Cdk2-pTpY/CycA substrate.     

The crystal structure of Clostridium perfringens SleM,a muramidase involved in cortical hydrolysis during spore germination

Clostridium perfringens spores employ two peptidoglycan lysins to degrade the spore cortex during germination. SleC initiates cortex hydrolysis to generate cortical fragments that are degraded further by the muramidase SleM. Here, we present the crystal structure of the C. perfringens S40 SleM protein at 1.8 Å. SleM comprises an N‐terminal catalytic domain that adopts an irregular α/β‐barrel fold that is common to GH25 family lysozymes, plus a C‐terminal fibronectin type III domain. The latter is involved in forming the SleM dimer that is evident in both the crystal structure and in solution. A truncated form of SleM that lacks the FnIII domain shows reduced activity against spore sacculi indicating that this domain may have a role in facilitating the position of substrate with respect to the enzyme's active site. Proteins 2016; 84:1681–1689. © 2016 Wiley Periodicals, Inc.     

Functional role of TRIM E3 ligase oligomerization and regulation of catalytic activity

TRIM E3 ubiquitin ligases regulate a wide variety of cellular processes and are particularly important during innate immune signalling events. They are characterized by a conserved tripartite motif in their N‐terminal portion which comprises a canonical RING domain, one or two B‐box domains and a coiled‐coil region that mediates ligase dimerization. Self‐association via the coiled‐coil has been suggested to be crucial for catalytic activity of TRIMs; however, the precise molecular mechanism underlying this observation remains elusive. Here, we provide a detailed characterization of the TRIM ligases TRIM25 and TRIM32 and show how their oligomeric state is linked to catalytic activity. The crystal structure of a complex between the TRIM25 RING domain and an ubiquitin‐loaded E2 identifies the structural and mechanistic features that promote a closed E2~Ub conformation to activate the thioester for ubiquitin transfer allowing us to propose a model for the regulation of activity in the full‐length protein. Our data reveal an unexpected diversity in the self‐association mechanism of TRIMs that might be crucial for their biological function.     

The dual specificity phosphatase Cdc25B, but not the closely related Cdc25C, is capable of inhibiting cellular proliferation in a manner dependent upon its catalytic activity

Cdc25B and Cdc25C are closely related dual specificity phosphatases that activate cyclin-dependent kinases by removal of inhibitory phosphorylations, thereby triggering entry into mitosis. Cdc25B, but not Cdc25C, has been implicated as an oncogene and been shown to be overexpressed in a variety of human tumors. Surprisingly, ectopic expression of Cdc25B, but not Cdc25C, inhibits cell proliferation in long term assays. Chimeric proteins generated from the two phosphatases show that the anti-proliferative activity is associated with the C-terminal end of Cdc25B. Indeed, the catalytic domain of Cdc25B is sufficient to suppress cell viability in a manner partially dependent upon its C-terminal 26 amino acids that is shown to influence substrate binding. Mutation analysis demonstrates that both the phosphatase activity of Cdc25B as well as its ability to interact with its substrates contribute to the inhibition of cell proliferation. These results demonstrate key differences in the biological activities of Cdc25B and Cdc25C caused by differential substrate affinity and recognition. This also argues that the antiproliferative activity of Cdc25B needs to be overcome for it to act as an oncogene during tumorigenesis.     

Structural and dynamical properties of KTS‐disintegrins: A comparison between Obtustatin and Lebestatin

Obtustatin and Lebestatin are lysine‐threonine‐serine (KTS)‐disintegrins, which are a family of low molecular weight polypeptides present in many viperidae venoms and are potent and specific inhibitors of collagen‐binding integrins. The integrin binding loop, harboring the ²¹KTS²³ motif, and the C‐terminal tail are known to be responsible for the selective binding to the α1β1 integrin. Despite a very high sequence homology (only two mutations are present in Lebestatin relative to Obtustatin, namely R24L and S38L), Lebestatin exhibits a higher inhibitory effect than Obtustatin on cell adhesion and cell migration to collagens I and IV. Here we show, by means of molecular dynamics simulations of the two polypeptides in aqueous solution, that Lebestatin possesses a higher flexibility of the C‐terminal tail and a greater solvent accessibility of the integrin binding loop than Obtustatin. It may be hypothesized that these properties may contribute to the higher binding‐affinity of Lebestatin to its biological partner. © 2012 Wiley Periodicals, Inc.     

Binding and inhibition of Cdc25 phosphatases by vitamin K analogues

A synthetic K vitamin analogue, 2-(2-mercaptothenol)-3-methyl-1,4-naphthoquinone or Cpd 5, was previously found to be a potent inhibitor of cell growth [Nishikawa et al., (1995) J. Biol. Chem. 270, 28304-28310]. The mechanisms of cell growth were hypothesized to include the inactivation of cellular protein tyrosine phosphatases, especially the Cdc25 family [Tamura et al. (2000) Cancer Res. 60, 1317-1325]. In this study, we synthesized PD 49, a new biotin containing Cpd 5 derivative, to search for evidence of direct interaction of these arylating analogues with Cdc25A, Cdc25B, and Cdc25C phosphatases. PD 49 was shown to directly bind to GST-Cdc25A, GST-Cdc25B, their catalytic fragments, and GST-Cdc25C. The binding could be competed with excess glutathione or Cpd 5, and a cysteine-to-serine mutation of the catalytic cysteine abolished binding. This was consistent with an involvement in binding of cysteine in the catalytic domain. This interaction between PD 49 and Cdc25 also occurred in lysates of treated cells. PD 49 also bound to protein phosphatases other than Cdc25. We found that the new analogue also inhibited Hep3B human hepatoma cell growth. This growth inhibition involved ERK1/2 phosphorylation and was inhibited by a MEK antagonist. The results demonstrate a direct interaction and binding between this growth-inhibiting K vitamin derivative with both purified as well as with cellular Cdc25A, Cdc25B, and Cdc25C.     

Membrane dynamics of γ‐secretase with the anterior pharynx‐defective 1B subunit

The four‐subunit protease complex γ‐secretase cleaves many single‐pass transmembrane (TM) substrates, including Notch and β‐amyloid precursor protein to generate amyloid‐β (Aβ), central to Alzheimer's disease. Two of the subunits anterior pharynx‐defective 1 (APH‐1) and presenilin (PS) exist in two homologous forms APH1‐A and APH1‐B, and PS1 and PS2. The consequences of these variations are poorly understood and could affect Aβ production and γ‐secretase medicine. Here, we developed the first complete structural model of the APH‐1B subunit using the published cryo‐electron microscopy (cryo‐EM) structures of APH1‐A (Protein Data Bank: 5FN2, 5A63, and 6IYC). We then performed all‐atom molecular dynamics simulations at 303 K in a realistic bilayer system to understand both APH‐1B alone and in γ‐secretase without and with substrate C83‐bound. We show that APH‐1B adopts a 7TM topology with a water channel topology similar to APH‐1A. We demonstrate direct transport of water through this channel, mainly via Glu84, Arg87, His170, and His196. The apo and holo states closely resemble the experimental cryo‐EM structures with APH‐1A, however with subtle differences: The substrate‐bound APH‐1B γ‐secretase was quite stable, but some TM helices of PS1 and APH‐1B rearranged in the membrane consistent with the disorder seen in the cryo‐EM data. This produces different accessibility of water molecules for the catalytic aspartates of PS1, critical for Aβ production. In particular, we find that the typical distance between the catalytic aspartates of PS1 and the C83 cleavage sites are shorter in APH‐1B, that is, it represents a more closed state, due to interactions with the C‐terminal fragment of PS1. Our structural‐dynamic model of APH‐1B alone and in γ‐secretase suggests generally similar topology but some notable differences in water accessibility which may be relevant to the protein's existence in two forms and their specific function and location.     

Is there nascent structure in the intrinsically disordered region of troponin I?

In striated muscle, the binding of calcium to troponin C (TnC) results in the removal of the C‐terminal region of the inhibitory protein troponin I (TnI) from actin. While structural studies of the muscle system have been successful in determining the overall organization of most of the components involved in force generation at the atomic level, the structure and dynamics of the C‐terminal region of TnI remains controversial. This domain of TnI is highly flexible, and it has been proposed that this intrinsically disordered region (IDR) regulates contraction via a “fly‐casting” mechanism. Different structures have been presented for this region using different methodologies: a single α‐helix, a “mobile domain” containing a small β‐sheet, an unstructured region, and a two helix segment. To investigate whether this IDR has in fact any nascent structure, we have constructed a skeletal TnC‐TnI chimera that contains the N‐domain of TnC (1–90), a short linker (GGAGG), and the C‐terminal region of TnI (97–182) and have acquired ¹⁵N NMR relaxation data for this chimera. We compare the experimental relaxation parameters with those calculated from molecular dynamic simulations using four models based upon the structural studies. Our experimental results suggest that the C‐terminal region of TnI does not contain any defined secondary structure, supporting the “fly‐casting” mechanism. We interpret the presence of a “plateau” in the ¹⁵N NMR relaxation data as being an intrinsic property of IDRs. We also identified a more rigid adjacent region of TnI that has implications for muscle performance under ischemic conditions. Proteins 2011. © 2011 Wiley‐Liss, Inc.     

The flexible C‐terminal arm of the Lassa arenavirus Z‐protein mediates interactions with multiple binding partners

The arenavirus genome encodes for a Z‐protein, which contains a RING domain that coordinates two zinc ions, and has been identified as having several functional roles at various stages of the virus life cycle. Z‐protein binds to multiple host proteins and has been directly implicated in the promotion of viral budding, repression of mRNA translation, and apoptosis of infected cells. Using homology models of the Z‐protein from Lassa strain arenavirus, replica exchange molecular dynamics (MD) was used to refine the structures, which were then subsequently clustered. Population‐weighted ensembles of low‐energy cluster representatives were predicted based upon optimal agreement of the chemical shifts computed with the SPARTA program with the experimental NMR chemical shifts. A member of the refined ensemble was indentified to be a potential binder of budding factor Tsg101 based on its correspondence to the structure of the HIV‐1 Gag late domain when bound to Tsg101. Members of these ensembles were docked against the crystal structure of human eIF4E translation initiation factor. Two plausible binding modes emerged based upon their agreement with experimental observation, favorable interaction energies and stability during MD trajectories. Mutations to Z are proposed that would either inhibit both binding mechanisms or selectively inhibit only one mode. The C‐terminal domain conformation of the most populated member of the representative ensemble shielded protein‐binding recognition motifs for Tsg101 and eIF4E and represents the most populated state free in solution. We propose that C‐terminal flexibility is key for mediating the different functional states of the Z‐protein. Proteins 2010. © 2010 Wiley‐Liss, Inc.     

Temperature‐accelerated molecular dynamics gives insights into globular conformations sampled in the free state of the AC catalytic domain

The catalytic domain of the adenyl cyclase (AC) toxin from Bordetella pertussis is activated by interaction with calmodulin (CaM), resulting in cAMP overproduction in the infected cell. In the X‐ray crystallographic structure of the complex between AC and the C terminal lobe of CaM, the toxin displays a markedly elongated shape. As for the structure of the isolated protein, experimental results support the hypothesis that more globular conformations are sampled, but information at atomic resolution is still lacking. Here, we use temperature‐accelerated molecular dynamics (TAMD) simulations to generate putative all‐atom models of globular conformations sampled by CaM‐free AC. As collective variables, we use centers of mass coordinates of groups of residues selected from the analysis of standard molecular dynamics (MD) simulations. Results show that TAMD allows extended conformational sampling and generates AC conformations that are more globular than in the complexed state. These structures are then refined via energy minimization and further unrestrained MD simulations to optimize inter‐domain packing interactions, thus resulting in the identification of a set of hydrogen bonds present in the globular conformations. Proteins 2014; 82:2483–2496. © 2014 Wiley Periodicals, Inc.     

Insights into the interaction of high potency inhibitor IRC‐083864 with phosphatase CDC25

CDC25 phosphatases play a crucial role in cell cycle regulation. They have been found to be over‐expressed in various human tumours and to be valuable targets for cancer treatment. Here, we report the first model of binding of the most potent CDC25 inhibitor to date, the bis‐quinone IRC‐083864, into CDC25B obtained by combining molecular modeling and NMR studies. Our study provides new insights into key interactions of the catalytic site inhibitor and CDC25B in the absence of any available experimental structure of CDC25 with a bound catalytic site inhibitor. The docking model reveals that IRC‐083864 occupies both the active site and the inhibitor binding pocket of the CDC25B catalytic domain. NMR saturation transfer difference and WaterLOGSY data indicate the binding zones of the inhibitor and support the docking model. Probing interactions of analogues of the two quinone units of IRC‐083864 with CDC25B demonstrate that IRC‐083864 competes with each monomer. Proteins 2017; 85:593–601. © 2016 Wiley Periodicals, Inc.     

Structure of dystrophia myotonica protein kinase

Dystrophia myotonica protein kinase (DMPK) is a serine/threonine kinase composed of a kinase domain and a coiled‐coil domain involved in the multimerization. The crystal structure of the kinase domain of DMPK bound to the inhibitor bisindolylmaleimide VIII (BIM‐8) revealed a dimeric enzyme associated by a conserved dimerization domain. The affinity of dimerisation suggested that the kinase domain alone is insufficient for dimerisation in vivo and that the coiled‐coil domains are required for stable dimer formation. The kinase domain is in an active conformation, with a fully‐ordered and correctly positioned αC helix, and catalytic residues in a conformation competent for catalysis. The conserved hydrophobic motif at the C‐terminal extension of the kinase domain is bound to the N‐terminal lobe of the kinase domain, despite being unphosphorylated. Differences in the arrangement of the C‐terminal extension compared to the closely related Rho‐associated kinases include an altered PXXP motif, a different conformation and binding arrangement for the turn motif, and a different location for the conserved NFD motif. The BIM‐8 inhibitor occupies the ATP site and has similar binding mode as observed in PDK1.