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

Binding patterns of inhibitors to different pockets of kinesin Eg5

The kinesin-5 family member, Eg5, plays very important role in the mitosis. As a mitotic protein, Eg5 is the target of various mitotic inhibitors. There are two targeting pockets in the motor domain of Eg5, which locates in the α2/L5/α3 region and the α4/α6 region respectively. We investigated the interactions between the different inhibitors and the two binding pockets of Eg5 by using all-atom molecular dynamics method. Combined the conformational analysis with the free-energy calculation, the binding patterns of inhibitors to the two binding pockets are shown. The α2/L5/α3 pocket can be divided into 4 regions. The structures and binding conformations of inhibitors in region 1 and 2 are highly conserved. The shape of α4/α6 pocket is alterable. The space of this pocket in ADP-binding state of Eg5 is larger than that in ADP·Pi-binding state due to the limitation of a hydrogen bond formed in the ADP·Pi-binding state. The results of this investigation provide the structural basis of the inhibitor-Eg5 interaction and offer a reference for the Eg5-targeted drug design.

Bound ion effects: Using machine learning method to study the kinesin Ncd's binding with microtubule

Drosophila Ncd proteins are motor proteins that play important roles in spindle organization. Ncd and the tubulin dimer are highly charged. Thus, it is crucial to investigate Ncd-tubulin dimer interactions in the presence of ions, especially ions that are bound or restricted at the Ncd-tubulin dimer binding interfaces. To consider the ion effects, widely used implicit solvent models treat ions implicitly in the continuous solvent environment without focusing on the individual ions' effects. But highly charged biomolecules such as the Ncd and tubulin dimer may capture some ions at highly charged regions as bound ions. Such bound ions are restricted to their binding sites; thus, they can be treated as part of the biomolecules. By applying multiscale computational methods, including the machine-learning-based Hybridizing Ions Treatment-2 program, molecular dynamics simulations, DelPhi, and DelPhiForce, we studied the interaction between the Ncd motor domain and the tubulin dimer using a hybrid solvent model, which considers the bound ions explicitly and the other ions implicitly in the solvent environment. To identify the importance of treating bound ions explicitly, we also performed calculations using the implicit solvent model without considering the individual bound ions. We found that the calculations of the electrostatic features differ significantly between those of the hybrid solvent model and the pure implicit solvent model. The analyses show that treating bound ions at highly charged regions explicitly is crucial for electrostatic calculations. This work proposes a machine-learning-based approach to handle the bound ions using the hybrid solvent model. Such an approach is not only capable of handling kinesin-tubulin complexes but is also appropriate for other highly charged biomolecules, such as DNA/RNA, viral capsid proteins, etc.

Modeling the motion of disease-associated KIF1A heterodimers

KIF1A is a member of the kinesin-3 motor protein family that transports synaptic vesicle precursors in axons. Mutations in the Kif1a gene cause neuronal diseases. Most patients are heterozygous and have both mutated and intact KIF1A alleles, suggesting that heterodimers composed of wild-type KIF1A and mutant KIF1A are likely involved in pathogenesis. In this study, we propose mathematical models to describe the motility of KIF1A heterodimers composed of wild-type KIF1A and mutant KIF1A. Our models precisely describe run length, run time, and velocity of KIF1A heterodimers using a few parameters obtained from two homodimers. The first model is a simple hand-over-hand model in which stepping and detachment rates from a microtubule of each head are identical to those in the respective homodimers. Although the velocities of heterodimers expected from this model were in good agreement with the experimental results, this model underestimated the run lengths and run times of some heterodimeric motors. To address this discrepancy, we propose the tethered-head affinity model, in which we hypothesize a tethered head, in addition to a microtubule-binding head, contributes to microtubule binding in a vulnerable one-head-bound state. The run lengths and run times of the KIF1A heterodimers predicted by the tethered-head affinity model matched well with experimental results, suggesting a possibility that the tethered head affects the microtubule binding of KIF1A. Our models provide insights into how each head contributes to the processive movement of KIF1A and can be used to estimate motile parameters of KIF1A heterodimers.

