“Snapshots” of Ispinesib-induced Conformational Changes in the Mitotic Kinesin Eg5

Kinesins comprise a superfamily of molecular motors that drive a wide variety of cellular physiologies, from cytoplasmic transport to formation of the bipolar spindle in mitosis. These differing roles are reflected in corresponding polymorphisms in key kinesin structural elements. One of these is a unique loop and stem motif found in all kinesins and referred to as loop 5 (L5). This loop is longest in the mitotic kinesin Eg5 and is the target for a number of small molecule inhibitors, including ispinesib, which is being used in clinical trials in patients with cancer. In this study, we have used x-ray crystallography to identify a new structure of an Eg5-ispinesib complex and have combined this with transient state kinetics to identify a plausible sequence of conformational changes that occur in response to ispinesib binding. Our results demonstrate that ispinesib-induced structural changes in L5 from Eg5 lead to subsequent changes in the conformation of the switch II loop and helix and in the neck linker. We conclude that L5 in Eg5 simultaneously regulates the structure of both the ATP binding site and the motor''s mechanical elements that generate force.     

Loop L5 Assumes Three Distinct Orientations during the ATPase Cycle of the Mitotic Kinesin Eg5: A TRANSIENT AND TIME-RESOLVED FLUORESCENCE STUDY*

Members of the kinesin superfamily of molecular motors differ in several key structural domains, which probably allows these molecular motors to serve the different physiologies required of them. One of the most variable of these is a stem-loop motif referred to as L5. This loop is longest in the mitotic kinesin Eg5, and previous structural studies have shown that it can assume different conformations in different nucleotide states. However, enzymatic domains often consist of a mixture of conformations whose distribution shifts in response to substrate binding or product release, and this information is not available from the “static” images that structural studies provide. We have addressed this issue in the case of Eg5 by attaching a fluorescent probe to L5 and examining its fluorescence, using both steady state and time-resolved methods. This reveals that L5 assumes an equilibrium mixture of three orientations that differ in their local environment and segmental mobility. Combining these studies with transient state kinetics demonstrates that there is a major shift in this distribution during transitions that interconvert weak and strong microtubule binding states. Finally, in conjunction with previous cryo-EM reconstructions of Eg5·microtubule complexes, these fluorescence studies suggest a model in which L5 regulates both nucleotide and microtubule binding through a set of reversible interactions with helix α3. We propose that these features facilitate the production of sustained opposing force by Eg5, which underlies its role in supporting formation of a bipolar spindle in mitosis.     

All-atom structural investigation of kinesin-microtubule complex constrained by high-quality cryo-electron-microscopy maps

In this study, we have performed a comprehensive structural investigation of three major biochemical states of a kinesin complexed with microtubule under the constraint of high-quality cryo-electron-microscopy (EM) maps. In addition to the ADP and ATP state which were captured by X-ray crystallography, we have also modeled the nucleotide-free or APO state for which no crystal structure is available. We have combined flexible fitting of EM maps with regular molecular dynamics simulations, hydrogen-bond analysis, and free energy calculation. Our APO-state models feature a subdomain rotation involving loop L2 and α6 helix of kinesin, and local structural changes in active site similar to a related motor protein, myosin. We have identified a list of hydrogen bonds involving key residues in the active site and the binding interface between kinesin and microtubule. Some of these hydrogen bonds may play an important role in coupling microtubule binding to ATPase activities in kinesin. We have validated our models by calculating the binding free energy between kinesin and microtubule, which quantitatively accounts for the observation of strong binding in the APO and ATP state and weak binding in the ADP state. This study will offer promising targets for future mutational and functional studies to investigate the mechanism of kinesin motors.     

Probing the structural and energetic basis of kinesin-microtubule binding using computational alanine-scanning mutagenesis

