Two TPX2-Dependent Switches Control the Activity of Aurora A

Aurora A is an important oncogenic kinase for mitotic spindle assembly and a potentially attractive target for human cancers. Its activation could be regulated by ATP cycle and its activator TPX2. To understand the activation mechanism of Aurora A, a series of 20 ns molecular dynamics (MD) simulations were performed on both the wild-type kinase and its mutants. Analyzing the three dynamic trajectories (Aurora A-ATP, Aurora A-ADP, and Aurora A-ADP-TPX2) at the residue level, for the first time we find two TPX2-dependent switches, i.e., switch-1 (Lys-143) and switch-2 (Arg-180), which are tightly associated with Aurora A activation. In the absence of TPX2, Lys-143 exhibits a “closed” state, and becomes hydrogen-bonded to ADP. Once TPX2 binding occurs, switch-1 is forced to “open” the binding site, thus pulling ADP away from Aurora A. Without facilitation of TPX2, switch-2 exits in an “open” conformation which accompanies the outward-flipping movement of P·Thr288 (in an inactive conformation), leading to the crucial phosphothreonine exposed and accessible for deactivation. However, with the binding of TPX2, switch-2 is forced to undergo a “closed” movement, thus capturing P·Thr288 into a buried position and locking its active conformation. Analysis of two Aurora A (K143A and R180A) mutants for the two switches further verifies their functionality and reliability in controlling Aurora activity. Our systems therefore suggest two switches determining Aurora A activation, which are important for the development of aurora kinase inhibitors

Anti-Fas mAb-induced apoptosis and cytolysis of airway tissue eosinophils aggravates rather than resolves established inflammation

BACKGROUND: Fas receptor-mediated eosinophil apoptosis is currently forwarded as a mechanism resolving asthma-like inflammation. This view is based on observations in vitro and in airway lumen with unknown translatability to airway tissues in vivo. In fact, apoptotic eosinophils have not been detected in human diseased airway tissues whereas cytolytic eosinophils abound and constitute a major mode of degranulation of these cells. Also, Fas receptor stimulation may bypass the apoptotic pathway and directly evoke cytolysis of non-apoptotic cells. We thus hypothesized that effects of anti-Fas mAb in vivo may include both apoptosis and cytolysis of eosinophils and, hence, that established eosinophilic inflammation may not resolve by this treatment. METHODS: Weeklong daily allergen challenges of sensitized mice were followed by airway administration of anti-Fas mAb. BAL was performed and airway-pulmonary tissues were examined using light and electron microscopy. Lung tissue analysis for CC-chemokines, apoptosis, mucus production and plasma exudation (fibrinogen) were performed. RESULTS: Anti-Fas mAb evoked apoptosis of 28% and cytolysis of 4% of eosinophils present in allergen-challenged airway tissues. Furthermore, a majority of the apoptotic eosinophils remained unengulfed and eventually exhibited secondary necrosis. A striking histopathology far beyond the allergic inflammation developed and included degranulated eosinophils, neutrophilia, epithelial derangement, plasma exudation, mucus-plasma plugs, and inducement of 6 CC-chemokines. In animals without eosinophilia anti-Fas evoked no inflammatory response. CONCLUSION: An efficient inducer of eosinophil apoptosis in airway tissues in vivo, anti-Fas mAb evoked unprecedented asthma-like inflammation in mouse allergic airways. This outcome may partly reflect the ability of anti-Fas to evoke direct cytolysis of non-apoptotic eosinophils in airway tissues. Additionally, since most apoptotic tissue eosinophils progressed into the pro-inflammatory cellular fate of secondary necrosis this may also explain the aggravated inflammation. Our data indicate that Fas receptor mediated eosinophil apoptosis in airway tissues in vivo may cause severe disease exacerbation due to direct cytolysis and secondary necrosis of eosinophils

Two TPX2-Dependent Switches Control the Activity of Aurora A

Aurora A is an important oncogenic kinase for mitotic spindle assembly and a potentially attractive target for human cancers. Its activation could be regulated by ATP cycle and its activator TPX2. To understand the activation mechanism of Aurora A, a series of 20 ns molecular dynamics (MD) simulations were performed on both the wild-type kinase and its mutants. Analyzing the three dynamic trajectories (Aurora A-ATP, Aurora A-ADP, and Aurora A-ADP-TPX2) at the residue level, for the first time we find two TPX2-dependent switches, i.e., switch-1 (Lys-143) and switch-2 (Arg-180), which are tightly associated with Aurora A activation. In the absence of TPX2, Lys-143 exhibits a “closed” state, and becomes hydrogen-bonded to ADP. Once TPX2 binding occurs, switch-1 is forced to “open” the binding site, thus pulling ADP away from Aurora A. Without facilitation of TPX2, switch-2 exits in an “open” conformation which accompanies the outward-flipping movement of P·Thr288 (in an inactive conformation), leading to the crucial phosphothreonine exposed and accessible for deactivation. However, with the binding of TPX2, switch-2 is forced to undergo a “closed” movement, thus capturing P·Thr288 into a buried position and locking its active conformation. Analysis of two Aurora A (K143A and R180A) mutants for the two switches further verifies their functionality and reliability in controlling Aurora activity. Our systems therefore suggest two switches determining Aurora A activation, which are important for the development of aurora kinase inhibitors

