Conformational Transition of Key Structural Features Involved in

mutations caused missense substitution of single amino acid in the kinase domain could promote constitutive ... provider but also as an allosteric mod...
0 downloads 3 Views 1MB Size
Subscriber access provided by University of Chicago Library

Conformational Transition of Key Structural Features Involved in Activation of ALK Induced by Two Neuroblastoma Mutations and ATP Binding: Insight from Accelerated Molecular Dynamics Simulations Muyang He, Weikang Li, Qingchuan Zheng, and Hong-Xing Zhang ACS Chem. Neurosci., Just Accepted Manuscript • DOI: 10.1021/acschemneuro.8b00105 • Publication Date (Web): 11 Apr 2018 Downloaded from http://pubs.acs.org on April 11, 2018

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 38 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Chemical Neuroscience

Conformational Transition of Key Structural Features Involved

in

Activation

of

ALK

Induced

by

two

Neuroblastoma Mutations and ATP Binding: Insight from Accelerated Molecular Dynamics Simulations

Mu-Yang He,2 Wei-Kang Li,2Qing-Chuan Zheng1,2*Hong-Xing Zhang2*

1. Key Laboratory for Molecular Enzymology and Engineering of the Ministry of Education, College of Life Science, Jilin University, Changchun 130012, People’s Republic of China

2. Laboratory of Theoretical and Computational Chemistry, Institute of Theoretical Chemistry, International Joint Research Laboratory of Nano-Micro Architecture Chemistry, Jilin University, Changchun 130023, People’s Republic of China

Correspondence to:

Qing-Chuan Zheng, Institute of Theoretical Chemistry, Jilin University, Changchun, 130023, People’s Republic of China.

Fax: (+86) 431-8849-8966, E-mail address: [email protected]

Hong-Xing Zhang, Institute of Theoretical Chemistry, Jilin University, Changchun, 130023, People’s Republic of China.

Fax:

(+86) 431-8849-8962, E-mail address: [email protected]

ACS Paragon Plus Environment

ACS Chemical Neuroscience 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ABSTRACT Deregulated kinase activity of Anaplastic Lymphoma Kinase (ALK) has been observed to be implicated in developing of tumor progress. The activation mechanism of ALK is proposed to be similar with other receptor tyrosine kinase (RTK), but the distinct static X-ray crystal conformation of ALK suggests its unique conformational transition. Herein, we have illustrated the dynamic conformational property of wild type ALK as well as the kinase activation equilibrium variation induced by two neuroblastoma mutations (R1275Q and Y1278S) and ATP binding by performed enhanced sampling accelerated Molecular Dynamics (aMD) simulations. The results suggested that the wild type ALK is mostly favored in the inactive state, whereas the mutations and ATP binding promoted clear shift toward the active-like conformation. The R1275Q mutant stabilize the active conformation by rigidified the αC-in conformation. The Y1278S mutant promotes the activation at the expense of a π-stacking hydrophobic cluster which plays critical role in the stabilization of inactive conformation of native ALK. ATP produces a more compact active site and thereby facilitates the activation of ALK. Taken together, these findings not only elucidate the diverse conformations in different ALKs, but also can shed light on new strategies for protein engineering and structural based drug design for ALK.

KEY WORDS: ALK, Kinase Activation, Neuroblastoma mutation, ATP, aMD

ACS Paragon Plus Environment

Page 2 of 38

Page 3 of 38 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Chemical Neuroscience

INTRODUCTION Anaplastic lymphoma kinase (ALK) is a member of the receptor tyrosine kinase (RTK) superfamily, which is postulated to play a pivotal role in normal development and function of the nervous system in embryos.

1, 2

Elevated activities of ALK are

associated with the development and progression of multiple human malignancies.2 Structurally, ALK is composed of an extracellular receptor, a transmembrane region, and an intracellular tyrosine kinase domain.2 The kinase domain of ALK adopts the canonical kinase-fold, which is built of a smaller N-terminal lobe and an α-helical C-terminal lobe.3 The ATP binding-site is at the interlobe cleft referred as the active site of ALK. A common activation mechanism is believed to exist for members of this superfamily which mainly controlled by conformational changes in two conserved structural motifs: the αC-helix and the activation loop (A-loop) (Figure 1). 4-8

On one hand, the active state is characterized by an “αC-in” conformation which is

proper for the catalytic salt bridge formation between conserved lysine and glutamic acid in N-terminal lobe. On the other hand, the A-loop assumes an extended conformation to make adequate space for peptide substrate binding. The X-ray crystal structures of dormant ALK reveal a unique partially inactive tyrosine kinase conformation which lack the complete repertoire of negative regulatory structural elements that exist in the inactive RTKs.3, 9, 10 In dormant ALK, αC-helix adopts an “in” conformation, and A-loop exhibits a semi-closed conformation, which is in an intermediate position between the inactive and active RTKs. It may highlight the fact that the activation of ALK is rather a unique mechanism which makes it both

ACS Paragon Plus Environment

ACS Chemical Neuroscience 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 4 of 38

interesting yet elusive. Substantial studies have demonstrated that mutations and ligand binding could play significant role as allosteric modulators in the activation of kinase.11-13 In ALK, mutations caused missense substitution of single amino acid in the kinase domain could promote constitutive, ligand independent tyrosine kinase activation.2, 14, 15 Many somatic and germline mutations of ALK arouse clinical notice on account of their activating or drug-resistance effects in many diseases, including neuroblastoma,16 non small cell lung cancer (NSCLC),17 and inflammatory myofibroblastic tumor (IMT).18 The hot spot mutation R1275Q is the single most common mutation

(approximately

43% of all mutations) detected in neuroblastoma,19 a pediatric disease that arises in the sympathetic nervous system, representing the most frequently diagnosed malignant tumor in the first year of life.16,

20

Point mutation is also observed in

Y1278S from the patients of neuroblastoma, although less frequently.21 Y1278 is the first tyrosine undergoes autocatalytic phosphorylation in A-loop, and its substitution also reported to promotes constitutive kinase activity of ALK.21,

22

To date, the

biological effects of these activating drivers are clear,21, 23, 24 whereas the molecular origins and modulating modes in the activation process of ALK have not been rationalized. In addition, recent studies state that ATP can serve not only as an energy provider but also as an allosteric modulator that regulates the process of kinase activation.25-28 In the ALK involved cellular signal transduction, ATP serves as substrate and provides phosphate group in a number of signal transduction cascades.23, 29-31

In contrast with the well understood role in phosphorylation, the allosteric

ACS Paragon Plus Environment

Page 5 of 38 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Chemical Neuroscience

behaviors of ATP in ALK are rarely illuminated. As ALK is highly expressed only in the nervous system during embryogenesis, it is difficult to implement experiments that pursue its mechanism of activation.2, 32, 33 Molecular dynamics (MD) simulations are increasingly being applied to study the structural biology of biomacromolecules.34-37 The computational study of ALK mainly focus on the drug resistance mechanisms and specific inhibitor design. Abundance of computational investigations have in-sighted both the first and second inhibitors drug resistance mechanisms at atomic-level induced by ALK mutations.38-41 Hou et.al have rationally designed a novel kind of inhibitors with excellent antiproliferation activity against ALK-positive cancer cells, 42-44 which have provided a new strategy in rational design of potent and selective ALK inhibitors to conquer drug resistance. Despite the plenty research about drug recognition of ALK, the activation properties of ALK have not been elaborated.

Thanks to the development of accelerated molecular dynamics