Serine 106 preserves the tertiary structure, function, and stability of a cyclophilin from Staphylococcus aureus

SaCyp, a staphylococcal cyclophilin involved in both protein folding and pathogenesis, has a Ser residue at position 106 and a Trp residue at position 136. While Ser 106 of SaCyp aligned with a cyclosporin A (CsA) binding Ala residue, its Trp 136 aligned with a Trp or a Phe residue of most other cyclophilins. To demonstrate the exact roles of Ser 106 and Trp 136 in SaCyp, we have elaborately studied rCyp[S106A] and rCyp[W136A], two-point mutants of a recombinant SaCyp (rCyp) harboring an Ala substitution at positions 106 and 136, respectively. Of the mutants, rCyp[W136A] showed the rCyp-like CsA binding affinity and peptidyl-prolyl cis-trans isomerase (PPIase) activity. Conversely, the PPIase activity, CsA binding affinity, stability, tertiary structure, surface hydrophobicity, and Trp accessibility of rCyp[S106A] notably differed from those of rCyp. The computational experiments also reveal that the structure, dimension, and fluctuation of SaCyp are not identical to those of SaCyp[S106A]. Furthermore, Ser at position 106 of SaCyp, compared to Ala at the same position, formed a higher number of non-covalent bonds with CsA. Collectively, Ser 106 is an indispensable residue for SaCyp that keeps its tertiary structure, function, and stability intact.Communicated by Ramaswamy H. Sarma.

Kinesin-3 motors are fine-tuned at the molecular level to endow distinct mechanical outputs

BACKGROUND

Kinesin-3 family motors drive diverse cellular processes and have significant clinical importance. The ATPase cycle is integral to the processive motility of kinesin motors to drive long-distance intracellular transport. Our previous work has demonstrated that kinesin-3 motors are fast and superprocessive with high microtubule affinity. However, chemomechanics of these motors remain poorly understood.

RESULTS

We purified kinesin-3 motors using the Sf9-baculovirus expression system and demonstrated that their motility properties are on par with the motors expressed in mammalian cells. Using biochemical analysis, we show for the first time that kinesin-3 motors exhibited high ATP turnover rates, which is 1.3- to threefold higher compared to the well-studied kinesin-1 motor. Remarkably, these ATPase rates correlate to their stepping rate, suggesting a tight coupling between chemical and mechanical cycles. Intriguingly, kinesin-3 velocities (KIF1A > KIF13A > KIF13B > KIF16B) show an inverse correlation with their microtubule-binding affinities (KIF1A < KIF13A < KIF13B < KIF16B). We demonstrate that this differential microtubule-binding affinity is largely contributed by the positively charged residues in loop8 of the kinesin-3 motor domain. Furthermore, microtubule gliding and cellular expression studies displayed significant microtubule bending that is influenced by the positively charged insert in the motor domain, K-loop, a hallmark of kinesin-3 family.

CONCLUSIONS

Together, we propose that a fine balance between the rate of ATP hydrolysis and microtubule affinity endows kinesin-3 motors with distinct mechanical outputs. The K-loop, a positively charged insert in the loop12 of the kinesin-3 motor domain promotes microtubule bending, an interesting phenomenon often observed in cells, which requires further investigation to understand its cellular and physiological significance.

Single-Molecule Studies on the Motion and Force Generation of the Kinesin-3 Motor KIF1A

KIF1A is a neuron-specific member of the kinesin-3 family of microtubule (MT) plus-end-directed motor proteins. It powers the migration of nuclei in differentiating brain stem cells and the transport of synaptic precursors and dense core vesicles in axons. Its dysfunction causes severe neurodevelopmental and neurodegenerative diseases termed KIF1A-associated neurological disorders (KAND). KAND mutations span the entirety of the KIF1A protein sequence, of which the majority are located within the motor domain and are thus predicted to affect the motor's motility and force-generating properties. Unfortunately, the molecular etiologies of KAND remain poorly understood, in part because KIF1A's molecular mechanism remains unclear. Here, we describe detailed methods for how to express a tail-truncated dimeric KIF1A in E. coli cells and provide step-by-step protocols for performing single-molecule studies with total internal reflection fluorescence microscopy and optical tweezers assays, which, when combined with structure-function studies, help to decipher KIF1A's molecular mechanism.