Kinesin-microtubule (MT) binding plays a critical role in facilitating and regulating the motor function of kinesins. To obtain a detailed structural and energetic picture of kinesin-MT binding, we performed large-scale computational alanine-scanning mutagenesis based on long-time molecular dynamics (MD) simulations of the kinesin-MT complex in both ADP and ATP states. First, we built three all-atom kinesin-MT models: human conventional kinesin bound to ADP and mouse KIF1A bound to ADP and ATP. Then, we performed 30 ns MD simulations followed by kinesin-MT binding free energy calculations for both the wild type and mutants obtained after substitution of each charged residue of kinesin with alanine. We found that the kinesin-MT binding free energy is dominated by van der Waals interactions and further enhanced by electrostatic interactions. The calculated mutational changes in kinesin-MT binding free energy are in excellent agreement with results of an experimental alanine-scanning study with a root-mean-square error of ~0.32 kcal/mol [Woehlke, G., et al. (1997) Cell 90, 207-216]. We identified a set of important charged residues involved in the tuning of kinesin-MT binding, which are clustered on several secondary structural elements of kinesin (including well-studied loops L7, L8, L11, and L12, helices α4, α5, and α6, and less-explored loop L2). In particular, we found several key residues that make different contributions to kinesin-MT binding in ADP and ATP states. The mutations of these residues are predicted to fine-tune the motility of kinesin by modulating the conformational transition between the ADP state and the ATP state of kinesin.     

A pathway of structural changes produced by monastrol binding to Eg5

Monastrol is a small molecule inhibitor that is specific for Eg5, a member of the kinesin 5 family of mitotic motors. Crystallographic models of Eg5 in the presence and absence of monastrol revealed that drug binding produces a variety of structural changes in the motor, including in loop L5 and the neck linker. What is not clear from static crystallographic models, however, is the sequence of structural changes produced by drug binding. Furthermore, because crystallographic structures can be influenced by the packing forces in the crystal, it also remains unclear whether these drug-induced changes occur in solution, at physiologically active concentrations of monastrol or of other drugs that target this site. We have addressed these issues by using a series of spectroscopic probes to monitor the structural consequences of drug binding. Our results demonstrated that the crystallographic model of an Eg5-ADP-monastrol ternary complex is consistent with several solution-based spectroscopic probes. Furthermore, the kinetics of these spectroscopic signal changes allowed us to determine the temporal sequence of drug-induced structural transitions. These results suggested that L5 may be an element in the pathway that links the state of the nucleotide-binding site to the neck linker in kinesin motors.     

Structural dynamics of the microtubule binding and regulatory elements in the kinesin-like calmodulin binding protein

Kinesins are molecular motors that power cell division and transport of various proteins and organelles. Their motor activity is driven by ATP hydrolysis and depends on interactions with microtubule tracks. Essential steps in kinesin movement rely on controlled alternate binding to and detaching from the microtubules. The conformational changes in the kinesin motors induced by nucleotide and microtubule binding are coordinated by structural elements within their motor domains. Loop L11 of the kinesin motor domain interacts with the microtubule and is implicated in both microtubule binding and sensing nucleotide bound to the active site of kinesin. Consistent with its proposed role as a microtubule sensor, loop L11 is rarely seen in crystal structures of unattached kinesins. Here, we report four structures of a regulated plant kinesin, the kinesin-like calmodulin binding protein (KCBP), determined by X-ray crystallography. Although all structures reveal the kinesin motor in the ATP-like conformation, its loop L11 is observed in different conformational states, both ordered and disordered. When structured, loop L11 adds three additional helical turns to the N-terminal part of the following helix α4. Although interactions with protein neighbors in the crystal support the ordering of loop L11, its observed conformation suggests the conformation for loop L11 in the microtubule-bound kinesin. Variations in the positions of other features of these kinesins were observed. A critical regulatory element of this kinesin, the calmodulin binding helix positioned at the C-terminus of the motor domain, is thought to confer negative regulation of KCBP. Calmodulin binds to this helix and inserts itself between the motor and the microtubule. Comparison of five independent structures of KCBP shows that the positioning of the calmodulin binding helix is not decided by crystal packing forces but is determined by the conformational state of the motor. The observed variations in the position of the calmodulin binding helix fit the regulatory mechanism previously proposed for this kinesin motor.     

Alpha7 Helix Plays an Important Role in the Conformational Stability of PTP1B

Abstract