(aMD) simulation, the states transition of protein can be clarified at atomic level.45-48 Thus, we performed aMD simulations on the apo wild type, R1275Q mutant, Y1278S mutant, and ATP-bound forms of ALK to explore the potential activation mechanisms in this study. Specifically, the introduction of free energy landscape method figures out the dynamical consequences of two neuroblastoma mutants and ATP binding on the conformational alterations of the key structural features associated with the activation of ALK. These insights on the ALK structural features and conformation dynamic can provide clues to understand the detailed activation mechanisms of ALK on atomistic level.

ACS Paragon Plus Environment

ACS Chemical Neuroscience 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

RESULTS AND DISCUSSION Conformational Free Energy Landscape of Apo-WT ALK. To study the conformational states and structural dynamics of ALKs, 50 ns conventional MD (cMD) followed by 500 ns aMD were performed for each system. As shown in Figure S1, it is clear that the RMSD values are convergent for both the cMD and aMD simulations. The protein shows no significant deviation from the initial structure in the preliminary 50 ns cMD simulations. In the subsequent aMD simulations, diverse conformations were sampled, thus provided the opportunity for the capture of activation process of ALK. Free energy profiles of the ALKs are explored utilizing aMD simulations in terms of the formation of catalytic salt bridge, and the A-loop extension, both of which have been noted previously to be critical for ALK activation. The free energy profiles were projected along two collective variables (CV1 and CV2) for the four systems. CV1 measures the deviation of A-loop from the

X-ray crystal structure (PDB code 3L9P), which is able to characterize

the open and closed conformation of A-loop. The deviation of A-loop is measured by the average deviation of the residues in the distal of A-loop, including Y1283 to M1290. These residues approaching to the active site would block the peptide binding sites, a conformation that regarded as inactive kinase state. On the contrary, the distal residues move outward would shape an open A-loop conformation which is favorable for the peptide substrate binding. The negative value denotes deviation of A-loop approaching the active site, and the positive value is on behalf of the outward

ACS Paragon Plus Environment

Page 6 of 38

Page 7 of 38 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Chemical Neuroscience

extension of A-loop. CV2 is defined as the distance of the salt-bridge-forming couple: K1150 and E1167, which is able to characterize the displacement of the functionally important component αC-helix of the N lobe in its transition between the “αC-in” and “αC-out” conformation. The value of CV2 less than 4 Å represents the formation of the salt bridge.49 The two dimensional free energy profiles of the Apo-WT ALK system is shown in Figure 2A. For the deviation of A-loop, the value sampled in this system mainly distributes between - 3Å to 5Å (with free energy less than 1 kcal/mol). Notably, no A-loop open conformation is sampled in Apo-WT ALK system. As shown in Figure S2, the conformation of A-loop with 5Å outward deviation is not full extended and is not able to provide sufficient space for the peptide binding. The 3Å inward deviation of A-loop obstructs the not only the peptide binding site but also the ATP binding site, which is very similar to the close conformation of A-loop observed in other inactive RTKs.50-52 These results demonstrate that the wild-type kinase domain on its own is indeed in an autoinhibited state. For the conformation of αC-helix, in the majority of simulation the distance drops below the value required for salt bridge to form between K1150 and E1167, representing an αC-in conformation. Remarkably, in the second low energy basin, both the distance below and more than 4 Å are sampled. The fairly little energy penalty suggests the toilless transition between αC-in and αC-out conformation in the Apo-WT ALK system. Apo-WT ALK Kinase Domain is Intrinsically Autoinhibited. The deepest free-energy minimum of Apo-WT corresponds to partial inactive

ACS Paragon Plus Environment

ACS Chemical Neuroscience 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

conformation with structural features similar to the X-ray crystal ALK structure (PDB code 3L9P),3 where the salt bridge between K1150 and E1167 is well maintained and the A-loop adopts a partial extended solvent exposed conformation (Figure 2A and 2B). In spite of the similarity, the aMD simulations provide more atomic details about sustain of the partial inactive conformation. As shown in Figure 3A, polar residues adorning the polar face of αC-helix and αAL (the first helix of A-loop) conduce to stabilize their relative position in the partial inactive conformation of ALK. The Aspartic acid in αC-helix (D1163) forms two hydrogen bonds with two positively charged residues in αAL (R1275 and R1279), which serve as the key factors that stabilize the αC-in conformation in this partial inactive state. For the conformation of A-loop, a series of hydrogen bonds set up the semi-closed architecture (listed in Table S1). As shown in Figure 3A, the hydrogen bond formed between Q1159 of αC-helix and the third A-loop tyrosine Y1283 constructs a sharp U-turn at the end of αAL that prevents A-loop shapes neither the full extension nor the closed conformation, which is regarded as the principle factor that stabilized the semi-closed A-loop conformation. The hydrogen bond formed between the main chain of D1276 and the first side chain nitrogen atom of R1284 further steady this U-turn. In addition, electrostatic and hydrophobic interactions among the distal part and the tail end of A-loop lend further stability to the semi-closed A-loop conformation. In the distal part of A-loop, a hydrogen bond between C1288 and L1291 is observed in this partial extended conformation, which may further restrain the adaptability of A-loop. The tail end residue of A-loop V1293 is fastened adjacent to αG-helix by a hydrophobic cluster

ACS Paragon Plus Environment

Page 8 of 38

Page 9 of 38 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Chemical Neuroscience

consisted of residue M1296, F1301, and L1339. As shown in Figure 3A, both of the hydrogen bond and the hydrophobic cluster mentioned above can contribute to the stability of a hydrophobic stem consisted of the residue M1273, L1291, and P1292. This hydrophobic stem has been observed in the crystal structure of ALK and considered to play a critical role in holding the semi-closed A-loop conformation. Besides the most populated partial inactive conformation, a conformation that matches the common criteria of inactive RTK is sampled. The ensembles in the second low energy basin are clustered into four clusters. One of which represents the value of CV1 is no more than -3 Å, and CV2 is more than 4 Å. The center structure of this cluster is chosen as the representative structure of inactive-like Apo-WT. In this inactive-like conformation, the αC-helix is rotated outward and displaced from the active site, and the A-loop adopts a closed conformation that obstruct both the ATP and the peptide binding sites (Figure 1D and Figure 3B). As expected, the interactions among these two activation relative motifs altered from the partial inactive conformation. In the inactive-like conformation, the charged side chain of D1163 in αC-helix makes a hydrogen bond with Y1278 who located in the end of αAL instead of with the two middle situated arginine. Another negative charged residue in αC-helix E1167 interacts with R1275 and R1279 leading to the deflection of the side chain of E1167. Regarding to the A-loop conformation, the hydrogen bond network are distinct with that emerged in the semi-closed state are listed in Table S2. One of the significant alteration is the disruption of hydrogen bond between Y1283 and Q1159. In the A-loop closed state, the side chain of Y1283 is hauled approaching the

ACS Paragon Plus Environment

ACS Chemical Neuroscience 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