Using a comprehensive approach to investigate the interaction between Kinesin-5/Eg5 and the microtubule

Kinesins are microtubule-based motor proteins that play important roles ranging from intracellular transport to cell division. Human Kinesin-5 (Eg5) is essential for mitotic spindle assembly during cell division. By combining molecular dynamics (MD) simulations with other multi-scale computational approaches, we systematically studied the interaction between Eg5 and the microtubule. We find the electrostatic feature on the motor domains of Eg5 provides attractive interactions to the microtubule. Additionally, the folding and binding energy analysis reveals that the Eg5 motor domain performs its functions best when in a weak acidic environment. Molecular dynamics analyses of hydrogen bonds and salt bridges demonstrate that, on the binding interfaces of Eg5 and the tubulin heterodimer, salt bridges play the most significant role in holding the complex. The salt bridge residues on the binding interface of Eg5 are mostly positive, while salt bridge residues on the binding interface of tubulin heterodimer are mostly negative. Such salt bridge residue distribution is consistent with electrostatic potential calculations. In contrast, the interface between α and β-tubulins is dominated by hydrogen bonds rather than salt bridges. Compared to the Eg5/α-tubulin interface, the Eg5/β-tubulin interface has a greater number of salt bridges and higher occupancy for salt bridges. This asymmetric salt bridge distribution may play a significant role in Eg5′s directionality. The residues involved in hydrogen bonds and salt bridges are identified in this work and may be helpful for anticancer drug design.

Monoallelic KIF1A-related disorders: a multicenter cross sectional study and systematic literature review

BACKGROUND

Monoallelic variants in the KIF1A gene are associated with a large set of clinical phenotypes including neurodevelopmental and neurodegenerative disorders, underpinned by a broad spectrum of central and peripheral nervous system involvement.

METHODS

In a multicenter study conducted in patients presenting spastic gait or complex neurodevelopmental disorders, we analyzed the clinical, genetic and neuroradiological features of 28 index cases harboring heterozygous variants in KIF1A. We conducted a literature systematic review with the aim to comparing our findings with previously reported KIF1A-related phenotypes.

RESULTS

Among 28 patients, we identified nine novel monoallelic variants, and one a copy number variation encompassing KIF1A. Mutations arose de novo in most patients and were prevalently located in the motor domain. Most patients presented features of a continuum ataxia-spasticity spectrum with only five cases showing a prevalently pure spastic phenotype and six presenting congenital ataxias. Seventeen mutations occurred in the motor domain of the Kinesin-1A protein, but location of mutation did not correlate with neurological and imaging presentations. When tested in 15 patients, muscle biopsy showed oxidative metabolism alterations (6 cases), impaired respiratory chain complexes II + III activity (3/6) and low CoQ10 levels (6/9). Ubiquinol supplementation (1gr/die) was used in 6 patients with subjective benefit.

CONCLUSIONS

This study broadened our clinical, genetic, and neuroimaging knowledge of KIF1A-related disorders. Although highly heterogeneous, it seems that manifestations of ataxia-spasticity spectrum disorders seem to occur in most patients. Some patients also present secondary impairment of oxidative metabolism; in this subset, ubiquinol supplementation therapy might be appropriate.

A conserved arginine residue in a staphylococcal anti-sigma factor is required to preserve its kinase activity, structure, and stability