The C-terminus of Protein Tyrosine Phosphatase 1B (PTP1B) includes an α-helix (α7), which forms an allosteric binding site 20 Å away from the active site. This helix is specific to PTP1B and its truncation decreases the catalytic activity significantly. Here, molecular dynamics (MD) simulations in the presence and absence of α7 were performed to investigate the role played by α7. The highly mobile α7 was found to maintain its contacts with loop 11 (L11)- α3 helix throughout the simulations. The interactions of Tyr152 on L11, Tyr176, Thr177 on the catalytically important WPD loop and Ser190 on α3 are important for the conformational stability and the concerted motions of the regions surrounding the WPD loop. In the absence of α7, L11 and WPD loop move away from their crystal structure conformations, resulting in the loss of the interactions in this region, and a decrease in the residue displacement correlations in the vicinity of WPD loop. Therefore, we suggest that one of the functionally important roles of α7 may be to limit the L11 and α3 motions, and, facilitate the WPD loop motions. Truncation of α7 in PTP1B is found to affect distant regions as well, such as the substrate recognition site and the phosphate binding-loop (P-loop), changing the conformations of these regions significantly. Our results show that the PTP1B specific α7 is important for the conformation and dynamics of the WPD loop, and also may play a role in ligand binding.     

Allosteric Drug Discrimination Is Coupled to Mechanochemical Changes in the Kinesin-5 Motor Core

Essential in mitosis, the human Kinesin-5 protein is a target for >80 classes of allosteric compounds that bind to a surface-exposed site formed by the L5 loop. Not established is why there are differing efficacies in drug inhibition. Here we compare the ligand-bound states of two L5-directed inhibitors against 15 Kinesin-5 mutants by ATPase assays and IR spectroscopy. Biochemical kinetics uncovers functional differences between individual residues at the N or C termini of the L5 loop. Infrared evaluation of solution structures and multivariate analysis of the vibrational spectra reveal that mutation and/or ligand binding not only can remodel the allosteric binding surface but also can transmit long range effects. Changes in L5-localized 3₁₀ helix and disordered content, regardless of substitution or drug potency, are experimentally detected. Principal component analysis couples these local structural events to two types of rearrangements in β-sheet hydrogen bonding. These transformations in β-sheet contacts are correlated with inhibitory drug response and are corroborated by wild type Kinesin-5 crystal structures. Despite considerable evolutionary divergence, our data directly support a theorized conserved element for long distance mechanochemical coupling in kinesin, myosin, and F₁-ATPase. These findings also suggest that these relatively rapid IR approaches can provide structural biomarkers for clinical determination of drug sensitivity and drug efficacy in nucleotide triphosphatases.     

Alpha7 helix plays an important role in the conformational stability of PTP1B

The C-terminus of Protein Tyrosine Phosphatase 1B (PTP1B) includes an α-helix α7), which forms an allosteric binding site 20 ? away from the active site. This helix is specific to PTP1B and its truncation decreases the catalytic activity significantly. Here, molecular dynamics (MD) simulations in the presence and absence of α7 were performed to investigate the role played by α7. The highly mobile α7 was found to maintain its contacts with loop 11 (L11)α3 helix throughout the simulations. The interactions of Tyr152 on L11, Tyr176, Thr177 on the catalytically important WPD loop and Ser190 on α3 are important for the conformational stability and the concerted motions of the regions surrounding the WPD loop. In the absence of α7, L11 and WPD loop move away from their crystal structure conformations, resulting in the loss of the interactions in this region, and a decrease in the residue displacement correlations in the vicinity of WPD loop. Therefore, we suggest that one of the functionally important roles of α7 may be to limit the L11 and α3 motions, and, facilitate the WPD loop motions. Truncation of α7 in PTP1B is found to affect distant regions as well, such as the substrate recognition site and the phosphate binding-loop (P-loop), changing the conformations of these regions significantly. Our results show that the PTP1B specific α7 is important for the conformation and dynamics of the WPD loop, and also may play a role in ligand binding.     

Proteolytic mapping of kinesin/ncd-microtubule interface: nucleotide-dependent conformational changes in the loops L8 and L12.