active site through the hydrogen bond chains consisted of E1167, R1279 and Y1283. Moreover, a second short helix after the αAL is observed in the forepart of A-loop. It is noteworthy that this second helical segment is commonly presented in the inactive conformations of other RTKs, such as EGFR,50 Src kinase,52 and is considered to be important in the formation of close A-loop conformation. In Apo-WT ALK, this helix is established by a network of hydrogen bonds between the main chains of residues S1281, Y1282, Y1283, R1284, and K1285. Indeed, the hydrogen bond network in this helix further reins the outward motion of residues in the middle part of A-loop. Notably, the residue M1273, L1291 and P1292 no longer construct the hydrophobic stem in this A-loop closed conformation. In contrast, M1290 packing against L1291 and P1292 forming a new hydrophobic stem that is close to the P-loop and αD-helix thus occupies both the ATP and peptide substrate binding pockets (Figure 2B). As shown in Figure 3C, L1291 and P1292 are held at the proline enrichment region by hydrophobic effect provided by P1329 and P1331, thus further consolidate the new hydrophobic stem. In addition, the terminal A-loop residue V1293 no longer stuck around αG-helix in the inactive state. As shown in Figure 3D, V1293 participates in a new hydrophobic cluster and sustains the closed conformation of A-loop. These hydrophobic interactions restrain the conformational dynamics of the A-loop and thus retain the active site in a “full-closed” state. In summary, the Apo-WT ALK kinase domain is intrinsically autoinhibited and the transition to its active state is hard to achieve. First of all, the negligible energy penalty for disruption of the K-E salt bridge indicates the unsteadiness of this catalytic

ACS Paragon Plus Environment

Page 10 of 38

Page 11 of 38 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Chemical Neuroscience

bridge in the Apo-WT system. In addition, the conformation of A-loop in the Apo-WT system can transform between the initial semi-closed and the closed conformations, both of which are incompatible with substrate binding. R1275Q Maintains the Stability of K-E Salt Bridge. In contrast to the Apo-WT, the mutants populate to some extent the active conformations, and generate new minima in different locations of the FEL. The conformational FEL of the R1275Q ALK system is shown in Figure 2A. The most stable state of R1275Q corresponds to an ensemble of structures which are in similar with the partial inactive conformation, where the A-loop is folded and assumes a semi-closed conformation and the αC-helix closely interacts with β3-strand (Figure 2C). For the distribution of CV1, the value sampled in this system mainly distributes between - 4Å to 3Å (with free energy less than 1 kcal/mol).The range of A-loop deviation in R1275Q represents the conformational variation between closed and semi-closed status, which is in the same way of the Apo-WT system. For the distribution of CV2, the main population is in the range less than 4 Å. In another word, the conformation adopted by αC-helix in R1275Q appears to result more from an energetic advantage for the “in” conformation rather than an energetic disadvantage for an alternate “out” one. Thus, it is reasonable to infer that R1275Q mutation is expected to destabilize the inactive conformation by maintain the formation of the K-E salt bridge thus keep the αC-in pre-active conformation. To unravel the molecular origins of this phenomenon, we analysis the structural variations between the native and R1275Q mutated ALKs. In the Apo-WT

ACS Paragon Plus Environment

ACS Chemical Neuroscience 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

system, the side chain orientation of glutamic acid 1167 who engages in the catalytic salt bridge formation exhibits high flexibility. In quite a few frames, E1167 turns outward and makes electrostatic interactions with the negative charged arginine R1275 in αAL (Figure 4A). The shift of E1167 directivity leads to the departure from its salt bridge formation partner K1150, which impacts the formation of the important salt bridge and thus drives αC-helix switching to the inactive conformations. On account of the hydrogen bond formed between E1167 and arginine in αAL, the side chain of E1167 is constraint around the charged R1275 and hence the disruption of the K-E salt bridge is hard to recover. Upon introduction of the R1275Q mutation with an uncharged and rather shorter side chain, the chance to hold E1167 to face outward the active site is reduced significantly. In another word, residue E1167 is more favor to make polar interaction with K1159 and maintains the αC-in conformation in the R1275Q system. Moreover, another hydrogen bond is observed between E1167 and F1269 (the phenylalanine belongs to the DFG motif) in the active site of R1275Q (Figure 4B). This hydrogen bond then produces a force, pulling αC-helix inward, which further prevents the outward shift of E1167 and thus assists the formation of the K-E salt bridge between αC-helix and β3-strand. In summary, the global effect of R1275Q mutation is to stabilize the K-E salt bridge and to elevate the free energy barriers for the αC-helix transition from the αC-in to αC-out conformation. R1275Q mutation decreases the outward shift of E1167 and enhancing its interaction with K1150 to keep a αC-in conformation of ALK, thereby protect the indispensable structural prerequisite for the kinase

ACS Paragon Plus Environment

Page 12 of 38

Page 13 of 38 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Chemical Neuroscience

activating process. Y1278S Shifts the Equilibrium to Activation State. In the case of the Y1278S mutant, the changes in conformational FEL is more visible and the equilibrium is shifted toward the active conformation (Figure 2A). The typical inactive conformation observed in Apo-WT has disappeared in Y1278S. The distance between the two atoms that forms the catalytic salt bridge sampled in Y1278S system is mainly populated less than 4Å. In analogue to the Apo-WT, the second populated energy minima incorporates the structures with both the K-E salt bridge formation and disruption conformations, indicating that the Y1278S mutation has little effect on the conformational switching of the αC-helix. For the conformation of A-loop, the free-energy landscape exhibits remarkable differences with that of Apo-WT. In the deepest free energy minimum, the deviation of A-loop in Y1278S distributes in the range between 3 Å to 10 Å, which is on behalf of the transition between the semi-closed and open conformations. When the representative structure of Y1278S is compared with the Apo-WT structures, it is observed the Y1278S mutation conduces the increased extension of A-loop (Figure 2B and 2D). In the most populated conformation of Apo-WT, the side chain aromatic ring of native Y1278 is favoring the establishment of a π-stacking hydrophobic cluster by packing against other three adjacent aromatic rings of phenylalanine (Figure 5A), including the phenylalanine F1174 at the base of αC-helix, the F1269 belonging the DFG motif, and F1245 situated right below the αC-helix. This hydrophobic cluster contributes the π-stacking interactions leading to an

ACS Paragon Plus Environment

ACS Chemical Neuroscience 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

extended aromatic box,which may take considerable effect on anchoring the relative positions of the aromatic side chains of these involving residues, thus maintains the relatively stationary position of αC-helix and αAL in Apo-WT system. Substitution of Y1278 with serine is expected to destabilize this hydrophobic cluster, as serine cannot be favorable accommodated in the π-stacking hydrophobic pocket occupied by Y1278. Evidently, the shorter side chain substitution S1278 obtains additional mobile space and thus increases the fluctuation of the end of αAL accordingly. The increased flexibility of the initial part of A-loop then gains the opportunity to leads to long-range process reorientation of the A-loop in Y1278S. It is worth noting that the increased flexibility induced by Y1278S mutation promotes the A-loop to a more open conformation rather than a more closed one. The underlying molecular mechanisms responsible for the directed conformational alteration of A-loop may origin from the structural unsuitability of the mutated S1278 in establishing the typical inactive-like conformation. In the typical inactive-like state of Apo-WT, residue Y1278 adopts the conformation by hydrogen bonded with the αC-helix aspartic acid D1163 (Figure 3B). The electrostatic interaction is considered to be one of the forceful factor that maintain the relative positions of αC-helix and initial part of A-loop and keeps their proper conformations accommodated in the inactive kinase conformation. Substituted S1278 with a shorter side chain is incompetent to form the hydrogen bond with D1163, thus will be unfavorable for stabilization of the inactive-like conformation. In addition, the steric flexible S1278 induces concomitant interaction changes which could further disrupt the mutual

ACS Paragon Plus Environment

Page 14 of 38

Page 15 of 38 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Chemical Neuroscience