A secondary RET mutation in the activation loop conferring resistance to vandetanib

Resistance to vandetanib, a type I RET kinase inhibitor, developed in a patient with metastatic lung adenocarcinoma harboring a CCDC6-RET fusion that initially exhibited a response to treatment. The resistant tumor acquired a secondary mutation resulting in a serine-to-phenylalanine substitution at codon 904 in the activation loop of the RET kinase domain. The S904F mutation confers resistance to vandetanib by increasing the ATP affinity and autophosphorylation activity of RET kinase. A reduced interaction with the drug is also observed in vitro for the S904F mutant by thermal shift assay. A crystal structure of the S904F mutant reveals a small hydrophobic core around F904 likely to enhance basal kinase activity by stabilizing an active conformer. Our findings indicate that missense mutations in the activation loop of the kinase domain are able to increase kinase activity and confer drug resistance through allosteric effects

ATP and its N6-substituted analogues: parameterization, molecular dynamics simulation and conformational analysis

In this work we used a combination of classical molecular dynamics and simulated annealing techniques to shed more light on the conformational flexibility of 12 adenosine triphosphate (ATP) analogues in a water environment. We present simulations in AMBER force field for ATP and 12 published analogues [Shah et al. (1997) Proc Natl Acad Sci USA 94: 3565–3570]. The calculations were carried out using the generalized Born (GB) solvation model in the presence of the cation Mg2+. The ion was placed at a close distance (2 Å) from the charged oxygen atoms of the beta and gamma phosphate groups of the −3 negatively charged ATP analogue molecules. Analysis of the results revealed the distribution of inter-proton distances H8–H1′ and H8–H2′ versus the torsion angle ψ (C4–N9-C1′–O4′) for all conformations of ATP analogues. There are two gaps in the distribution of torsion angle ψ values: the first is between −30 and 30 degrees and is described by cis-conformation; and the second is between 90 and 175 degrees, which mostly covers a region of anti conformation. Our results compare favorably with results obtained in experimental assays [Jiang and Mao (2002) Polyhedron 21:435–438]

Ras Conformational Switching: Simulating Nucleotide-Dependent Conformational Transitions with Accelerated Molecular Dynamics

Ras mediates signaling pathways controlling cell proliferation and development by cycling between GTP- and GDP-bound active and inactive conformational states. Understanding the complete reaction path of this conformational change and its intermediary structures is critical to understanding Ras signaling. We characterize nucleotide-dependent conformational transition using multiple-barrier-crossing accelerated molecular dynamics (aMD) simulations. These transitions, achieved for the first time for wild-type Ras, are impossible to observe with classical molecular dynamics (cMD) simulations due to the large energetic barrier between end states. Mapping the reaction path onto a conformer plot describing the distribution of the crystallographic structures enabled identification of highly populated intermediate structures. These structures have unique switch orientations (residues 25–40 and 57–75) intermediate between GTP and GDP states, or distinct loop3 (46–49), loop7 (105–110), and α5 C-terminus (159–166) conformations distal from the nucleotide-binding site. In addition, these barrier-crossing trajectories predict novel nucleotide-dependent correlated motions, including correlations of α2 (residues 66–74) with α3-loop7 (93–110), loop2 (26–37) with loop10 (145–151), and loop3 (46–49) with α5 (152–167). The interconversion between newly identified Ras conformations revealed by this study advances our mechanistic understanding of Ras function. In addition, the pattern of correlated motions provides new evidence for a dynamic linkage between the nucleotide-binding site and the membrane interacting C-terminus critical for the signaling function of Ras. Furthermore, normal mode analysis indicates that the dominant collective motion that occurs during nucleotide-dependent conformational exchange, and captured in aMD (but absent in cMD) simulations, is a low-frequency motion intrinsic to the structure

Nucleotide Binding Switches the Information Flow in Ras GTPases

The Ras superfamily comprises many guanine nucleotide-binding proteins (G proteins) that are essential to intracellular signal transduction. The guanine nucleotide-dependent intrinsic flexibility patterns of five G proteins were investigated in atomic detail through Molecular Dynamics simulations of the GDP- and GTP-bound states (SGDP and SGTP, respectively). For all the considered systems, the intrinsic flexibility of SGDP was higher than that of SGTP, suggesting that Guanine Exchange Factor (GEF) recognition and nucleotide switch require higher amplitude motions than effector recognition or GTP hydrolysis. Functional mode, dynamic domain, and interaction energy correlation analyses highlighted significant differences in the dynamics of small G proteins and Gα proteins, especially in the inactive state. Indeed, SGDP of Gαt, is characterized by a more extensive energy coupling between nucleotide binding site and distal regions involved in GEF recognition compared to small G proteins, which attenuates in the active state. Moreover, mechanically distinct domains implicated in nucleotide switch could be detected in the presence of GDP but not in the presence of GTP. Finally, in small G proteins, functional modes are more detectable in the inactive state than in the active one and involve changes in solvent exposure of two highly conserved amino acids in switches I and II involved in GEF recognition. The average solvent exposure of these amino acids correlates in turn with the rate of GDP release, suggesting for them either direct or indirect roles in the process of nucleotide switch. Collectively, nucleotide binding changes the information flow through the conserved Ras-like domain, where GDP enhances the flexibility of mechanically distinct portions involved in nucleotide switch, and favors long distance allosteric communication (in Gα proteins), compared to GTP