We used a battery of proteases to probe the footprint of microtubules on kinesin and ncd, and to search for nucleotide-induced conformational changes in these two oppositely-directed yet homologous molecular motors. Proteolytic cleavage sites were identified by N-terminal microsequencing and electrospray mass spectrometry, and then mapped onto the recently-determined atomic structures of ncd and kinesin. In both kinesin and ncd, microtubule binding shields a set of cleavage sites within or immediately flanking the loops L12, L8 and L11 and, in ncd, the loop L2. Even in the absence of microtubules, exchange of ADP for AMPPNP in the motor active site drives conformational shifts involving these loops. In ncd, a chymotryptic cleavage at Y622 in L12 is protected in the strong binding AMPPNP conformation, but cleaved in the weak binding ADP conformation. In kinesin, a thermolysin cleavage at L154 in L8 is protected in AMPPNP but cleaved in ADP. We speculate that ATP turnover in the active site governs microtubule binding by cyclically retracting or displaying the loops L8 and L12. Curiously, the retracted state of the loops corresponds to microtubule strong binding. Conceivably, nucleotide-dependent display of loops works as a reversible block on strong binding.     

Structural insights into the catalytic active site and activity of human Nit2/ω-amidase: kinetic assay and molecular dynamics simulation

Human nitrilase-like protein 2 (hNit2) is a putative tumor suppressor, recently identified as ω-amidase. hNit2/ω-amidase plays a crucial metabolic role by catalyzing the hydrolysis of α-ketoglutaramate (the α-keto analog of glutamine) and α-ketosuccinamate (the α-keto analog of asparagine), yielding α-ketoglutarate and oxaloacetate, respectively. Transamination between glutamine and α-keto-γ-methiolbutyrate closes the methionine salvage pathway. Thus, hNit2/ω-amidase links sulfur metabolism to the tricarboxylic acid cycle. To elucidate the catalytic specificity of hNit2/ω-amidase, we performed molecular dynamics simulations on the wild type enzyme and its mutants to investigate enzyme-substrate interactions. Binding free energies were computed to characterize factors contributing to the substrate specificity. The predictions resulting from these computations were verified by kinetic analyses and mutational studies. The activity of hNit2/ω-amidase was determined with α-ketoglutaramate and succinamate as substrates. We constructed three catalytic triad mutants (E43A, K112A, and C153A) and a mutant with a loop 116-128 deletion to validate the role of key residues and the 116-128 loop region in substrate binding and turnover. The molecular dynamics simulations successfully verified the experimental trends in the binding specificity of hNit2/ω-amidase toward various substrates. Our findings have revealed novel structural insights into the binding of substrates to hNit2/ω-amidase. A catalytic triad and the loop residues 116-128 of hNit2 play an essential role in supporting the stability of the enzyme-substrate complex, resulting in the generation of the catalytic products. These observations are predicted to be of benefit in the design of new inhibitors or activators for research involving cancer and hyperammonemic diseases.     

The Structural Basis of Force Generation by the Mitotic Motor Kinesin-5

Kinesin-5 is required for forming the bipolar spindle during mitosis. Its motor domain, which contains nucleotide and microtubule binding sites and mechanical elements to generate force, has evolved distinct properties for its spindle-based functions. In this study, we report subnanometer resolution cryoelectron microscopy reconstructions of microtubule-bound human kinesin-5 before and after nucleotide binding and combine this information with studies of the kinetics of nucleotide-induced neck linker and cover strand movement. These studies reveal coupled, nucleotide-dependent conformational changes that explain many of this motor''s properties. We find that ATP binding induces a ratchet-like docking of the neck linker and simultaneous, parallel docking of the N-terminal cover strand. Loop L5, the binding site for allosteric inhibitors of kinesin-5, also undergoes a dramatic reorientation when ATP binds, suggesting that it is directly involved in controlling nucleotide binding. Our structures indicate that allosteric inhibitors of human kinesin-5, which are being developed as anti-cancer therapeutics, bind to a motor conformation that occurs in the course of normal function. However, due to evolutionarily defined sequence variations in L5, this conformation is not adopted by invertebrate kinesin-5s, explaining their resistance to drug inhibition. Together, our data reveal the precision with which the molecular mechanism of kinesin-5 motors has evolved for force generation.     

Multiple conformations of the nucleotide site of Kinesin family motors in the triphosphate state