stabilization of the closed A-loop conformation. As shown in Figure 5B, the Y1278S mutation promotes the fluctuation of the neighboring “RAS” sequence motif, which has been predicted as a controlling unit for the activation process of ALK.53 In the A-loop closed conformation, R1279, A1280, and S1281 hydrogen bonded with residues situated distal part of A-loop, including Y1283, R1284, and K1285, respectively (Figure 5C). These hydrogen bonds help construct the second A-loop helical segment and thus restraint these distal residues of A-loop in the position obstructed the peptide substrate binding pocket. In Y1278S, the set of hydrogen bonds have dismantled with the fluctuation of “RAS” motif (Figure 5D). Without the electrostatic attraction, the distal residues are able to shift outward by little energy penalty, thus leading to an extended A-loop conformation. In addition, it is rationed to speculate that the less restraint S1278 is more easily to exposure in the solvent, which could increase the possibility of the phosphorylation of S1278. Hence, the activation of ALK may be more facile. Overall, the Y1278S substitution is readily accommodated in the active form of the kinase. The substitution of tyrosine with serine could disrupt these connections favored for the ALK inactive state, thereby promotes the activation balance to shift toward active state. Allosteric Effects of ATP Promotes Activation of ALK. Many allosteric effects on the kinase activation on account of ATP binding has been reported.27, 28 Herein, the ATP induced allosteric effects are observed directly at an atomic level by the accelerated MD simulations. As shown in Figure 2A, two energy

ACS Paragon Plus Environment

ACS Chemical Neuroscience 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

basins can be found in the free energy landscape of ATP-bound system. The first minima represents the ensembles that correspond to the partial inactive conformation which is similar with that emerged in Apo-WT system. The second minima which is the also the most populated one represents the ensembles subjected to the allosteric effect induced by ATP, which is regarded as the activate-like conformation of ALK with the αC-in and A-loop open conformation (Figure 2E). Compared with the representative structures of Apo-WT, we can suggest that ATP binding may bring the allosteric effect on promoting the connection between the active site and the two activated related elements, αC-helix and A-loop, thus establish a more compact active site in the ALK kinase domain. The allosteric details are mainly reflected in two aspect, one is the e αC-helix approaching to the active site, the other is the 4 Å shift of αAL approaching to the active site. As shown in Figure 6A, the triphosphate moiety of ATP engages in the electrostatic interactions with K1150 who participates in the catalytic salt bridge. The favorable electrostatic attraction anchors the side chain of K1150 orientated to the salt bridge partner E1167, thus protects the stability of the catalytic salt bridge. Another crucial connection is the hydrogen bond formed between the triphosphate moiety of ATP and the αAL arginine R1275. As mentioned above, the positive charged R1275 plays important role in the connection between αC-helix and A-loop by forming hydrogen bond with the αC-helix aspartic acid D1063. Herein, the interactions among ATP, R1275, and D1063 construct a hydrogen bonds chain (ATP-R1275-D1063), which further pulls the αC-helix approaching to the active site.

ACS Paragon Plus Environment

Page 16 of 38

Page 17 of 38 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Chemical Neuroscience

In addition to the αC-helix, aMD simulations shows that A-loop undergoes a conformational rearrangement when ATP binds to ALK. In the ATP-bound system, the A-loop closed conformation has no longer been sampled, whereas the A-loop prefer to assume an open conformation where the sharp U-turn at the end of αAL has vanished (Figure 2E). In this open conformation, A-loop is fully extended to make space favor for the peptide substrate binding, which is regarded as a determinant step for the kinase activation.8,

10

This observation also suggested that the outward motion of

A-loop could be coupled to the intrinsic interactional changes in the regions of the active site and the A-loop (Figure 6A and 6B). In the Apo-WT, D1276 perfectly construct the hydrogen bond with the first side chain nitrogen atom of R1284, thereby generates the sharp U-turn of the A-loop, a structural controlling instrument for shaping the closed or semi-closed A-loop conformations (Figure 7A and 7B). By contrast, the αAL moves nearly 4 Å forward due to the pulling of ATP through the hydrogen bond between ATP and R1275 (Figure 6B). Due to spatial mismatch, main chain oxygen atom of D1276 losses the opportunity to make the crucial hydrogen bond with the first side chain nitrogen atom of R1284, which leads to the destruction of constrained U-turn and long-range process of reorientation of the A-loop into an open conformation (Figure 7C). In addition to the disruption of the U-turn, some structural features in favor of the open conformation of A-loop are observed in the ATP-bound system. As shown in Figure 7D, the A-loop distal residues C1288 and K1285 form a hydrogen bond, thus moor the distal residues in a position that no longer blocks the substrate binding site in ATP-bound ALK. Besides, residue M1273

ACS Paragon Plus Environment

ACS Chemical Neuroscience 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

at the initial part of A-loop packs against with P1292 and V1293 to form a hydrophobic stem in ATP-bound ALK (Figure 7D), which is similar to the hydrophobic stem in Apo-WT. Notably, the residues participate in this hydrophobic stem are more departure from the active site than the residues involve in the hydrophobic stem in Apo-WT system. This structural feature has been observed in active RTKs to firming the open A-loop conformation.8 As shown in Figure 7D, P1292 and V1293 participate hydrophobic interactions with L1291, P1298, F1301, and L1339, thereupon render the terminal part of A-loop adjacent to αG helix to consolidate the open conformation. Altogether, the allosteric performance of ATP in ALK is the induction of a more compact active site. This allosteric effect of ATP on the activation of ALK is twofold: one is that it stabilizes the αC-in conformation; the other is that it destabilizes the U-turn structure of A-loop and promote it to an open conformation. CONCLUSIONS Anaplastic Lymphoma Kinase (ALK) has been established as a therapeutic molecular target in neuroblastoma

19, 22

. It has been suggested that alterations in the equilibrium

between its inactive and active conformations are coupled with its oncogenic potential.19 To understand the mechanism of activation, herein, the structural features and conformation dynamics of the different forms of ALK are explored through accelerate MD (aMD) enhanced sampling. Our results indicate that the active state is rarely visited in the wild type ALK (Apo-WT). This phenomenon is supposed to attribute to the low catalytic activity,54

ACS Paragon Plus Environment

Page 18 of 38

Page 19 of 38 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Chemical Neuroscience

which is similar with the behavior of EGFR monomer.55, 56 The energy minimum state identified form aMD of Apo-WT is partially inactive which is in similar with the X-ray structure of ALK.3 Additionally, the simulation of Apo-WT also revealed new low-energy state corresponding to the typical receptor tyrosine kinase (RTK) inactive conformation. Regarding to the conformational stabilization in each conformation, the hydrogen bond networks play a crucial role in the construction of the conformational framework, while the hydrophobic interactions further cemented the respective structure. In the case of mutations and ATP-bound ALK, the equilibriums are shifted to more active state to different extents. Substitution of R1275 with glutamine is expected to secure αC-helix assuming a “αC-in” conformation, which is favored for the formation of structural important K-E salt bridge. For Y1278S, the foremost variation is the disruption of the π-stacking hydrophobic box. The substitution of tyrosine with serine releases the terminal of αAL helix, and therefore disrupts the U-turn which is critical to maintain the closed and semi-closed conformation of A-loop, and hence shift the A-loop to an open conformation. For the ATP-bound ALK, the conformation equilibrium is shifted to the active state by the allosteric effects induced by ATP. Interacting with R1275 situated in αAL, ATP promotes a more compact active site which is favoring for the αC-helix assuming the “in” conformation. In addition, binding of ATP may force considerable migration of αAL, which is accommodated in the reorientation of A-loop to an open conformation. The complex changes seen in the conformational transition of the mutants and ATP-bound systems are attributed to both the local alterations in the network of hydrogen bonding