RsbW, σ^B, and RsbV, encoded by Staphylococcus aureus and related bacteria, act as an anti-sigma factor, an sigma factor, and an anti-anti-sigma factor, respectively. The interaction between RsbW and σ^B blocks the transcription initiation activity of the latter protein. RsbW also functions as a serine kinase and phosphorylates RsbV in the presence of ATP. Our modeling study indicates that the RsbW-RsbV complex is stabilized by twenty-four intermolecular non-covalent bonds. Of the bond-forming RsbW residues, Arg 23, and Glu 49 are conserved residues. To understand the roles of Arg 23 in RsbW, rRsbW[R23A], a recombinant S. aureus RsbW (rRsbW) harboring Arg to Ala change at position 23, was investigated using various probes. The results reveal that rRsbW[R23A], like rRsbW, exists as the dimers in the aqueous solution. However, rRsbW[R23A], unlike rRsbW, neither interacted with a chimeric RsbV (rRsbV) nor formed the phosphorylated rRsbV in the presence of ATP. Furthermore, the tertiary structure and hydrophobic surface area of rRsbW[R23A] matched little with those of rRsbW. Conversely, both rRsbW[R23A] and rRsbW showed interaction with a recombinant σ^B (rσ^B). rRsbW and rRsbW[R23A] were also unfolded via the formation of at least one intermediate in the presence of urea. However, the thermodynamic stability of rRsbW significantly differed from that of rRsbW[R23A]. Our molecular dynamics (MD) simulation study also reveals the substantial change of structure, dimension, and stability of RsbW due to the above mutation. The ways side chain of critical Arg 23 contributes to maintaining the tertiary structure, and stability of RsbW was elaborately discussed.Communicated by Ramaswamy H. Sarma.

Molecular dynamics simulation of tropomyosin bound to actins/myosin in the closed and open states

Tropomyosin (Tpm) is a dimeric coiled-coil protein that binds to filamentous actin, and regulates actin-myosin interaction by moving between three positions corresponding to the blocked, closed, and open states. To elucidate how Tpm undergoes transitions between these functional states, we have built structural models and conducted extensive molecular dynamics simulations of the Tpm-actins/myosin complex in the closed and open states (total simulation time >1.4 μs). Based on the simulation trajectories, we have analyzed the dynamics and energetics of a truncated Tpm interacting with actins/myosin under the physiological conditions. Our simulations have shown distinct dynamics of four Tpm periods (P3-P6), featuring pronounced biased fluctuations of P4 and P5 toward the open position in the closed state, which is consistent with a conformational selection mechanism for Tpm-regulated myosin binding. Additionally, we have identified key residues of Tpm specifically binding to actins/myosin in the closed and open state. Some of them were validated as functionally important in comparison with past functional/clinical studies, and the rest will make promising targets for future mutational experiments.

Computational Approaches to Prioritize Cancer Driver Missense Mutations

Cancer is a complex disease that is driven by genetic alterations. There has been a rapid development of genome-wide techniques during the last decade along with a significant lowering of the cost of gene sequencing, which has generated widely available cancer genomic data. However, the interpretation of genomic data and the prediction of the association of genetic variations with cancer and disease phenotypes still requires significant improvement. Missense mutations, which can render proteins non-functional and provide a selective growth advantage to cancer cells, are frequently detected in cancer. Effects caused by missense mutations can be pinpointed by in silico modeling, which makes it more feasible to find a treatment and reverse the effect. Specific human phenotypes are largely determined by stability, activity, and interactions between proteins and other biomolecules that work together to execute specific cellular functions. Therefore, analysis of missense mutations’ effects on proteins and their complexes would provide important clues for identifying functionally important missense mutations, understanding the molecular mechanisms of cancer progression and facilitating treatment and prevention. Herein, we summarize the major computational approaches and tools that provide not only the classification of missense mutations as cancer drivers or passengers but also the molecular mechanisms induced by driver mutations. This review focuses on the discussion of annotation and prediction methods based on structural and biophysical data, analysis of somatic cancer missense mutations in 3D structures of proteins and their complexes, predictions of the effects of missense mutations on protein stability, protein-protein and protein-nucleic acid interactions, and assessment of conformational changes in protein conformations induced by mutations.