Identifying conformational changes in kinesin family motors associated with nucleotide and microtubule (MT) binding is essential to determining an atomic-level model for force production and motion by the motors. Using the mobility of nucleotide analog spin probes bound at the active sites of kinesin family motors to monitor conformational changes, we previously demonstrated that, in the ADP state, the open nucleotide site closes upon MT binding [Naber, N., Minehardt, T. J., Rice, S., Chen, X., Grammer, J., Matuska, M., et al. (2003). Closing of the nucleotide pocket of kinesin family motors upon binding to microtubules. Science, 300, 798-801]. We now extend these studies to kinesin-1 (K) and ncd (nonclaret disjunctional protein) motors in ATP and ATP-analog states. Our results reveal structural differences between several triphosphate and transition-state analogs bound to both kinesin and ncd in solution. The spectra of kinesin/ncd in the presence of SLADP•AlFx/BeFx and kinesin, with the mutation E236A (K-E236A; does not hydrolyze ATP) bound to ATP, show an open conformation of the nucleotide pocket similar to that seen in the kinesin/ncd•ADP states. In contrast, the triphosphate analogs K•SLAMPPNP and K-E236A•SLAMPPNP induce a more immobilized component of the electron paramagnetic resonance spectrum, implying closing of the nucleotide site. The MT-bound states of all of the triphosphate analogs reveal two novel spectral components. The equilibrium between these two components is only weakly dependent on temperature. Both components have more restricted mobility than observed in MT-bound diphosphate states. Thus, the closing of the nucleotide pocket when the diphosphate state binds to MTs is amplified in the triphosphate state, perhaps promoting accelerated ATP hydrolysis. Consistent with this idea, molecular dynamics simulations show a good correlation between our spectroscopic data, X-ray crystallography, and the electron microscopy of MT-bound triphosphate-analog states.     

Molecular mechanism of ADP-ribose hydrolysis by human NUDT5 from structural and kinetic studies

Human NUDT5 (hNUDT5) is an ADP-ribose (ADPR) pyrophosphatase (ADPRase) that plays important roles in controlling the intracellular levels of ADPR and preventing non-enzymatic ADP-ribosylation of proteins by hydrolyzing ADPR to AMP and ribose 5′-phosphate. We report the crystal structure of hNUDT5 in complex with a non-hydrolyzable ADPR analogue, α,β-methyleneadenosine diphosphoribose, and three Mg^2 + ions representing the transition state of the enzyme during catalysis. Analysis of this structure and comparison with previously reported hNUDT5 structures identify key residues involved in substrate binding and catalysis. In the transition-state structure, three metal ions are bound at the active site and are coordinated by surrounding residues and water molecules. A conserved water molecule is at an ideal position for nucleophilic attack on the α-phosphate of ADPR. The side chain of Glu166 on loop L9 changes its conformation to interact with the conserved water molecule compared with that in the substrate-bound structure and appears to function as a catalytic base. Mutagenesis and kinetic studies show that Trp28 and Trp46 are important for the substrate binding; Arg51 is involved in both the substrate binding and the catalysis; and Glu112 and Glu116 of the Nudix motif, Glu166 on loop L9, and Arg111 are critical for the catalysis. The structural and biochemical data together reveal the molecular basis of the catalytic mechanism of ADPR hydrolysis by hNUDT5. Specifically, Glu166 functions as a catalytic base to deprotonate a conserved water molecule that acts as a nucleophile to attack the α-phosphate of ADPR, and three Mg^2 + ions are involved in the activation of the nucleophile and the binding of the substrate. Structural comparison of different ADPRases also suggests that most dimeric ADPRases may share a similar catalytic mechanism of ADPR hydrolysis.     

Mapping the Structural and Dynamical Features of Kinesin Motor Domains