ACS Paragon Plus Environment

ACS Chemical Neuroscience 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

and hydrophobic interactions and to long-range allosteric couplings. Collectively, the aMD simulations present novel ALK conformations which is corresponding to the typical inactive and active state of RTKs, and gave rational explanations for the activation molecular origins induced by R1275Q, Y1275S mutations, as well as by ATP binding. The gained atomistic insight may provide valuable information for further atomistic and mechanistic studies of the oncogenic properties of ALK and as a background in designing more effective inhibitors as potential drugs.

ACS Paragon Plus Environment

Page 20 of 38

Page 21 of 38 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Chemical Neuroscience

METHODS Simulation Systems Preparation. The initial structure for native ALK kinase domain was obtained from the Protein Data Bank (PDB), with PDB code 3L9P.3 Based on the WT structure, we generated the starting structures of R1275Q and Y1278S mutants by Discovery studio (DS) (Discovery Studio 4.1, Accelrys, San Diego, CA). The starting point for ATP-bound system was obtained by superposition of a homologous ATP-containing protein kinase (insulin receptor kinase in complex with Mg2+ and ATP; PDB code 3BU5 )57 on the ADP-bound ALK structure (PDB code 3LCT).3 Thus, the ATP molecule coupling with Mg2+ was optimally located in the nucleotide-binding site of ALK. All missing hydrogen atoms were added to ALKs using the Leap program form AMBER16 package.58 Sodium ions (Na+) were added to each system based on a coulomb potential grip in order to keep the whole system neutral. Each system was then solvated with the TIP3P water box model in truncated octahedron box in which protein atom was at least a 10 Å distance away from the nearest edge of the box. Molecular Dynamics Simulation. Four independent conventional MD (cMD) simulations were initially performed for each starting system, wild type apo form ALK (Apo-WT), R1275Q mutant, Y1278S mutant, and bound by ATP (ATP-bound), using AMBER16 software package employed the AMBER ff14SB force field.59 The particle Mesh Ewald (PME) algorithm was applied for long-range electrostatics interactions calculation, and a non-bonded cutoff of 10 Å was used for real space interactions. The SHAKE

ACS Paragon Plus Environment

ACS Chemical Neuroscience 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

algorithm was utilized to restrain all covalent bonds involving hydrogen atoms. The time step was set to 2 fs. Firstly, two steps of minimizations were exerted prior to cMD simulations, including a 5000 steps of the steepest descent (SD) method and a further 5000 steps of conjugate gradient minimization. Secondly, each system was gently heated from 0 to 300 K in 300 ps by Langevin dynamics and then followed by 500 ps of constant pressure equilibration with everything released in NVT conditions. Finally, 50 ns cMD production runs for all four systems were then performed under NPT ensemble condition with periodic boundary conditions. Accelerated Molecular Dynamics. Dual-boost aMD is a good enhanced sampling approach in studying the conformational transition of proteins, which has the capability to facilitate more sampling of the intermediate states and possibly even conversion between inactive and active states.46, 60 In this study, four dual- boost aMD simulations of ALK were initiate from the final structure of the corresponding cMD simulations. In dual-boost aMD, two sets of parameters are need to be specified are Etotal (the inverse strength boost factor for the total potential energy), αtotal (average total potential energy threshold), Edihed (average dihedral energy threshold); αdihed (inverse strength boost factor for the dihedral energy). The input parameters for aMD are as follows:61 Etotal=Etotal_avg + Natoms × 0.16

(1)

αtotal= Natoms × 0.16

(2)

Edihed=Edihed_avg+Nresdues ×4

(3)

ACS Paragon Plus Environment

Page 22 of 38

Page 23 of 38 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Chemical Neuroscience

αdihed=Nresdues ×4/5

(4)

where Natoms and Nresdues are the total number of atoms and residues, respectively, Etotal_avg and Edihed_avg are the average total and dihedral energies calculated form the conventional MD simulations, respectively. In all of the four ALK systems, 500 ns aMD simulations were performed after the cMD simulations. Potential of Mean Force Calculation. The potential of mean force (PMF)

62

was used to generate the free energy landscape

(FEL) to study the free energetic change that accompanied the conformational change of ALK. The free energy landscape is calculated from the simulation trajectories using the following equation: ∆G (x, y) = kBT ln g(x, y) Where kB is Boltzmann constant, T is the simulation temperature, and g(x, y) is the normalized joint probability distribution. The minimum energy is set to zero. In this work, we use K1150-E1167 distance and A-loop deviation as the two reaction coordinates. The free energy landscapes were generated using 50 bins in both x and y directions. The free energy landscape includes the entire aMD trajectory data set. Clustering method isperformed to obtain a representative structure of each energy basin. In the most populated cluster of each basin, the central structure, i.e., the structure with smallest distance to all the other members of the cluster, was selected as representative of the basin.

ACS Paragon Plus Environment

ACS Chemical Neuroscience 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ASSOCIATED CONTENT Supporting Information Figure S1: Comparison of active structure of insulin receptor (A) and inactive structure of insulin receptor (B) and partial inactive structure of ALK (C). Figure S2: Deviation of A-loop in Apo-WT ALK system. Table S1-S2: summary of key hydrogen bonds in partial inactive and typical inactive-like Apo-WT ALK. (PDF)

AUTHOR INFORMATION Corresponding Author [email protected] Notes The authors declare no competing financial interest.

Author Contributions Qing-Chuan Zheng and Hong-Xing Zhang conceived the project. Mu-Yang He and Wei-Kang Li designed and carried out the experiments. All of the authors analyzed the data and discussed the results. The manuscript was written through contributions of all of the authors. All of the authors read and approved the final manuscript.

ACKNOWLEDGMENTS This work is supported by Natural Science Foundation of China (Grant Nos.21273095

ACS Paragon Plus Environment

Page 24 of 38

Page 25 of 38 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Chemical Neuroscience

and Grant Nos.21773084).

ACS Paragon Plus Environment

ACS Chemical Neuroscience 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