Forces and Disease: Electrostatic force differences caused by mutations in kinesin motor domains can distinguish between disease-causing and non-disease-causing mutations

The ability to predict if a given mutation is disease-causing or not has enormous potential to impact human health. Typically, these predictions are made by assessing the effects of mutation on macromolecular stability and amino acid conservation. Here we report a novel feature: the electrostatic component of the force acting between a kinesin motor domain and tubulin. We demonstrate that changes in the electrostatic component of the binding force are able to discriminate between disease-causing and non-disease-causing mutations found in human kinesin motor domains using the receiver operating characteristic (ROC). Because diseases may originate from multiple effects not related to kinesin-microtubule binding, the prediction rate of 0.843 area under the ROC plot due to the change in magnitude of the electrostatic force alone is remarkable. These results reflect the dependence of kinesin’s function on motility along the microtubule, which suggests a precise balance of microtubule binding forces is required.

Investigation of alanine mutations affecting insulin-like growth factor (IGF) I binding to IGF binding proteins

Binding properties of wild type (WT) and six single amino acid substituted variants (E3A, E9A, D12A, D20A, F23A, and E58A) of insulin-like growth factor I (IGF-I) were analyzed with respect to their binding details to IGF binding proteins (IGFBPs) by molecular dynamics (MD) simulations. The binding sites and binding interactions on IGF-I and IGFBPs are screened and compared with the static X-ray structure. Electrostatic interaction is the primary driving force of the interaction between IGF-I and IGFBPs. Mutation may cause the rearrangement of binding sites, however, the unfolding of protein induced by mutation is not obvious in this work. We also provide the detailed picture of binding factors. And the results show that, whether the unfolding of helix occurs or not, the Ala mutation will change the molecular atmosphere of the binding interface by the rearrangement of conformation, and further affects the binding residues and binding interactions.

Decrypting the structural, dynamic, and energetic basis of a monomeric kinesin interacting with a tubulin dimer in three ATPase states by all-atom molecular dynamics simulation

We have employed molecular dynamics (MD) simulation to investigate, with atomic details, the structural dynamics and energetics of three major ATPase states (ADP, APO, and ATP state) of a human kinesin-1 monomer in complex with a tubulin dimer. Starting from a recently solved crystal structure of ATP-like kinesin-tubulin complex by the Knossow lab, we have used flexible fitting of cryo-electron-microscopy maps to construct new structural models of the kinesin-tubulin complex in APO and ATP state, and then conducted extensive MD simulations (total 400 ns for each state), followed by flexibility analysis, principal component analysis, hydrogen bond analysis, and binding free energy analysis. Our modeling and simulation have revealed key nucleotide-dependent changes in the structure and flexibility of the nucleotide-binding pocket (featuring a highly flexible and open switch I in APO state) and the tubulin-binding site, and allosterically coupled motions driving the APO to ATP transition. In addition, our binding free energy analysis has identified a set of key residues involved in kinesin-tubulin binding. On the basis of our simulation, we have attempted to address several outstanding issues in kinesin study, including the possible roles of β-sheet twist and neck linker docking in regulating nucleotide release and binding, the structural mechanism of ADP release, and possible extension and shortening of α4 helix during the ATPase cycle. This study has provided a comprehensive structural and dynamic picture of kinesin's major ATPase states, and offered promising targets for future mutational and functional studies to investigate the molecular mechanism of kinesin motors.

The family-specific K-loop influences the microtubule on-rate but not the superprocessivity of kinesin-3 motors

The kinesin-3 family–specific, positively charged insert, the K-loop, in loop 12 of the motor domain plays a critical role in cargo transport by enhancing the initial interaction of cargo-bound dimeric motors with the microtubule. The replacement of the K-loop, however, does not abolish the superprocessive motion of this class of kinesin motors.