The Distinct Conformational Dynamics of K-Ras and H-Ras A59G

Ras proteins regulate signaling cascades crucial for cell proliferation and differentiation by switching between GTP- and GDP-bound conformations. Distinct Ras isoforms have unique physiological functions with individual isoforms associated with different cancers and developmental diseases. Given the small structural differences among isoforms and mutants, it is currently unclear how these functional differences and aberrant properties arise. Here we investigate whether the subtle differences among isoforms and mutants are associated with detectable dynamical differences. Extensive molecular dynamics simulations reveal that wild-type K-Ras and mutant H-Ras A59G are intrinsically more dynamic than wild-type H-Ras. The crucial switch 1 and switch 2 regions along with loop 3, helix 3, and loop 7 contribute to this enhanced flexibility. Removing the gamma-phosphate of the bound GTP from the structure of A59G led to a spontaneous GTP-to-GDP conformational transition in a 20-ns unbiased simulation. The switch 1 and 2 regions exhibit enhanced flexibility and correlated motion when compared to non-transitioning wild-type H-Ras over a similar timeframe. Correlated motions between loop 3 and helix 5 of wild-type H-Ras are absent in the mutant A59G reflecting the enhanced dynamics of the loop 3 region. Taken together with earlier findings, these results suggest the existence of a lower energetic barrier between GTP and GDP states of the mutant. Molecular dynamics simulations combined with principal component analysis of available Ras crystallographic structures can be used to discriminate ligand- and sequence-based dynamic perturbations with potential functional implications. Furthermore, the identification of specific conformations associated with distinct Ras isoforms and mutants provides useful information for efforts that attempt to selectively interfere with the aberrant functions of these species

Profiling the immune landscape in mucinous ovarian carcinoma

Objective: Mucinous ovarian carcinoma (MOC) is a rare histotype of ovarian cancer, with low response rates to standard chemotherapy, and very poor survival for patients diagnosed at advanced stage. There is a limited understanding of the MOC immune landscape, and consequently whether immune checkpoint inhibitors could be considered for a subset of patients. Methods: We performed multicolor immunohistochemistry (IHC) and immunofluorescence (IF) on tissue microarrays in a cohort of 126 MOC patients. Cell densities were calculated in the epithelial and stromal components for tumor-associated macrophages (CD68+/PD-L1+, CD68+/PD-L1-), T cells (CD3+/CD8-, CD3+/CD8+), putative T-regulatory cells (Tregs, FOXP3+), B cells (CD20+/CD79A+), plasma cells (CD20-/CD79a+), and PD-L1+ and PD-1+ cells, and compared these values with clinical factors. Univariate and multivariable Cox Proportional Hazards assessed overall survival. Unsupervised k-means clustering identified patient subsets with common patterns of immune cell infiltration. Results: Mean densities of PD1+ cells, PD-L1- macrophages, CD4+ and CD8+ T cells, and FOXP3+ Tregs were higher in the stroma compared to the epithelium. Tumors from advanced (Stage III/IV) MOC had greater epithelial infiltration of PD-L1- macrophages, and fewer PD-L1+ macrophages compared with Stage I/II cancers (p = 0.004 and p = 0.014 respectively). Patients with high epithelial density of FOXP3+ cells, CD8+/FOXP3+ cells, or PD-L1- macrophages, had poorer survival, and high epithelial CD79a + plasma cells conferred better survival, all upon univariate analysis only. Clustering showed that most MOC (86%) had an immune depleted (cold) phenotype, with only a small proportion (11/76,14%) considered immune inflamed (hot) based on T cell and PD-L1 infiltrates. Conclusion: In summary, MOCs are mostly immunogenically ‘cold’, suggesting they may have limited response to current immunotherapies

Elucidation of the ATP7B N-Domain Mg2+-ATP Coordination Site and Its Allosteric Regulation

The diagnostic of orphan genetic disease is often a puzzling task as less attention is paid to the elucidation of the pathophysiology of these rare disorders at the molecular level. We present here a multidisciplinary approach using molecular modeling tools and surface plasmonic resonance to study the function of the ATP7B protein, which is impaired in the Wilson disease. Experimentally validated in silico models allow the elucidation in the Nucleotide binding domain (N-domain) of the Mg2+-ATP coordination site and answer to the controversial role of the Mg2+ ion in the nucleotide binding process. The analysis of protein motions revealed a substantial effect on a long flexible loop branched to the N-domain protein core. We demonstrated the capacity of the loop to disrupt the interaction between Mg2+-ATP complex and the N-domain and propose a role for this loop in the allosteric regulation of the nucleotide binding process