REFERENCE 1. Iwahara, T., Fujimoto, J., Wen, D., Cupples, R., Bucay, N., Arakawa, T., Mori, S., Ratzkin, B., and Yamamoto, T. (1997) Molecular characterization of ALK, a receptor tyrosine kinase expressed specifically in the nervous system, Oncogene 14, 439-449. 2. Roskoski, R., Jr. (2013) Anaplastic lymphoma kinase (ALK): Structure, oncogenic activation, and pharmacological inhibition, Pharmacological Research 68, 68-94. 3. Lee, C. C., Jia, Y., Li, N., Sun, X., Ng, K., Ambing, E., Gao, M.-Y., Hua, S., Chen, C., Kim, S., Michellys, P.-Y., Lesley, S. A., Harris, J. L., and Spraggon, G. (2010) Crystal structure of the ALK (anaplastic lymphoma kinase) catalytic domain, Biochemical Journal 430, 425-437. 4. Hubbard, S. R. (1997) Crystal structure of the activated insulin receptor tyrosine kinase in complex with peptide substrate and ATP analog, The EMBO journal 16, 5572-5581. 5. Favelyukis, S., Till, J. H., Hubbard, S. R., and Miller, W. T. (2001) Structure and autoregulation of the insulin-like growth factor 1 receptor kinase, Nature Structural and Molecular Biology 8, 1058. 6. Munshi, S., Hall, D. L., Kornienko, M., Darke, P. L., and Kuo, L. C. (2003) Structure of apo, unactivated insulin-like growth factor-1 receptor kinase at 1.5 Å resolution, Acta Crystallographica Section D: Biological Crystallography 59, 1725-1730. 7. Hubbard, S. R., Wei, L., and Hendrickson, W. A. (1994) Crystal structure of the tyrosine kinase domain of the human insulin receptor, Nature 372, 746. 8. Hubbard, S. R. (1997) Crystal structure of the activated insulin receptor tyrosine kinase in complex with peptide substrate and ATP analog, The EMBO journal 16, 5572-5581. 9. Hubbard, S. R., Wei, L., Ellis, L., and Hendrickson, W. A. (1994) Crystal structure of the tyrosine kinase domain of the human insulin receptor, Nature 372, 746-754. 10. Hanks, S. K., Quinn, A. M., and Hunter, T. (1988) The protein kinase family: conserved features and deduced phylogeny of the catalytic domains, Science (New York, N.Y.) 241, 42-52. 11. Shen, Q., Cheng, F., Song, H., Lu, W., Zhao, J., An, X., Liu, M., Chen, G., Zhao, Z., and Zhang, J. (2017) Proteome-Scale Investigation of Protein Allosteric Regulation Perturbed by Somatic Mutations in 7,000 Cancer Genomes, American Journal of Human Genetics 100, 5-20. 12. Jiang, H., Deng, R., Yang, X., Shang, J., Lu, S., Zhao, Y., Song, K., Liu, X., Zhang, Q., Chen, Y., Chinn, Y. E., Wu, G., Li, J., Chen, G., Yu, J., and Zhang, J. (2017) Peptidomimetic inhibitors of APC-Asef interaction block colorectal cancer migration, Nature Chemical Biology 13, 994-+. 13. Huang, W., Wang, G., Shen, Q., Liu, X., Lu, S., Geng, L., Huang, Z., and Zhang, J. (2015) ASBench: benchmarking sets for allosteric discovery, Bioinformatics 31, 2598-2600. 14. Mesaros, E. F., Ott, G. R., and Dorsey, B. D. (2014) Anaplastic lymphoma kinase inhibitors as anticancer therapeutics: a patent review, Expert Opinion on Therapeutic Patents 24, 417-442.

ACS Paragon Plus Environment

Page 26 of 38

Page 27 of 38 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Chemical Neuroscience

15. Palmer, R. H., Vernersson, E., Grabbe, C., and Hallberg, B. (2009) Anaplastic lymphoma kinase: signalling in development and disease, Biochemical Journal 420, 345-361. 16. Mosse, Y. P., Laudenslager, M., Longo, L., Cole, K. A., Wood, A., Attiyeh, E. F., Laquaglia, M. J., Sennett, R., Lynch, J. E., Perri, P., Laureys, G., Speleman, F., Kim, C., Hou, C., Hakonarson, H., Torkamani, A., Schork, N. J., Brodeur, G. M., Tonini, G. P., Rappaport, E., Devoto, M., and Maris, J. M. (2008) Identification of ALK as a major familial neuroblastoma predisposition gene, Nature 455, 930-U922. 17. Soda, M., Choi, Y. L., Enomoto, M., Takada, S., Yamashita, Y., Ishikawa, S., Fujiwara, S.-i., Watanabe, H., Kurashina, K., Hatanaka, H., Bando, M., Ohno, S., Ishikawa, Y., Aburatani, H., Niki, T., Sohara, Y., Sugiyama, Y., and Mano, H. (2007) Identification of the transforming EML4-ALK fusion gene in non-small-cell lung cancer, Nature 448, 561-U563. 18. Griffin, C. A., Hawkins, A. L., Dvorak, C., Henkle, C., Ellingham, T., and Perlman, E. J. (1999) Recurrent involvement of 2p23 in inflammatory myofibroblastic tumors, Cancer Research 59, 2776-2780. 19. Bresler, S. C., Weiser, D. A., Huwe, P. J., Park, J. H., Krytska, K., Ryles, H., Laudenslager, M., Rappaport, E. F., Wood, A. C., McGrady, P. W., Hogarty, M. D., London, W. B., Radhakrishnan, R., Lemmon, M. A., and Mosse, Y. P. (2014) ALK Mutations Confer Differential Oncogenic Activation and Sensitivity to ALK Inhibition Therapy in Neuroblastoma, Cancer Cell 26, 682-694. 20. Cheung, N.-K. V., and Dyer, M. A. (2013) Neuroblastoma: developmental biology, cancer genomics and immunotherapy, Nature Reviews Cancer 13, 397-411. 21. Janoueix-Lerosey, I., Lequin, D., Brugieres, L., Ribeiro, A., de Pontual, L., Combaret, V., Raynal, V., Puisieux, A., Schleiermacher, G., Pierron, G., Valteau-Couanet, D., Frebourg, T., Michon, J., Lyonnet, S., Amiel, J., and Delattre, O. (2008) Somatic and germline activating mutations of the ALK kinase receptor in neuroblastoma, Nature 455, 967-U951. 22. Schleiermacher, G., Javanmardi, N., Bernard, V., Leroy, Q., Cappo, J., Rio-Frio, T., Pierron, G., Combaret, V., Speleman, F., Djos, A., Ora, I., Hedborg, F., Holmqvist, B.-M., Abrahamsson, J., Peuchmaur, M., Michon, J., Janoueix-Lerosey, I., Kogner, P., Delattre, O., and Martinsson, T. (2014) Emergence of new ALK mutations at relapse of neuroblastoma, Journal of Clinical Oncology 32. 23. Johnson, L. N., and Barford, D. (1993) The effects of phosphorylation on the structure and function of proteins, Annual review of biophysics and biomolecular structure 22, 199-232. 24. Chen, Y., Takita, J., Choi, Y. L., Kato, M., Ohira, M., Sanada, M., Wang, L., Soda, M., Kikuchi, A., Igarashi, T., Nakagawara, A., Hayashi, Y., Mano, H., and Ogawa, S. (2008) Oncogenic mutations of ALK kinase in neuroblastoma, Nature 455, 971-U956. 25. Huang, Z., Zhu, L., Cao, Y., Wu, G., Liu, X., Chen, Y., Wang, Q., Shi, T., Zhao, Y., Wang, Y., Li, W., Li, Y., Chen, H., Chen, G., and Zhang, J. (2011) ASD: a comprehensive database of allosteric proteins and modulators, Nucleic Acids Research 39, D663-D669. 26. Huang, Z., Mou, L., Shen, Q., Lu, S., Li, C., Liu, X., Wang, G., Li, S., Geng, L., Liu, Y.,

ACS Paragon Plus Environment