The kinesin-3 family (KIF) is one of the largest among the kinesin superfamily and an important driver of a variety of cellular transport events. Whereas all kinesins contain the highly conserved kinesin motor domain, different families have evolved unique motor features that enable different mechanical and functional outputs. A defining feature of kinesin-3 motors is the presence of a positively charged insert, the K-loop, in loop 12 of their motor domains. However, the mechanical and functional output of the K-loop with respect to processive motility of dimeric kinesin-3 motors is unknown. We find that, surprisingly, the K-loop plays no role in generating the superprocessive motion of dimeric kinesin-3 motors (KIF1, KIF13, and KIF16). Instead, we find that the K-loop provides kinesin-3 motors with a high microtubule affinity in the motor's ADP-bound state, a state that for other kinesins binds only weakly to the microtubule surface. A high microtubule affinity results in a high landing rate of processive kinesin-3 motors on the microtubule surface. We propose that the family-specific K-loop contributes to efficient kinesin-3 cargo transport by enhancing the initial interaction of dimeric motors with the microtubule track.

Probing the flexibility of tropomyosin and its binding to filamentous actin using molecular dynamics simulations

Tropomyosin (Tm) is a coiled-coil protein that binds to filamentous actin (F-actin) and regulates its interactions with actin-binding proteins like myosin by moving between three positions on F-actin (the blocked, closed, and open positions). To elucidate the molecular details of Tm flexibility in relation to its binding to F-actin, we conducted extensive molecular dynamics simulations for both Tm alone and Tm-F-actin complex in the presence of explicit solvent (total simulation time >400 ns). Based on the simulations, we systematically analyzed the local flexibility of the Tm coiled coil using multiple parameters. We found a good correlation between the regions with high local flexibility and a number of destabilizing regions in Tm, including six clusters of core alanines. Despite the stabilization by F-actin binding, the distribution of local flexibility in Tm is largely unchanged in the absence and presence of F-actin. Our simulations showed variable fluctuations of individual Tm periods from the closed position toward the open position. In addition, we performed Tm-F-actin binding calculations based on the simulation trajectories, which support the importance of Tm flexibility to Tm-F-actin binding. We identified key residues of Tm involved in its dynamic interactions with F-actin, many of which have been found in recent mutational studies to be functionally important, and the rest of which will make promising targets for future mutational experiments.

Predicting the Impact of Missense Mutations on Protein-Protein Binding Affinity

The crucial prerequisite for proper biological function is the protein's ability to establish highly selective interactions with macromolecular partners. A missense mutation that alters the protein binding affinity may cause significant perturbations or complete abolishment of the function, potentially leading to diseases. The availability of computational methods to evaluate the impact of mutations on protein-protein binding is critical for a wide range of biomedical applications. Here, we report an efficient computational approach for predicting the effect of single and multiple missense mutations on protein-protein binding affinity. It is based on a well-tested simulation protocol for structure minimization, modified MM-PBSA and statistical scoring energy functions with parameters optimized on experimental sets of several thousands of mutations. Our simulation protocol yields very good agreement between predicted and experimental values with Pearson correlation coefficients of 0.69 and 0.63 and root-mean-square errors of 1.20 and 1.90 kcal mol^-1 for single and multiple mutations, respectively. Compared with other available methods, our approach achieves high speed and prediction accuracy and can be applied to large datasets generated by modern genomics initiatives. In addition, we report a crucial role of water model and the polar solvation energy in estimating the changes in binding affinity. Our analysis also reveals that prediction accuracy and effect of mutations on binding strongly depends on the type of mutation and its location in a protein complex.

Mechanical amplification mechanism of kinesin's β-domain

Conventional kinesin's force generation process always takes place on the leading head and the generated force is transmitted to the trailing head through two neck linkers. To guarantee a strong force to be transmitted to the trailing head so that it can be detached from microtubule surface, the neck linker of the leading head must have a large enough forward displacement, which is proposed to be achieved by the amplifying function of the β-domain. However, the experimental result shows that the forward displacement of the β-domain itself appears too small. To elucidate the function of the β-domain, we make a detailed analysis of the mechanical relationship between the two motor heads and, based on the results of molecular dynamics simulation and mechanical analysis, we calculate the forward displacement of the neck linker of the leading head during the ATP binding induced motor head rotation. We show that β-domain achieves its amplifying function together with β0, so that neck linker can have a forward displacement during motor head rotation. This displacement of neck linker is large enough to cause detachment of the trailing head. Based on these results, a possible initiation mechanism of neck linker docking is proposed.