Kinesin motor proteins drive intracellular transport by coupling ATP hydrolysis to conformational changes that mediate directed movement along microtubules. Characterizing these distinct conformations and their interconversion mechanism is essential to determining an atomic-level model of kinesin action. Here we report a comprehensive principal component analysis of 114 experimental structures along with the results of conventional and accelerated molecular dynamics simulations that together map the structural dynamics of the kinesin motor domain. All experimental structures were found to reside in one of three distinct conformational clusters (ATP-like, ADP-like and Eg5 inhibitor-bound). These groups differ in the orientation of key functional elements, most notably the microtubule binding α4–α5, loop8 subdomain and α2b-β4-β6-β7 motor domain tip. Group membership was found not to correlate with the nature of the bound nucleotide in a given structure. However, groupings were coincident with distinct neck-linker orientations. Accelerated molecular dynamics simulations of ATP, ADP and nucleotide free Eg5 indicate that all three nucleotide states could sample the major crystallographically observed conformations. Differences in the dynamic coupling of distal sites were also evident. In multiple ATP bound simulations, the neck-linker, loop8 and the α4–α5 subdomain display correlated motions that are absent in ADP bound simulations. Further dissection of these couplings provides evidence for a network of dynamic communication between the active site, microtubule-binding interface and neck-linker via loop7 and loop13. Additional simulations indicate that the mutations G325A and G326A in loop13 reduce the flexibility of these regions and disrupt their couplings. Our combined results indicate that the reported ATP and ADP-like conformations of kinesin are intrinsically accessible regardless of nucleotide state and support a model where neck-linker docking leads to a tighter coupling of the microtubule and nucleotide binding regions. Furthermore, simulations highlight sites critical for large-scale conformational changes and the allosteric coupling between distal functional sites.     

Identification of novel scaffolds to inhibit human mitotic kinesin Eg5 targeting the second allosteric binding site using in silico methods

Human mitotic kinesins are potential anticancer drug targets because of their essential role in mitotic cell division. The kinesin Eg5 (Kinesin-5, kif11) has gained much attention in this regard and has many inhibitors in different phases of clinical trials. All drug candidates considered for Eg5 so far binds to the binding site (Site 1) formed by the loop L5, helices α2 and α3 and are uncompetitive to ATP/ADP. Recently, it has been reported that Eg5 also has a second binding site (Site 2) formed by helices α4 and α6. In the current work, we have screened the compounds in the diversity set-III from National Cancer Institute (NCI) and Zinc database to identify potential inhibitors for Eg5 that specifically binds to the site 2. The compounds were ranked based on the glide extra precision docking scores and the top ranked compounds were found to have pyridazine scaffold. The top five compounds were further evaluated for other drug like properties. Stability of protein-ligand complexes were analyzed using molecular dynamic simulations. Our studies suggest that pyridazine analogs have good MDCK, permeability properties and high binding affinity to the human Eg5.     

Dynamics of the active site loops in catalyzing aminoacylation reaction in seryl and histidyl tRNA synthetases

Aminoacylation reaction is the first step of protein biosynthesis. The catalytic reorganization at the active site of aminoacyl tRNA synthetases (aaRSs) is driven by the loop motions. There remain lacunae of understanding concerning the catalytic loop dynamics in aaRSs. We analyzed the functional loop dynamics in seryl tRNA synthetase from Methanopyrus kandleri (^mkSerRS) and histidyl tRNA synthetases from Thermus thermophilus (^ttHisRS), respectively, using molecular dynamics. Results confirm that the motif 2 loop and other active site loops are flexible spots within the catalytic domain. Catalytic residues of the loops form a network of interaction with the substrates to form a reactive state. The loops undergo transitions between closed state and open state and the relaxation of the constituent residues occurs in femtosecond to nanosecond time scale. Order parameters are higher for constituent catalytic residues which form a specific network of interaction with the substrates to form a reactive state compared to the Gly residues within the loop. The development of interaction is supported from mutation studies where the catalytic domain with mutated loop exhibits unfavorable binding energy with the substrates. During the open-close motion of the loops, the catalytic residues make relaxation by ultrafast librational motion as well as fast diffusive motion and subsequently relax rather slowly via slower diffusive motion. The Gly residues act as a hinge to facilitate the loop closing and opening by their faster relaxation behavior. The role of bound water is analyzed by comparing implicit solvent-based and explicit solvent-based simulations. Loops fail to form catalytically competent geometry in absence of water. The present result, for the first time reveals the nature of the active site loop dynamics in aaRS and their influence on catalysis.     

Identification of the protein binding region of S-trityl-L-cysteine, a new potent inhibitor of the mitotic kinesin Eg5