ACS Chemical Neuroscience 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Wu, J., Chen, G., and Zhang, J. (2014) ASD v2.0: updated content and novel features focusing on allosteric regulation, Nucleic Acids Research 42, D510-D516. 27. Hattori, M., and Gouaux, E. (2012) Molecular mechanism of ATP binding and ion channel activation in P2X receptors, Nature 485, 207-U291. 28. Lu, S., Deng, R., Jiang, H., Song, H., Li, S., Shen, Q., Huang, W., Nussinov, R., Yu, J., and Zhang, J. (2015) The Mechanism of ATP-Dependent Allosteric Protection of Akt Kinase Phosphorylation, Structure 23, 1725-1734. 29. Bilodeau, M. T., Cunningham, A. M., Hungate, R. W., and Koester, T. J. (2000) Tyrosine kinase inhibitors, Google Patents. 30. Palmer, R. H., Vernersson, E., Grabbe, C., and Hallberg, B. (2009) Anaplastic lymphoma kinase: signalling in development and disease, Biochemical Journal 420, 345-361. 31. Shaw, A. T., and Solomon, B. (2011) Targeting anaplastic lymphoma kinase in lung cancer, Clinical Cancer Research 17, 2081-2086. 32. Iwahara, T., Fujimoto, J., Wen, D., Cupples, R., Bucay, N., Arakawa, T., Mori, S., Ratzkin, B., and Yamamoto, T. (1997) Molecular characterization of ALK, a receptor tyrosine kinase expressed specifically in the nervous system, Oncogene 14, 439. 33. Morris, S. W., Kirstein, M. N., Valentine, M. B., Dittmer, K. G., Shapiro, D. N., Saltman, D. L., and Look, A. T. (1994) Fusion of a kinase gene, ALK, to a nucleolar protein gene, NPM, in non-Hodgkin's lymphoma, Science 263, 1281-1284. 34. Sun, H., Li, Y., Tian, S., Wang, J., and Hou, T. (2014) P-loop Conformation Governed Crizotinib Resistance in G2032R-Mutated ROS1 Tyrosine Kinase: Clues from Free Energy Landscape, Plos Computational Biology 10. 35. Khoddami, M., Nadri, H., Moradi, A., and Sakhteman, A. (2015) Homology modeling, molecular dynamic simulation, and docking based binding site analysis of human dopamine (D4) receptor, Journal of Molecular Modeling 21. 36. Cui, Y.-L., Zheng, Q.-C., Zhang, J.-L., Xue, Q., Wang, Y., and Zhang, H.-X. (2013) Molecular Dynamic Investigations of the Mutational Effects on Structural Characteristics and Tunnel Geometry in CYP17A1, Journal of Chemical Information and Modeling 53, 3308-3317. 37. Park, J., McDonald, J. J., Petter, R. C., and Houk, K. N. (2016) Molecular Dynamics Analysis of Binding of Kinase Inhibitors to WT EGFR and the T790M Mutant, Journal of Chemical Theory and Computation 12, 2066-2078. 38. Chen, J., Wang, J., and Zhu, W. (2017) Mutation L1196M-induced conformational changes and the drug resistant mechanism of anaplastic lymphoma kinase studied by free energy perturbation and umbrella sampling, Physical Chemistry Chemical Physics 19, 30239-30248. 39. Ni, Z., and Zhang, T.-C. (2015) Computationally unraveling how ceritinib overcomes drug-resistance mutations in ALK-rearranged lung cancer, Journal of Molecular Modeling 21. 40. Nagarajan, N., Alphonse, C. R. W., Gnana, P. V. S., and Rajaretinam, R. K. (2017) Molecular Dynamics Validation of Crizotinib Resistance to ALK Mutations (L1196M and G1269A) and Identification of Specific Inhibitors, J. Cell. Biochem. 118, 3462-3471.

ACS Paragon Plus Environment

Page 28 of 38

Page 29 of 38 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Chemical Neuroscience

41. Sun, H., Li, Y., Li, D., and Hou, T. (2013) Insight into Crizotinib Resistance Mechanisms Caused by Three Mutations in ALK Tyrosine Kinase using Free Energy Calculation Approaches, Journal of Chemical Information and Modeling 53, 2376-2389. 42. Pan, P., Yu, H., Liu, Q., Kong, X., Chen, H., Chen, J., Liu, Q., Li, D., Kang, Y., Sun, H., Zhou, W., Tian, S., Cui, S., Zhu, F., Li, Y., Huang, Y., and Hou, T. (2017) Combating Drug-Resistant Mutants of Anaplastic Lymphoma Kinase with Potent and Selective Type-I-1/2 Inhibitors by Stabilizing Unique DFG-Shifted Loop Conformation, Acs Central Science 3, 1208-1220. 43. Kong, X., Sun, H., Pan, P., Zhu, F., Chang, S., Xu, L., Li, Y., and Hou, T. (2018) Importance of protein flexibility in molecular recognition: a case study on Type-I1/2 inhibitors of ALK, Physical Chemistry Chemical Physics 20, 4851-4863. 44. Kong, X., Pan, P., Li, D., Tian, S., Li, Y., and Hou, T. (2015) Importance of protein flexibility in ranking inhibitor affinities: modeling the binding mechanisms of piperidine carboxamides as Type I1/2 ALK inhibitors, Physical Chemistry Chemical Physics 17, 6098-6113. 45. Markwick, P. R. L., and McCammon, J. A. (2011) Studying functional dynamics in bio-molecules using accelerated molecular dynamics, Physical Chemistry Chemical Physics 13, 20053-20065. 46. Hamelberg, D., Mongan, J., and McCammon, J. A. (2004) Accelerated molecular dynamics: A promising and efficient simulation method for biomolecules, Journal of Chemical Physics 120, 11919-11929. 47. Markwick, P. R. L., Pierce, L. C. T., Goodin, D. B., and McCammon, J. A. (2011) Adaptive Accelerated Molecular Dynamics (Ad-AMD) Revealing the Molecular Plasticity of P450cam, Journal of Physical Chemistry Letters 2, 158-164. 48. Wereszczynski, J., and McCammon, J. A. (2012) Nucleotide-dependent mechanism of Get3 as elucidated from free energy calculations, Proceedings of the National Academy of Sciences of the United States of America 109, 7759-7764. 49. Kumar, S., and Nussinov, R. (2002) Close-range electrostatic interactions in proteins, Chembiochem 3, 604-617. 50. Zhang, X., Gureasko, J., Shen, K., Cole, P. A., and Kuriyan, J. (2006) An allosteric mechanism for activation of the kinase domain of epidermal growth factor receptor, Cell 125, 1137-1149. 51. Raaf, J., Issinger, O.-G., and Niefind, K. (2009) First inactive conformation of CK2α, the catalytic subunit of protein kinase CK2, Journal of molecular biology 386, 1212-1221. 52. Xu, W., Doshi, A., Lei, M., Eck, M. J., and Harrison, S. C. (1999) Crystal structures of c-Src reveal features of its autoinhibitory mechanism, Molecular cell 3, 629-638. 53. Donella-Deana, A., Marin, O., Cesaro, L., Gunby, R. H., Ferrarese, A., Coluccia, A. M., Tartari, C. J., Mologni, L., Scapozza, L., and Gambacorti-Passerini, C. (2005) Unique substrate specificity of anaplastic lymphoma kinase (ALK): development of phosphoacceptor peptides for the assay of ALK activity, Biochemistry 44, 8533-8542. 54. Donella-Deana, A., Marin, O., Cesaro, L., Gunby, R. H., Ferrarese, A., Coluccia, A. M. L., Tartari, C. J., Mologni, L., Scapozza, L., Gambacorti-Passerini, C., and Pinna, L. A.

ACS Paragon Plus Environment

ACS Chemical Neuroscience 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