All-atom molecular dynamics simulations of actin-myosin interactions: a comparative study of cardiac α myosin, β myosin, and fast skeletal muscle myosin

Myosins are a superfamily of actin-binding motor proteins with significant variations in kinetic properties (such as actin binding affinity) between different isoforms. It remains unknown how such kinetic variations arise from the structural and dynamic tuning of the actin-myosin interface at the amino acid residue level. To address this key issue, we have employed molecular modeling and simulations to investigate, with atomistic details, the isoform dependence of actin-myosin interactions in the rigor state. By combining electron microscopy-based docking with homology modeling, we have constructed three all-atom models for human cardiac α and β and rabbit fast skeletal muscle myosin in complex with three actin subunits in the rigor state. Starting from these models, we have performed extensive all-atom molecular dynamics (MD) simulations (total of 100 ns per system) and then used the MD trajectories to calculate actin-myosin binding free energies with contributions from both electrostatic and nonpolar forces. Our binding calculations are in good agreement with the experimental finding of isoform-dependent differences in actin binding affinity between these myosin isoforms. Such differences are traced to changes in actin-myosin electrostatic interactions (i.e., hydrogen bonds and salt bridges) that are highly dynamic and involve several flexible actin-binding loops. By partitioning the actin-myosin binding free energy to individual myosin residues, we have also identified key myosin residues involved in the actin-myosin interactions, some of which were previously validated experimentally or implicated in cardiomyopathy mutations, and the rest make promising targets for future mutational experiments.

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.

Kinesin motor proteins transport cargo along microtubule tracks to support essential cellular functions including cell growth, axonal signaling and the separation of chromosomes during cell division. All kinesins contain one or more conserved motor domains that modulate binding and movement along microtubules via cycles of ATP hydrolysis. Important conformational transitions occurring during this cycle have been characterized with extensive crystallographic studies. However, the link between the observed conformations and the mechanisms involved in conformational change and microtubule interaction modulation remain unclear. Here we describe a comprehensive principal component analysis of available motor domain crystallographic structures supplemented with extensive unbiased conventional and accelerated molecular dynamics simulations that together characterize the response of kinesin motor domains to ATP binding and hydrolysis. Our studies reveal atomic details of conformational transitions, as well as novel nucleotide-dependent dynamical couplings, of distal regions and residues potentially important for the allosteric link between nucleotide and microtubule binding sites.

Motor domain phosphorylation modulates kinesin-1 transport

Disruptions in microtubule motor transport are associated with a variety of neurodegenerative diseases. Post-translational modification of the cargo-binding domain of the light and heavy chains of kinesin has been shown to regulate transport, but less is known about how modifications of the motor domain affect transport. Here we report on the effects of phosphorylation of a mammalian kinesin motor domain by the kinase JNK3 at a conserved serine residue (Ser-175 in the B isoform and Ser-176 in the A and C isoforms). Phosphorylation of this residue has been implicated in Huntington disease, but the mechanism by which Ser-175 phosphorylation affects transport is unclear. The ATPase, microtubule-binding affinity, and processivity are unchanged between a phosphomimetic S175D and a nonphosphorylatable S175A construct. However, we find that application of force differentiates between the two. Placement of negative charge at Ser-175, through phosphorylation or mutation, leads to a lower stall force and decreased velocity under a load of 1 piconewton or greater. Sedimentation velocity experiments also show that addition of a negative charge at Ser-175 favors the autoinhibited conformation of kinesin. These observations imply that when cargo is transported by both dynein and phosphorylated kinesin, a common occurrence in the cell, there may be a bias that favors motion toward the minus-end of microtubules. Such bias could be used to tune transport in healthy cells when properly regulated but contribute to a disease state when misregulated.