Human Eg5, a mitotic motor of the kinesin superfamily, is involved in the formation and maintenance of the mitotic spindle. The recent discovery of small molecules that inhibit HsEg5 by binding to its catalytic motor domain leading to mitotic arrest has attracted more interest in Eg5 as a potential anticancer drug target. We have used hydrogen-deuterium exchange mass spectrometry and directed mutagenesis to identify the secondary structure elements that form the binding sites of new Eg5 inhibitors, in particular for S-trityl-l-cysteine, a potent inhibitor of Eg5 activity in vitro and in cell-based assays. The binding of this inhibitor modifies the deuterium incorporation rate of eight peptides that define two areas within the motor domain: Tyr125-Glu145 and Ile202-Leu227. Replacement of the Tyr125-Glu145 region with the equivalent region in the Neurospora crassa conventional kinesin heavy chain prevents the inhibition of the Eg5 ATPase activity by S-trityl-l-cysteine. We show here that S-trityl-l-cysteine and monastrol both bind to the same region on Eg5 by induced fit in a pocket formed by helix alpha3-strand beta5 and loop L5-helix alpha2, and both inhibitors trigger similar local conformational changes within the interaction site. It is likely that S-trityl-l-cysteine and monastrol inhibit HsEg5 by a similar mechanism. The common inhibitor binding region appears to represent a "hot spot" for HsEg5 that could be exploited for further inhibitor screening.     

Deprotonation states of the two active site water molecules regulate the binding of protein phosphatase 5 with its substrate: A molecular dynamics study

Protein phosphatase 5 (PP5), mainly localized in human brain, can dephosphorylate tau protein whose high level of phosphorylation is related to Alzheimer's disease. Similar to other protein phosphatases, PP5 has a conserved motif in the catalytic domain that contains two binding sites for manganese (Mn²⁺) ions. Structural data indicate that two active site water molecules, one bridging the two Mn²⁺ ions and the other terminally coordinated with one of the Mn²⁺ ions (Mn1), are involved in catalysis. Recently, a density functional theory study revealed that the two water molecules can be both deprotonated to keep a neutral active site for catalysis. The theoretical study gives us an insight into the catalytic mechanism of PP5, but the knowledge of how the deprotonation states of the two water molecules affect the binding of PP5 with its substrate is still lacking. To approach this problem, molecular dynamics simulations were performed to model the four possible deprotonation states. Through structural, dynamical and energetic analyses, the results demonstrate that the deprotonation states of the two water molecules affect the structure of the active site including the distance between the two Mn²⁺ ions and their coordination, impact the interaction energy of residues R275, R400 and H304 which directly interact with the substrate phosphoserine, and mediate the dynamics of helix αJ which is involved in regulation of the enzyme's activity. Furthermore, the deprotonation state that is preferable for PP5 binding of its substrate has been identified. These findings could provide new design strategy for PP5 inhibitor.     

The role of loop ZA and Pro371 in the function of yeast Gcn5p bromodomain revealed through molecular dynamics and experiment

Biological experiments were combined with molecular dynamics simulations to understand the importance of amino acidic residues present in the bromodomain of the yeast histone acetyltransferase Gcn5p. It was found that residue Pro371 plays an important role in the molecular recognition of the acetylated histone H4 tail by Gcn5p bromodomain. Crystallographic analysis of the complex showed that this residue does not directly interact with the histone substrate. It has been demonstrated that a double mutation Pro371Thr and Met372Ala in the Gcn5p bromodomain impairs chromatin remodeling activity. It is demonstrated here that, in this double mutant and in the fully deleted bromodomain strain, there is lower growth under amino acid deprivation conditions. By in vitro surface plasmon resonance (Biacore) experiments it is shown that the binding affinity of the double mutation to acetyl lysine 16 histone H4 peptide decreases. Molecular dynamics simulations were used to explain this loss in acetyl lysine-Gcn5p bromodomain affinity, in the double mutant. By comparing nanosecond molecular dynamics trajectories of the native as well as the single and doubly mutated bromodomain, it is concluded that the presence of Pro371 is important to the functionality of the Gcn5p bromodomain. In the simulation a point mutation involving this highly conserved residue induced an increase in the flexibility of the ZA loop, which in turn modulated the exposure of the binding pocket to the acetyl lysine. The combined double mutations (Pro371Thr-Met372Ala) not only markedly perturb the motion of the ZA loop but also destabilize the entire structure of the bromodomain.