(2005) Unique substrate specificity of anaplastic lymphoma kinase (ALK): Development of phosphoacceptor peptides for the assay of ALK activity, Biochemistry 44, 8533-8542. 55. Shan, Y., Arkhipov, A., Kim, E. T., Pan, A. C., and Shaw, D. E. (2013) Transitions to catalytically inactive conformations in EGFR kinase, Proceedings of the National Academy of Sciences of the United States of America 110, 7270-7275. 56. Sutto, L., and Luigi Gervasio, F. (2013) Effects of oncogenic mutations on the conformational free-energy landscape of EGFR kinase, Proceedings of the National Academy of Sciences of the United States of America 110, 10616-10621. 57. Wu, J., Tseng, Y. D., Xu, C.-F., Neubert, T. A., White, M. F., and Hubbard, S. R. (2008) Structural and biochemical characterization of the KRLB region in insulin receptor substrate-2, Nature Structural & Molecular Biology 15, 251-258. 58. D.A. Case, D. S. C., T.E. Cheatham, III, T.A. Darden, R.E. Duke, T.J. Giese, H. Gohlke, A.W. Goetz, D. Greene, N. Homeyer, S. Izadi, A. Kovalenko, T.S. Lee, S. LeGrand, P. Li, C. Lin, J. Liu, T. Luchko, R. Luo, D. Mermelstein, K.M. Merz, G. Monard, H. Nguyen, I. Omelyan, A. Onufriev, F. Pan, R. Qi, D.R. Roe, A. Roitberg, C. Sagui, C.L. Simmerling, W.M. Botello-Smith, J. Swails, R.C. Walker, J. Wang, R.M. Wolf, X. Wu, L. Xiao, D.M. York and P.A. Kollman. (2017) AMBER 2017,, In University of California, San Francisco. 59. Wang, J. M., Wolf, R. M., Caldwell, J. W., Kollman, P. A., and Case, D. A. (2005) Development and testing of a general amber force field (vol 25, pg 1157, 2004), Journal of Computational Chemistry 26, 114-114. 60. Hamelberg, D., de Oliveira, C. A. F., and McCammon, J. A. (2007) Sampling of slow diffusive conformational transitions with accelerated molecular dynamics, Journal of Chemical Physics 127. 61. Pierce, L. C. T., Salomon-Ferrer, R., de Oliveira, C. A. F., McCammon, J. A., and Walker, R. C. (2012) Routine Access to Millisecond Time Scale Events with Accelerated Molecular Dynamics, Journal of Chemical Theory and Computation 8, 2997-3002. 62. Roux, B. (1995) The calculation of the potential of mean force using computer simulations, Computer physics communications 91, 275-282.

ACS Paragon Plus Environment

Page 30 of 38

Page 31 of 38 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Chemical Neuroscience

For Table of Contents Use Only

Conformational Transition of Key Structural Features Involved

in

Activation

of

ALK

Induced

by

two

Neuroblastoma Mutations and ATP Binding: Insight from Accelerated Molecular Dynamics Simulations

Mu-Yang He,2 Wei-Kang Li,2Qing-Chuan Zheng1,2*Hong-Xing Zhang2*

ACS Paragon Plus Environment

ACS Chemical Neuroscience 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Figure 1. Comparison of active structure of insulin receptor (A) and inactive structure of insulin receptor (B) and partial inactive structure of ALK (C).Key structural elements are colored in yellow (αC-helix) and purple (A-loop). Conserved residues form the catalytic salt bridge are shown in sticks. (D) ATP and peptide substrate binding site in the kinase domain of ALK. ATP binding site is shown as red mesh, peptide substrate binding site is shown in slate mesh. 165x52mm (300 x 300 DPI)

ACS Paragon Plus Environment

Page 32 of 38

Page 33 of 38 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Chemical Neuroscience

Figure 2. Free energy landscape and representative structure of Apo-WT ALK and two mutants and ATPbound ALK. (A) FEL of four system as a function of CV1 and CV2 in Å. (B) Representative structure of ApoWT. (C) Representative structure of R1275Q. (D) Representative structure of Y1278S. (D) Representative structure of ATP-bound ALK. In the representative structures shown below the FEL the αC-helix is colored in yellow, the A-loop is colored in purple. 165x83mm (300 x 300 DPI)

ACS Paragon Plus Environment

ACS Chemical Neuroscience 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Figure 3. Representative structures of Apo-WT system. Residues involve in conformation stationary interactions are shown in sticks. The αC-helix is colored in yellow and the A-loop is colored in purple. Hydrogen bonds are shown in yellow dash. Residues involve in the hydrophobic stem is colored in orange. Residues involve in the local hydrophobic cluster is colored in cyan. (A) Internal interactions of the partial inactive conformation of Apo-WT. (B) Internal interactions of the typical inactive-like conformation of ApoWT. (C) Hydrophobic cluster around L1291 and P1292 in the typical inactive-like conformation of Apo-WT. (D) Hydrophobic cluster around V1239 in the typical inactive-like conformation of Apo-WT. 165x80mm (300 x 300 DPI)

ACS Paragon Plus Environment

Page 34 of 38

Page 35 of 38 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Chemical Neuroscience

Figure 4. K-E salt bridge comparison among the typical Apo-WT ALK and R1275Q mutant. The hydrogen bonds are shown in yellow dash. The salt bridge is shown in red dash. The αC-helix is colored in yellow. (A) Interaction between E1167 and R1275 pulls the side chain of E1167 departure from K1150 in typical ApoWT ALK. (B) Interaction between E1167 and F1269 pulls the side chain of E1167 approaching K1150 in R1275Q. 165x60mm (300 x 300 DPI)

ACS Paragon Plus Environment

ACS Chemical Neuroscience 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Figure 5. Conformational variation induced by Y1278S. In all of the representative structures, the αC-helix is colored in yellow and the A-loop is colored in purple. (A) The π-stacking hydrophobic cluster in representative structure of Apo-WT system. The residues involving this cluster are shown in sticks and colored in green. (B) Deviation of “RAS” sequence motif in Y1278S. “RAS” residues in typical inactive-like Apo-WT and Y1278S are colored in white and slate, respectively. (C) Interactions between ”RAS“ with the residues in the distal part of A-loop in typical inactive-like Apo-WT. Hydrogen bonds are shown in yellow dash. (D) Shift of the A-loop distal part residues in Y1278S. In all of the representative structures above, the αC-helix is colored in yellow and the A-loop is colored in purple. 82x71mm (300 x 300 DPI)

ACS Paragon Plus Environment

Page 36 of 38

Page 37 of 38 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Chemical Neuroscience

Figure 6. ATP induces a more compact active site. In ATP-bound ALK representative structure, the αC-helix is colored in yellow and the A-loop is colored in purple. The representative structure of typical inactive-like Apo-WT are colored in white. Hydrogen bonds are shown in yellow dash. (A) Inward shift of αC-helix in ATPbound system. (B) Migration of the αAL in ATP-bound ALK. 82x36mm (300 x 300 DPI)

ACS Paragon Plus Environment

ACS Chemical Neuroscience 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Figure 7. Disruption of the U-turn and extension of A-loop in the ATP-bound ALK. A-loop is colored in purple in each structures. Hydrogen bonds are shown in yellow dash. Residues participate in structure stabilization are shown in white sticks. (A) Formation of U-turn in the partial inactive Apo-WT. (B) Formation of U-turn in the typical inactive-like Apo-WT. (C) Disruption of the U-turn in the ATP-bound ALK. (D) Open A-loop conformation of ATP-bound ALK. Residues involve in hydrophobic stem are colored in orange. 82x56mm (300 x 300 DPI)

ACS Paragon Plus Environment

Page 38 of 38