Mutations in DNA-binding loop of NFAT5 transcription factor produce unique outcomes on protein-DNA binding and dynamics

The nuclear factor of activated T cells 5 (NFAT5 or TonEBP) is a Rel family transcriptional activator and is activated by hypertonic conditions. Several studies point to a possible connection between nuclear translocation and DNA binding; however, the mechanism of NFAT5 nuclear translocation and the effect of DNA binding on retaining NFAT5 in the nucleus are largely unknown. Recent experiments showed that different mutations introduced in the DNA-binding loop and dimerization interface were important for DNA binding and some of them decreased the nuclear–cytoplasm ratio of NFAT5. To understand the mechanisms of these mutations, we model their effect on protein dynamics and DNA binding. We show that the NFAT5 complex without DNA is much more flexible than the complex with DNA. Moreover, DNA binding considerably stabilizes the overall dimeric complex and the NFAT5 dimer is only marginally stable in the absence of DNA. Two sets of NFAT5 mutations from the same DNA-binding loop are found to have different mechanisms of specific and nonspecific binding to DNA. The R217A/E223A/R226A (R293A/E299A/R302A using isoform c numbering) mutant is characterized by significantly compromised binding to DNA and higher complex flexibility. On the contrary, the T222D (T298D in isoform c) mutation, a potential phosphomimetic mutation, makes the overall complex more rigid and does not significantly affect the DNA binding. Therefore, the reduced nuclear–cytoplasm ratio of NFAT5 can be attributed to reduced binding to DNA for the triple mutant, while the T222D mutant suggests an additional mechanism at work.

Structure-based molecular simulations reveal the enhancement of biased Brownian motions in single-headed kinesin

Kinesin is a family of molecular motors that move unidirectionally along microtubules (MT) using ATP hydrolysis free energy. In the family, the conventional two-headed kinesin was experimentally characterized to move unidirectionally through “walking” in a hand-over-hand fashion by coordinated motions of the two heads. Interestingly a single-headed kinesin, a truncated KIF1A, still can generate a biased Brownian movement along MT, as observed by in vitro single molecule experiments. Thus, KIF1A must use a different mechanism from the conventional kinesin to achieve the unidirectional motions. Based on the energy landscape view of proteins, for the first time, we conducted a set of molecular simulations of the truncated KIF1A movements over an ATP hydrolysis cycle and found a mechanism exhibiting and enhancing stochastic forward-biased movements in a similar way to those in experiments. First, simulating stand-alone KIF1A, we did not find any biased movements, while we found that KIF1A with a large friction cargo-analog attached to the C-terminus can generate clearly biased Brownian movements upon an ATP hydrolysis cycle. The linked cargo-analog enhanced the detachment of the KIF1A from MT. Once detached, diffusion of the KIF1A head was restricted around the large cargo which was located in front of the head at the time of detachment, thus generating a forward bias of the diffusion. The cargo plays the role of a diffusional anchor, or cane, in KIF1A “walking.”

It is one of the major issues in biophysics how molecular motors such as conventional two-headed kinesin convert the chemical energy released at ATP hydrolysis into mechanical work. While most molecular motors move with more than one catalytic domain working in coordinated fashions, there are some motors that can move with only a single catalytic domain, which provides us a possibly simpler case to understand. A single-headed kinesin, KIF1A, with only one catalytic domain, has been characterized by in vitro single-molecule assay to generate a biased Brownian movement along the microtubule. Here, we conducted a set of structure-based coarse-grained molecular simulations for KIF1A system over an ATP hydrolysis cycle for the first time to our knowledge. Without cargo the simulated stand-alone KIF1A could not generate any directional movement, while a large-friction cargo-analog linked to the C-terminus of KIF1A clearly enhanced the forward-biased Brownian movement of KIF1A significantly. Interestingly, the cargo-analog here is not merely load but an important promoter to enable efficient movements of KIF1A.

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.